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

Databases of ligand-binding pockets and protein-ligand interactions

Many research groups and institutions have created a variety of databases curating experimental and predicted data related to protein-ligand binding. The landscape of available databases is dynamic, with new databases emerging and established databases becoming defunct. Here, we review the current state of databases that contain binding pockets and protein-ligand binding interactions. We have compiled a list of such databases, fifty-three of which are currently available for use. We discuss variation in how binding pockets are defined and summarize pocket-finding methods. We organize the fifty-three databases into subgroups based on goals and contents, and describe standard use cases. We also illustrate that pockets within the same protein are characterized differently across different databases. Finally, we assess critical issues of sustainability, accessibility and redundancy.

Copper-binding proteins and exonic splicing enhancers and silencers

Eukaryotic DNA codes not only for proteins but contains a wealth of information required for accurate splicing of messenger RNA precursors and inclusion of constitutively or alternatively spliced exons in mature transcripts. This "auxiliary" splicing code has been characterized as exonic splicing enhancers and silencers (ESE and ESS). The exact interplay between protein and splicing codes is, however, poorly understood. Here, we show that exons encoding copper-coordinating amino acids in human cuproproteins lack ESEs and/or have an excess of ESSs, yet RNA sequencing and expressed sequence tags data show that they are more efficiently included in mature transcripts by the splicing machinery than average exons. Their largely constitutive inclusion in messenger RNA is facilitated by stronger splice sites, including polypyrimidine tracts, consistent with an important role of the surrounding intron architecture in ensuring high expression of metal-binding residues during evolution. ESE/ESS profiles of codons and entire exons that code for copper-coordinating residues were very similar to those encoding residues that coordinate zinc but markedly different from those that coordinate calcium. Together, these results reveal how the traditional and auxiliary splicing motifs responded to constraints of metal coordination in proteins.

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.

MetalHawk: Enhanced Classification of Metal Coordination Geometries by Artificial Neural Networks

The chemical properties of metal complexes are strongly dependent on the number and geometrical arrangement of ligands coordinated to the metal center. Existing methods for determining either coordination number or geometry rely on a trade-off between accuracy and computational costs, which hinders their application to the study of large structure data sets. Here, we propose MetalHawk (https://github.com/vrettasm/MetalHawk), a machine learning-based approach to perform simultaneous classification of metal site coordination number and geometry through artificial neural networks (ANNs), which were trained using the Cambridge Structural Database (CSD) and Metal Protein Data Bank (MetalPDB). We demonstrate that the CSD-trained model can be used to classify sites belonging to the most common coordination numbers and geometry classes with balanced accuracy equal to 96.51% for CSD-deposited metal sites. The CSD-trained model was also found to be capable of classifying bioinorganic metal sites from the MetalPDB database, with balanced accuracy equal to 84.29% on the whole PDB data set and to 91.66% on manually reviewed sites in the PDB validation set. Moreover, we report evidence that the output vectors of the CSD-trained model can be considered as a proxy indicator of metal-site distortions, showing that these can be interpreted as a low-dimensional representation of subtle geometrical features present in metal site structures.

The Peptide Antibiotic Corramycin Adopts a β-Hairpin-like Structure and Is Inactivated by the Kinase ComG

The rapid development of antibiotic resistance, especially among difficult-to-treat Gram-negative bacteria, is recognized as a serious and urgent threat to public health. The detection and characterization of novel resistance mechanisms are essential to better predict the spread and evolution of antibiotic resistance. Corramycin is a novel and modified peptidic antibiotic with activity against several Gram-negative pathogens. We demonstrate that the kinase ComG, part of the corramycin biosynthetic gene cluster, phosphorylates and thereby inactivates corramycin, leading to the resistance of the host. Remarkably, we found that the closest structural homologues of ComG are aminoglycoside phosphotransferases; however, ComG shows no activity toward this class of antibiotics. The crystal structure of ComG in complex with corramycin reveals that corramycin adopts a β-hairpin-like structure and allowed us to define the changes leading to a switch in substrate from sugar to peptide. Bioinformatic analyses suggest a limited occurrence of ComG-like proteins, which along with the absence of cross-resistance to clinically used drugs positions corramycin as an attractive antibiotic for further development.

On the hunt for metalloenzyme inhibitors: Investigating the presence of metal-coordinating compounds in screening libraries and chemical spaces

Metalloenzymes play vital roles in various biological processes, requiring the search for inhibitors to develop treatment options for diverse diseases. While compound library screening is a conventional approach, the exploration of virtual chemical spaces housing trillions of compounds has emerged as an alternative strategy. In this study, we investigated the suitability of selected screening libraries and chemical spaces for discovering inhibitors of metalloenzymes featuring common ions (Mg2+, Mn2+, and Zn2+). First, metal-coordinating groups from ligands interacting with ions in the Protein Data Bank were extracted. Subsequently, the prevalence of these groups in two focused screening libraries (Life Chemicals' chelator library, comprising 6,428 compounds, and Otava's chelator fragment library, with 1,784 fragments) as well as two chemical spaces (GalaXi and REAL space, containing billions of virtual products) was investigated. In total, 1,223 metal-coordinating groups were identified, with about a quarter of these groups found within the examined libraries and spaces. Our results indicate that these can serve as valuable starting points for drug discovery targeting metalloenzymes. In addition, this study suggests ways to improve libraries and spaces for better success in finding potential inhibitors for metalloenzymes.

A comprehensive review of protein-centric predictors for biomolecular interactions: from proteins to nucleic acids and beyond

Proteins interact with diverse ligands to perform a large number of biological functions, such as gene expression and signal transduction. Accurate identification of these protein-ligand interactions is crucial to the understanding of molecular mechanisms and the development of new drugs. However, traditional biological experiments are time-consuming and expensive. With the development of high-throughput technologies, an increasing amount of protein data is available. In the past decades, many computational methods have been developed to predict protein-ligand interactions. Here, we review a comprehensive set of over 160 protein-ligand interaction predictors, which cover protein-protein, protein-nucleic acid, protein-peptide and protein-other ligands (nucleotide, heme, ion) interactions. We have carried out a comprehensive analysis of the above four types of predictors from several significant perspectives, including their inputs, feature profiles, models, availability, etc. The current methods primarily rely on protein sequences, especially utilizing evolutionary information. The significant improvement in predictions is attributed to deep learning methods. Additionally, sequence-based pretrained models and structure-based approaches are emerging as new trends.

Early transcriptional changes of heavy metal resistance and multiple efflux genes in Xanthomonas campestris pv. campestris under copper and heavy metal ion stress

BACKGROUND

Copper-induced gene expression in Xanthomonas campestris pv. campestris (Xcc) is typically evaluated using targeted approaches involving qPCR. The global response to copper stress in Xcc and resistance to metal induced damage is not well understood. However, homologs of heavy metal efflux genes from the related Stenotrophomonas genus are found in Xanthomonas which suggests that metal related efflux may also be present.

METHODS AND RESULTS

Gene expression in Xcc strain BrA1 exposed to 0.8 mM CuSO4.5H2O for 15 minutes was captured using RNA-seq analysis. Changes in expression was noted for genes related to general stress responses and oxidoreductases, biofilm formation, protein folding chaperones, heat-shock proteins, membrane lipid profile, multiple drug and efflux (MDR) transporters, and DNA repair were documented. At this timepoint only the cohL (copper homeostasis/tolerance) gene was upregulated as well as a chromosomal czcCBA efflux operon. An additional screen up to 4 hrs using qPCR was conducted using a wider range of heavy metals. Target genes included a cop-containing heavy metal resistance island and putative metal efflux genes. Several efflux pumps, including a copper resistance associated homolog from S. maltophilia, were upregulated under toxic copper stress. However, these pumps were also upregulated in response to other toxic heavy metals. Additionally, the temporal expression of the coh and cop operons was also observed, demonstrating co-expression of tolerance responses and later activation of part of the cop operon.

CONCLUSIONS

Overall, initial transcriptional responses focused on combating oxidative stress, mitigating protein damage and potentially increasing resistance to heavy metals and other biocides. A putative copper responsive efflux gene and others which might play a role in broader heavy metal resistance were also identified. Furthermore, the expression patterns of the cop operon in conjunction with other copper responsive genes allowed for a better understanding of the fate of copper ions in Xanthomonas. This work provides useful evidence for further evaluating MDR and other efflux pumps in metal-specific homeostasis and tolerance phenotypes in the Xanthomonas genus. Furthermore, non-canonical copper tolerance and resistance efflux pumps were potentially identified. These findings have implications for interpreting MIC differences among strains with homologous copLAB resistance genes, understanding survival under copper stress, and resistance in disease management.

Protein3D: Enabling analysis and extraction of metal-containing sites from the Protein Data Bank with molSimplify

Metalloenzymes catalyze a wide range of chemical transformations, with the active site residues playing a key role in modulating chemical reactivity and selectivity. Unlike smaller synthetic catalysts, a metalloenzyme active site is embedded in a larger protein, which makes interrogation of electronic properties and geometric features with quantum mechanical calculations challenging. Here we implement the ability to fetch crystallographic structures from the Protein Data Bank and analyze the metal binding sites in the program molSimplify. We show the usefulness of the newly created protein3D class to extract the local environment around non-heme iron enzymes containing a two histidine motif and prepare 372 structures for quantum mechanical calculations. Our implementation of protein3D serves to expand the range of systems molSimplify can be used to analyze and will enable high-throughput study of metal-containing active sites in proteins.

Differences in heavy metal binding to cysteine-containing coiled-coil peptides

One third of all structurally characterised proteins contain a metal; however, the interplay between metal-binding and peptide/protein folding has yet to be fully elucidated. To better understand how metal binding affects peptide folding, a range of metals should be studied within a specific scaffold. To this end, we modified a histidine-containing coiled-coil peptide to create a cysteine-containing scaffold, named CX3C, which was designed to bind heavy metal ions. In addition, we generated a peptide named CX2C, which contains a binding site more commonly found in natural proteins. Using a combination of analytical techniques including circular dichroism (CD) spectroscopy, UV-Vis spectroscopy and size-exclusion chromatography coupled to multi-angle light scattering (SEC-MALS), we examined the differences in the metal-binding properties of the two peptides. Both peptides are largely unfolded in the apo state due to the disruption of the hydrophobic core by inclusion of the polar cysteine residues. However, this unfolding is overcome by the addition of Cd(II), Pb(II) and Hg(II), and helical assemblies are formed. Both peptides have differing affinities for these metal ions, a fact likely attributed to the differing sizes of the ions. We also show that the oligomerisation state of the peptide complexes and the coordination geometries of the metal ions differ between the two peptide scaffolds. These findings highlight that subtle changes in the primary structure of a peptide can have considerable implications for metal binding.

Discovery of Novel Metalloenzyme Inhibitors Based on Property Characterization: Strategy and Application for HDAC1 Inhibitors

Metalloenzymes are ubiquitously present in the human body and are relevant to a variety of diseases. However, the development of metalloenzyme inhibitors is limited by low specificity and poor drug-likeness associated with metal-binding fragments (MBFs). A generalized drug discovery strategy was established, which is characterized by the property characterization of zinc-dependent metalloenzyme inhibitors (ZnMIs). Fifteen potential Zn2+-binding fragments (ZnBFs) were identified, and a customized pharmacophore feature was defined based on these ZnBFs. The customized feature was set as a required feature and applied to a search for novel inhibitors for histone deacetylase 1 (HDAC1). Ten potential HDAC1 inhibitors were recognized, and one of them (compound 9) was a known potent HDAC1 inhibitor. The results demonstrated the effectiveness of our strategy to identify novel inhibitors for zinc-dependent metalloenzymes.

Exonic splicing code and coordination of divalent metals in proteins

Exonic sequences contain both protein-coding and RNA splicing information but the interplay of the protein and splicing code is complex and poorly understood. Here, we have studied traditional and auxiliary splicing codes of human exons that encode residues coordinating two essential divalent metals at the opposite ends of the Irving-Williams series, a universal order of relative stabilities of metal-organic complexes. We show that exons encoding Zn2+-coordinating amino acids are supported much less by the auxiliary splicing motifs than exons coordinating Ca2+. The handicap of the former is compensated by stronger splice sites and uridine-richer polypyrimidine tracts, except for position -3 relative to 3' splice junctions. However, both Ca2+ and Zn2+ exons exhibit close-to-constitutive splicing in multiple tissues, consistent with their critical importance for metalloprotein function and a relatively small fraction of expendable, alternatively spliced exons. These results indicate that constraints imposed by metal coordination spheres on RNA splicing have been efficiently overcome by the plasticity of exon-intron architecture to ensure adequate metalloprotein expression.

Juggling cadmium detoxification and zinc homeostasis: A division of labour between the two C. elegans metallothioneins

The chemical properties of toxic cadmium and essential zinc are very similar, and organisms require intricate mechanisms that drive selective handling of metals. Previously regarded as unspecific "metal sponges", metallothioneins (MTLs) are emerging as metal selectivity filters. By utilizing C. elegans mtl-1 and mtl-2 knockout strains, metal accumulation in single worms, single copy fluorescent-tagged transgenes, isoform specific qPCR and lifespan studies it was possible to demonstrate that the handling of cadmium and zinc by the two C. elegans metallothioneins differs fundamentally: the MTL-2 protein can handle both zinc and cadmium, but when it becomes unavailable, either via a knockout or by elevated cadmium exposure, MTL-1 takes over zinc handling, leaving MTL-2 to sequester cadmium. This division of labour is reflected in the folding behaviour of the proteins: MTL-1 folded well in presence of zinc but not cadmium, the reverse was the case for MTL-2. These differences are in part mediated by a zinc-specific mononuclear His3Cys site in the C-terminal insertion of MTL-1; its removal affected the entire C-terminal domain and may shift its metal selectivity towards zinc. Overall, we uncover how metallothionein isoform-specific responses and protein properties allow C. elegans to differentiate between toxic cadmium and essential zinc.

MESPEUS: a database of metal coordination groups in proteins

MESPEUS is a freely accessible database which uses carefully selected metal coordination groups found in metalloprotein structures taken from the Protein Data Bank. The database contains geometrical information of metal sites within proteins, including 40 metal types. In order to completely determine the metal coordination, the symmetry-related units of a given protein structure are taken into account and are generated using the appropriate space group symmetry operations. This permits a more complete description of the metal coordination geometry by including all coordinating atoms. The user-friendly web interface allows users to directly search for a metal site of interest using several useful options, including searching for metal elements, metal-donor distances, coordination number, donor residue group, and structural resolution. These searches can be carried out singly or in combination. The details of a metal site and the metal site(s) in the whole structure can be graphically displayed using the interactive web interface. MESPEUS is automatically updated monthly by synchronizing with the PDB database. An investigation for the Mg-ATP interaction is given to demonstrate how MESPEUS can be used to extract information about metal sites by selecting structure and coordination features. MESPEUS is available at http://mespeus.nchu.edu.tw/.

Metal-binding proteins and proteases in RNA viruses: unravelling functional diversity and expanding therapeutic horizons

IMPORTANCE

Metal-binding proteins are pivotal components with diverse functions in organisms, including viruses. Despite their significance, many metalloproteins in viruses remain uncharacterized, posing challenges to understanding viral systems. This study addresses this knowledge gap by identifying and analyzing metal-binding proteins and proteases in RNA viruses. The findings emphasize the prevalence of these proteins as essential functional classes within viruses and shed light on the role of metal ions and metalloproteins in viral replication and pathogenesis. Moreover, this research serves as a crucial foundation for further investigations in this field, offering the potential for developing innovative antiviral strategies. Additionally, the study enhances our understanding of the distribution and evolutionary patterns of metal-binding proteases in major human viruses. Continually exploring metal-binding proteomes across diverse viruses will deepen our knowledge of metal-dependent biological processes and provide valuable insights for combating viral infections, including respiratory viruses and other life-threatening diseases.

Fe(II), Mn(II), and Zn(II) Binding to the C-Terminal Region of FeoB Protein: An Insight into the Coordination Chemistry and Specificity of the Escherichia coli Fe(II) Transporter

The interactions between two peptide ligands [Ac763CCAASTTGDCH773 (P1) and Ac743RRARSRVDIELLATRKSVSSCCAASTTGDCH773 (P2)] derived from the cytoplasmic C-terminal region of Eschericha coli FeoB protein and Fe(II), Mn(II), and Zn(II) ions were investigated. The Feo system is regarded as the most important bacterial Fe(II) acquisition system, being one of the key virulence factors, especially in anaerobic conditions. Located in the inner membrane of Gram-negative bacteria, FeoB protein transports Fe(II) from the periplasm to the cytoplasm. Despite its crucial role in bacterial pathogenicity, the mechanism in which the metal ion is trafficked through the membrane is not yet elucidated. In the gammaproteobacteria class, the cytoplasmic C-terminal part of FeoB contains conserved cysteine, histidine, and glutamic and aspartic acid residues, which could play a vital role in Fe(II) binding in the cytoplasm, receiving the metal ion from the transmembrane helices. In this work, we characterized the complexes formed between the whole cytosolic C-terminal sequence of E. coli FeoB (P2) and its key polycysteine region (P1) with Fe(II), Mn(II), and Zn(II) ions, exploring the specificity of the C-terminal region of FeoB. With the help of a variety of potentiometric, spectroscopic (electron paramagnetic resonance and NMR), and spectrometric (electrospray ionization mass spectrometry) techniques and molecular dynamics, we propose the metal-binding modes of the ligands, compare their affinities toward the metal ions, and discuss the possible physiological role of the C-terminal region of E. coli FeoB.

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.

Characterizing metal-biomolecule interactions by mass spectrometry

Metal micronutrients are essential for life and exist in a delicate balance to maintain an organism's health. The labile nature of metal-biomolecule interactions clouds the understanding of metal binders and metal-mediated conformational changes that are influential to health and disease. Mass spectrometry (MS)-based methods and technologies have been developed to better understand metal micronutrient dynamics in the intra- and extracellular environment. In this review, we describe the challenges associated with studying labile metals in human biology and highlight MS-based methods for the discovery and study of metal-biomolecule interactions.

Effective assessment of lanthanide ion delivery into live cells by paramagnetic NMR spectroscopy

We report an effective assessment of lanthanide ion (Ln3+) delivery into live cells by paramagnetic NMR spectroscopy. Free Ln3+ ions are toxic to live cells resulting in a gradual leakage of target proteins to the extracellular media. The citrate-Ln3+ complex is an efficient and mild reagent over the free Ln3+ form for live cell delivery.

A systematic in silico report on iron and zinc proteome of Zea mays

Zea mays is an essential staple food crop across the globe. Maize contains macro and micronutrients but is limited in essential mineral micronutrients such as Fe and Zn. Worldwide, serious health concerns have risen due to the deficiencies of essential nutrients in human diets, which rigorously jeopardizes economic development. In the present study, the systematic in silico approach has been used to predict Fe and Zn binding proteins from the whole proteome of maize. A total of 356 and 546 putative proteins have been predicted, which contain sequence and structural motifs for Fe and Zn ions, respectively. Furthermore, the functional annotation of these predicted proteins, based on their domains, subcellular localization, gene ontology, and literature support, showed their roles in distinct cellular and biological processes, such as metabolism, gene expression and regulation, transport, stress response, protein folding, and proteolysis. The versatile roles of these shortlisted putative Fe and Zn binding proteins of maize could be used to manipulate many facets of maize physiology. Moreover, in the future, the predicted Fe and Zn binding proteins may act as relevant, novel, and economical markers for various crop improvement programs.

Elucidating the zinc-binding proteome of Fusarium oxysporum f. sp. lycopersici with particular emphasis on zinc-binding effector proteins

Fusarium oxysporum f. sp. lycopersici is a soil-borne phytopathogenic species which causes vascular wilt disease in the Solanum lycopersicum (tomato). Due to the continuous competition for zinc usage by Fusarium and its host during infection makes zinc-binding proteins a hotspot for focused investigation. Zinc-binding effector proteins are pivotal during the infection process, working in conjunction with other essential proteins crucial for its biological activities. This work aims at identifying and analysing zinc-binding proteins and zinc-binding proteins effector candidates of Fusarium. We have identified three hundred forty-six putative zinc-binding proteins; among these proteins, we got two hundred and thirty zinc-binding proteins effector candidates. The functional annotation, subcellular localization, and Gene Ontology analysis of these putative zinc-binding proteins revealed their probable role in wide range of cellular and biological processes such as metabolism, gene expression, gene expression regulation, protein biosynthesis, protein folding, cell signalling, DNA repair, and RNA processing. Sixteen proteins were found to be putatively secretory in nature. Eleven of these were putative zinc-binding protein effector candidates may be involved in pathogen-host interaction during infection. The information obtained here may enhance our understanding to design, screen, and apply the zinc-metal ion-based antifungal agents to protect the S. lycopersicum and control the vascular wilt caused by F. oxysporum.

Mutations observed in somatic evolution reveal underlying gene mechanisms

Highly sensitive DNA sequencing techniques have allowed the discovery of large numbers of somatic mutations in normal tissues. Some mutations confer a competitive advantage over wild-type cells, generating expanding clones that spread through the tissue. Competition between mutant clones leads to selection. This process can be considered a large scale, in vivo screen for mutations increasing cell fitness. It follows that somatic missense mutations may offer new insights into the relationship between protein structure, function and cell fitness. We present a flexible statistical method for exploring the selection of structural features in data sets of somatic mutants. We show how this approach can evidence selection of specific structural features in key drivers in aged tissues. Finally, we show how drivers may be classified as fitness-enhancing and fitness-suppressing through different patterns of mutation enrichment. This method offers a route to understanding the mechanism of protein function through in vivo mutant selection.

Atomically Accurate Design of Metalloproteins with Predefined Coordination Geometries

We report a new computational protein design method for the construction of oligomeric protein assemblies around metal centers with predefined coordination geometries. We apply this method to design two homotrimeric assemblies, Tet4 and TP1, with tetrahedral and trigonal-pyramidal tris(histidine) metal coordination geometries, respectively, and demonstrate that both assemblies form the targeted metal centers with ≤0.2 Å accuracy. Although Tet4 and TP1 are constructed from the same parent protein building block, they are distinct in terms of their overall architectures, the environment surrounding the metal centers, and their metal-based reactivities, illustrating the versatility of our approach.

Targeting Metalloenzymes: The "Achilles' Heel" of Viruses and Parasites

Metalloenzymes are central to the regulation of a wide range of essential viral and parasitic functions, including protein degradation, nucleic acid modification, and many others. Given the impact of infectious diseases on human health, inhibiting metalloenzymes offers an attractive approach to disease therapy. Metal-chelating agents have been expansively studied as antivirals and antiparasitics, resulting in important classes of metal-dependent enzyme inhibitors. This review provides the recent advances in targeting the metalloenzymes of viruses and parasites that impose a significant burden on global public health, including influenza A and B, hepatitis B and C, and human immunodeficiency viruses as well as Trypanosoma brucei and Trypanosoma cruzi.

A Comparison of Bonded and Nonbonded Zinc(II) Force Fields with NMR Data

Classical molecular dynamics (MD) simulations are widely used to inspect the behavior of zinc(II)-proteins at the atomic level, hence the need to properly model the zinc(II) ion and the interaction with its ligands. Different approaches have been developed to represent zinc(II) sites, with the bonded and nonbonded models being the most used. In the present work, we tested the well-known zinc AMBER force field (ZAFF) and a recently developed nonbonded force field (NBFF) to assess how accurately they reproduce the dynamic behavior of zinc(II)-proteins. For this, we selected as benchmark six zinc-fingers. This superfamily is extremely heterogenous in terms of architecture, binding mode, function, and reactivity. From repeated MD simulations, we computed the order parameter (S2) of all backbone N-H bond vectors in each system. These data were superimposed to heteronuclear Overhauser effect measurements taken by NMR spectroscopy. This provides a quantitative estimate of the accuracy of the FFs in reproducing protein dynamics, leveraging the information about the protein backbone mobility contained in the NMR data. The correlation between the MD-computed S2 and the experimental data indicated that both tested FFs reproduce well the dynamic behavior of zinc(II)-proteins, with comparable accuracy. Thus, along with ZAFF, NBFF represents a useful tool to simulate metalloproteins with the advantage of being extensible to diverse systems such as those bearing dinuclear metal sites.

Structural and biochemical insights into Zn2+-bound EF-hand proteins, EFhd1 and EFhd2

EF-hand proteins, which contain a Ca²⁺-binding EF-hand motif, are involved in regulating diverse cellular functions. Ca²⁺ binding induces conformational changes that modulate the activities of EF-hand proteins. Moreover, these proteins occasionally modify their activities by coordinating metals other than Ca²⁺, including Mg²⁺, Pb²⁺ and Zn²⁺, within their EF-hands. EFhd1 and EFhd2 are homologous EF-hand proteins with similar structures. Although separately localized within cells, both are actin-binding proteins that modulate F-actin rearrangement through Ca²⁺-independent actin-binding and Ca²⁺-dependent actin-bundling activity. Although Ca²⁺ is known to affect the activities of EFhd1 and EFhd2, it is not known whether their actin-related activities are affected by other metals. Here, the crystal structures of the EFhd1 and EFhd2 core domains coordinating Zn²⁺ ions within their EF-hands are reported. The presence of Zn²⁺ within EFhd1 and EFhd2 was confirmed by analyzing anomalous signals and the difference between anomalous signals using data collected at the peak positions as well as low-energy remote positions at the Zn K-edge. EFhd1 and EFhd2 were also found to exhibit Zn²⁺-independent actin-binding and Zn²⁺-dependent actin-bundling activity. This suggests the actin-related activities of EFhd1 and EFhd2 could be regulated by Zn²⁺ as well as Ca²⁺.

Global Molecular Response of Paracoccidioides brasiliensis to Zinc Deprivation: Analyses at Transcript, Protein and MicroRNA Levels

Zinc is one of the main micronutrients for all organisms. One of the defense mechanisms used by the host includes the sequestration of metals used in fungal metabolism, such as iron and zinc. There are several mechanisms that maintain the balance in the intracellular zinc supply. MicroRNAs are effector molecules of responses between the pathogen and host, favoring or preventing infection in many microorganisms. Fungi of the Paracoccidioides genus are thermodimorphic and the etiological agents of paracoccidioidomycosis (PCM). In the current pandemic scenario world mycosis studies continue to be highly important since a significant number of patients with COVID-19 developed systemic mycoses, co-infections that complicated their clinical condition. The objective was to identify transcriptomic and proteomic adaptations in Paracoccidioides brasiliensis during zinc deprivation. Nineteen microRNAs were identified, three of which were differentially regulated. Target genes regulated by those microRNAs are elements of zinc homeostasis such as ZRT1, ZRT3 and COT1 transporters. Transcription factors that have zinc in their structure are also targets of those miRNAs. Transcriptional and proteomic data suggest that P. brasiliensis undergoes metabolic remodeling to survive zinc deprivation and that miRNAs may be part of the regulatory process.

Metalloproteome of human-infective RNA viruses: a study towards understanding the role of metal ions in virology

Metalloproteins and metal-based inhibitors have been shown to effectively combat infectious diseases, particularly those caused by RNA viruses. In this study, a diverse set of bioinformatics methods was employed to identify metal-binding proteins of human RNA viruses. Seventy-three viral proteins with a high probability of being metal-binding proteins were identified. These proteins included 40 zinc-, 47 magnesium- and 14 manganese-binding proteins belonging to 29 viral species and eight significant viral families, including Coronaviridae, Flaviviridae and Retroviridae. Further functional characterization has revealed that these proteins play a critical role in several viral processes, including viral replication, fusion and host viral entry. They fall under the essential categories of viral proteins, including polymerase and protease enzymes. Magnesium ion is abundantly predicted to interact with these viral enzymes, followed by zinc. In addition, this study also examined the evolutionary aspects of predicted viral metalloproteins, offering essential insights into the metal utilization patterns among different viral species. The analysis indicates that the metal utilization patterns are conserved within the functional classes of the proteins. In conclusion, the findings of this study provide significant knowledge on viral metalloproteins that can serve as a valuable foundation for future research in this area.

The putative metal-binding proteome of the Coronaviridae family

Metalloproteins are well-known for playing various physicochemical processes in all life forms, including viruses. Some life-threatening viruses (such as some members of the Coronaviridae family of viruses) are emerged and remerged frequently and are rapidly transmitted throughout the globe. This study aims to identify and characterize the metal-binding proteins (MBPs) of the Coronaviridae family of viruses and further provides insight into the MBP's role in sustaining and propagating viruses inside a host cell and in the outer environment. In this study, the available proteome of the Coronaviridae family was exploited. Identified potential MBPs were analyzed for their functional domains, structural aspects, and subcellular localization. We also demonstrate phylogenetic aspects of all predicted MBPs among other Coronaviridae family members to understand the evolutionary trend among their respective hosts. A total of 256 proteins from 51 different species of coronaviruses are predicted as MBPs. These MBPs perform various key roles in the replication and survival of viruses within the host cell. Cysteine, aspartic acid, threonine, and glutamine are key amino acid residues interacting with respective metal ions. Our observations also indicate that the metalloproteins of this family of viruses circulated and evolved in different hosts, which supports the zoonotic nature of coronaviruses. The comprehensive information on MBPs of the Coronaviridae family may be further helpful in designing novel therapeutic metalloprotein targets. Moreover, the study of viral MBPs can also help to understand the roles of MBPs in virus pathogenesis and virus-host interactions.

UniProt: the Universal Protein Knowledgebase in 2023

The aim of the UniProt Knowledgebase is to provide users with a comprehensive, high-quality and freely accessible set of protein sequences annotated with functional information. In this publication we describe enhancements made to our data processing pipeline and to our website to adapt to an ever-increasing information content. The number of sequences in UniProtKB has risen to over 227 million and we are working towards including a reference proteome for each taxonomic group. We continue to extract detailed annotations from the literature to update or create reviewed entries, while unreviewed entries are supplemented with annotations provided by automated systems using a variety of machine-learning techniques. In addition, the scientific community continues their contributions of publications and annotations to UniProt entries of their interest. Finally, we describe our new website (https://www.uniprot.org/), designed to enhance our users' experience and make our data easily accessible to the research community. This interface includes access to AlphaFold structures for more than 85% of all entries as well as improved visualisations for subcellular localisation of proteins.

Annotation of biologically relevant ligands in UniProtKB using ChEBI

MOTIVATION

To provide high quality, computationally tractable annotation of binding sites for biologically relevant (cognate) ligands in UniProtKB using the chemical ontology ChEBI (Chemical Entities of Biological Interest), to better support efforts to study and predict functionally relevant interactions between protein sequences and structures and small molecule ligands.

RESULTS

We structured the data model for cognate ligand binding site annotations in UniProtKB and performed a complete reannotation of all cognate ligand binding sites using stable unique identifiers from ChEBI, which we now use as the reference vocabulary for all such annotations. We developed improved search and query facilities for cognate ligands in the UniProt website, REST API and SPARQL endpoint that leverage the chemical structure data, nomenclature and classification that ChEBI provides.

AVAILABILITY AND IMPLEMENTATION

Binding site annotations for cognate ligands described using ChEBI are available for UniProtKB protein sequence records in several formats (text, XML and RDF) and are freely available to query and download through the UniProt website (www.uniprot.org), REST API (www.uniprot.org/help/api), SPARQL endpoint (sparql.uniprot.org/) and FTP site (https://ftp.uniprot.org/pub/databases/uniprot/).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Metal trafficking in the cell: Combining atomic resolution with cellular dimension

Metals are widely present in biological systems as simple ions or complex cofactors, and are involved in a variety of processes essential for life. Their transport inside cells and insertion into the binding sites of the proteins that need metals to function occur through complex and selective pathways involving dedicated multiprotein machineries specifically and transiently interacting with each other, often sharing the coordination of metal ions and/or cofactors. The understanding of these machineries requires integrated approaches, ranging from bioinformatics to experimental investigations, possibly in the cellular context. In this review, we report two case studies where the use of integrated in vitro and in cellulo approaches is necessary to clarify at atomic resolution essential aspects of metal trafficking in cells.

Structural modeling of protein ensembles between E3 RING ligases and SARS-CoV-2: The role of zinc binding domains

BACKGROUND

The ubiquitin system is a modification process with many different cellular functions including immune signaling and antiviral functions. E3 ubiquitin ligases are enzymes that recruit an E2 ubiquitin-conjugating enzyme bound to ubiquitin in order to catalyze the transfer of ubiquitin from the E2 to a protein substrate. The RING E3s, the most abundant type of ubiquitin ligases, are characterized by a zinc (II)-binding domain called RING (Really Interesting New Gene). Viral replication requires modifying and hijacking key cellular pathways within host cells such as cellular ubiquitination. There are well-established examples where a viral proteins bind to RING E3s, redirecting them to degrade otherwise long-lived host proteins or inhibiting E3's ubiquitination activity. Recently, three binary interactions between SARS-CoV-2 proteins and innate human immune signaling Ε3 RING ligases: NSP15-RNF41, ORF3a-TRIM59 and NSP9-MIB1 have been experimentally established.

METHODS

In this work, we have investigated the mode of the previous experimentally supported NSP15-RNF41, ORF3a,-TRIM59 and NSP9-MIB1 binary interactions by in silico methodologies intending to provide structural insights of E3-virus interplay that can help identify potential inhibitors that could block SARS-CoV-2 infection of immune cells.

CONCLUSION

In silico methodologies have shown that the above human E3 ligases interact with viral partners through their Zn(II) binding domains. This RING mediated formation of stable SARS-CoV-2-E3 complexes indicates a critical structural role of RING domains in immune system disruption by SARS-CoV-2-infection.

DATA AVAILABILITY

The data used to support the findings of this research are included within the article and are labeled with references.

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.

Biological Effects of Human Exposure to Environmental Cadmium

Cadmium (Cd) is a toxic metal for the human organism and for all ecosystems. Cd is naturally found at low levels; however, higher amounts of Cd in the environment result from human activities as it spreads into the air and water in the form of micropollutants as a consequence of industrial processes, pollution, waste incineration, and electronic waste recycling. The human body has a limited ability to respond to Cd exposure since the metal does not undergo metabolic degradation into less toxic species and is only poorly excreted. The extremely long biological half-life of Cd essentially makes it a cumulative toxin; chronic exposure causes harmful effects from the metal stored in the organs. The present paper considers exposure and potential health concerns due to environmental cadmium. Exposure to Cd compounds is primarily associated with an elevated risk of lung, kidney, prostate, and pancreatic cancer. Cd has also been linked to cancers of the breast, urinary system, and bladder. The multiple mechanisms of Cd-induced carcinogenesis include oxidative stress with the inhibition of antioxidant enzymes, the promotion of lipid peroxidation, and interference with DNA repair systems. Cd²⁺ can also replace essential metal ions, including redox-active ones. A total of 12 cancer types associated with specific genes coding for the Cd-metalloproteome were identified in this work. In addition, we summarize the proper treatments of Cd poisoning, based on the use of selected Cd detoxifying agents and chelators, and the potential for preventive approaches to counteract its chronic exposure.

Studying Peptide-Metal Ion Complex Structures by Solution-State NMR

Metal chelation can provide structural stability and form reactive centers in metalloproteins. Approximately one third of known protein structures are metalloproteins, and metal binding, or the lack thereof, is often implicated in disease, making it necessary to be able to study these systems in detail. Peptide-metal complexes are both present in nature and can provide a means to focus on the binding region of a protein and control experimental variables to a high degree. Structural studies of peptide complexes with metal ions by nuclear magnetic resonance (NMR) were surveyed for all the essential metal complexes and many non-essential metal complexes. The various methods used to study each metal ion are presented together with examples of recent research. Many of these metal systems have been individually reviewed and this current overview of NMR studies of metallopeptide complexes aims to provide a basis for inspiration from structural studies and methodology applied in the field.

Metal ion availability and homeostasis as drivers of metabolic evolution and enzyme function

Metal ions are potent catalysts and have been available for cellular biochemistry at all stages of evolution. Growing evidence suggests that metal catalysis was critical for the origin of the very first metabolic reactions. With approximately 80% of modern metabolic pathways being dependent on metal ions, metallocatalysis and homeostasis continue to be essential for intracellular metabolic networks and physiology. However, the genetic network that controls metal ion homeostasis and the impact of metal availability on metabolism is poorly understood. Here, we review recent work on gene and protein evolution relevant for better understanding metal ion biology and its role in metabolism. We highlight the importance of analysing the origin and evolution of enzyme catalysis in the context of catalytically relevant metal ions, summarise unanswered questions essential for developing a comprehensive understanding of metal ion homeostasis and advocate for the consideration of metal ion properties and availability in the design and directed evolution of novel enzymes and pathways.

Alignment-free metal ion-binding site prediction from protein sequence through pretrained language model and multi-task learning

More than one-third of the proteins contain metal ions in the Protein Data Bank. Correct identification of metal ion-binding residues is important for understanding protein functions and designing novel drugs. Due to the small size and high versatility of metal ions, it remains challenging to computationally predict their binding sites from protein sequence. Existing sequence-based methods are of low accuracy due to the lack of structural information, and time-consuming owing to the usage of multi-sequence alignment. Here, we propose LMetalSite, an alignment-free sequence-based predictor for binding sites of the four most frequently seen metal ions in BioLiP (Zn2+, Ca2+, Mg2+ and Mn2+). LMetalSite leverages the pretrained language model to rapidly generate informative sequence representations and employs transformer to capture long-range dependencies. Multi-task learning is adopted to compensate for the scarcity of training data and capture the intrinsic similarities between different metal ions. LMetalSite was shown to surpass state-of-the-art structure-based methods by more than 19.7, 14.4, 36.8 and 12.6% in area under the precision recall on the four independent tests, respectively. Further analyses indicated that the self-attention modules are effective to learn the structural contexts of residues from protein sequence. We provide the data sets, source codes and trained models of LMetalSite at https://github.com/biomed-AI/LMetalSite.

Glaucoma-TrEl: A web-based interactive database to build evidence-based hypotheses on the role of trace elements in glaucoma

Objective

Glaucoma is a chronic neurological disease that is associated with high intraocular pressure (IOP), causes gradual damage to retinal ganglion cells, and often culminates in vision loss. Recent research suggests that glaucoma is a complex multifactorial disease in which multiple interlinked genes and pathways play a role during onset and development. Also, differential availability of trace elements seems to play a role in glaucoma pathophysiology, although their mechanism of action is unknown. The aim of this work is to disseminate a web-based repository on interactions between trace elements and protein-coding genes linked to glaucoma pathophysiology.

Results

In this study, we present Glaucoma-TrEl, a web database containing information about interactions between trace elements and protein-coding genes that are linked to glaucoma. In the database, we include interactions between 437 unique genes and eight trace elements. Our analysis found a large number of interactions between trace elements and protein-coding genes mutated or linked to the pathophysiology of glaucoma. We associated genes interacting with multiple trace elements to pathways known to play a role in glaucoma. The web-based platform provides an easy-to-use and interactive tool, which serves as an information hub facilitating future research work on trace elements in glaucoma.

A cold‐active esterase enhances mesophilic properties through Mn 2+ binding

A key aspect of adaptation to cold environments is the production of cold-active enzymes by psychrophilic organisms. These enzymes not only have high activity at low temperatures, but also exhibit remarkable structural flexibility and thermolability. In this context, the role of metal ions has been little explored, and the few available studies seem to suggest that metal binding counteracts structural flexibility. This article reports an investigation into the role of the binding of manganese ion (Mn²⁺ ) in the thermal adaptation of an esterase (M-Est) of the GDSx family, identified in the genome of the Antarctic bacterium Marinomonas sp. ef1. M-Est is specific for esters containing acetate groups and turned out to be a highly thermolabile cold-active enzyme, with a catalysis optimum temperature of 5 °C and a melting temperature of 31.7 °C. A combination of biochemical and computational analyses, including molecular dynamics simulations, revealed that M-Est binds Mn²⁺ ions via a single binding site located on the surface of the enzyme, close to the active site. Although the interaction between M-Est and Mn²⁺ induces only local conformational changes involving the active site, quite surprisingly they trigger an improvement in both thermal stability and catalytic efficiency under mild temperature conditions. These results, together with the conservation of the Mn²⁺ binding site among psychrophilic and psychrotolerant homologues, suggest that Mn²⁺ binding may be a useful, albeit atypical, strategy to mitigate the detrimental effects of temperature on true cold-active enzymes.

Design of a Flexible, Zn-Selective Protein Scaffold that Displays Anti-Irving-Williams Behavior

Selective metal binding is a key requirement not only for the functions of natural metalloproteins but also for the potential applications of artificial metalloproteins in heterogeneous environments such as cells and environmental samples. The selection of transition-metal ions through protein design can, in principle, be achieved through the appropriate choice and the precise positioning of amino acids that comprise the primary metal coordination sphere. However, this task is made difficult by the intrinsic flexibility of proteins and the fact that protein design approaches generally lack the sub-Å precision required for the steric selection of metal ions. We recently introduced a flexible/probabilistic protein design strategy (MASCoT) that allows metal ions to search for optimal coordination geometry within a flexible, yet covalently constrained dimer interface. In an earlier proof-of-principle study, we used MASCoT to generate an artificial metalloprotein dimer, (AB)₂, which selectively bound Co^II and Ni^II over Cu^II (as well as other first-row transition-metal ions) through the imposition of a rigid octahedral coordination geometry, thus countering the Irving-Williams trend. In this study, we set out to redesign (AB)₂ to examine the applicability of MASCoT to the selective binding of other metal ions. We report here the design and characterization of a new flexible protein dimer, B₂, which displays Zn^II selectivity over all other tested metal ions including Cu^II both in vitro and in cellulo. Selective, anti-Irving-Williams Zn^II binding by B₂ is achieved through the formation of a unique trinuclear Zn coordination motif in which His and Glu residues are rigidly placed in a tetrahedral geometry. These results highlight the utility of protein flexibility in the design and discovery of selective binding motifs.

Calcium dysregulation potentiates wild-type myocilin misfolding: implications for glaucoma pathogenesis

Myocilin is secreted from trabecular meshwork cells to an eponymous extracellular matrix that is critical for maintaining intraocular pressure. Missense mutations found in the myocilin olfactomedin domain (OLF) lead to intracellular myocilin misfolding and are causative for the heritable form of early-onset glaucoma. The OLF domain contains a unique internal, hetero-dinuclear calcium site. Here, we tested the hypothesis that calcium dysregulation causes wild-type (WT) myocilin misfolding reminiscent of that observed for disease variants. Using two cellular models expressing WT myocilin, we show that the Ca2+ ATPase channel blocker thapsigargin inhibits WT myocilin secretion. Intracellular WT myocilin is at least partly insoluble and aggregated in the endoplasmic reticulum (ER), and stains positively with an amyloid dye. By comparing the effect of thapsigargin on WT myocilin to that on a de novo secretion-competent Ca2+-free variant D478S, we discern that non-secretion of WT myocilin is due initially to calcium dysregulation, and is potentiated further by resultant ER stress. In E. coli, depletion of calcium leads to recombinant expression of misfolded isolated WT OLF but the D478S variant is still produced as a folded monomer. Treatment of cells expressing a double mutant composed of D478S and either disease variants P370L or Y437H with thapsigargin promotes its misfolding and aggregation, demonstrating the limits of D478S to correct secretion defects. Taken together, the heterodinuclear calcium site is a liability for proper folding of myocilin. Our study suggests a molecular mechanism by which WT myocilin misfolding may contribute broadly to glaucoma-associated ER stress. This study explores the effect of calcium depletion on myocilin olfactomedin domain folding.

The zinc proteome of SARS-CoV-2

Zinc is an essential element for human health. Among its many functions, zinc(II) modulates the immune response to infections and, at high concentrations or in the presence of ionophores, inhibits the replication of various RNA viruses. Structural biology studies on severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) revealed that zinc(II) is the most common metal ion that binds to viral proteins. However, the number of zinc(II)-binding sites identified by experimental methods is far from exhaustive, as metal ions may be lost during protein purification protocols. To better define the zinc(II)-binding proteome of coronavirus, we leveraged the wealth of deposited structural data and state-of-the-art bioinformatics methods. Through this in silico approach, 15 experimental zinc(II) sites were identified and a further 22 were predicted in Spike, open reading frame (ORF)3a/d, ORF8, and several nonstructural proteins, highlighting an essential role of zinc(II) in viral replication. Furthermore, the structural relationships between viral and eukaryotic sites (typically zinc fingers) indicate that SARS-CoV-2 can compete with human proteins for zinc(II) binding. Given the double-edged effect of zinc(II) ions, both essential and toxic to coronavirus, only the complete elucidation of the structural and regulatory zinc(II)-binding sites can guide selective antiviral strategies based on zinc supplementation.

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.

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.

Impacts of targeting different hydration free energy references on the development of ion potentials

Hydration free energy (HFE) as the most important solvation parameter is often targeted in ion model development, even though the reported values differ by dozens of kcal mol^-1 mainly due to the experimentally undetermined HFE of the proton ΔG°(H⁺). The choice of ΔG°(H⁺) obviously affects the hydration of single ions and the relative HFE between the ions with different (magnitude or sign) charges, and the impacts of targeted HFEs on the ion solvation and ion-ion interactions are largely unrevealed. Here we designed point charge models of K⁺, Mg²⁺, Al³⁺, and Cl^- ions targeting a variety of HFE references and then investigated the HFE influences on the simulations of dilute and concentrated ion solutions and of the salt ion pairs in gas, liquid, and solid phases. Targeting one more property of ion-water oxygen distances (IOD) leaves the ion-water binding distance invariant, while the binding strength increases with the decreasing (more negative) HFE of ions as a result of a decrease in ΔG°(H⁺) for the cation and an increase in ΔG°(H⁺) for the anion. The increase in ΔG°(H⁺) leads to strengthened cation-anion interactions and thus to close ion-ion contacts, low osmotic pressures, and small activity derivatives in concentrated ion solutions as well as too stable ion pairs of the salts in different phases. The ion diffusivity and water exchange rates around the ions are simply not HFE dependent but rather more complex. Targeting both the aqueous IOD and salt crystal properties of KCl was also attempted and the comparison between different models indicates the complexity and challenge in obtaining a balanced performance between different phases using classical force fields. Our results also support that a real ΔG°(H⁺) value of -259.8 kcal mol^-1 recommended by Hünenberger and Reif guides ion models to reproduce ion-water and ion-ion interactions reasonably at relatively low salt concentrations. Simulations of a metalloprotein show that a relatively more positive ΔG°(H⁺) for Mg²⁺ model is better for a reasonable description of the metal binding network.