E2-E3 ubiquitin enzyme pairing - partnership in provoking or mitigating cancers

The ubiquitin-proteasome system (UPS) modulates carcinogenesis through ubiquitination of cancer-related target proteins, leading to their degradation in the proteasome. This may deactivate tumor suppressors or activate tumor promoters- either way causing homeostatic imbalance. As major components of the UPS, the E2 and E3 enzymes are recognized as pivotal determinants of substrate recognition and ubiquitination. Identification of E2-E3 pairing selectivity is particularly pertinent to early diagnosis and potential development of targeted cancer therapeutics. This review is motivated by recent findings and new insights into the molecular dynamics of ubiquitination triggered by specific E2-E3 pairing, leading to cancer initiation and progression if cancer suppressors are degraded or cancer suppression (if cancer promoters are degraded), respectively. We provide an overview of strategies employed in screening for E2-E3 interactions based on up-to-date studies focusing on the E2-E3 interface motifs. Of considerable recent interest is how E2 and E3 might switch their functional partnerships via UBE2O, which suggests an emerging significance on how UBE2O might influence E2-E3 pairing. Thus, a reflection on the role of UBE2O is included. Finally, we deliberate on the rational and cautious development of anti-cancer cocktail drugs which specifically target E2-E3 interacting residues for precision in cancer-killing with minimal side-effects. To this end, a list of potential future research is proposed.     

Dynamics Govern Specificity of a Protein-Protein Interface: Substrate Recognition by Thrombin

Biomolecular recognition is crucial in cellular signal transduction. Signaling is mediated through molecular interactions at protein-protein interfaces. Still, specificity and promiscuity of protein-protein interfaces cannot be explained using simplistic static binding models. Our study rationalizes specificity of the prototypic protein-protein interface between thrombin and its peptide substrates relying solely on binding site dynamics derived from molecular dynamics simulations. We find conformational selection and thus dynamic contributions to be a key player in biomolecular recognition. Arising entropic contributions complement chemical intuition primarily reflecting enthalpic interaction patterns. The paradigm “dynamics govern specificity” might provide direct guidance for the identification of specific anchor points in biomolecular recognition processes and structure-based drug design.     

Intrinsic α‐helical and β‐sheet conformational preferences: A computational case study of alanine

A fundamental question in protein science is what is the intrinsic propensity for an amino acid to be in an α-helix, β-sheet, or other backbone dihedral angle (-ψ) conformation. This question has been hotly debated for many years because including all protein crystal structures from the protein database, increases the probabilities for α-helical structures, while experiments on small peptides observe that β-sheet-like conformations predominate. We perform molecular dynamics (MD) simulations of a hard-sphere model for Ala dipeptide mimetics that includes steric interactions between nonbonded atoms and bond length and angle constraints with the goal of evaluating the role of steric interactions in determining protein backbone conformational preferences. We find four key results. For the hard-sphere MD simulations, we show that (1) β-sheet structures are roughly three and half times more probable than α-helical structures, (2) transitions between α-helix and β-sheet structures only occur when the backbone bond angle τ (N–C_α–C) is greater than 110°, and (3) the probability distribution of τ for Ala conformations in the “bridge” region of-ψ space is shifted to larger angles compared to other regions. In contrast, (4) the distributions obtained from Amber and CHARMM MD simulations in the bridge regions are broader and have increased τ compared to those for hard sphere simulations and from high-resolution protein crystal structures. Our results emphasize the importance of hard-sphere interactions and local stereochemical constraints that yield strong correlations between -ψ conformations and τ.     

Entropy-driven mechanism of an E3 ligase

Ubiquitin-like modifications are macromolecular chemistry for which our understanding of the enzymatic mechanisms is lacking. Most E3 ligases in ubiquitin-like modifications do not directly participate in chemistry but are thought to confer allosteric effects; however, the nature of the allosteric effects has been elusive. Recent molecular dynamics simulations suggested that an E3 binding enhances the population of the conformational states of the E2·SUMO thioester that favor reactions. In this study, we conducted the first temperature-dependent enzyme kinetic analysis to investigate the role of an E3 on activation entropy and enthalpy. The small ubiquitin-like modifier (SUMO) E3, RanBP2, confers unusually large, favorable activation entropy to lower the activation energy of the reaction. Mutants of RanBP2, designed to alter the flexibilities of the E2·SUMO thioester, showed a direct correlation of their favorable entropic effects with their ability to restrict the conformational flexibility of the E2·SUMO thioester. While the more favorable activation entropy is consistent with the previously suggested role of E3 in conformational selection, the large positive entropy suggests a significant role of solvent in catalysis. Indeed, molecular dynamics simulations in explicit water revealed that the more stable E2·SUMO thioester upon E3 binding results in stabilization of a large number of bound water molecules. Liberating such structured water at the transition state can result in large favorable activation entropy but unfavorable activation enthalpy. The entropy-driven mechanism of the E3 is consistent with the lack of structural conservation among E3s despite their similar functions. This study also illustrates how proteins that bind both SUMO and E2 can function as E3s and how intrinsically unstructured proteins can enhance macromolecular chemistry in addition to their known advantages in protein--protein interactions.     

Dynamic diagnosis of familial prion diseases supports the β2-α2 loop as a universal interference target

Background

Mutations in the cellular prion protein associated to familial prion disorders severely increase the likelihood of its misfolding into pathogenic conformers. Despite their postulation as incompatible elements with the native fold, these mutations rarely modify the native state structure. However they variably have impact on the thermodynamic stability and metabolism of PrP^C and on the properties of PrP^Sc aggregates. To investigate whether the pathogenic mutations affect the dynamic properties of the HuPrP(125-229) α-fold and find possible common patterns of effects that could help in prophylaxis we performed a dynamic diagnosis of ten point substitutions.

Methodology/Principal Findings

Using all-atom molecular dynamics simulations and novel analytical tools we have explored the effect of D178N, V180I, T183A, T188K, E196K, F198S, E200K, R208H, V210I and E211Q mutations on the dynamics of HuPrP(125-228) α-fold. We have found that while preserving the native state, all mutations produce dynamic changes which perturb the coordination of the α2-α3 hairpin to the rest of the molecule and cause the reorganization of the patches for intermolecular recognition, as the disappearance of those for conversion inhibitors and the emergence of an interaction site at the β2-α2 loop region.

Conclusions/Significance

Our results suggest that pathogenic mutations share a common pattern of dynamical alterations that converge to the conversion of the β2-α2 loop into an interacting region that can be used as target for interference treatments in genetic diseases.     

Bistable cell fate specification as a result of stochastic fluctuations and collective spatial cell behaviour

Background

In culture, isogenic mammalian cells typically display enduring phenotypic heterogeneity that arises from fluctuations of gene expression and other intracellular processes. This diversity is not just simple noise but has biological relevance by generating plasticity. Noise driven plasticity was suggested to be a stem cell-specific feature.

Results

Here we show that the phenotypes of proliferating tissue progenitor cells such as primary mononuclear muscle cells can also spontaneously fluctuate between different states characterized by the either high or low expression of the muscle-specific cell surface molecule CD56 and by the corresponding high or low capacity to form myotubes. Although this capacity is a cell-intrinsic property, the cells switch their phenotype under the constraints imposed by the highly heterogeneous microenvironment created by their own collective movement. The resulting heterogeneous cell population is characterized by a dynamic equilibrium between “high CD56” and “low CD56” phenotype cells with distinct spatial distribution. Computer simulations reveal that this complex dynamic is consistent with a context-dependent noise driven bistable model where local microenvironment acts on the cellular state by encouraging the cell to fluctuate between the phenotypes until the low noise state is found.

Conclusions

These observations suggest that phenotypic fluctuations may be a general feature of any non-terminally differentiated cell. The cellular microenvironment created by the cells themselves contributes actively and continuously to the generation of fluctuations depending on their phenotype. As a result, the cell phenotype is determined by the joint action of the cell-intrinsic fluctuations and by collective cell-to-cell interactions.     

Insights into diterpene cyclization from structure of bifunctional abietadiene synthase from Abies grandis

Abietadiene synthase from Abies grandis (AgAS) is a model system for diterpene synthase activity, catalyzing class I (ionization-initiated) and class II (protonation-initiated) cyclization reactions. Reported here is the crystal structure of AgAS at 2.3 Å resolution and molecular dynamics simulations of that structure with and without active site ligands. AgAS has three domains (α, β, and γ). The class I active site is within the C-terminal α domain, and the class II active site is between the N-terminal γ and β domains. The domain organization resembles that of monofunctional diterpene synthases and is consistent with proposed evolutionary origins of terpene synthases. Molecular dynamics simulations were carried out to determine the effect of substrate binding on enzymatic structure. Although such studies of the class I active site do lead to an enclosed substrate-Mg²⁺ complex similar to that observed in crystal structures of related plant enzymes, it does not enforce a single substrate conformation consistent with the known product stereochemistry. Simulations of the class II active site were more informative, with observation of a well ordered external loop migration. This “loop-in” conformation not only limits solvent access but also greatly increases the number of conformational states accessible to the substrate while destabilizing the nonproductive substrate conformation present in the “loop-out” conformation. Moreover, these conformational changes at the class II active site drive the substrate toward the proposed transition state. Docked substrate complexes were further assessed with regard to the effects of site-directed mutations on class I and II activities.     

Structural analysis of P450 AmphL from Streptomyces nodosus provides insights into substrate selectivity of polyene macrolide antibiotic biosynthetic P450s

AmphL is a cytochrome P450 enzyme that catalyzes the C8 oxidation of 8-deoxyamphotericin B to the polyene macrolide antibiotic, amphotericin B. To understand this substrate selectivity, we solved the crystal structure of AmphL to a resolution of 2.0 Å in complex with amphotericin B and performed molecular dynamics (MD) simulations. A detailed comparison with the closely related P450, PimD, which catalyzes the epoxidation of 4,5-desepoxypimaricin to the macrolide antibiotic, pimaricin, reveals key catalytic structural features responsible for stereo- and regio-selective oxidation. Both P450s have a similar access channel that runs parallel to the active site I helix over the surface of the heme. Molecular dynamics simulations of substrate binding reveal PimD can “pull” substrates further into the P450 access channel owing to additional electrostatic interactions between the protein and the carboxyl group attached to the hemiketal ring of 4,5-desepoxypimaricin. This substrate interaction is absent in AmphL although the additional substrate -OH groups in 8-deoxyamphotericin B help to correctly position the substrate for C8 oxidation. Simulations of the oxy-complex indicates that these -OH groups may also participate in a proton relay network required for O₂ activation as has been suggested for two other macrolide P450s, PimD and P450eryF. These findings provide experimentally testable models that can potentially contribute to a new generation of novel macrolide antibiotics with enhanced antifungal and/or antiprotozoal efficacy.     

Multiscale modeling of nucleosome dynamics

Nucleosomes form the fundamental building blocks of chromatin. Subtle modifications of the constituent histone tails mediate chromatin stability and regulate gene expression. For this reason, it is important to understand structural dynamics of nucleosomes at atomic levels. We report a novel multiscale model of the fundamental chromatin unit, a nucleosome, using a simplified model for rapid discrete molecular dynamics simulations and an all-atom model for detailed structural investigation. Using a simplified structural model, we perform equilibrium simulations of a single nucleosome at various temperatures. We further reconstruct all-atom nucleosome structures from simulation trajectories. We find that histone tails bind to nucleosomal DNA via strong salt-bridge interactions over a wide range of temperatures, suggesting a mechanism of chromatin structural organization whereby histone tails regulate inter- and intranucleosomal assemblies via binding with nucleosomal DNA. We identify specific regions of the histone core H2A/H2B-H4/H3-H3/H4-H2B/H2A, termed “cold sites”, which retain a significant fraction of contacts with adjoining residues throughout the simulation, indicating their functional role in nucleosome organization. Cold sites are clustered around H3-H3, H2A-H4 and H4-H2A interhistone interfaces, indicating the necessity of these contacts for nucleosome stability. Essential dynamics analysis of simulation trajectories shows that bending across the H3-H3 is a prominent mode of intranucleosomal dynamics. We postulate that effects of salts on mononucleosomes can be modeled in discrete molecular dynamics by modulating histone-DNA interaction potentials. Local fluctuations in nucleosomal DNA vary significantly along the DNA sequence, suggesting that only a fraction of histone-DNA contacts make strong interactions dominating mononucleosomal dynamics. Our findings suggest that histone tails have a direct functional role in stabilizing higher-order chromatin structure, mediated by salt-bridge interactions with adjacent DNA.     

Structural Basis for Influence of Viral Glycans on Ligand Binding by Influenza Hemagglutinin

Binding of cell surface glycans by influenza hemagglutinin controls viral attachment and infection of host cells. This binding is a three-way interaction between viral proteins, host glycans, and viral glycans; many structural details of this interaction have been difficult to resolve. Here, we use a series of 100-ns molecular dynamics simulations to further analyze available crystallographic data on hemagglutinin-ligand interactions. Based on our simulations, we predict that the viral glycans contact the host glycans within 1-2 residues of the ligand-binding site. We also predict that the glycan-glycan interactions contain both stabilizing and destabilizing components. These predictions suggest a structural means to explain why changes to viral glycosylation alter the efficiency and selectivity of ligand binding. We also predict that the proximity of these interactions to the ligand-binding pocket will impact the binding affinity of small glycomimetic ligands analogous to the influenza neuraminidase inhibitors currently in clinical use.     

Crystal structures of bacterial peptidoglycan amidase AmpD and an unprecedented activation mechanism

AmpD is a cytoplasmic peptidoglycan (PG) amidase involved in bacterial cell-wall recycling and in induction of β-lactamase, a key enzyme of β-lactam antibiotic resistance. AmpD belongs to the amidase_2 family that includes zinc-dependent amidases and the peptidoglycan-recognition proteins (PGRPs), highly conserved pattern-recognition molecules of the immune system. Crystal structures of Citrobacter freundii AmpD were solved in this study for the apoenzyme, for the holoenzyme at two different pH values, and for the complex with the reaction products, providing insights into the PG recognition and the catalytic process. These structures are significantly different compared with the previously reported NMR structure for the same protein. The NMR structure does not possess an accessible active site and shows the protein in what is proposed herein as an inactive “closed” conformation. The transition of the protein from this inactive conformation to the active “open” conformation, as seen in the x-ray structures, was studied by targeted molecular dynamics simulations, which revealed large conformational rearrangements (as much as 17 Å) in four specific regions representing one-third of the entire protein. It is proposed that the large conformational change that would take the inactive NMR structure to the active x-ray structure represents an unprecedented mechanism for activation of AmpD. Analysis is presented to argue that this activation mechanism might be representative of a regulatory process for other intracellular members of the bacterial amidase_2 family of enzymes.     

Recognition of Histone H3K4 Trimethylation by the Plant Homeodomain of PHF2 Modulates Histone Demethylation

Distinct lysine methylation marks on histones create dynamic signatures deciphered by the “effector” modules, although the underlying mechanisms remain unclear. We identified the plant homeodomain- and Jumonji C domain-containing protein PHF2 as a novel histone H3K9 demethylase. We show in biochemical and crystallographic analyses that PHF2 recognizes histone H3K4 trimethylation through its plant homeodomain finger and that this interaction is essential for PHF2 occupancy and H3K9 demethylation at rDNA promoters. Our study provides molecular insights into the mechanism by which distinct effector domains within a protein cooperatively modulate the “cross-talk” of histone modifications.     

Integrin Activation Dynamics between the RGD-binding Site and the Headpiece Hinge

Integrins form mechanical links between the extracellular matrix and the cytoskeleton. Although integrin activation is known to be regulated by an allosteric conformational change, which can be induced from the extracellular or intracellular end of the molecule, little is known regarding the sequence of structural events by which signals propagate between distant sites. Here, we reveal with molecular dynamics simulations of the FnIII₁₀-bound α_Vβ₃ integrin headpiece how the binding pocket and interdomain βA/hybrid domain hinge on the distal end of the βA domain are allosterically linked via a hydrophobic T-junction between the middle of the α1 helix and top of the α7 helix. The key results of this study are: 1) that this T-junction is induced by ligand binding and hinge opening, and thus displays bidirectionality; 2) that formation of this junction can be accelerated by ligand-mediated force; and 3) how formation of this junction is inhibited by Ca²⁺ in place of Mg²⁺ at the site adjacent to the metal ion-dependent adhesion site (“ADMIDAS”). Together with recent experimental evidence that integrin complexes can form catch bonds (i.e. become strengthened under force), as well as earlier evidence that Ca²⁺ at the ADMIDAS results in lower binding affinity, these simulations provide a common structural model for the dynamic process by which integrins become activated.     

Molecular Probing of the HPV-16 E6 Protein Alpha Helix Binding Groove with Small Molecule Inhibitors

The human papillomavirus (HPV) HPV E6 protein has emerged as a central oncoprotein in HPV-associated cancers in which sustained expression is required for tumor progression. A majority of the E6 protein interactions within the human proteome use an alpha-helix groove interface for binding. The UBE3A/E6AP HECT domain ubiquitin ligase binds E6 at this helix-groove interface. This enables formation of a trimeric complex with p53, resulting in destruction of this tumor suppressor. While recent x-ray crystal structures are useful, examples of small molecule probes that can modulate protein interactions at this interface are limited. To develop insights useful for potential structure-based design of ligands for HPV E6, a series of 2,6-disubstituted benzopyranones were prepared and tested as competitive antagonists of E6-E6AP helix-groove interactions. These small molecule probes were used in both binding and functional assays to evaluate recognition features of the E6 protein. Evidence for an ionic functional group interaction within the helix groove was implicated by the structure-activity among the highest affinity ligands. The molecular topographies of these protein-ligand interactions were evaluated by comparing the binding and activities of single amino acid E6 mutants with the results of molecular dynamic simulations. A group of arginine residues that form a rim-cap over the E6 helix groove offer compensatory roles in binding and recognition of the small molecule probes. The flexibility and impact on the overall helix-groove shape dictated by these residues offer new insights for structure-based targeting of HPV E6.     

Folded conformation of an immunostimulating tetrapeptide rigin: high temperature molecular dynamics simulation study

Employing high temperature quenched molecular dynamics (QMD) simulations the conformational energy space of an immunostimulating tetrapeptide rigin: H-Gly³⁴¹-Gln-Pro-Arg³⁴⁴-OH, is explored. Using distance dependent dielectric (=r_ij) 31 different low energy starting structures with identical sequence were computed for their conformational preferences. According to the hypothesis of O'Connors et al. [J. Med. Chem. 35 (1992), 2870], 83 low-energy conformers resulted from unrestrained molecular dynamics (MD) simulations, could be classified into two energy minimized families: A and B, comprised of 64 (Pro C^γ-endo orientation) and 19 (Pro C^γ-exo orientation) structures, respectively. An examination of these families revealed the existence of a remarkably similar folded backbone conformation: torsion angles being φ_i+1 ≈−65°, ψ_i+1 ≈−65°, φ_i+2 ≈−65°, ψ_i+2 ≈−60°, characterizing a distorted type III β-turn structure across the central Gln-Pro segment. The folded conformation of rigin is devoid of a classical 1 ← 4 intra-molecular hydrogen bond nevertheless, the conformation is stabilized by an effective ‘salt-bridge’, i.e., Gly H₃N⁺… COO⁻ Arg interaction. Surprisingly, in both the families the unusual folded side-chain dispositions of the Gln residue favor the formation of a unique intra-residue ‘main-chain to side-chain’ H-bond, i.e., N–H…N^ε interaction, encompassing a seven-membered ring motif. The conformational attributes may be valuable in de novo construction of structure-based drug candidates having sufficient stimulating activity.     

Deciphering the mechanisms of HPV E6 mutations in the destabilization of E6/E6AP/p53 complex

In epithelial tumors, oncoprotein E6 binds with the ubiquitin ligase E6AP to form E6/E6AP heterodimer; then this heterodimer recruits p53 to form E6/E6AP/p53 heterotrimer and induces p53 degradation. Recent experiments demonstrated that three E6 single-site mutants (F47R, R102A, and L50E) can inhibit the E6/E6AP/p53 heterotrimer formation and rescue p53 from the degradation pathway. However, the molecular mechanism underlying mutation-induced heterotrimer inhibition remains largely elusive. Herein, we performed extensive molecular dynamics simulations (totally ～13 μs) on both heterodimer and heterotrimer to elucidate at an atomic level how each p53-degradation-defective HPV16 E6 mutant reduces the structural stabilities of the two complexes. Our simulations reveal that the three E6 mutations destabilize the structure of E6/E6AP/p53 complex through distinct mechanisms. Although F47R^E6 mutation has no effect on the structure of E6/E6AP heterodimer, it results in an electrostatic repulsion between R47^E6 and R290^p53, which is unfavorable for E6-p53 binding. R102A^E6 mutation destabilizes the structure of E6/E6AP heterodimer and significantly disrupts hydrophobic and cation-π interactions between F47^E6 and E286^p53/L298^p53/R290^p53. L50E^E6 mutation impairs both E6 interdomain interactions (especially F47-K108 cation-π interaction) and E6-E6AP intermolecular interactions important for the stabilization of E6/E6AP heterodimer. This study identifies the intra- and intermolecular interactions crucial for the complex stability, which may provide mechanistic insights into the inhibition of complex formation by the three HPV16 E6 mutations.     

Global self-organization of the cellular metabolic structure

Background

Over many years, it has been assumed that enzymes work either in an isolated way, or organized in small catalytic groups. Several studies performed using “metabolic networks models” are helping to understand the degree of functional complexity that characterizes enzymatic dynamic systems. In a previous work, we used “dissipative metabolic networks” (DMNs) to show that enzymes can present a self-organized global functional structure, in which several sets of enzymes are always in an active state, whereas the rest of molecular catalytic sets exhibit dynamics of on-off changing states. We suggested that this kind of global metabolic dynamics might be a genuine and universal functional configuration of the cellular metabolic structure, common to all living cells. Later, a different group has shown experimentally that this kind of functional structure does, indeed, exist in several microorganisms.

Methodology/Principal Findings

Here we have analyzed around 2.500.000 different DMNs in order to investigate the underlying mechanism of this dynamic global configuration. The numerical analyses that we have performed show that this global configuration is an emergent property inherent to the cellular metabolic dynamics. Concretely, we have found that the existence of a high number of enzymatic subsystems belonging to the DMNs is the fundamental element for the spontaneous emergence of a functional reactive structure characterized by a metabolic core formed by several sets of enzymes always in an active state.

Conclusions/Significance

This self-organized dynamic structure seems to be an intrinsic characteristic of metabolism, common to all living cellular organisms. To better understand cellular functionality, it will be crucial to structurally characterize these enzymatic self-organized global structures.     

Molecular and Mechanical Causes of Microtubule Catastrophe and Aging

Tubulin polymers, microtubules, can switch abruptly from the assembly to shortening. These infrequent transitions, termed “catastrophes”, affect numerous cellular processes but the underlying mechanisms are elusive. We approached this complex stochastic system using advanced coarse-grained molecular dynamics modeling of tubulin-tubulin interactions. Unlike in previous simplified models of dynamic microtubules, the catastrophes in this model arise owing to fluctuations in the composition and conformation of a growing microtubule tip, most notably in the number of protofilament curls. In our model, dynamic evolution of the stochastic microtubule tip configurations over a long timescale, known as the system’s “aging”, gives rise to the nonexponential distribution of microtubule lifetimes, consistent with experiment. We show that aging takes place in the absence of visible changes in the microtubule wall or tip, as this complex molecular-mechanical system evolves slowly and asymptotically toward the steady-state level of the catastrophe-promoting configurations. This new, to our knowledge, theoretical basis will assist detailed mechanistic investigations of the mechanisms of action of different microtubule-binding proteins and drugs, thereby enabling accurate control over the microtubule dynamics to treat various pathologies.