Searching for Frataxin Function: Exploring the Analogy with Nqo15, the Frataxin-like Protein of Respiratory Complex I from Thermus thermophilus.

Nqo15 is a subunit of respiratory complex I of the bacterium Thermus thermophilus, with strong structural similarity to human frataxin (FXN), a protein involved in the mitochondrial disease Friedreich's ataxia (FRDA). Recently, we showed that the expression of recombinant Nqo15 can ameliorate the respiratory phenotype of FRDA patients' cells, and this prompted us to further characterize both the Nqo15 solution's behavior and its potential functional overlap with FXN, using a combination of in silico and in vitro techniques. We studied the analogy of Nqo15 and FXN by performing extensive database searches based on sequence and structure. Nqo15's folding and flexibility were investigated by combining nuclear magnetic resonance (NMR), circular dichroism, and coarse-grained molecular dynamics simulations. Nqo15's iron-binding properties were studied using NMR, fluorescence, and specific assays and its desulfurase activation by biochemical assays. We found that the recombinant Nqo15 isolated from complex I is monomeric, stable, folded in solution, and highly dynamic. Nqo15 does not share the iron-binding properties of FXN or its desulfurase activation function.

A STRP-ed definition of Structured Tandem Repeats in Proteins.

Tandem Repeat Proteins (TRPs) are a class of proteins with repetitive amino acid sequences that have been studied extensively for over two decades. Different features at the level of sequence, structure, function and evolution have been attributed to them by various authors. And yet many of its salient features appear only when looking at specific subclasses of protein tandem repeats. Here, we attempt to rationalize the existing knowledge on Tandem Repeat Proteins (TRPs) by pointing out several dichotomies. The emerging picture is more nuanced than generally assumed and allows us to draw some boundaries of what is not a "proper" TRP. We conclude with an operational definition of a specific subset, which we have denominated STRPs (Structural Tandem Repeat Proteins), which separates a subclass of tandem repeats with distinctive features from several other less well-defined types of repeats. We believe that this definition will help researchers in the field to better characterize the biological meaning of this large yet largely understudied group of proteins.

CardIAP: calcium transients confocal image analysis tool.

One of the main topics of cardiovascular research is the study of calcium (Ca2+) handling, as even small changes in Ca2+ concentration can alter cell functionality (Bers, Annu Rev Physiol, 2014, 76, 107-127). Ionic calcium (Ca2+) plays the role of a second messenger in eukaryotic cells, associated with cellular functions such as cell cycle regulation, transport, motility, gene expression, and regulation. The use of fluorometric techniques in isolated cells loaded with Ca2+-sensitive fluorescent probes allows quantitative measurement of dynamic events occurring in living, functioning cells. The Cardiomyocytes Images Analyzer Python (CardIAP) application addresses the need to analyze and retrieve information from confocal microscopy images systematically, accurately, and rapidly. Here we present CardIAP, an open-source tool developed entirely in Python, freely available and useable in an interactive web application. In addition, CardIAP can be used as a standalone Python library and freely installed via PIP, making it easy to integrate into biomedical imaging pipelines. The images that can be generated in the study of the heart have the particularity of requiring both spatial and temporal analysis. CardIAP aims to open the field of cardiomyocytes and intact hearts image processing. The improvement in the extraction of information from the images will allow optimizing the usage of resources and animals. With CardIAP, users can run the analysis to both, the complete image, and portions of it in an easy way, and replicate it on a series of images. This analysis provides users with information on the spatial and temporal changes in calcium releases and characterizes them. The web application also allows users to extract calcium dynamics data in downloadable tables, simplifying the calculation of alternation and discordance indices and their classification. CardIAP aims to provide a tool that could assist biomedical researchers in studying the underlying mechanisms of anomalous calcium release phenomena.

Combining Protein Conformational Diversity and Phylogenetic Information Using CoDNaS and CoDNaS-Q.

CoDNaS (http://ufq.unq.edu.ar/codnas/) and CoDNaS-Q (http://ufq.unq.edu.ar/codnasq) are repositories of proteins with different degrees of conformational diversity. Following the ensemble nature of the native state, conformational diversity represents the structural differences between the conformers in the ensemble. Each entry in CoDNaS and CoDNaS-Q contains a redundant collection of experimentally determined conformers obtained under different conditions. These conformers represent snapshots of the protein dynamism. While CoDNaS contains examples of conformational diversity at the tertiary level, a recent development, CoDNaS-Q, contains examples at the quaternary level. In the emerging age of accurate protein structure prediction by machine learning approaches, many questions remain open regarding the characterization of protein dynamism. In this context, most bioinformatics resources take advantage of distinct features derived from protein alignments, however, the complexity and heterogeneity of information makes it difficult to recover reliable biological signatures. Here we present five protocols to explore tertiary and quaternary conformational diversity at the individual protein level as well as for the characterization of the distribution of conformational diversity at the protein family level in a phylogenetic context. These protocols can provide curated protein families with experimentally known conformational diversity, facilitating the exploration of sequence determinants of protein dynamism. © 2023 Wiley Periodicals LLC. Basic Protocol 1: Assessing conformational diversity with CoDNaS Alternate Protocol 1: Assessing conformational diversity at the quaternary level with CoDNaS-Q Basic Protocol 2: Exploring conformational diversity in a protein family Alternate Protocol 2: Exploring quaternary conformational diversity in a protein family Basic Protocol 3: Representing conformational diversity in a phylogenetic context.

Impact of protein conformational diversity on AlphaFold predictions.

MOTIVATION: After the outstanding breakthrough of AlphaFold in predicting protein 3D models, new questions appeared and remain unanswered. The ensemble nature of proteins, for example, challenges the structural prediction methods because the models should represent a set of conformers instead of single structures. The evolutionary and structural features captured by effective deep learning techniques may unveil the information to generate several diverse conformations from a single sequence. Here, we address the performance of AlphaFold2 predictions obtained through ColabFold under this ensemble paradigm. RESULTS: Using a curated collection of apo-holo pairs of conformers, we found that AlphaFold2 predicts the holo form of a protein in ∼70% of the cases, being unable to reproduce the observed conformational diversity with the same error for both conformers. More importantly, we found that AlphaFold2's performance worsens with the increasing conformational diversity of the studied protein. This impairment is related to the heterogeneity in the degree of conformational diversity found between different members of the homologous family of the protein under study. Finally, we found that main-chain flexibility associated with apo-holo pairs of conformers negatively correlates with the predicted local model quality score plDDT, indicating that plDDT values in a single 3D model could be used to infer local conformational changes linked to ligand binding transitions. AVAILABILITY AND IMPLEMENTATION: Data and code used in this manuscript are publicly available at https://gitlab.com/sbgunq/publications/af2confdiv-oct2021. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

CoDNaS-RNA: a database of conformational diversity in the native state of RNA.

SUMMARY: Conformational changes in RNA native ensembles are central to fulfill many of their biological roles. Systematic knowledge of the extent and possible modulators of this conformational diversity is desirable to better understand the relationship between RNA dynamics and function. We have developed CoDNaS-RNA as the first database of conformational diversity in RNA molecules. Known RNA structures are retrieved and clustered to identify alternative conformers of each molecule. Pairwise structural comparisons between all conformers within each cluster allows to measure the variability of the molecule. Additional annotations about structural features, molecular interactions and biological function are provided. All data in CoDNaS-RNA is free to download and available as a public website that can be of interest for researchers in computational biology and other life science disciplines. AVAILABILITY AND IMPLEMENTATION: The data underlying this article are available at http://ufq.unq.edu.ar/codnasrna or https://codnas-rna.bioinformatica.org/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

"Protein" no longer means what it used to.

Every biologist knows that the word protein describes a group of macromolecules essential to sustain life on Earth. As biologists, we are invariably trained under a protein paradigm established since the early twentieth century. However, in recent years, the term protein unveiled itself as an euphemism to describe the overwhelming heterogeneity of these compounds. Most of our current studies are targeted on carefully selected subsets of proteins, but we tend to think and write about these as representative of the whole population. Here we discuss how seeking for universal definitions and general rules in any arbitrarily segmented study would be misleading about the conclusions. Of course, it is not our purpose to discourage the use of the word protein. Instead, we suggest to embrace the extended universe of proteins to reach a deeper understanding of their full potential, realizing that the term encompasses a group of molecules very heterogeneous in terms of size, shape, chemistry and functions, i.e. the term protein no longer means what it used to.

RepeatsDB in 2021: improved data and extended classification for protein tandem repeat structures.

The RepeatsDB database (URL: https://repeatsdb.org/) provides annotations and classification for protein tandem repeat structures from the Protein Data Bank (PDB). Protein tandem repeats are ubiquitous in all branches of the tree of life. The accumulation of solved repeat structures provides new possibilities for classification and detection, but also increasing the need for annotation. Here we present RepeatsDB 3.0, which addresses these challenges and presents an extended classification scheme. The major conceptual change compared to the previous version is the hierarchical classification combining top levels based solely on structural similarity (Class > Topology > Fold) with two new levels (Clan > Family) requiring sequence similarity and describing repeat motifs in collaboration with Pfam. Data growth has been addressed with improved mechanisms for browsing the classification hierarchy. A new UniProt-centric view unifies the increasingly frequent annotation of structures from identical or similar sequences. This update of RepeatsDB aligns with our commitment to develop a resource that extracts, organizes and distributes specialized information on tandem repeat protein structures.

Revenant: a database of resurrected proteins.

Revenant is a database of resurrected proteins coming from extinct organisms. Currently, it contains a manually curated collection of 84 resurrected proteins derived from bibliographic data. Each protein is extensively annotated, including structural, biochemical and biophysical information. Revenant contains a browse capability designed as a timeline from where the different proteins can be accessed. The oldest Revenant entries are between 4200 and 3500 million years ago, while the younger entries are between 8.8 and 6.3 million years ago. These proteins have been resurrected using computational tools called ancestral sequence reconstruction techniques combined with wet-laboratory synthesis and expression. Resurrected proteins are commonly used, with a noticeable increase during the past years, to explore and test different evolutionary hypotheses such as protein stability, to explore the origin of new functions, to get biochemical insights into past metabolisms and to explore specificity and promiscuous behaviour of ancient proteins.

Ensembles from Ordered and Disordered Proteins Reveal Similar Structural Constraints during Evolution.

The conformations accessible to proteins are determined by the inter-residue interactions between amino acid residues. During evolution, structural constraints that are required for protein function providing biologically relevant information can exist. Here, we studied the proportion of sites evolving under structural constraints in two very different types of ensembles, those coming from ordered and disordered proteins. Using a structurally constrained model of protein evolution, we found that both types of ensembles show comparable, near 40%, number of positions evolving under structural constraints. Among these sites, ~68% are in disordered regions and ~57% of them show long-range inter-residue contacts. Also, we found that disordered ensembles are redundant in reference to their structurally constrained evolutionary information and could be described on average with ~11 conformers. Despite the different complexity of the studied ensembles and proteins, the similar constraints reveal a comparable level of selective pressure to maintain their biological functions. These results highlight the importance of the evolutionary information to recover meaningful biological information to further characterize conformational ensembles.

Bioinformatics calls the school: Use of smartphones to introduce Python for bioinformatics in high schools.

The dynamic nature of technological developments invites us to rethink the learning spaces. In this context, science education can be enriched by the contribution of new computational resources, making the educational process more up-to-date, challenging, and attractive. Bioinformatics is a key interdisciplinary field, contributing to the understanding of biological processes that is often underrated in secondary schools. As a useful resource in learning activities, bioinformatics could help in engaging students to integrate multiple fields of knowledge (logical-mathematical, biological, computational, etc.) and generate an enriched and long-lasting learning environment. Here, we report our recent project in which high school students learned basic concepts of programming applied to solving biological problems. The students were taught the Python syntax, and they coded simple tools to answer biological questions using resources at hand. Notably, these were built mostly on the students' own smartphones, which proved to be capable, readily available, and relevant complementary tools for teaching. This project resulted in an empowering and inclusive experience that challenged differences in social background and technological accessibility.

On the dynamical incompleteness of the Protein Data Bank.

Major scientific challenges that are beyond the capability of individuals need to be addressed by multi-disciplinary and multi-institutional consortia. Examples of these endeavours include the Human Genome Project, and more recently, the Structural Genomics (SG) initiative. The SG initiative pursues the expansion of structural coverage to include at least one structural representative for each protein family to derive the remaining structures using homology modelling. However, biological function is inherently connected with protein dynamics that can be studied by knowing different structures of the same protein. This ensemble of structures provides snapshots of protein conformational diversity under native conditions. Thus, sequence redundancy in the Protein Data Bank (PDB) (i.e. crystallization of the same protein under different conditions) is therefore an essential input contributing to experimentally based studies of protein dynamics and providing insights into protein function. In this work, we show that sequence redundancy, a key concept for exploring protein dynamics, is highly biased and fundamentally incomplete in the PDB. Additionally, our results show that dynamical behaviour of proteins cannot be inferred using homologous proteins. Minor to moderate changes in sequence can produce great differences in dynamical behaviour. Nonetheless, the structural and dynamical incompleteness of the PDB is apparently unrelated concepts in SG. While the first could be reversed by promoting the extension of the structural coverage, we would like to emphasize that further focused efforts will be needed to amend the incompleteness of the PDB in terms of dynamical information content, essential to fully understand protein function.

Exploring Protein Conformational Diversity.

The native state of proteins is composed of conformers in dynamical equilibrium. In this chapter, different issues related to conformational diversity are explored using a curated and experimentally based database called CoDNaS (Conformational Diversity in the Native State). This database is a collection of redundant structures for the same sequence. CoDNaS estimates the degree of conformational diversity using different global and local structural similarity measures. It allows the user to explore how structural differences among conformers change as a function of several structural features providing further biological information. This chapter explores the measurement of conformational diversity and its relationship with sequence divergence. Also, it discusses how proteins with high conformational diversity could affect homology modeling techniques.

Dynamics fingerprints of active conformers of epidermal growth factor receptor kinase.

Epidermal growth factor receptor (EGFR) is a prototypical cell-surface receptor that plays a key role in the regulation of cellular signaling, proliferation and differentiation. Mutations of its kinase domain have been associated with the development of a variety of cancers and, therefore, it has been the target of drug design. Single amino acid substitutions (SASs) in this domain have been proven to alter the equilibrium of pre-existing conformer populations. Despite the advances in structural descriptions of its so-called active and inactive conformations, the associated dynamics aspects that characterize them have not been thoroughly studied yet. As the dynamic behaviors and molecular motions of proteins are important for a complete understanding of their structure-function relationships we present a novel procedure, using (or based on) normal mode analysis, to identify the collective dynamics shared among different conformers in EGFR kinase. The method allows the comparison of patterns of low-frequency vibrational modes defining representative directions of motions. Our procedure is able to emphasize the main similarities and differences between the collective dynamics of different conformers. In the case of EGFR kinase, two representative directions of motions have been found as dynamics fingerprints of the active conformers. Protein motion along both directions reveals to have a significant impact on the cavity volume of the main pocket of the active site. Otherwise, the inactive conformers exhibit a more heterogeneous distribution of collective motions. © 2018 Wiley Periodicals, Inc.

Correction: Pockets as structural descriptors of EGFR kinase conformations.

[This corrects the article DOI: 10.1371/journal.pone.0189147.].

Two novel unstable hemoglobin variants due to in-frame deletions of key amino acids in the ß-globin chain.

Hemoglobinopathies are the most common autosomal recessive disorders and are mostly inherited in a recessive manner. However, certain mutations can affect the globin chain stability, leading to dominant forms of thalassemia. The aim of this work was the molecular and structural characterization of two heterozygous in-frame deletions, leading to ß-globin variants in pediatric patients in Argentina. The HBB gene of the probands and their parents was sequenced, and other markers of globin chain imbalance were analyzed. Several structural analyses were performed, and the effect of the mutations on the globin chain stability was analyzed. In Hb JC-Paz, HBB:c.29_37delCTGCCGTTA (p.Ala10_Thr12del), detected in an Argentinean boy, one α-helix turn is expected to be lost. In Hb Tavapy, HBB:c.182_187delTGAAGG (p.Val60_Lys61del), the deleted residues are close to distal histidine (His63) in the heme pocket. Both mutations are predicted to have a destabilizing effect. The development of computational structural models and bioinformatics algorithms is expected to become a useful tool to understand the impact of the mutations leading to dominant thalassemia.

Pockets as structural descriptors of EGFR kinase conformations.

Epidermal Growth Factor Receptor (EGFR), a tyrosine kinase receptor, is one of the main tumor markers in different types of cancers. The kinase native state is mainly composed of two populations of conformers: active and inactive. Several sequence variations in EGFR kinase region promote the differential enrichment of conformers with higher activity. Some structural characteristics have been proposed to differentiate kinase conformations, but these considerations could lead to ambiguous classifications. We present a structural characterisation of EGFR kinase conformers, focused on active site pocket comparisons, and the mapping of known pathological sequence variations. A structural based clustering of this pocket accurately discriminates active from inactive, well-characterised conformations. Furthermore, this main pocket contains, or is in close contact with, ≈65% of cancer-related variation positions. Although the relevance of protein dynamics to explain biological function has been extensively recognised, the usage of the ensemble of conformations in dynamic equilibrium to represent the functional state of proteins and the importance of pockets, cavities and/or tunnels was often neglected in previous studies. These functional structures and the equilibrium between them could be structurally analysed in wild type as well as in sequence variants. Our results indicate that biologically important pockets, as well as their shape and dynamics, are central to understanding protein function in wild-type, polymorphic or disease-related variations.

Conformational diversity analysis reveals three functional mechanisms in proteins.

Protein motions are a key feature to understand biological function. Recently, a large-scale analysis of protein conformational diversity showed a positively skewed distribution with a peak at 0.5 Å C-alpha root-mean-square-deviation (RMSD). To understand this distribution in terms of structure-function relationships, we studied a well curated and large dataset of ~5,000 proteins with experimentally determined conformational diversity. We searched for global behaviour patterns studying how structure-based features change among the available conformer population for each protein. This procedure allowed us to describe the RMSD distribution in terms of three main protein classes sharing given properties. The largest of these protein subsets (~60%), which we call "rigid" (average RMSD = 0.83 Å), has no disordered regions, shows low conformational diversity, the largest tunnels and smaller and buried cavities. The two additional subsets contain disordered regions, but with differential sequence composition and behaviour. Partially disordered proteins have on average 67% of their conformers with disordered regions, average RMSD = 1.1 Å, the highest number of hinges and the longest disordered regions. In contrast, malleable proteins have on average only 25% of disordered conformers and average RMSD = 1.3 Å, flexible cavities affected in size by the presence of disordered regions and show the highest diversity of cognate ligands. Proteins in each set are mostly non-homologous to each other, share no given fold class, nor functional similarity but do share features derived from their conformer population. These shared features could represent conformational mechanisms related with biological functions.

Addressing the Role of Conformational Diversity in Protein Structure Prediction.

Computational modeling of tertiary structures has become of standard use to study proteins that lack experimental characterization. Unfortunately, 3D structure prediction methods and model quality assessment programs often overlook that an ensemble of conformers in equilibrium populates the native state of proteins. In this work we collected sets of publicly available protein models and the corresponding target structures experimentally solved and studied how they describe the conformational diversity of the protein. For each protein, we assessed the quality of the models against known conformers by several standard measures and identified those models ranked best. We found that model rankings are defined by both the selected target conformer and the similarity measure used. 70% of the proteins in our datasets show that different models are structurally closest to different conformers of the same protein target. We observed that model building protocols such as template-based or ab initio approaches describe in similar ways the conformational diversity of the protein, although for template-based methods this description may depend on the sequence similarity between target and template sequences. Taken together, our results support the idea that protein structure modeling could help to identify members of the native ensemble, highlight the importance of considering conformational diversity in protein 3D quality evaluations and endorse the study of the variability of the native structure for a meaningful biological analysis.

Disorder transitions and conformational diversity cooperatively modulate biological function in proteins.

Structural differences between conformers sustain protein biological function. Here, we studied in a large dataset of 745 intrinsically disordered proteins, how ordered-disordered transitions modulate structural differences between conformers as derived from crystallographic data. We found that almost 50% of the proteins studied show no transitions and have low conformational diversity while the rest show transitions and a higher conformational diversity. In this last subset, 60% of the proteins become more ordered after ligand binding, while 40% more disordered. As protein conformational diversity is inherently connected with protein function our analysis suggests differences in structure-function relationships related to order-disorder transitions.