Human cytochrome C natural variants: Studying the membrane binding properties of G41S and Y48H by fluorescence energy transfer and molecular dynamics.

Cytochrome C (cyt C), the protein involved in oxidative phosphorylation, plays several other crucial roles necessary for both cell life and death. Studying natural variants of cyt C offers the possibility to better characterize the structure-to-function relationship that modulates the different activities of this protein. Naturally mutations in human cyt C (G41S and Y48H) occur in the protein central Ω-loop and cause thrombocytopenia 4. In this study, we have investigated the binding of such variants and of wild type (wt) cyt C to synthetic cardiolipin-containing vesicles. The mutants have a lower propensity in membrane binding, displaying higher dissociation constants with respect to the wt protein. Compressibility measurements reveal that both variants are more flexible than the wt, suggesting that the native central Ω-loop is important for the interaction with membranes. Such hypothesis is supported by molecular dynamics simulations. A minimal distance analysis indicates that in the presence of cardiolipin the central Ω-loop of the mutants is no more in contact with the membrane, as it happens instead in the case of wt cyt C. Such finding might provide a hint for the reduced membrane binding capacity of the variants and their enhanced peroxidase activity in vivo.

Investigation of the structure of protein-polymer conjugates in solution reveals the impact of protein deuteration and the size of the polymer on its thermal stability.

The conjugation of proteins with polymers offers immense biotechnological potential by creating novel macromolecules. This article presents experimental findings on the structural properties of maltose-binding protein (MBP) conjugated with linear biodegradable polyphosphoester polymers with different molecular weights. We studied isotopic effects on both proteins and polymers. Circular dichroism and fluorescence spectroscopy and small-angle neutron scattering reveal that the conjugation process destabilizes the protein, affecting the secondary more than the tertiary structure, even at room temperature, and that the presence of two domains in the MBP may contribute to its observed instability. Notably, unfolding temperatures differ between native MBP and the conjugates. In particular, this study sheds light on the complex interplay of factors such as the deuteration influencing protein stability and conformational changes in the conjugation processes. The perdeuteration influences the hydrogen bond network and hydrophobic interactions in the case of the MBP protein. The perdeuteration of the protein influences the hydrogen bond network and hydrophobic interactions. This is evident in the decreased thermal stability of deuterated MBP protein, in the conjugate, especially with high-molecular-mass polymers.

Conformational Dynamics of Lipoxygenases and Their Interaction with Biological Membranes.

Lipoxygenases (LOXs) are a family of enzymes that includes different fatty acid oxygenases with a common tridimensional structure. The main functions of LOXs are the production of signaling compounds and the structural modifications of biological membranes. These features of LOXs, their widespread presence in all living organisms, and their involvement in human diseases have attracted the attention of the scientific community over the last decades, leading to several studies mainly focused on understanding their catalytic mechanism and designing effective inhibitors. The aim of this review is to discuss the state-of-the-art of a different, much less explored aspect of LOXs, that is, their interaction with lipid bilayers. To this end, the general architecture of six relevant LOXs (namely human 5-, 12-, and 15-LOX, rabbit 12/15-LOX, coral 8-LOX, and soybean 15-LOX), with different specificity towards the fatty acid substrates, is analyzed through the available crystallographic models. Then, their putative interface with a model membrane is examined in the frame of the conformational flexibility of LOXs, that is due to their peculiar tertiary structure. Finally, the possible future developments that emerge from the available data are discussed.

Head or tail? A molecular dynamics approach to the complex structure of TNF-associated factor TRAF2.

Tumor necrosis factor receptor-associated factor proteins (TRAFs) are trimeric proteins that play a fundamental role in signaling, acting as intermediaries between the tumor necrosis factor (TNF) receptors and the proteins that transmit the downstream signal. The monomeric subunits of all the TRAF family members share a common tridimensional structure: a C-terminal globular domain and a long coiled-coil tail characterizing the N-terminal section. In this study, the dependence of the TRAF2 dynamics on the length of its tail was analyzed in silico. In particular, we used the available crystallographic structure of a C-terminal fragment of TRAF2 (168 out of 501 a.a.), TRAF2-C, and that of a longer construct, addressed as TRAF2-plus, that we have re-constructed using the AlphaFold2 code. The results indicate that the longer N-terminal tail of TRAF2-plus has a strong influence on the dynamics of the globular regions in the protein C-terminal head. In fact, the quaternary interactions among the TRAF2-C subunits change asymmetrically in time, while the movements of TRAF2-plus monomers are rather limited and more ordered than those of the shorter construct. Such findings shed a new light on the dynamics of TRAF subunits and on the protein mechanism in vivo, since TRAF monomer-trimer equilibrium is crucial for several reasons (receptor recognition, membrane binding, hetero-oligomerization).

Revisited role of TRAF2 and TRAF2 C-terminal domain in endoplasmic reticulum stress-induced autophagy in HAP1 leukemia cells.

The scaffold protein Tumor Necrosis Factor Receptor-Associated Factor 2 (TRAF2) has been reported to play a key role in the endoplasmic reticulum (ER) stress-induced activation of c-Jun N-terminal Kinase (JNK) and hence autophagy. Autophagy is a highly conserved catabolic process, whose dysregulation is involved in the pathogenesis of various human diseases, including cancer. We investigated the involvement of TRAF2 in autophagy regulation in the human leukemic HAP1 cell line, under both basal and ER stress conditions. In TRAF2-knockout HAP1 cell line (KO), the basal autophagic flux was higher than in the parental cell line (WT). Moreover, tunicamycin-induced ER stress stimulated JNK activation and autophagy both in WT and KO HAP1. On the other hand, re-expression of a TRAF2 C-terminal fragment (residues ,310-501), in a TRAF2-KO cellular background, rendered HAP1 cells unable to activate both JNK and autophagy upon ER stress induction. Of note, this apparent dominant negative effect of the C-terminal fragment was observed even in the absence of the endogenous, full-length TRAF2 molecule. Furthermore, the expression of the C-terminal fragment resulted in both protein kinase B (AKT) pathway activation and increased resistance to the toxic effects induced by prolonged ER stress conditions. These findings indicate that TRAF2 is dispensable for the activation of both JNK and autophagy in HAP1 cells, while the TRAF2 C-terminal domain may play an autonomous role in regulating the cellular response to ER stress.

The Odd Faces of Oligomers: The Case of TRAF2-C, A Trimeric C-Terminal Domain of TNF Receptor-Associated Factor.

TNF Receptor Associated Factor 2 (TRAF2) is a trimeric protein that belongs to the TNF receptor associated factor family (TRAFs). The TRAF2 oligomeric state is crucial for receptor binding and for its interaction with other proteins involved in the TNFR signaling. The monomer-trimer equilibrium of a C- terminal domain truncated form of TRAF2 (TRAF2-C), plays also a relevant role in binding the membrane, causing inward vesiculation. In this study, we have investigated the conformational dynamics of TRAF2-C through circular dichroism, fluorescence, and dynamic light scattering, performing temperature-dependent measurements. The data indicate that the protein retains its oligomeric state and most of its secondary structure, while displaying a significative increase in the heterogeneity of the tyrosines signal, increasing the temperature from ≈15 to ≈35 °C. The peculiar crowding of tyrosine residues (12 out of 18) at the three subunit interfaces and the strong dependence on the trimer concentration indicate that such conformational changes mainly involve the contact areas between each pair of monomers, affecting the oligomeric state. Molecular dynamic simulations in this temperature range suggest that the interfaces heterogeneity is an intrinsic property of the trimer that arises from the continuous, asymmetric approaching and distancing of its subunits. Such dynamics affect the results of molecular docking on the external protein surface using receptor peptides, indicating that the TRAF2-receptor interaction in the solution might not involve three subunits at the same time, as suggested by the static analysis obtainable from the crystal structure. These findings shed new light on the role that the TRAF2 oligomeric state might have in regulating the protein binding activity in vivo.

Polymorphism on human aromatase affects protein dynamics and substrate binding: spectroscopic evidence.

Human aromatase is a member of the cytochrome P450 superfamily, involved in steroid hormones biosynthesis. In particular, it converts androgen into estrogens being therefore responsible for the correct sex steroids balance. Due to its capacity in producing estrogens it has also been considered as a promising target for breast cancer therapy. Two single-nucleotide polymorphisms (R264C and R264H) have been shown to alter aromatase activity and they have been associated to an increased or decreased risk for estrogen-dependent pathologies. Here, the effect of these mutations on the protein dynamics is investigated by UV/FTIR and time resolved fluorescence spectroscopy. H/D exchange rates were measured by FTIR for the three proteins in the ligand-free, substrate- and inhibitor-bound forms and the data indicate that the wild-type enzyme undergoes a conformational change leading to a more compact tertiary structure upon substrate or inhibitor binding. Indeed, the H/D exchange rates are decreased when a ligand is present. In the variants, the exchange rates in the ligand-free and -bound forms are similar, indicating that a structural change is lacking, despite the single amino acid substitution is located in the peripheral shell of the protein molecule. Moreover, the fluorescence lifetimes data show that the quenching effect on tryptophan-224 observed upon ligand binding in the wild-type, is absent in both variants. Since this residue is located in the catalytic pocket, these findings suggest that substrate entrance and/or retention in the active site is partially compromised in both mutants. A contact network analysis demonstrates that the protein structure is organized in two main clusters, whose connectivity is altered by ligand binding, especially in correspondence of helix-G, where the amino acid substitutions occur. Our findings demonstrate that SNPs resulting in mutations on aromatase surface modify the protein flexibility that is required for substrate binding and catalysis. The cluster analysis provides a rationale for such effect, suggesting helix G as a possible target for aromatase inhibition.

Conformational Heterogeneity and Cooperative Effects of Mammalian ALOX15.

Arachidonic acid lipoxygenases (ALOXs) have been suggested to function as monomeric enzymes, but more recent data on rabbit ALOX15 indicated that there is a dynamic monomer-dimer equilibrium in aqueous solution. In the presence of an active site ligand (the ALOX15 inhibitor RS7) rabbit ALOX15 was crystalized as heterodimer and the X-ray coordinates of the two monomers within the dimer exhibit subtle structural differences. Using native polyacrylamide electrophoresis, we here observed that highly purified and predominantly monomeric rabbit ALOX15 and human ALOX15B are present in two conformers with distinct electrophoretic mobilities. In silico docking studies, molecular dynamics simulations, site directed mutagenesis experiments and kinetic measurements suggested that in aqueous solutions the two enzymes exhibit motional flexibility, which may impact the enzymatic properties.

The Puzzling Problem of Cardiolipin Membrane-Cytochrome c Interactions: A Combined Infrared and Fluorescence Study.

The interaction of cytochrome c (cyt c) with natural and synthetic membranes is known to be a complex phenomenon, involving both protein and lipid conformational changes. In this paper, we combined infrared and fluorescence spectroscopy to study the structural transformation occurring to the lipid network of cardiolipin-containing large unilamellar vesicles (LUVs). The data, collected at increasing protein/lipid ratio, demonstrate the existence of a multi-phase process, which is characterized by: (i) the interaction of cyt c with the lipid polar heads; (ii) the lipid anchorage of the protein on the membrane surface; and (iii) a long-distance order/disorder transition of the cardiolipin acyl chains. Such effects have been quantitatively interpreted introducing specific order parameters and discussed in the frame of the models on cyt c activity reported in literature.

Non-symmetrical structural behavior of a symmetric protein: the case of homo-trimeric TRAF2 (tumor necrosis factor-receptor associated factor 2).

The oligomeric state of TRAF2 (tumor necrosis factor-receptor associated factor 2), a TNF (tumor necrosis factor) receptor-associated factor, is crucial for membrane binding and probably plays a fundamental role in regulating the protein function in vivo. In this study we have combined molecular dynamics with the protein contact network approach to characterize the interaction of the three identical subunits of TRAF2. The average structure obtained after a 225 ns simulation reveals that two clusters of different size are formed, one of which includes almost completely two subunits, while the third monomer appears to be more independent. This picture is also confirmed by the estimated average number of inter-subunit contacts and by the comparison of side chains mobility in each monomer. The analysis of equilibrium pressure-induced dissociation measurements supports such findings, indicating that the dimeric-monomeric (2 + 1) might be prevalent with respect to the trimeric configuration, especially in the case of more diluted samples. These findings suggest that the formation of monomeric species, which is crucial for the formation of intra-luminal vesicles, might depend on preferential asymmetric interactions among the three subunits.Communicated by Ramaswamy H. Sarma.

Disclosing Allostery Through Protein Contact Networks.

Proteins are located in the twilight zone between chemistry and biology, where a peculiar kind of complexity starts. Proteins are the smallest 'devices' showing a sensible adaptation to their environment by the production of appropriate behavior when facing a specific stimulus. This fact qualifies (from the 'effector' side) proteins as nanomachines working as catalysts, motors, or switches. However (from the sensor side), the need to single out the 'specific stimulus' out of thermal noise qualifies proteins as information processing devices. Allostery corresponds to the modification of the configuration (in a broad sense) of the protein molecule in response to a specific stimulus in a non-strictly local way, thereby connecting the sensor and effector sides of the nanomachine. This is why the 'disclosing' of allostery phenomenon is at the very heart of protein function; in this chapter, we will demonstrate how a network-based representation of protein structure in terms of nodes (aminoacid residues) and edges (effective contacts between residues) is the natural language for getting rid of allosteric phenomena and, more in general, of protein structure/function relationships.

Protein Assembly: Defining the Strength of Protein-Protein Interactions Coupling the Theory with Experiments.

In this paper we report a procedure to analyze protein homodimer interfaces.We approached the problem by means of a topological methodology. In particular, we analyzed the subunits interface of about 50 homodimers and we have defined a few parameters that allow to organize these proteins in six different classes. The main characteristics of each class of homodimers have been discussed also taking into account their stabilization energy, as reported in the literature from the experimental measurements. A paradigmatic example for each class has been reported and a graphical representation proposed in order to better explain the meaning of the parameters chosen.

A role of Gln596 in fine-tuning mammalian ALOX15 specificity, protein stability and allosteric properties.

His596 of human ALOX12 has been suggested to interact with the COO--group of arachidonic acid during ALOX catalysis. In mammalian ALOX15 orthologs Gln596 occupies this position and this amino acid exchange might contribute to the functional differences between the two ALOX-isoforms. To explore the role of Gln596 for ALOX15 functionality we mutated this amino acid to different residues in rabbit and human ALOX15 and investigated the impact of these mutations on structural, catalytic and allosteric enzyme properties. To shed light on the molecular basis of the observed functional alterations we performed in silico substrate docking studies and molecular dynamics simulations and also explored the impact of Gln596 exchange on the protein structure. The combined theoretical and experimental data suggest that Gln596 may not directly interact with the COO--group of arachidonic acid. In contrast, mutations at Gln596 destabilize the secondary and tertiary structure of ALOX15 orthologs, which may be related to a disturbance of the electrostatic interaction network with other amino acids in the immediate surrounding. Moreover, our MD-simulations suggest that the geometry of the dimer interface depends on the structure of substrate bound inside the substrate-binding pocket and that Gln596Ala exchange impairs the allosteric properties of the enzyme. Taken together, these data indicate the structural and functional importance of Gln596 for ALOX15 catalysis.

Studying the TRAF2 binding to model membranes: The role of subunits dissociation.

The ability of a C-terminal truncated form of TRAF2 to bind synthetic vesicles has been quantitatively studied by steady-state fluorescence energy transfer from the protein to large unilamellar vesicles (LUVs) prepared with different lipid mixtures. The dissociation constants, the free energy of binding, and the average number of phospholipids interacting with truncated TRAF2 have been evaluated from the corresponding binding curves. The results indicate that the protein strongly interacts with the lipid bilayer, preferentially in the monomeric state. These findings have been discussed in terms of their possible role in the activity of TRAF2 in vivo.

New insight into the interaction of TRAF2 C-terminal domain with lipid raft microdomains.

In this study we provide the first evidence of the interaction of a truncated-TRAF2 with lipid raft microdomains. We have analyzed this interaction by measuring the diffusion coefficient of the protein in large and giant unilamellar vesicles (LUVs and GUVs, respectively) obtained both from synthetic lipid mixtures and from natural extracts. Steady-state fluorescence measurements performed with synthetic vesicles indicate that this truncated form of TRAF2 displays a tighter binding to raft-like LUVs with respect to the control (POPC-containing LUVs), and that this process depends on the protein oligomeric state. Generalized Polarization measurements and spectral phasor analysis revealed that truncated-TRAF2 affects the membrane fluidity, especially when vesicles are heated up at physiological temperature. The addition of nanomolar concentration of TRAF2 in GUVs also seems to exert a mechanical action, as demonstrated by the formation of intraluminal vesicles, a process in which ganglioside GM1 plays a crucial role.

Exploring the stability of dimers through protein structure topology.

Protein homodimers pose some intriguing questions about the relation between structure and stability. We approached the problem by means of a topological methodology based on protein contact networks. We correlated local interface descriptors with structure and energy global properties of the systems under analysis. We demonstrated that the graph energy, formerly applied to the analysis of unconjugated hydrocarbons structures, is the bridge between the topological and energetic description of protein complexes. This is a first step for the generation of a "protein structural formula", analogous to the molecular graphs in organic chemistry.

TNFR-Associated Factor-2 (TRAF2): Not Only a Trimer.

TNF receptor-associated factors (TRAFs) are characterized by an oligomeric structure that plays a fundamental role in the binding process with membrane receptors. In this work, we studied the trimer-to-monomer (T ↔ 3M) equilibrium transition of the TRAF2 C-terminal domain using both chemical (dilution/guanidinium hydrochloride) and mechanical stress (high pressure) to induce the dissociation of the native protein into subunits. The experimental results and computer simulations indicate that stable monomers exist and that their population accounts for 15% of the total TRAF2 molecules already at a physiological intracellular concentration (≈1 µM), being instead the predominant species in the nanomolar concentration range. Because the total amount of TRAF2 changes during a cell cycle, the monomer-trimer equilibrium can be crucial for regulating the activities of TRAF2 in vivo.

Liposome-Templated Hydrogel Nanoparticles as Vehicles for Enzyme-Based Therapies.

Several diseases are related to the lack or to the defective activity of a particular enzyme; therefore, these proteins potentially represent a very interesting class of therapeutics. However, their application is hampered by their rapid degradation and immunogenic side effects. Most attempts to increase the bioavailability of therapeutic enzymes are based on formulations in which the protein is entrapped within a scaffold structure but needs to be released to exert its activity. In this work, an alternative method will be described, designed to keep the enzyme in its active form inside a nanoparticle (NP) without the need to release it, thus maintaining the protective action of the nanoscaffold during the entire period of administration. In this approach, liposomes were used as nanotemplates for the synthesis of polyacrylamide hydrogel NPs under nondenaturing conditions, optimizing the polymer properties to obtain a mesh size small enough to limit the enzyme release while allowing the free diffusion of its substrates and products. The enzyme Cu, Zn-superoxide dismutase was chosen as a test case for this study, but our results indicate that the approach is generalizable to other enzymes. Biocompatible, size-tunable nanoparticles have been obtained, with a good encapsulation efficiency (37%), in which the enzyme maintains its activity. This system represents a promising tool for enzyme-based therapy, which would protect the protein from antibodies and degradation while allowing it to exert its catalytic activity.

Highly specific ubiquitin-competing molecules effectively promote frataxin accumulation and partially rescue the aconitase defect in Friedreich ataxia cells.

Friedreich ataxia is an inherited neurodegenerative disease that leads to progressive disability. There is currently no effective treatment and patients die prematurely. The underlying genetic defect leads to reduced expression of the mitochondrial protein frataxin. Frataxin insufficiency causes mitochondrial dysfunction and ultimately cell death, particularly in peripheral sensory ganglia. There is an inverse correlation between the amount of residual frataxin and the severity of disease progression; therefore, therapeutic approaches aiming at increasing frataxin levels are expected to improve patients' conditions. We previously discovered that a significant amount of frataxin precursor is degraded by the ubiquitin/proteasome system before its functional mitochondrial maturation. We also provided evidence for the therapeutic potential of small molecules that increase frataxin levels by docking on the frataxin ubiquitination site, thus preventing frataxin ubiquitination and degradation. We called these compounds ubiquitin-competing molecules (UCM). By extending our search for effective UCM, we identified a set of new and more potent compounds that more efficiently promote frataxin accumulation. Here we show that these compounds directly interact with frataxin and prevent its ubiquitination. Interestingly, these UCM are not effective on the ubiquitin-resistant frataxin mutant, indicating their specific action on preventing frataxin ubiquitination. Most importantly, these compounds are able to promote frataxin accumulation and aconitase rescue in cells derived from patients, strongly supporting their therapeutic potential.

Modules identification in protein structures: the topological and geometrical solutions.

The identification of modules in protein structures has major relevance in structural biology, with consequences in protein stability and functional classification, adding new perspectives in drug design. In this work, we present the comparison between a topological (spectral clustering) and a geometrical (k-means) approach to module identification, in the frame of a multiscale analysis of the protein architecture principles. The global consistency of an adjacency matrix based technique (spectral clustering) and a method based on full rank geometrical information (k-means) give a proof-of-concept of the relevance of protein contact networks in structure determination. The peculiar "small-world" character of protein contact graphs is established as well, pointing to average shortest path as a mesoscopic crucial variable to maximize the efficiency of within-molecule signal transmission. The specific nature of protein architecture indicates topological approach as the most proper one to highlight protein functional domains, and two new representations linking sequence and topological role of aminoacids are demonstrated to be of use for protein structural analysis. Here we present a case study regarding azurin, a small copper protein implied in the Pseudomonas aeruginosa respiratory chain. Its pocket molecular shape and its electron transfer function have challenged the method, highlighting its potentiality to catch jointly the structure and function features of protein structures through their decomposition into modules.