A comprehensive analysis of SARS-CoV-2 missense mutations indicates that all possible amino acid replacements in the viral proteins occurred within the first two-and-a-half years of the pandemic.

The surveillance of COVID-19 pandemic has led to the determination of millions of genome sequences of the SARS-CoV-2 virus, with the accumulation of a wealth of information never collected before for an infectious disease. Exploring the information retrieved from the GISAID database reporting at that time >13 million genome sequences, we classified the 141,639 unique missense mutations detected in the first two-and-a-half years (up to October 2022) of the pandemic. Notably, our analysis indicates that 98.2 % of all possible conservative amino acid replacements occurred. Even non-conservative mutations were highly represented (73.9 %). For a significant number of residues (3 %), all possible replacements with the other nineteen amino acids have been observed. These observations strongly indicate that, in this time interval, the virus explored all possible alternatives in terms of missense mutations for all sites of its polypeptide chain and that those that are not observed severely affect SARS-CoV-2 integrity. The implications of the present findings go well beyond the structural biology of SARS-CoV-2 as the huge amount of information here collected and classified may be valuable for the elucidation of the sequence-structure-function relationships in proteins.

A Comprehensive Analysis of the Structural Recognition between KCTD Proteins and Cullin 3.

KCTD ((K)potassium Channel Tetramerization Domain-containing) proteins constitute an emerging class of proteins involved in fundamental physio-pathological processes. In these proteins, the BTB domain, which represents the defining element of the family, may have the dual role of promoting oligomerization and favoring functionally important partnerships with different interactors. Here, by exploiting the potential of recently developed methodologies for protein structure prediction, we report a comprehensive analysis of the interactions of all KCTD proteins with their most common partner Cullin 3 (Cul3). The data here presented demonstrate the impressive ability of this approach to discriminate between KCTDs that interact with Cul3 and those that do not. Indeed, reliable and stable models of the complexes were only obtained for the 15 members of the family that are known to interact with Cul3. The generation of three-dimensional models for all KCTD-Cul3 complexes provides interesting clues on the determinants of the structural basis of this partnership as clear structural differences emerged between KCTDs that bind or do not bind Cul3. Finally, the availability of accurate three-dimensional models for KCTD-Cul3 interactions may be valuable for the ad hoc design and development of compounds targeting specific KCTDs that are involved in several common diseases.

Structural and Dynamic-Based Characterization of the Recognition Patterns of E7 and TRP-2 Epitopes by MHC Class I Receptors through Computational Approaches.

A detailed comprehension of MHC-epitope recognition is essential for the design and development of new antigens that could be effectively used in immunotherapy. Yet, the high variability of the peptide together with the large abundance of MHC variants binding makes the process highly specific and large-scale characterizations extremely challenging by standard experimental techniques. Taking advantage of the striking predictive accuracy of AlphaFold, we report a structural and dynamic-based strategy to gain insights into the molecular basis that drives the recognition and interaction of MHC class I in the immune response triggered by pathogens and/or tumor-derived peptides. Here, we investigated at the atomic level the recognition of E7 and TRP-2 epitopes to their known receptors, thus offering a structural explanation for the different binding preferences of the studied receptors for specific residues in certain positions of the antigen sequences. Moreover, our analysis provides clues on the determinants that dictate the affinity of the same epitope with different receptors. Collectively, the data here presented indicate the reliability of the approach that can be straightforwardly extended to a large number of related systems.

Self-Assembled Materials Based on Fully Aromatic Peptides: The Impact of Tryptophan, Tyrosine, and Dopa Residues.

Peptides are able to self-organize in structural elements including cross-ß structures. Taking advantage of this tendency, in the last decades, peptides have been scrutinized as molecular elements for the development of multivalent supramolecular architectures. In this context, different classes of peptides, also with completely aromatic sequences, were proposed. Our previous studies highlighted that the (FY)3 peptide, which alternates hydrophobic phenylalanine and more hydrophilic tyrosine residues, is able to self-assemble, thanks to the formation of both polar and apolar interfaces. It was observed that the replacement of Phe and Tyr residues with other noncoded aromatic amino acids like 2-naphthylalanine (Nal) and Dopa affects the interactions among peptides with consequences on the supramolecular organization. Herein, we have investigated the self-assembling behavior of two novel (FY)3 analogues with Trp and Dopa residues in place of the Phe and Tyr ones, respectively. Additionally, PEGylation of the N-terminus was analyzed too. The supramolecular organization, morphology, and capability to gel were evaluated using complementary techniques, including fluorescence, Fourier transform infrared spectroscopy, and scanning electron microscopy. Structural periodicities along and perpendicular to the fiber axis were detected by grazing incidence wide-angle X-ray scattering. Finally, molecular dynamics studies provided interesting insights into the atomic structure of the cross-ß that constitutes the basic motif of the assemblies formed by these novel peptide systems.

Structural Insights into Protein-Aptamer Recognitions Emerged from Experimental and Computational Studies.

Aptamers are synthetic nucleic acids that are developed to target with high affinity and specificity chemical entities ranging from single ions to macromolecules and present a wide range of chemical and physical properties. Their ability to selectively bind proteins has made these compounds very attractive and versatile tools, in both basic and applied sciences, to such an extent that they are considered an appealing alternative to antibodies. Here, by exhaustively surveying the content of the Protein Data Bank (PDB), we review the structural aspects of the protein-aptamer recognition process. As a result of three decades of structural studies, we identified 144 PDB entries containing atomic-level information on protein-aptamer complexes. Interestingly, we found a remarkable increase in the number of determined structures in the last two years as a consequence of the effective application of the cryo-electron microscopy technique to these systems. In the present paper, particular attention is devoted to the articulated architectures that protein-aptamer complexes may exhibit. Moreover, the molecular mechanism of the binding process was analyzed by collecting all available information on the structural transitions that aptamers undergo, from their protein-unbound to the protein-bound state. The contribution of computational approaches in this area is also highlighted.

New insight into the traditional model of the coagulation cascade and its regulation: illustrated review of a three-dimensional view.

The coagulation process relies on an intricate network of three-dimensional structural interactions and subtle biological regulations. In the present review, we illustrate the state of the art of the structural biology of the coagulation cascade by surveying the Protein Data Bank and the EBI AlphaFold databases. Investigations performed in the last decade have provided structural information on essentially all players involved in the process. Indeed, the initial characterization of specific and rather canonical domains has been progressively extended to complicated multidomain proteins. Recently, the application of cryogenic electron microscopy techniques has unraveled the structural features of highly complex coagulation factors, which has led to enhanced understanding. This review initially focuses on the structure of the individual factors as a function of their involvement in intrinsic, extrinsic, and common pathways. A specific emphasis is given to what is known or unknown on the structural basis of each step of the cascade. Available data providing clues on the structural recognition of the factors involved in the functional partnerships of the pathways are illustrated. Recent structures of important complexes formed by these proteins with regulators are described, focusing on the drugs used as anticoagulants and on their reversal agents. Finally, we highlight the different roles that innovative biomolecules such as aptamers may have in the regulation of the cascade.

A Comprehensive Analysis of the Intrinsic Visible Fluorescence Emitted by Peptide/Protein Amyloid-like Assemblies.

Amyloid aggregation is a widespread process that involves proteins and peptides with different molecular complexity and amino acid composition. The structural motif (cross-ß) underlying this supramolecular organization generates aggregates endowed with special mechanical and spectroscopic properties with huge implications in biomedical and technological fields, including emerging precision medicine. The puzzling ability of these assemblies to emit intrinsic and label-free fluorescence in regions of the electromagnetic spectrum, such as visible and even infrared, usually considered to be forbidden in the polypeptide chain, has attracted interest for its many implications in both basic and applied science. Despite the interest in this phenomenon, the physical basis of its origin is still poorly understood. To gain a global view of the available information on this phenomenon, we here provide an exhaustive survey of the current literature in which original data on this fluorescence have been reported. The emitting systems have been classified in terms of their molecular complexity, amino acid composition, and physical state. Information about the wavelength of the radiation used for the excitation as well as the emission range/peak has also been retrieved. The data collected here provide a picture of the complexity of this multifaceted phenomenon that could be helpful for future studies aimed at defining its structural and electronic basis and/or stimulating new applications.

A Structure-Based Mechanism for the Denaturing Action of Urea, Guanidinium Ion and Thiocyanate Ion.

An exhaustive analysis of all the protein structures deposited in the Protein Data Bank, here performed, has allowed the identification of hundredths of protein-bound urea molecules and the structural characterization of such binding sites. It emerged that, even though urea molecules are largely involved in hydrogen bonds with both backbone and side chains, they are also able to make van der Waals contacts with nonpolar moieties. As similar findings have also been previously reported for guanidinium and thiocyanate, this observation suggests that promiscuity is a general property of protein denaturants. Present data provide strong support for a mechanism based on the protein-denaturant direct interactions with a denaturant binding model to equal and independent sites. In this general framework, our investigations also highlight some interesting insights into the different denaturing power of urea compared to guanidinium/thiocyanate.

Alphafold Predictions Provide Insights into the Structural Features of the Functional Oligomers of All Members of the KCTD Family.

Oligomerization endows proteins with some key properties such as extra-stabilization, long-range allosteric regulation(s), and partnerships not accessible to their monomeric counterparts. How oligomerization is achieved and preserved during evolution is a subject of remarkable scientific relevance. By exploiting the abilities of the machine-learning algorithms implemented in AlphaFold (AF) in predicting protein structures, herein, we report a comprehensive analysis of the structural states of functional oligomers of all members of the KCTD protein family. Interestingly, our approach led to the identification of reliable three-dimensional models for the pentameric states of KCNRG, KCTD6, KCTD4, KCTD7, KCTD9, and KCTD14 and possibly for KCTD11 and KCTD21 that are involved in key biological processes and that were previously uncharacterized from a structural point of view. Although for most of these proteins, the CTD domains lack any sequence similarity, they share some important structural features, such as a propeller-like structure with a central cavity delimited by five exposed and regular ß-strands. Moreover, the structure of the related proteins KCTD7 and KCTD14, although pentameric, appears to be characterized by a different organization of the CTD region, with the five chains forming a circle-like structure with a large cavity. Our predictions also suggest that other members of the family, such as KCTD10, KCTD13, and TNFAIP1, present a strong propensity to assume dimeric states. Although the structures of the functional oligomers reported herein represent models that require additional validations, they provide a consistent and global view of KCTD protein oligomerization.

Local Backbone Geometry Plays a Critical Role in Determining Conformational Preferences of Amino Acid Residues in Proteins.

The definition of the structural basis of the conformational preferences of the genetically encoded amino acid residues is an important yet unresolved issue of structural biology. In order to gain insights into this intricate topic, we here determined and compared the amino acid propensity scales for different (φ, ψ) regions of the Ramachandran plot and for different secondary structure elements. These propensities were calculated using the Chou-Fasman approach on a database of non-redundant protein chains retrieved from the Protein Data Bank. Similarities between propensity scales were evaluated by linear regression analyses. One of the most striking and unexpected findings is that distant regions of the Ramachandran plot may exhibit significantly similar propensity scales. On the other hand, contiguous regions of the Ramachandran plot may present anticorrelated propensities. In order to provide an interpretative background to these results, we evaluated the role that the local variability of protein backbone geometry plays in this context. Our analysis indicates that (dis)similarities of propensity scales between different regions of the Ramachandran plot are coupled with (dis)similarities in the local geometry. The concept that similarities of the propensity scales are dictated by the similarity of the NCαC angle and not necessarily by the similarity of the (φ, ψ) conformation may have far-reaching implications in the field.

Atomic-Level View of the Functional Transition in Vertebrate Hemoglobins: The Case of Antarctic Fish Hbs.

Tetrameric hemoglobins (Hbs) are prototypal systems for studies aimed at unveiling basic structure-function relationships as well as investigating the molecular/structural basis of adaptation of living organisms to extreme conditions. However, a chronological analysis of decade-long studies conducted on Hbs is illuminating on the difficulties associated with the attempts of gaining functional insights from static structures. Here, we applied molecular dynamics (MD) simulations to explore the functional transition from the T to the R state of the hemoglobin of the Antarctic fish Trematomus bernacchii (HbTb). Our study clearly demonstrates the ability of the MD technique to accurately describe the transition of HbTb from the T to R-like states, as shown by a number of global and local structural indicators. A comparative analysis of the structural states that HbTb assumes in the simulations with those detected in previous MD analyses conducted on HbA (human Hb) highlights interesting analogies (similarity of the transition pathway) and differences (distinct population of intermediate states). In particular, the ability of HbTb to significantly populate intermediate states along the functional pathway explains the observed propensity of this protein to assume these structures in the crystalline state. It also explains some functional data reported on the protein that indicate the occurrence of other functional states in addition to the canonical R and T ones. These findings are in line with the emerging idea that the classical two-state view underlying tetrameric Hb functionality is probably an oversimplification and that other structural states play important roles in these proteins. The ability of MD simulations to accurately describe the functional pathway in tetrameric Hbs suggests that this approach may be effectively applied to unravel the molecular and structural basis of Hbs exhibiting peculiar functional properties as a consequence of the environmental adaptation of the host organism.

A Protein Data Bank survey of multimodal binding of thiocyanate to proteins: Evidence for thiocyanate promiscuity.

Over the last one and half century, a myriad of studies has demonstrated that Hofmeister ions have a major impact on protein stability and solubility. Nevertheless, the definition of the physico-chemical basis of their activity has proved to be highly challenging and controversial. Here, by exploiting the enormous information content of the Protein Data Bank, we explored the binding to proteins of thiocyanate, the anion of the series exerting the highest solubilization/destabilization effects. The survey, which led to the identification and characterization of 712 thiocyanate binding sites, provides a comprehensive and atomic-level view of the varied interactions that the ion forms with proteins. The inspection of these sites highlights a limited tendency of thiocyanate to interact with structured water molecules, in line with the reported poor hydration of the ion. On the other hand, the thiocyanate makes interactions with protein nonpolar moieties, especially with the backbone Cα atom. In as many as 104 cases, the ion exclusively makes nonpolar contacts. In conclusion, these findings suggest that the ability of thiocyanate to bind all types of protein exposed patches may lead to the formation of a negatively charged electrostatic barrier that could prevent protein-protein aggregation and promote protein solubility. Moreover, the denaturing action of thiocyanate may be ascribed to its ability to establish multiple attractive interactions with protein surfaces.

The temporal correlation between positive testing and death in Italy: from the first phase to the later evolution of the COVID-19 pandemic.

BACKGROUND AND AIM: After the global spread of the novel coronavirus disease 2019 (COVID-19), research has concentrated its efforts on several aspects of the epidemiological burden of pandemic. In this frame, the presented study follows a previous analysis of the temporal link between cases and deaths during the first epidemic wave (Phase 1) in Italy (March-June 2020). METHODS: We here analyze the COVID-19 epidemic in the time span from March 2020 to June 2021. RESULTS: The elaboration of the curves of cases and deaths allows identifying the temporal shift between the positive testing and the fatal event, which corresponds to one week from W2 to W33, two weeks from W34 to W41, and three weeks from W42 to W67. Based on this finding, we calculate the Weekly Lethality Rate (WLR). The WLR was grossly overestimated (~13.5%) in Phase 1, while a mean value of 2.6% was observed in most of Phase 2 (starting from October 2020), with a drop to 1.4% in the last investigated weeks. CONCLUSIONS: Overall, these findings offer an interesting insight into the magnitude and time evolution of the lethality burden attributable to COVID-19 during the entire pandemic period in Italy. In particular, the analysis highlighted the impact of the effectiveness of public health and social measures, of changes in disease management, and of preventive strategies over time.  (www.actabiomedica.it).

AlphaFold-Predicted Structures of KCTD Proteins Unravel Previously Undetected Relationships among the Members of the Family.

One of the most striking features of KCTD proteins is their involvement in apparently unrelated yet fundamental physio-pathological processes. Unfortunately, comprehensive structure-function relationships for this protein family have been hampered by the scarcity of the structural data available. This scenario is rapidly changing due to the release of the protein three-dimensional models predicted by AlphaFold (AF). Here, we exploited the structural information contained in the AF database to gain insights into the relationships among the members of the KCTD family with the aim of facilitating the definition of the structural and molecular basis of key roles that these proteins play in many biological processes. The most important finding that emerged from this investigation is the discovery that, in addition to the BTB domain, the vast majority of these proteins also share a structurally similar domain in the C-terminal region despite the absence of general sequence similarities detectable in this region. Using this domain as reference, we generated a novel and comprehensive structure-based pseudo-phylogenetic tree that unraveled previously undetected similarities among the protein family. In particular, we generated a new clustering of the KCTD proteins that will represent a solid ground for interpreting their many functions.

A global analysis of conservative and non-conservative mutations in SARS-CoV-2 detected in the first year of the COVID-19 world-wide diffusion.

The ability of SARS-CoV-2 to rapidly mutate represents a remarkable complicancy. Quantitative evaluations of the effects that these mutations have on the virus structure/function is of great relevance and the availability of a large number of SARS-CoV-2 sequences since the early phases of the pandemic represents a unique opportunity to follow the adaptation of the virus to humans. Here, we evaluated the SARS-CoV-2 amino acid mutations and their progression by analyzing publicly available viral genomes at three stages of the pandemic (2020 March 15th and October 7th, 2021 February 7th). Mutations were classified in conservative and non-conservative based on the probability to be accepted during the evolution according to the Point Accepted Mutation substitution matrices and on the similarity of the encoding codons. We found that the most frequent substitutions are T > I, L > F, and A > V and we observe accumulation of hydrophobic residues. These findings are consistent among the three stages analyzed. We also found that non-conservative mutations are less frequent than conservative ones. This finding may be ascribed to a progressive adaptation of the virus to the host. In conclusion, the present study provides indications of the early evolution of the virus and tools for the global and genome-specific evaluation of the possible impact of mutations on the structure/function of SARS-CoV-2 variants.

Exosite Binding in Thrombin: A Global Structural/Dynamic Overview of Complexes with Aptamers and Other Ligands.

Thrombin is the key enzyme of the entire hemostatic process since it is able to exert both procoagulant and anticoagulant functions; therefore, it represents an attractive target for the developments of biomolecules with therapeutic potential. Thrombin can perform its many functional activities because of its ability to recognize a wide variety of substrates, inhibitors, and cofactors. These molecules frequently are bound to positively charged regions on the surface of protein called exosites. In this review, we carried out extensive analyses of the structural determinants of thrombin partnerships by surveying literature data as well as the structural content of the Protein Data Bank (PDB). In particular, we used the information collected on functional, natural, and synthetic molecular ligands to define the anatomy of the exosites and to quantify the interface area between thrombin and exosite ligands. In this framework, we reviewed in detail the specificity of thrombin binding to aptamers, a class of compounds with intriguing pharmaceutical properties. Although these compounds anchor to protein using conservative patterns on its surface, the present analysis highlights some interesting peculiarities. Moreover, the impact of thrombin binding aptamers in the elucidation of the cross-talk between the two distant exosites is illustrated. Collectively, the data and the work here reviewed may provide insights into the design of novel thrombin inhibitors.

The Introduction of a Cysteine Residue Modulates The Mechanical Properties of Aromatic-Based Solid Aggregates and Self-Supporting Hydrogels.

Peptide-based hydrogels, originated by multiscale self-assembling phenomenon, have been proposed as multivalent tools in different technological areas. Structural studies and molecular dynamics simulations pointed out the capability of completely aromatic peptides to gelificate if hydrophilic and hydrophobic forces are opportunely balanced. Here, the effect produced by the introduction of a Cys residue in the heteroaromatic sequence of (FY)3 and in its PEGylated variant was evaluated. The physicochemical characterization indicates that both FYFCFYF and PEG8-FYFCFYF are able to self-assemble in supramolecular nanostructures whose basic cross-ß motif resembles the one detected in the ancestor (FY)3 assemblies. However, gelification occurs only for FYFCFYF at a concentration of 1.5 wt%. After cross-linking of cysteine residues, the hydrogel undergoes to an improvement of the rigidity compared to the parent (FY)3 assemblies as suggested by the storage modulus (G') that increases from 970 to 3360 Pa. The mechanical properties of FYFCFYF are compatible with its potential application in bone tissue regeneration. Moreover, the avalaibility of a Cys residue in the middle of the peptide sequence could allow the hydrogel derivatization with targeting moieties or with biologically relevant molecules.

Quaternary Structure Transitions of Human Hemoglobin: An Atomic-Level View of the Functional Intermediate States.

Human hemoglobin (HbA) is one of the prototypal systems used to investigate structure-function relationships in proteins. Indeed, HbA has been used to develop the basic concepts of protein allostery, although the atomic-level mechanism underlying the HbA functionality is still highly debated. This is due to the fact that most of the three-dimensional structural information collected over the decades refers to the endpoints of HbA functional transition with little data available for the intermediate states. Here, we report molecular dynamics (MD) simulations by focusing on the relevance of the intermediate states of the protein functional transition unraveled by the crystallographic studies carried out on vertebrate Hbs. Fully atomistic simulations of the HbA T-state indicate that the protein undergoes a spontaneous transition toward the R-state. The inspection of the trajectory structures indicates that the protein significantly populates the intermediate HL-(C) state previously unraveled by crystallography. In the structural transition, it also assumes the intermediate states crystallographically detected in Antarctic fish Hbs. This finding suggests that HbA and Antarctic fish Hbs, in addition to the endpoints of the transitions, also share a similar deoxygenation pathway despite a distace of hundreds of millions of years in the evolution scale. Finally, using the essential dynamic sampling methodology, we gained some insights into the reverse R to T transition that is not spontaneously observed in classic MD simulations.

Analysis of the time evolution of COVID-19 lethality during the first epidemic wave in Italy.

BACKGROUND AND AIM: While the entire world is still experiencing the dramatic emergency due to SARS-CoV-2, Italy has a prominent position since it has been the locus of the first major outbreak among Western countries. The aim of this study is the evaluation of temporal connection between SARS-CoV-2 positive tests (cases) and deaths in Italy in the first wave of the epidemic. METHODS: A temporal link between cases and deaths was determined by comparing their daily/weekly trends using surveillance data of the period March 2-June 2020. RESULTS: The monitoring of the cases/deaths evolution during the first wave of the outbreak highlights a striking correlation between infections of a certain week and deaths of the following one. We defined a weekly lethality rate that is virtually unchanged over the entire months of April and May until the first week of June (≈13.6%). Due to the rather low number of cases/deaths, this parameter starts to fluctuate in the following three weeks. CONCLUSIONS: The analysis indicates that the weekly lethality rate is virtually unchanged over the entire first wave of the epidemic, despite the progressive increase of the testing. As observed for the overall lethality, this parameter uniformly presents rather high values. The definition of a temporal link between cases and deaths will likely represent a useful tool for highlighting analogies and differences between the first and the second wave of the pandemic and for evaluating the effectiveness, even if partial, of the strategies applied during the ongoing outbreak. (www.actabiomedica.it).

Structural and functional analysis of the simultaneous binding of two duplex/quadruplex aptamers to human α-thrombin.

The long-range communication between the two exosites of human α-thrombin (thrombin) tightly modulates the protein-effector interactions. Duplex/quadruplex aptamers represent an emerging class of very effective binders of thrombin. Among them, NU172 and HD22 aptamers are at the forefront of exosite I and II recognition, respectively. The present study investigates the simultaneous binding of these two aptamers by combining a structural and dynamics approach. The crystal structure of the ternary complex formed by the thrombin with NU172 and HD22_27mer provides a detailed view of the simultaneous binding of these aptamers to the protein, inspiring the design of novel bivalent thrombin inhibitors. The crystal structure represents the starting model for molecular dynamics studies, which point out the cooperation between the binding at the two exosites. In particular, the binding of an aptamer to its exosite reduces the intrinsic flexibility of the other exosite, that preferentially assumes conformations similar to those observed in the bound state, suggesting a predisposition to interact with the other aptamer. This behaviour is reflected in a significant increase of the anticoagulant activity of NU172 when the inactive HD22_27mer is bound to exosite II, providing a clear evidence of the synergic action of the two aptamers.