Hunting down zinc(II)-binding sites in proteins with distance matrices.

MOTIVATION: In recent years, high-throughput sequencing technologies have made available the genome sequences of a huge variety of organisms. However, the functional annotation of the encoded proteins often still relies on low-throughput and costly experimental studies. Bioinformatics approaches offer a promising alternative to accelerate this process. In this work, we focus on the binding of zinc(II) ions, which is needed for 5%-10% of any organism's proteins to achieve their physiologically relevant form. RESULTS: To implement a predictor of zinc(II)-binding sites in the 3D structures of proteins, we used a neural network, followed by a filter of the network output against the local structure of all known sites. The latter was implemented as a function comparing the distance matrices of the Cα and Cß atoms of the sites. We called the resulting tool Master of Metals (MOM). The structural models for the entire proteome of an organism generated by AlphaFold can be used as input to our tool in order to achieve annotation at the whole organism level within a few hours. To demonstrate this, we applied MOM to the yeast proteome, obtaining a precision of about 76%, based on data for homologous proteins. AVAILABILITY AND IMPLEMENTATION: Master of Metals has been implemented in Python and is available at https://github.com/cerm-cirmmp/Master-of-metals.

Upgrading of the general AMBER force field 2 for fluorinated alcohol biosolvents: A validation for water solutions and melittin solvation.

The standard GAFF2 force field parameterization has been refined for the fluorinated alcohols 2,2,2-trifluoroethanol (TFE), 1,1,1,3,3,3-hexafluoro-2-propanol (HFIP), and 1,1,1,3,3,3-hexafluoropropan-2-one (HFA), which are commonly used to study proteins and peptides in biomimetic media. The structural and dynamic properties of both proteins and peptides are significantly influenced by the biomimetic environment created by the presence of these cosolvents in aqueous solutions. Quantum mechanical calculations on stable conformers were used to parameterize the atomic charges. Different systems, such as pure liquids, aqueous solutions, and systems formed by melittin protein and cosolvent/water solutions, have been used to validate the new models. The calculated macroscopic and structural properties are in agreement with experimental findings, supporting the validity of the newly proposed models.

A New Paradigm of Multiheme Cytochrome Evolution by Grafting and Pruning Protein Modules.

Multiheme cytochromes play key roles in diverse biogeochemical cycles, but understanding the origin and evolution of these proteins is a challenge due to their ancient origin and complex structure. Up until now, the evolution of multiheme cytochromes composed by multiple redox modules in a single polypeptide chain was proposed to occur by gene fusion events. In this context, the pentaheme nitrite reductase NrfA and the tetraheme cytochrome c554 were previously proposed to be at the origin of the extant octa- and nonaheme cytochrome c involved in metabolic pathways that contribute to the nitrogen, sulfur, and iron biogeochemical cycles by a gene fusion event. Here, we combine structural and character-based phylogenetic analysis with an unbiased root placement method to refine the evolutionary relationships between these multiheme cytochromes. The evidence show that NrfA and cytochrome c554 belong to different clades, which suggests that these two multiheme cytochromes are products of truncation of ancestral octaheme cytochromes related to extant octaheme nitrite reductase and MccA, respectively. From our phylogenetic analysis, the last common ancestor is predicted to be an octaheme cytochrome with nitrite reduction ability. Evolution from this octaheme framework led to the great diversity of extant multiheme cytochromes analyzed here by pruning and grafting of protein modules and hemes. By shedding light into the evolution of multiheme cytochromes that intervene in different biogeochemical cycles, this work contributes to our understanding about the interplay between biology and geochemistry across large time scales in the history of Earth.

Ruthenium(II) Polypyridyl Complexes and Metronidazole Derivatives: A Powerful Combination in the Design of Photoresponsive Antibacterial Agents Effective under Hypoxic Conditions.

Ruthenium(II) polypyridyl complexes (RPCs) are gaining momentum in photoactivated chemotherapy (PACT), thanks to the possibility of overcoming the classical reliance on molecular oxygen of photodynamic therapy while preserving the selective drug activation by using light. However, notwithstanding the intriguing perspectives, the translation of such an approach in the development of new antimicrobials has been only barely considered. Herein, MTZH-1 and MTZH-2, two novel analogues of metronidazole (MTZ), a mainstay drug in the treatment of anaerobic bacterial infections, were designed and inserted in the strained ruthenium complexes [Ru(tpy)(dmp)(MTZ-1)]PF6 (Ru2) and [Ru(tpy)(dmp)(MTZ-2)]PF6 (Ru3) (tpy = terpyridine, dmp = 2,9-dimethyl-1,10-phenanthroline) (Chart 1). Analogously to the parental compound [Ru(tpy)(dmp)(5NIM)]PF6 (Ru1) (5-nitroimidazolate), the Ru(II)-imidazolate coordination of MTZ derivatives resulted in promising Ru(II) photocages, capable to easily unleash the bioactive ligands upon light irradiation and increase the antibacterial activity against Bacillus subtilis, which was chosen as a model of Gram-positive bacteria. The photoreleased 5-nitroimidazole-based ligands led to remarkable phototoxicities under hypoxic conditions (<1% O2), with the lead compound Ru3 that exhibited the highest potency across the series, being comparable to the one of the clinical drug MTZ. Besides, the chemical architectures of MTZ derivatives made their interaction with NimAunfavorable, being NimA a model of reductases responsible for bacterial resistance against 5-nitroimidazole-based antibiotics, thus hinting at their possible use to combat antimicrobial resistance. This work may therefore provide fundamental knowledge in the design of novel photoresponsive tools to be used in the fight against infectious diseases. For the first time, the effectiveness of the "photorelease antimicrobial therapy" under therapeutically relevant hypoxic conditions was demonstrated.

Learning to Identify Physiological and Adventitious Metal-Binding Sites in the Three-Dimensional Structures of Proteins by Following the Hints of a Deep Neural Network.

Thirty-eight percent of protein structures in the Protein Data Bank contain at least one metal ion. However, not all these metal sites are biologically relevant. Cations present as impurities during sample preparation or in the crystallization buffer can cause the formation of protein-metal complexes that do not exist in vivo. We implemented a deep learning approach to build a classifier able to distinguish between physiological and adventitious zinc-binding sites in the 3D structures of metalloproteins. We trained the classifier using manually annotated sites extracted from the MetalPDB database. Using a 10-fold cross validation procedure, the classifier achieved an accuracy of about 90%. The same neural classifier could predict the physiological relevance of non-heme mononuclear iron sites with an accuracy of nearly 80%, suggesting that the rules learned on zinc sites have general relevance. By quantifying the relative importance of the features describing the input zinc sites from the network perspective and by analyzing the characteristics of the MetalPDB datasets, we inferred some common principles. Physiological sites present a low solvent accessibility of the aminoacids forming coordination bonds with the metal ion (the metal ligands), a relatively large number of residues in the metal environment (≥20), and a distinct pattern of conservation of Cys and His residues in the site. Adventitious sites, on the other hand, tend to have a low number of donor atoms from the polypeptide chain (often one or two). These observations support the evaluation of the physiological relevance of novel metal-binding sites in protein structures.

Structural Bioinformatics and Deep Learning of Metalloproteins: Recent Advances and Applications.

All living organisms require metal ions for their energy production and metabolic and biosynthetic processes. Within cells, the metal ions involved in the formation of adducts interact with metabolites and macromolecules (proteins and nucleic acids). The proteins that require binding to one or more metal ions in order to be able to carry out their physiological function are called metalloproteins. About one third of all protein structures in the Protein Data Bank involve metalloproteins. Over the past few years there has been tremendous progress in the number of computational tools and techniques making use of 3D structural information to support the investigation of metalloproteins. This trend has been boosted by the successful applications of neural networks and machine/deep learning approaches in molecular and structural biology at large. In this review, we discuss recent advances in the development and availability of resources dealing with metalloproteins from a structure-based perspective. We start by addressing tools for the prediction of metal-binding sites (MBSs) using structural information on apo-proteins. Then, we provide an overview of the methods for and lessons learned from the structural comparison of MBSs in a fold-independent manner. We then move to describing databases of metalloprotein/MBS structures. Finally, we summarizing recent ML/DL applications enhancing the functional interpretation of metalloprotein structures.

MetalPDB in 2018: a database of metal sites in biological macromolecular structures.

MetalPDB (http://metalweb.cerm.unifi.it/) is a database providing information on metal-binding sites detected in the three-dimensional (3D) structures of biological macromolecules. MetalPDB represents such sites as 3D templates, called Minimal Functional Sites (MFSs), which describe the local environment around the metal(s) independently of the larger context of the macromolecular structure. The 2018 update of MetalPDB includes new contents and tools. A major extension is the inclusion of proteins whose structures do not contain metal ions although their sequences potentially contain a known MFS. In addition, MetalPDB now provides extensive statistical analyses addressing several aspects of general metal usage within the PDB, across protein families and in catalysis. Users can also query MetalPDB to extract statistical information on structural aspects associated with individual metals, such as preferred coordination geometries or aminoacidic environment. A further major improvement is the functional annotation of MFSs; the annotation is manually performed via a password-protected annotator interface. At present, ∼50% of all MFSs have such a functional annotation. Other noteworthy improvements are bulk query functionality, through the upload of a list of PDB identifiers, and ftp access to MetalPDB contents, allowing users to carry out in-depth analyses on their own computational infrastructure.

HIV-1 Tat Protein Enters Dysfunctional Endothelial Cells via Integrins and Renders Them Permissive to Virus Replication.

Previous work has shown that the Tat protein of Human Immunodeficiency Virus (HIV)-1 is released by acutely infected cells in a biologically active form and enters dendritic cells upon the binding of its arginine-glycine-aspartic acid (RGD) domain to the α5ß1, αvß3, and αvß5 integrins. The up-regulation/activation of these integrins occurs in endothelial cells exposed to inflammatory cytokines that are increased in HIV-infected individuals, leading to endothelial cell dysfunction. Here, we show that inflammatory cytokine-activated endothelial cells selectively bind and rapidly take up nano-micromolar concentrations of Tat, as determined by flow cytometry. Protein oxidation and low temperatures reduce Tat entry, suggesting a conformation- and energy-dependent process. Consistently, Tat entry is competed out by RGD-Tat peptides or integrin natural ligands, and it is blocked by anti-α5ß1, -αvß3, and -αvß5 antibodies. Moreover, modelling-docking calculations identify a low-energy Tat-αvß3 integrin complex in which Tat makes contacts with both the αv and ß3 chains. It is noteworthy that internalized Tat induces HIV replication in inflammatory cytokine-treated, but not untreated, endothelial cells. Thus, endothelial cell dysfunction driven by inflammatory cytokines renders the vascular system a target of Tat, which makes endothelial cells permissive to HIV replication, adding a further layer of complexity to functionally cure and/or eradicate HIV infection.

Upgrading and Validation of the AMBER Force Field for Histidine and Cysteine Zinc(II)-Binding Residues in Sites with Four Protein Ligands.

We developed and validated a novel force field in the context of the AMBER parameterization for the simulation of zinc(II)-binding proteins. The proposed force field assumes nonbonded spherical interactions between the central zinc(II) and the coordinating residues. A crucial innovative aspect of our approach is to account for the polarization effects of the cation by redefining the atomic charges of the coordinating residues and an adjustment of Lennard-Jones parameters of Zn-interacting atoms to reproduce mean distance distributions. The optimal transferable parametrization was obtained by performing accurate quantum mechanical calculations on a training set of high-quality protein structures, encompassing the most common folds of zinc(II) sites. The addressed sites contain a zinc(II) ion tetra-coordinated by histidine and cysteine residues and represent about 70% of all physiologically relevant zinc(II) sites in the Protein Data Bank. Molecular dynamics simulations with explicit solvent, carried out on several zinc(II)-binding proteins not included in the training set, show that our model for zinc(II) sites preserves the tetra-coordination of the metal site with remarkable stability, yielding zinc(II)-X mean distances similar to experimental data. Finally, the model was tested by evaluating the zinc(II)-binding affinities, using the alchemical free energy perturbation approach. The calculated dissociation constants correlate satisfactorily with the experimental counterpart demonstrating the validity and transferability of the proposed parameterization for zinc(II)-binding proteins.

MetalPredator: a web server to predict iron-sulfur cluster binding proteomes.

MOTIVATION: The prediction of the iron-sulfur proteome is highly desirable for biomedical and biological research but a freely available tool to predict iron-sulfur proteins has not been developed yet. RESULTS: We developed a web server to predict iron-sulfur proteins from protein sequence(s). This tool, called MetalPredator, is able to process complete proteomes rapidly with high recall and precision. AVAILABILITY AND IMPLEMENTATION: The web server is freely available at: http://metalweb.cerm.unifi.it/tools/metalpredator/ CONTACT: andreini@cerm.unifi.it SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Identification of the zinc, copper and cadmium metalloproteome of the protozoon Tetrahymena thermophila by systematic bioinformatics.

Tetrahymena thermophila (T. thermophila) is a ciliated protozoon that can detect freshwater pollution by heavy metals ("whole-cell biosensor"). This work employed a systematic bioinformatics approach to predict and analyze the metalloproteome of T. thermophila for the essential Zn, Cu and the non-essential Cd. 3784 metal-binding proteins were identified compared to the 456 annotated so far in UniProt. The localization, functional classification, and the functionally enriched network of the newly identified metalloproteome are presented. Cd toxicity could be explained in terms of the metal replacing Cu and especially Zn in MAPKs, transporters and antioxidant enzymes. The predicted results for Cd toxicity and responses reflect those observed experimentally in different organisms after their exposure to Cd.

Exploiting Bacterial Operons To Illuminate Human Iron-Sulfur Proteins.

Organisms from all kingdoms of life use iron-sulfur proteins (FeS-Ps) in a multitude of functional processes. We applied a bioinformatics approach to investigate the human portfolio of FeS-Ps. Sixty-one percent of human FeS-Ps bind Fe4S4 clusters, whereas 39% bind Fe2S2 clusters. However, this relative ratio varies significantly depending on the specific cellular compartment. We compared the portfolio of human FeS-Ps to 12 other eukaryotes and to about 700 prokaryotes. The comparative analysis of the organization of the prokaryotic homologues of human FeS-Ps within operons allowed us to reconstruct the human functional networks involving the conserved FeS-Ps common to prokaryotes and eukaryotes. These functional networks have been maintained during evolution and thus presumably represent fundamental cellular processes. The respiratory chain and the ISC machinery for FeS-P biogenesis are the two conserved processes that involve the majority of human FeS-Ps. Purine metabolism is another process including several FeS-Ps, in which BOLA proteins possibly have a regulatory role. The analysis of the co-occurrence of human FeS-Ps with other proteins highlighted numerous links between the iron-sulfur cluster machinery and the response mechanisms to cell damage, from repair to apoptosis. This relationship probably relates to the production of reactive oxygen species within the biogenesis and degradation of FeS-Ps.

Minimal Functional Sites in Metalloproteins and Their Usage in Structural Bioinformatics.

Metal ions play a functional role in numerous biochemical processes and cellular pathways. Indeed, about 40% of all enzymes of known 3D structure require a metal ion to be able to perform catalysis. The interactions of the metals with the macromolecular framework determine their chemical properties and reactivity. The relevant interactions involve both the coordination sphere of the metal ion and the more distant interactions of the so-called second sphere, i.e., the non-bonded interactions between the macromolecule and the residues coordinating the metal (metal ligands). The metal ligands and the residues in their close spatial proximity define what we call a minimal functional site (MFS). MFSs can be automatically extracted from the 3D structures of metal-binding biological macromolecules deposited in the Protein Data Bank (PDB). They are 3D templates that describe the local environment around a metal ion or metal cofactor and do not depend on the overall macromolecular structure. MFSs provide a different view on metal-binding proteins and nucleic acids, completely focused on the metal. Here we present different protocols and tools based upon the concept of MFS to obtain deeper insight into the structural and functional properties of metal-binding macromolecules. We also show that structure conservation of MFSs in metalloproteins relates to local sequence similarity more strongly than to overall protein similarity.

MetalPDB: a database of metal sites in biological macromolecular structures.

We present here MetalPDB (freely accessible at http://metalweb.cerm.unifi.it), a novel resource aimed at conveying the information available on the three-dimensional (3D) structures of metal-binding biological macromolecules in a consistent and effective manner. This is achieved through the systematic and automated representation of metal-binding sites in proteins and nucleic acids by way of Minimal Functional Sites (MFSs). MFSs are 3D templates that describe the local environment around the metal(s) independently of the larger context of the macromolecular structure embedding the site(s), and are the central objects of MetalPDB design. MFSs are grouped into equistructural (broadly defined as sites found in corresponding positions in similar structures) and equivalent sites (equistructural sites that contain the same metals), allowing users to easily analyse similarities and variations in metal-macromolecule interactions, and to link them to functional information. The web interface of MetalPDB allows access to a comprehensive overview of metal-containing biological structures, providing a basis to investigate the basic principles governing the properties of these systems. MetalPDB is updated monthly in an automated manner.

MetalS(3), a database-mining tool for the identification of structurally similar metal sites.

We have developed a database search tool to identify metal sites having structural similarity to a query metal site structure within the MetalPDB database of minimal functional sites (MFSs) contained in metal-binding biological macromolecules. MFSs describe the local environment around the metal(s) independently of the larger context of the macromolecular structure. Such a local environment has a determinant role in tuning the chemical reactivity of the metal, ultimately contributing to the functional properties of the whole system. The database search tool, which we called MetalS(3) (Metal Sites Similarity Search), can be accessed through a Web interface at http://metalweb.cerm.unifi.it/tools/metals3/ . MetalS(3) uses a suitably adapted version of an algorithm that we previously developed to systematically compare the structure of the query metal site with each MFS in MetalPDB. For each MFS, the best superposition is kept. All these superpositions are then ranked according to the MetalS(3) scoring function and are presented to the user in tabular form. The user can interact with the output Web page to visualize the structural alignment or the sequence alignment derived from it. Options to filter the results are available. Test calculations show that the MetalS(3) output correlates well with expectations from protein homology considerations. Furthermore, we describe some usage scenarios that highlight the usefulness of MetalS(3) to obtain mechanistic and functional hints regardless of homology.

MACiE: exploring the diversity of biochemical reactions.

MACiE (which stands for Mechanism, Annotation and Classification in Enzymes) is a database of enzyme reaction mechanisms, and can be accessed from http://www.ebi.ac.uk/thornton-srv/databases/MACiE/. This article presents the release of Version 3 of MACiE, which not only extends the dataset to 335 entries, covering 182 of the EC sub-subclasses with a crystal structure available (~90%), but also incorporates greater chemical and structural detail. This version of MACiE represents a shift in emphasis for new entries, from non-homologous representatives covering EC reaction space to enzymes with mechanisms of interest to our users and collaborators with a view to exploring the chemical diversity of life. We present new tools for exploring the data in MACiE and comparing entries as well as new analyses of the data and new searches, many of which can now be accessed via dedicated Perl scripts.

A database overview of metal-coordination distances in metalloproteins.

Metalloproteins are ubiquitous in all living organisms and take part in a very wide range of biological processes. For this reason, their experimental characterization is crucial to obtain improved knowledge of their structure and biological functions. The three-dimensional structure represents highly relevant information since it provides insight into the interaction between the metal ion(s) and the protein fold. Such interactions determine the chemical reactivity of the bound metal. The available PDB structures can contain errors due to experimental factors such as poor resolution and radiation damage. A lack of use of distance restraints during the refinement and validation process also impacts the structure quality. Here, the aim was to obtain a thorough overview of the distribution of the distances between metal ions and their donor atoms through the statistical analysis of a data set based on more than 115 000 metal-binding sites in proteins. This analysis not only produced reference data that can be used by experimentalists to support the structure-determination process, for example as refinement restraints, but also resulted in an improved insight into how protein coordination occurs for different metals and the nature of their binding interactions. In particular, the features of carboxylate coordination were inspected, which is the only type of interaction that is commonly present for nearly all metals.

FindGeo: a tool for determining metal coordination geometry.

UNLABELLED: Metals are essential for the structure and function of many proteins and nucleic acids. The geometrical arrangement of the atoms that coordinate a metal in a biological macromolecule is an important determinant of the specificity and role of that metal. At present, however, this information can be retrieved only from the literature, which sometimes contains an improper or incorrect description of the geometry, and often lacks it altogether. Thus, we developed FindGeo to quickly and easily determine the coordination geometry of selected, or all, metals in a given structure. FindGeo works by superimposing the metal-coordinating atoms in the input structure to a library of templates with alternative ideal geometries, which are ranked by RMSD to identify the best geometry assignment. AVAILABILITY: FindGeo is freely available as a web service and as a stand-alone program at http://metalweb.cerm.unifi.it/tools/findgeo/.

MetalS2: a tool for the structural alignment of minimal functional sites in metal-binding proteins and nucleic acids.

We developed a new software tool, MetalS(2), for the structural alignment of Minimal Functional Sites (MFSs) in metal-binding biological macromolecules. MFSs are 3D templates that describe the local environment around the metal(s) independently of the larger context of the macromolecular structure. Such local environment has a determinant role in tuning the chemical reactivity of the metal, ultimately contributing to the functional properties of the whole system. On our example data sets, MetalS(2) unveiled structural similarities that other programs for protein structure comparison do not consistently point out and overall identified a larger number of structurally similar MFSs. MetalS(2) supports the comparison of MFSs harboring different metals and/or with different nuclearity and is available both as a stand-alone program and a Web tool ( http://metalweb.cerm.unifi.it/tools/metals2/).

Metal-induced structural variability of mononuclear metal-binding sites from a database perspective.

Metalloproteins are ubiquitous in all kingdoms of life. Their role and function are tightly related to the local structure of the metal-binding site. In this regard, the MetalPDB database is an invaluable tool since it stores the 3D structure of metal-binding sites and of their corresponding apo forms. In this work, we exploited MetalPDB to compute extensive statistics over >3000 clusters of mononuclear sites about the rearrangements occurring upon change in metalation state. For each cluster, we matched the holo and apo sites so that it was possible to average the distances between all possible pairs of Cα and donor atoms and thus quantitatively assess structural variations by computing the Δ values (mean apo distance - mean holo distance). For most of the structures the backbone is rigid with little to no rearrangement, while donor atoms experience significant changes of their relative position when the metal is removed. Sodium and potassium sites are an exception to this general observation. This is most likely caused by their preference for coordination by the main-chain oxygen atoms, making the rearrangement of donor atoms superimposable to that of the backbone. Magnesium and calcium show a different behavior, despite their chemical similarity: calcium sites undergo a larger reorganization upon metalation although both metals have similar percentage of backbone oxygen as donor atoms. We ascribe this observation to the structural and energetic factors regulating the selectivity for calcium over magnesium.