Human Frataxin Folds Via an Intermediate State. Role of the C-Terminal Region

Fifty years of biophysics in Argentina

In 1972, a group of young Argentinean scientists nucleated in the so-called Membrane Club constituted the Biophysical Society of Argentina (SAB). Over the years, this Society has grown and embraced new areas of research and emerging technologies. In this commentary, we provide an overview of the early stages of biophysics development in Argentina and highlight some of the notable achievements made during the past five decades. The SAB Annual Meetings have been a platform for intense scientific discussions, and the Society has fostered numerous international connections, becoming a hallmark of SAB activities over these 50 years. Initially centered on membrane biophysics, SAB focus has since expanded to encompass diverse fields such as molecular, cellular, and systems biophysics.

Selection of synthetic proteins to modulate the human frataxin function

Frataxin is a kinetic activator of the mitochondrial supercomplex for iron-sulfur cluster assembly. Low frataxin expression or a decrease in its functionality results in Friedreich's Ataxia (FRDA). With the aim of creating new molecular tools to study this metabolic pathway, and ultimately, to explore new therapeutic strategies, we have investigated the possibility of obtaining small proteins exhibiting a high affinity for frataxin. In this study, we applied the ribosome display approach, using human frataxin as the target. We focused on Affi_224, one of the proteins that we were able to select after five rounds of selection. We have studied the interaction between both proteins and discussed some applications of this specific molecular tutor, concerning the modulation of the supercomplex activity. Affi_224 and frataxin showed a K_D value in the nanomolar range, as judged by surface plasmon resonance analysis. Most likely, it binds to the frataxin acidic ridge, as suggested by the analysis of chemical shift perturbations (nuclear magnetic resonance) and computational simulations. Affi_224 was able to increase Cys NFS1 desulfurase activation exerted by the FRDA frataxin variant G130V. Importantly, Affi_224 interacts with frataxin in a human cellular model. Our results suggest quaternary addition may be a new tool to modulate frataxin function in vivo. Nevertheless, more functional experiments under physiological conditions should be carried out to evaluate Affi_224 effectiveness in FRDA cell models.

Natural flavonoids effectively block the CD81 receptor of hepatocytes and inhibit HCV infection: a computational drug development approach

Hepatitis C virus (HCV) infection is a major public health concern, and almost two million people are infected per year globally. This is occurred by the diverse spectrum of viral genotypes, which are directly associated with chronic liver disease (fibrosis, and cirrhosis). Indeed, the viral genome encodes three principal proteins as sequentially core, E1, and E2. Both E1 and E2 proteins play a crucial role in the attachment of the host system, but E2 plays a more fundamental role in attachment. The researchers have found the "E2-CD81 complex" at the entry site, and therefore, CD81 is the key receptor for HCV entrance in both humans, and chimpanzees. So, the researchers are trying to block the host CD81 receptor and halt the virus entry within the cellular system via plant-derived compounds. Perhaps that is why the current research protocol is designed to perform an in silico analysis of the flavonoid compounds for targeting the tetraspanin CD81 receptor of hepatocytes. To find out the best flavonoid compounds from our library, web-based tools (Swiss ADME, pKCSM), as well as computerized tools like the PyRx, PyMOL, BIOVIA Discovery Studio Visualizer, Ligplot+ V2.2, and YASARA were employed. For molecular docking studies, the flavonoid compounds docked with the targeted CD81 protein, and herein, the best-outperformed compounds are Taxifolin, Myricetin, Puerarin, Quercetin, and (-)-Epicatechin, and outstanding binding affinities are sequentially - 7.5, - 7.9, - 8.2, - 8.4, and - 8.5 kcal/mol, respectively. These compounds have possessed more interactions with the targeted protein. To validate the post docking data, we analyzed both 100 ns molecular dynamic simulation, and MM-PBSA via the YASARA simulator, and finally finds the more significant outcomes. It is concluded that in the future, these compounds may become one of the most important alternative antiviral agents in the fight against HCV infection. It is suggested that further in vivo, and in vitro research studies should be done to support the conclusions of this in silico research workflow.

Modelling Protein Plasticity: The Example of Frataxin and Its Variants

Frataxin (FXN) is a protein involved in storage and delivery of iron in the mitochondria. Single-point mutations in the FXN gene lead to reduced production of functional frataxin, with the consequent dyshomeostasis of iron. FXN variants are at the basis of neurological impairment (the Friedreich's ataxia) and several types of cancer. By using altruistic metadynamics in conjunction with the maximal constrained entropy principle, we estimate the change of free energy in the protein unfolding of frataxin and of some of its pathological mutants. The sampled configurations highlight differences between the wild-type and mutated sequences in the stability of the folded state. In partial agreement with thermodynamic experiments, where most of the analyzed variants are characterized by lower thermal stability compared to wild type, the D104G variant is found with a stability comparable to the wild-type sequence and a lower water-accessible surface area. These observations, obtained with the new approach we propose in our work, point to a functional switch, affected by single-point mutations, of frataxin from iron storage to iron release. The method is suitable to investigate wide structural changes in proteins in general, after a proper tuning of the chosen collective variable used to perform the transition.

Protein Mutations and Stability, a Link with Disease: The Case Study of Frataxin

Protein mutations may lead to pathologies by causing protein misfunction or propensity to degradation. For this reason, several studies have been performed over the years to determine the capability of proteins to retain their native conformation under stress condition as well as factors to explain protein stabilization and the mechanisms behind unfolding. In this review, we explore the paradigmatic example of frataxin, an iron binding protein involved in Fe–S cluster biogenesis, and whose impairment causes a neurodegenerative disease called Friedreich’s Ataxia (FRDA). We summarize what is known about most common point mutations identified so far in heterozygous FRDA patients, their effects on frataxin structure and function and the consequences of its binding with partners.

Molecular Dynamics Simulation to Uncover the Mechanisms of Protein Instability During Freezing

Freezing is a common process applied in the pharmaceutical industry to store and transport biotherapeutics. Herewith, multi-scale molecular dynamics simulations of Lactate dehydrogenase (LDH) protein in phosphate buffer with/without ice formation performed to uncover the still poorly understood mechanisms and molecular details of protein destabilization upon freezing. Both fast and slow ice growing conditions were simulated at 243 K from one or two-side of the simulation box, respectively. The rate of ice formation at all-atom simulations was crucial to LDH stability, as faster freezing rates resulted in enhanced structural stability maintained by a higher number of intramolecular hydrogen bonds, less flexible protein's residues, lower solvent accessibility and greater structural compactness. Further, protein aggregation investigated by coarse-grained simulations was verified to be initiated by extended protein structures and retained by electrostatic interactions of the salt bridges between charged residues and hydrogen bonds between polar residues of the protein. Lastly, the study of free energy of dissociation through steered molecular dynamics simulation revealed LDH was destabilized by the solvation of the hydrophobic core and the loss of hydrophobic interactions. For the first time, experimentally validated molecular simulations revealed the detailed mechanisms of LDH destabilization upon ice formation and cryoconcentration of solutes.

Effects of Fe2+/Fe3+ Binding to Human Frataxin and Its D122Y Variant, as Revealed by Site-Directed Spin Labeling (SDSL) EPR Complemented by Fluorescence and Circular Dichroism Spectroscopies

Frataxin is a highly conserved protein whose deficiency results in the neurodegenerative disease Friederich's ataxia. Frataxin's actual physiological function has been debated for a long time without reaching a general agreement; however, it is commonly accepted that the protein is involved in the biosynthetic iron-sulphur cluster (ISC) machinery, and several authors have pointed out that it also participates in iron homeostasis. In this work, we use site-directed spin labeling coupled to electron paramagnetic resonance (SDSL EPR) to add new information on the effects of ferric and ferrous iron binding on the properties of human frataxin in vitro. Using SDSL EPR and relating the results to fluorescence experiments commonly performed to study iron binding to FXN, we produced evidence that ferric iron causes reversible aggregation without preferred interfaces in a concentration-dependent fashion, starting at relatively low concentrations (micromolar range), whereas ferrous iron binds without inducing aggregation. Moreover, our experiments show that the ferrous binding does not lead to changes of protein conformation. The data reported in this study reveal that the currently reported binding stoichiometries should be taken with caution. The use of a spin label resistant to reduction, as well as the comparison of the binding effect of Fe2+ in wild type and in the pathological D122Y variant of frataxin, allowed us to characterize the Fe2+ binding properties of different protein sites and highlight the effect of the D122Y substitution on the surrounding residues. We suggest that both Fe2+ and Fe3+ might play a relevant role in the context of the proposed FXN physiological functions.

Relationship between activity and stability: Design and characterization of stable variants of human frataxin

The relationships between conformational dynamics, stability and protein function are not obvious. Frataxin (FXN) is an essential protein that forms part of a supercomplex dedicated to the iron-sulfur (Fe-S) cluster assembly within the mitochondrial matrix. In humans, the loss of FXN expression or a decrease in its functionality results in Friedreich's Ataxia, a cardio-neurodegenerative disease. Recently, the way in which FXN interacts with the rest of the subunits of the supercomplex was uncovered. This opens a window to explore relationships between structural dynamics and function. In this study, we prepared a set of FXN variants spanning a broad range of conformational stabilities. Variants S160I, S160M and A204R were more stable than the wild-type and showed similar biological activity. Additionally, we prepared SILCAR, a variant that combines S160I, L203C and A204R mutations. SILCAR was 2.4 kcal mol-1 more stable and equally active. Some of the variants were significantly more resistant to proteolysis than the wild-type FXN. SILCAR showed the highest resistance, suggesting a more rigid structure. It was corroborated by means of molecular dynamics simulations. Relaxation dispersion NMR experiments comparing SILCAR and wild-type variants suggested similar internal motions in the microsecond to millisecond timescale. Instead, variant S157I showed higher denaturation resistance but a significant lower function, similarly to that observed for the FRDA variant N146K. We concluded that the contribution of particular side chains to the conformational stability of FXN might be highly subordinated to their impact on both the protein function and the stability of the functional supercomplex.

A Highly Conserved Iron-Sulfur Cluster Assembly Machinery between Humans and Amoeba Dictyostelium discoideum: The Characterization of Frataxin

Several biological activities depend on iron-sulfur clusters ([Fe-S]). Even though they are well-known in several organisms their function and metabolic pathway were poorly understood in the majority of the organisms. We propose to use the amoeba Dictyostelium discoideum, as a biological model to study the biosynthesis of [Fe-S] at the molecular, cellular and organism levels. First, we have explored the D. discoideum genome looking for genes corresponding to the subunits that constitute the molecular machinery for Fe-S cluster assembly and, based on the structure of the mammalian supercomplex and amino acid conservation profiles, we inferred the full functionality of the amoeba machinery. After that, we expressed the recombinant mature form of D. discoideum frataxin protein (DdFXN), the kinetic activator of this pathway. We characterized the protein and its conformational stability. DdFXN is monomeric and compact. The analysis of the secondary structure content, calculated using the far-UV CD spectra, was compatible with the data expected for the FXN fold, and near-UV CD spectra were compatible with the data corresponding to a folded protein. In addition, Tryptophan fluorescence indicated that the emission occurs from an apolar environment. However, the conformation of DdFXN is significantly less stable than that of the human FXN, (4.0 vs. 9.0 kcal mol-1, respectively). Based on a sequence analysis and structural models of DdFXN, we investigated key residues involved in the interaction of DdFXN with the supercomplex and the effect of point mutations on the energetics of the DdFXN tertiary structure. More than 10 residues involved in Friedreich's Ataxia are conserved between the human and DdFXN forms, and a good correlation between mutational effect on the energetics of both proteins were found, suggesting the existence of similar sequence/function/stability relationships. Finally, we integrated this information in an evolutionary context which highlights particular variation patterns between amoeba and humans that may reflect a functional importance of specific protein positions. Moreover, the complete pathway obtained forms a piece of evidence in favor of the hypothesis of a shared and highly conserved [Fe-S] assembly machinery between Human and D. discoideum.

Identification of novel inhibitor of protein tyrosine phosphatases delta: structure-based pharmacophore modeling, virtual screening, flexible docking, molecular dynamics simulation, and post-molecular dynamics analysis

Owing to their unique functions in regulating the synapse activity of protein tyrosine phosphatases delta (PTPδ) that has drawn special attention for developing drugs to autism spectrum disorders (ASDs). In this study, the PTPδ pharmacophore was first established by the structure-based pharmacophore method. Subsequently, 10 compounds contented Lipinski's rule of five was acquired by the virtual screening of the PTPδ pharmacophore against ZINC and PubChem databases. Then, the 10 identified molecules were discovered that had better binding affinity than a known PTPδ inhibitors compound SCHEMBL16375396. Two compounds SCHEMBL16375408 and ZINC19796658 with high binding score, low toxicity were gained. They were observed by docking analysis and molecular dynamics simulations that the novel potential inhibitors not only possessed the same function as SCHEMBL16375396 did in inhibiting PTPδ, but also had more favorable conformation to bind with the catalytic active regions. This study provides a new method for identify PTPδ inhibitor for the treatment of ASDs disease.Communicated by Ramaswamy H. Sarma.

Folding and Dynamics Are Strongly pH-Dependent in a Psychrophile Frataxin

Protein dynamics, folding, and thermodynamics represent a central aspect of biophysical chemistry. pH, temperature, and denaturant perturbations inform our understanding of diverse contributors to stability and rates. In this work, we performed a thermodynamic analysis using a combined experimental and computational approach to gain insights into the role of electrostatics in the folding reaction of a psychrophile frataxin variant from Psychromonas ingrahamii. This folding reaction is strongly modulated by pH with a single, narrow, and well-defined transition state with ∼80% compactness, ∼70% electrostatic interactions, and ∼60% hydration shell compared to the native state (αD = 0.82, αH = 0.67, and αΔCp = 0.59). Our results are best explained by a two-proton/two-state model with very different pKa values of the native and denatured states (∼5.5 and ∼8.0, respectively). As a consequence, the stability strongly increases from pH 8.0 to 6.0 (|ΔΔG°| = 5.2 kcal mol-1), mainly because of a decrease in the TΔS°. Variation of ΔH° and ΔS° at pH below 7.0 is dominated by a change in ΔHf⧧ and ΔSf⧧, while at pH above 7.0, it is governed by ΔHu⧧ and ΔSu⧧. Molecular dynamics simulations showed that these pH modulations could be explained by the fluctuations of two regions, rich in electrostatic contacts, whose dynamics are pH-dependent and motions are strongly correlated. Results presented herein contribute to the understanding of the stability and dynamics of this frataxin variant, pointing to an intrinsic feature of the family topology to support different folding mechanisms.

Evaluating the predictions of the protein stability change upon single amino acid substitutions for the FXN CAGI5 challenge

Frataxin (FXN) is a highly conserved protein found in prokaryotes and eukaryotes that is required for efficient regulation of cellular iron homeostasis. Experimental evidence associates amino acid substitutions of the FXN to Friedreich Ataxia, a neurodegenerative disorder. Recently, new thermodynamic experiments have been performed to study the impact of somatic variations identified in cancer tissues on protein stability. The Critical Assessment of Genome Interpretation (CAGI) data provider at the University of Rome measured the unfolding free energy of a set of variants (FXN challenge data set) with far-UV circular dichroism and intrinsic fluorescence spectra. These values have been used to calculate the change in unfolding free energy between the variant and wild-type proteins at zero concentration of denaturant ( Δ Δ G H 2 O ) . The FXN challenge data set, composed of eight amino acid substitutions, was used to evaluate the performance of the current computational methods for predicting the Δ Δ G H 2 O value associated with the variants and to classify them as destabilizing and not destabilizing. For the fifth edition of CAGI, six independent research groups from Asia, Australia, Europe, and North America submitted 12 sets of predictions from different approaches. In this paper, we report the results of our assessment and discuss the limitations of the tested algorithms.

Global Implications of Local Unfolding Phenomena, Probed by Cysteine Reactivity in Human Frataxin

Local events that affect specific regions of proteins are of utmost relevance for stability and function. The aim of this study is to quantitatively assess the importance of locally-focused dynamics by means of a simple chemical modification procedure. Taking human Frataxin as a working model, we investigated local fluctuations of the C-terminal region (the last 16 residues of the protein) by means of three L → C replacement mutants: L98C, L200C and L203C. The conformation and thermodynamic stability of each variant was assessed. All the variants exhibited native features and high stabilities: 9.1 (wild type), 8.1 (L198C), 7.0 (L200C) and 10.0 kcal mol^-1 (L203C). In addition, kinetic rates of Cys chemical modification by DTNB and DTDPy were measured, conformational dynamics data were extracted and free energy for the local unfolding of the C-terminal region was estimated. The analysis of these results indicates that the conformation of the C-terminal region fluctuates with partial independence from global unfolding events. Additionally, numerical fittings of the kinetic model of the process suggest that the local transition occurs in the seconds to minutes timescale. In fact, standard free energy differences for local unfolding were found to be significantly lower than those of the global unfolding reaction, showing that chemical modification results may not be explained in terms of the global unfolding reaction alone. These results provide unequivocal experimental evidence of local phenomena with global effects and contribute to understanding how global and local stability are linked to protein dynamics.

PFDB: A standardized protein folding database with temperature correction

We constructed a standardized protein folding kinetics database (PFDB) in which the logarithmic rate constants of all listed proteins are calculated at the standard temperature (25 °C). A temperature correction based on the Eyring–Kramers equation was introduced for proteins whose folding kinetics were originally measured at temperatures other than 25 °C. We verified the temperature correction by comparing the logarithmic rate constants predicted and experimentally observed at 25 °C for 14 different proteins, and the results demonstrated improvement of the quality of the database. PFDB consists of 141 (89 two-state and 52 non-two-state) single-domain globular proteins, which has the largest number among the currently available databases of protein folding kinetics. PFDB is thus intended to be used as a standard for developing and testing future predictive and theoretical studies of protein folding. PFDB can be accessed from the following link: http://lee.kias.re.kr/~bala/PFDB.

Frataxin Structure and Function

Mammalian frataxin is a small mitochondrial protein involved in iron sulfur cluster assembly. Frataxin deficiency causes the neurodegenerative disease Friedreich's Ataxia. Valuable knowledge has been gained on the structural dynamics of frataxin, metal-ion-protein interactions, as well as on the effect of mutations on protein conformation, stability and internal motions. Additionally, laborious studies concerning the enzymatic reactions involved have allowed for understanding the capability of frataxin to modulate Fe-S cluster assembly function. Remarkably, frataxin biological function depends on its interaction with some proteins to form a supercomplex, among them NFS1 desulfurase and ISCU, the scaffolding protein. By combining multiple experimental tools including high resolution techniques like NMR and X-ray, but also SAXS, crosslinking and mass-spectrometry, it was possible to build a reliable model of the structure of the desulfurase supercomplex NFS1/ACP-ISD11/ISCU/frataxin. In this chapter, we explore these issues showing how the scientific view concerning frataxin structure-function relationships has evolved over the last years.

Biophysical characterisation of the recombinant human frataxin precursor

Friedreich's ataxia is a disease caused by a decrease in the levels of expression or loss of functionality of the mitochondrial protein frataxin (FXN). The development of an active and stable recombinant variant of FXN is important for protein replacement therapy. Although valuable data about the mature form FXN81-210 has been collected, not enough information is available about the conformation of the frataxin precursor (FXN1-210). We investigated the conformation, stability and function of a recombinant precursor variant (His6-TAT-FXN1-210), which includes a TAT peptide in the N-terminal region to assist with transport across cell membranes. His6-TAT-FXN1-210 was expressed in Escherichia coli and conditions were found for purifying folded protein free of aggregation, oxidation or degradation, even after freezing and thawing. The protein was found to be stable and monomeric, with the N-terminal stretch (residues 1-89) mostly unstructured and the C-terminal domain properly folded. The experimental data suggest a complex picture for the folding process of full-length frataxin in vitro: the presence of the N-terminal region increased the tendency of FXN to aggregate at high temperatures but this could be avoided by the addition of low concentrations of GdmCl. The purified precursor was translocated through cell membranes. In addition, immune response against His6-TAT-FXN1-210 was measured, suggesting that the C-terminal fragment was not immunogenic at the assayed protein concentrations. Finally, the recognition of recombinant FXN by cellular proteins was studied to evaluate its functionality. In this regard, cysteine desulfurase NFS1/ISD11/ISCU was activated in vitro by His6-TAT-FXN1-210. Moreover, the results showed that His6-TAT-FXN1-210 can be ubiquitinated in vitro by the recently identified frataxin E3 ligase RNF126, in a similar way as the FXN1-210, suggesting that the His6-TAT extension does not interfere with the ubiquitination machinery.

Insights on the conformational dynamics of human frataxin through modifications of loop-1

Human frataxin (FXN) is a highly conserved mitochondrial protein involved in iron homeostasis and activation of the iron-sulfur cluster assembly. FXN deficiency causes the neurodegenerative disease Friedreich's Ataxia. Here, we investigated the effect of alterations in loop-1, a stretch presumably essential for FXN function, on the conformational stability and dynamics of the native state. We generated four loop-1 variants, carrying substitutions, insertions and deletions. All of them were stable and well-folded proteins. Fast local motions (ps-ns) and slower long-range conformational dynamics (μs-ms) were altered in some mutants as judged by NMR. Particularly, loop-1 modifications impact on the dynamics of a distant region that includes residues from the β-sheet, helix α1 and the C-terminal. Remarkably, all the mutants retain the ability to activate cysteine desulfurase, even when two of them exhibit a strong decrease in iron binding, revealing a differential sensitivity of these functional features to loop-1 perturbation. Consequently, we found that even for a small and relatively rigid protein, engineering a loop segment enables to alter conformational dynamics through a long-range effect, preserving the native-state structure and important aspects of function.

How Molecular Size Impacts RMSD Applications in Molecular Dynamics Simulations

The root-mean-square deviation (RMSD) is a similarity measure widely used in analysis of macromolecular structures and dynamics. As increasingly larger macromolecular systems are being studied, dimensionality effects such as the "curse of dimensionality" (a diminishing ability to discriminate pairwise differences between conformations with increasing system size) may exist and significantly impact RMSD-based analyses. For such large bimolecular systems, whether the RMSD or other alternative similarity measures might suffer from this "curse" and lose the ability to discriminate different macromolecular structures had not been explicitly addressed. Here, we show such dimensionality effects for both weighted and nonweighted RMSD schemes. We also provide a mechanism for the emergence of the "curse of dimensionality" for RMSD from the law of large numbers by showing that the conformational distributions from which RMSDs are calculated become increasingly similar as the system size increases. Our findings suggest the use of weighted RMSD schemes for small proteins (less than 200 residues) and nonweighted RMSD for larger proteins when analyzing molecular dynamics trajectories.