Evaluation of potential miRNA sponge effects of SARS genomes in human.

To date the coronavirus family is composed of seven different viruses which were commonly known as cold viruses until the appearance of the severe acute respiratory coronavirus (SARS-CoV) in 2002, the middle east respiratory syndrome coronavirus (MERS) in 2012 and the severe acute respiratory coronavirus 2 (SARS-CoV-2) which caused the COVID-19 global pandemic in 2019. Using bioinformatic approaches we tested the potential interactions of human miRNAs, expressed in pulmonary epithelial cells, with the available coronavirus genomes. Putative miRNA binding sites were then compared between pathogenic and non pathogenic virus groups. The pathogenic group shares 6 miRNA binding sites that can be potentially involved in the sequestration of miRNAs already known to be associated with deep vein thrombosis. We then analysed ∼100k SARS-CoV-2 variant genomes for their potential interaction with human miRNAs and this study highlighted a group of 97 miRNA binding sites which is present in all the analysed genomes. Among these, we identified 6 miRNA binding sites specific for SARS-CoV-2 and the other two pathogenic viruses whose down-regulation has been seen associated with deep vein thrombosis and cardiovascular diseases. Interestingly, one of these miRNAs, namely miR-20a-5p, whose expression decreases with advancing age, is involved in cytokine signaling, cell differentiation and/or proliferation. We hypothesize that depletion of poorly expressed miRNA could be related with disease severity.

A novel structure-based encoding for machine-learning applied to the inference of SH3 domain specificity.

MOTIVATION: Unravelling the rules underlying protein-protein and protein-ligand interactions is a crucial step in understanding cell machinery. Peptide recognition modules (PRMs) are globular protein domains which focus their binding targets on short protein sequences and play a key role in the frame of protein-protein interactions. High-throughput techniques permit the whole proteome scanning of each domain, but they are characterized by a high incidence of false positives. In this context, there is a pressing need for the development of in silico experiments to validate experimental results and of computational tools for the inference of domain-peptide interactions. RESULTS: We focused on the SH3 domain family and developed a machine-learning approach for inferring interaction specificity. SH3 domains are well-studied PRMs which typically bind proline-rich short sequences characterized by the PxxP consensus. The binding information is known to be held in the conformation of the domain surface and in the short sequence of the peptide. Our method relies on interaction data from high-throughput techniques and benefits from the integration of sequence and structure data of the interacting partners. Here, we propose a novel encoding technique aimed at representing binding information on the basis of the domain-peptide contact residues in complexes of known structure. Remarkably, the new encoding requires few variables to represent an interaction, thus avoiding the 'curse of dimension'. Our results display an accuracy >90% in detecting new binders of known SH3 domains, thus outperforming neural models on standard binary encodings, profile methods and recent statistical predictors. The method, moreover, shows a generalization capability, inferring specificity of unknown SH3 domains displaying some degree of similarity with the known data.

Three-dimensional profiles: a new tool to identify protein surface similarities.

We report a procedure for the description and comparison of protein surfaces, which is based on a three-dimensional (3D) transposition of the profile method for sensitive protein homology sequence searches. Although the principle of the method can be applied to detect similarities to a single protein surface, the possibility of extending this approach to protein families displaying common structural and/or functional properties, makes it a more powerful tool. In analogy to profiles derived from the multiple alignment of protein sequences, we derive a 3D surface profile from a protein structure or from a multiple structure alignment of several proteins. The 3D profile is used to screen the protein structure database, searching for similar protein surfaces. The application of the procedure to SH2 and SH3 binding pockets and to the nucleotide binding pocket associated with the p-loop structural motif is described. The SH2 and SH3 3D profiles can identify all the SH2 and SH3 binding regions present in the test dataset; the p-loop 3D profile is able to recognize all the p-loop-containing proteins present in the test dataset. Analysis of the p-loop 3D profile allowed the identification of a positive charge whose position is conserved in space but not in sequence. The best ranking non-p-loop-containing protein is an ADP-forming succinyl coenzyme A synthetase, whose nucleotide-binding region has not yet been identified.

Correlated mutations contain information about protein-protein interaction.

Many proteins have evolved to form specific molecular complexes and the specificity of this interaction is essential for their function. The network of the necessary inter-residue contacts must consequently constrain the protein sequences to some extent. In other words, the sequence of an interacting protein must reflect the consequence of this process of adaptation. It is reasonable to assume that the sequence changes accumulated during the evolution of one of the interacting proteins must be compensated by changes in the other. Here we apply a method for detecting correlated changes in multiple sequence alignments to a set of interacting protein domains and show that positions where changes occur in a correlated fashion in the two interacting molecules tend to be close to the protein-protein interfaces. This leads to the possibility of developing a method for predicting contacting pairs of residues from the sequence alone. Such a method would not need the knowledge of the structure of the interacting proteins, and hence would be both radically different and more widely applicable than traditional docking methods. We indeed demonstrate here that the information about correlated sequence changes is sufficient to single out the right inter-domain docking solution amongst many wrong alternatives of two-domain proteins. The same approach is also used here in one case (haemoglobin) where we attempt to predict the interface of two different proteins rather than two protein domains. Finally, we report here a prediction about the inter-domain contact regions of the heat- shock protein Hsc70 based only on sequence information.

ESCHER: a new docking procedure applied to the reconstruction of protein tertiary structure.

Evaluation of Surface Complementarity, Hydrogen bonding, and Electrostatic interaction in molecular Recognition (ESCHER) is a new docking procedure consisting of three modules that work in series. The first module evaluates the geometric complementarity and produces a set of rough solutions for the docking problem. The second module identifies molecular collisions within those solutions, and the third evaluates their electrostatic complementarity. We describe the algorithm and its application to the docking of cocrystallized protein domains and unbound components of protein-protein complexes. Furthermore, ESCHER has been applied to the reassociation of secondary and supersecondary structure elements. The possibility of applying a docking method to the problem of protein structure prediction is discussed.