Validation and correction of Zn-CysxHisy complexes

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.

The CCP4 suite: integrative software for macromolecular crystallography

The Collaborative Computational Project No. 4 (CCP4) is a UK-led international collective with a mission to develop, test, distribute and promote software for macromolecular crystallography. The CCP4 suite is a multiplatform collection of programs brought together by familiar execution routines, a set of common libraries and graphical interfaces. The CCP4 suite has experienced several considerable changes since its last reference article, involving new infrastructure, original programs and graphical interfaces. This article, which is intended as a general literature citation for the use of the CCP4 software suite in structure determination, will guide the reader through such transformations, offering a general overview of the new features and outlining future developments. As such, it aims to highlight the individual programs that comprise the suite and to provide the latest references to them for perusal by crystallographers around the world.

AlphaFill: enriching AlphaFold models with ligands and cofactors

Artificial intelligence-based protein structure prediction approaches have had a transformative effect on biomolecular sciences. The predicted protein models in the AlphaFold protein structure database, however, all lack coordinates for small molecules, essential for molecular structure or function: hemoglobin lacks bound heme; zinc-finger motifs lack zinc ions essential for structural integrity and metalloproteases lack metal ions needed for catalysis. Ligands important for biological function are absent too; no ADP or ATP is bound to any of the ATPases or kinases. Here we present AlphaFill, an algorithm that uses sequence and structure similarity to 'transplant' such 'missing' small molecules and ions from experimentally determined structures to predicted protein models. The algorithm was successfully validated against experimental structures. A total of 12,029,789 transplants were performed on 995,411 AlphaFold models and are available together with associated validation metrics in the alphafill.eu databank, a resource to help scientists make new hypotheses and design targeted experiments.

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.

Simplified quality assessment for small-molecule ligands in the Protein Data Bank

More than 70% of the experimentally determined macromolecular structures in the Protein Data Bank (PDB) contain small-molecule ligands. Quality indicators of ∼643,000 ligands present in ∼106,000 PDB X-ray crystal structures have been analyzed. Ligand quality varies greatly with regard to goodness of fit between ligand structure and experimental data, deviations in bond lengths and angles from known chemical structures, and inappropriate interatomic clashes between the ligand and its surroundings. Based on principal component analysis, correlated quality indicators of ligand structure have been aggregated into two largely orthogonal composite indicators measuring goodness of fit to experimental data and deviation from ideal chemical structure. Ranking of the composite quality indicators across the PDB archive enabled construction of uniformly distributed composite ranking score. This score is implemented at RCSB.org to compare chemically identical ligands in distinct PDB structures with easy-to-interpret two-dimensional ligand quality plots, allowing PDB users to quickly assess ligand structure quality and select the best exemplars.

The missing link: covalent linkages in structural models

Covalent linkages between constituent blocks of macromolecules and ligands have been subject to inconsistent treatment during the model-building, refinement and deposition process. This may stem from a number of sources, including difficulties with initially detecting the covalent linkage, identifying the correct chemistry, obtaining an appropriate restraint dictionary and ensuring its correct application. The analysis presented herein assesses the extent of problems involving covalent linkages in the Protein Data Bank (PDB). Not only will this facilitate the remediation of existing models, but also, more importantly, it will inform and thus improve the quality of future linkages. By considering linkages of known type in the CCP4 Monomer Library (CCP4-ML), failure to model a covalent linkage is identified to result in inaccurate (systematically longer) interatomic distances. Scanning the PDB for proximal atom pairs that do not have a corresponding type in the CCP4-ML reveals a large number of commonly occurring types of unannotated potential linkages; in general, these may or may not be covalently linked. Manual consideration of the most commonly occurring cases identifies a number of genuine classes of covalent linkages. The recent expansion of the CCP4-ML is discussed, which has involved the addition of over 16 000 and the replacement of over 11 000 component dictionaries using AceDRG. As part of this effort, the CCP4-ML has also been extended using AceDRG link dictionaries for the aforementioned linkage types identified in this analysis. This will facilitate the identification of such linkage types in future modelling efforts, whilst concurrently easing the process involved in their application. The need for a universal standard for maintaining link records corresponding to covalent linkages, and references to the associated dictionaries used during modelling and refinement, following deposition to the PDB is emphasized. The importance of correctly modelling covalent linkages is demonstrated using a case study, which involves the covalent linkage of an inhibitor to the main protease in various viral species, including SARS-CoV-2. This example demonstrates the importance of properly modelling covalent linkages using a comprehensive restraint dictionary, as opposed to just using a single interatomic distance restraint or failing to model the covalent linkage at all.

Modelling covalent linkages in CCP4

In this contribution, the current protocols for modelling covalent linkages within the CCP4 suite are considered. The mechanism used for modelling covalent linkages is reviewed: the use of dictionaries for describing changes to stereochemistry as a result of the covalent linkage and the application of link-annotation records to structural models to ensure the correct treatment of individual instances of covalent linkages. Previously, linkage descriptions were lacking in quality compared with those of contemporary component dictionaries. Consequently, AceDRG has been adapted for the generation of link dictionaries of the same quality as for individual components. The approach adopted by AceDRG for the generation of link dictionaries is outlined, which includes associated modifications to the linked components. A number of tools to facilitate the practical modelling of covalent linkages available within the CCP4 suite are described, including a new restraint-dictionary accumulator, the Make Covalent Link tool and AceDRG interface in Coot, the 3D graphical editor JLigand and the mechanisms for dealing with covalent linkages in the CCP4i2 and CCP4 Cloud environments. These integrated solutions streamline and ease the covalent-linkage modelling workflow, seamlessly transferring relevant information between programs. Current recommended practice is elucidated by means of instructive practical examples. By summarizing the different approaches to modelling linkages that are available within the CCP4 suite, limitations and potential pitfalls that may be encountered are highlighted in order to raise awareness, with the intention of improving the quality of future modelled covalent linkages in macromolecular complexes.

Exploring the Role of Glycans in the Interaction of SARS-CoV-2 RBD and Human Receptor ACE2

COVID-19 is a highly infectious respiratory disease caused by the novel coronavirus SARS-CoV-2. It has become a global pandemic and its frequent mutations may pose new challenges for vaccine design. During viral infection, the Spike RBD of SARS-CoV-2 binds the human host cell receptor ACE2, enabling the virus to enter the host cell. Both the Spike and ACE2 are densely glycosylated, and it is unclear how distinctive glycan types may modulate the interaction of RBD and ACE2. Detailed understanding of these determinants is key for the development of novel therapeutic strategies. To this end, we perform extensive all-atom simulations of the (i) RBD-ACE2 complex without glycans, (ii) RBD-ACE2 with oligomannose MAN9 glycans in ACE2, and (iii) RBD-ACE2 with complex FA2 glycans in ACE2. These simulations identify the key residues at the RBD-ACE2 interface that form contacts with higher probabilities, thus providing a quantitative evaluation that complements recent structural studies. Notably, we find that this RBD-ACE2 contact signature is not altered by the presence of different glycoforms, suggesting that RBD-ACE2 interaction is robust. Applying our simulated results, we illustrate how the recently prevalent N501Y mutation may alter specific interactions with host ACE2 that facilitate the virus-host binding. Furthermore, our simulations reveal how the glycan on Asn90 of ACE2 can play a distinct role in the binding and unbinding of RBD. Finally, an energetics analysis shows that MAN9 glycans on ACE2 decrease RBD-ACE2 affinity, while FA2 glycans lead to enhanced binding of the complex. Together, our results provide a more comprehensive picture of the detailed interplay between virus and human receptor, which is much needed for the discovery of effective treatments that aim at modulating the physical-chemical properties of this virus.

Quantitative trace metal determinations in cell culture media using LS-APGD-MS and ICP-OES with free/bound species differentiation following polymer fiber separations

Liquid sampling-atmospheric pressure glow discharge-mass spectrometry (LS-APGD-MS) and inductively coupled plasma-optical emission spectroscopy (ICP-OES) were employed for the quantification of trace metals in cell culture media and their capabilities compared. The LS-APGD is interfaced here to a compact mass spectrometer (Advion CMS) towards the development of an at-bioreactor process monitoring strategy. Both techniques have been previously employed for the quantification of trace metals in samples of various complexities, making them a natural choice for this application. They have also demonstrated comparable analytical figures of merit including limits of detection (LOD), matrix tolerance, etc. While cell culture media is a complex sample, the ICP-OES technique was unaffected by the matrix. However, the LS-APGD-MS suffered from increases in spectral background. Despite this, both techniques achieved appropriate LODs for all metals analyzed in this work (Cu, Fe, Zn, Co, Mn, Ni; LOD < 100 ng mL^-1), except for Mn and Ni via LS-APGD-MS. To overcome the increased background seen on the LS-APGD-MS, a capillary channeled polymer (C-CP) polypropylene (PPY) fiber stationary phase was employed as an on-line separation for the removal of organic components prior to sample introduction into the plasma. It was further determined that Ni was retained on the column, preventing the detection of this element via LS-APGD-MS, and insights into metal speciation were discussed. Following implementation of this on-line separation strategy, the agreement between the techniques was acceptable for all analytes, and was excellent for Cu, Fe, and Zn.

LAHMA: structure analysis through local annotation of homology-matched amino acids

Comparison of homologous structure models is a key step in analyzing protein structure. With a wealth of homologous structures, comparison becomes a tedious process, and often only a small (user-biased) selection of data is used. A multitude of structural superposition algorithms are then typically used to visualize the structures together in 3D and to compare them. Here, the Local Annotation of Homology-Matched Amino acids (LAHMA) website (https://lahma.pdb-redo.eu) is presented, which compares any structure model with all of its close homologs from the PDB-REDO databank. LAHMA displays structural features in sequence space, allowing users to uncover differences between homologous structure models that can be analyzed for their relevance to chemistry or biology. LAHMA visualizes numerous structural features, also allowing one-click comparison of structure-quality plots (for example the Ramachandran plot) and `in-browser' structural visualization of 3D models.

Local and global analysis of macromolecular atomic displacement parameters

This paper describes the global and local analysis of atomic displacement parameters (ADPs) of macromolecules in X-ray crystallography. The distribution of ADPs is shown to follow the shifted inverse-gamma distribution or a mixture of these distributions. The mixture parameters are estimated using the expectation-maximization algorithm. In addition, a method for the resolution- and individual ADP-dependent local analysis of neighbouring atoms has been designed. This method facilitates the detection of mismodelled atoms, heavy-metal atoms and disordered and/or incorrectly modelled ligands. Both global and local analyses can be used to detect errors in atomic models, thus helping in the (re)building, refinement and validation of macromolecular structures. This method can also serve as an additional validation tool during PDB deposition.

Non-native fold of the putative VPS39 zinc finger domain

Background: The multi-subunit homotypic fusion and vacuole protein sorting (HOPS) membrane-tethering complex is involved in regulating the fusion of late endosomes and autophagosomes with lysosomes in eukaryotes. The C-terminal regions of several HOPS components have been shown to be required for correct complex assembly, including the C-terminal really interesting new gene (RING) zinc finger domains of HOPS components VPS18 and VPS41. We sought to structurally characterise the putative C-terminal zinc finger domain of VPS39, which we hypothesised may be important for binding of VPS39 to cellular partners or to other HOPS components.

Methods: We recombinantly expressed, purified and solved the crystal structure of the proposed zinc-binding region of VPS39.

Results: In the structure, this region forms an anti-parallel β-hairpin that is incorporated into a homotetrameric eight-stranded β-barrel. However, the fold is stabilised by coordination of zinc ions by residues from the purification tag and an intramolecular disulphide bond between two predicted zinc ligands.

Conclusions: We solved the structure of the VPS39 C-terminal domain adopting a non-native fold. Our work highlights the risk of non-native folds when purifying small zinc-containing domains with hexahistidine tags. However, the non-native structure we observe may have implications for rational protein design.

Non-native fold of the putative VPS39 zinc finger domain

Background: The multi-subunit homotypic fusion and vacuole protein sorting (HOPS) membrane-tethering complex is involved in regulating the fusion of late endosomes and autophagosomes with lysosomes in eukaryotes. The C-terminal regions of several HOPS components have been shown to be required for correct complex assembly, including the C-terminal really interesting new gene (RING) zinc finger domains of HOPS components VPS18 and VPS41. We sought to structurally characterise the putative C-terminal zinc finger domain of VPS39, which we hypothesised may be important for binding of VPS39 to cellular partners or to other HOPS components. Methods: We recombinantly expressed, purified and solved the crystal structure of the proposed zinc-binding region of VPS39. Results: In the structure, this region forms an anti-parallel β-hairpin that is incorporated into a homotetrameric eight-stranded β-barrel. However, the fold is stabilised by coordination of zinc ions by residues from the purification tag and an intramolecular disulphide bond between two predicted zinc ligands. Conclusions: We solved the structure of the VPS39 C-terminal domain adopting a non-native fold. Our work highlights the risk of non-native folds when purifying small zinc-containing domains with hexahistidine tags. However, the non-native structure we observe may have implications for rational protein design.

Structures of virus-like capsids formed by the Drosophila neuronal Arc proteins

Arc, a neuronal gene critical for synaptic plasticity, originated through domestication of retrotransposon Gag genes and mediates intercellular mRNA transfer. We report high-resolution structures of retrovirus-like capsids formed by Drosophila dArc1 and dArc2 that have surface spikes and putative internal RNA-binding domains. These data demonstrate that virus-like capsid-forming properties of Arc are evolutionarily conserved and provide a structural basis for understanding their function in intercellular communication.

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.

Structural and functional characterization of the transcriptional regulator Rv3488 of Mycobacterium tuberculosis H37Rv

Rv3488 of Mycobacterium tuberculosis H37Rv has been assigned to the phenolic acid decarboxylase repressor (PadR) family of transcriptional regulators that play key roles in multidrug resistance and virulence of prokaryotes. The binding of cadmium, zinc, and several other metals to Rv3488 was discovered and characterized by isothermal titration calorimetery to be an exothermic process. Crystal structures of apo-Rv3488 and Rv3488 in complex with cadmium or zinc ions were determined by X-ray crystallography. The structure of Rv3488 revealed a dimeric protein with N-terminal winged-helix-turn-helix DNA-binding domains composed of helices α1, α2, α3, and strands β1 and β2, with the dimerization interface being formed of helices α4 and α1. The overall fold of Rv3488 was similar to PadR-s2 and metal sensor transcriptional regulators. In the crystal structure of Rv3488–Cd complex, two octahedrally coordinated Cd2+ ions were present, one for each subunit. The same sites were occupied by zinc ions in the structure of Rv3488–Zn, with two additional zinc ions complexed in one monomer. EMSA studies showed specific binding of Rv3488 with its own 30-bp promoter DNA. The functional role of Rv3488 was characterized by expressing the rv3488 gene under the control of hsp60 promoter in Mycobacterium smegmatis. Expression of Rv3488 increased the intracellular survival of recombinant M. smegmatis in murine macrophage cell line J774A.1 and also augmented its tolerance to Cd2+ ions. Overall, the studies show that Rv3488 may have transcription regulation and metal-detoxifying functions and its expression in M. smegmatis increases intracellular survival, perhaps by counteracting toxic metal stress.

Structure of the Cdc48 ATPase with its ubiquitin-binding cofactor Ufd1-Npl4

Many poly-ubiquitinated proteins are extracted from membranes or complexes by a conserved ATPase, called Cdc48 in yeast and p97/VCP in mammals, before proteasomal degradation¹. Each Cdc48 hexamer contains two stacked ATPase rings (D1 and D2) and six N-terminal (N) domains². Cdc48 binds various cofactors, including a heterodimer of Ufd1 and Npl4³. Here, we report structures of the Cdc48-Ufd1-Npl4 ATPase complex from Chaetomium thermophilum. Npl4 interacts through its UBX-like domain with a Cdc48 N domain, and uses two Zn²⁺-finger domains to anchor an enzymatically inactive Mpr1/Pad1 N-terminal (MPN) domain, homologous to domains found in several isopeptidases, to the top of the D1 ATPase ring. The MPN domain of Npl4 is located above Cdc48’s central pore, similarly to the MPN of the de-ubiquitinase Rpn11 in the proteasome⁴. Our results indicate that Npl4 is unique among Cdc48 cofactors, and suggest a mechanism for how poly-ubiquitinated substrates bind to and translocate into the ATPase.

ISOLDE: a physically realistic environment for model building into low-resolution electron-density maps

This paper introduces ISOLDE, a new software package designed to provide an intuitive environment for high-fidelity interactive remodelling/refinement of macromolecular models into electron-density maps. ISOLDE combines interactive molecular-dynamics flexible fitting with modern molecular-graphics visualization and established structural biology libraries to provide an immersive interface wherein the model constantly acts to maintain physically realistic conformations as the user interacts with it by directly tugging atoms with a mouse or haptic interface or applying/removing restraints. In addition, common validation tasks are accelerated and visualized in real time. Using the recently described 3.8 Å resolution cryo-EM structure of the eukaryotic minichromosome maintenance (MCM) helicase complex as a case study, it is demonstrated how ISOLDE can be used alongside other modern refinement tools to avoid common pitfalls of low-resolution modelling and improve the quality of the final model. A detailed analysis of changes between the initial and final model provides a somewhat sobering insight into the dangers of relying on a small number of validation metrics to judge the quality of a low-resolution model.

An allosteric binding site of the α7 nicotinic acetylcholine receptor revealed in a humanized acetylcholine-binding protein

Nicotinic acetylcholine receptors (nAChRs) belong to the family of pentameric ligand-gated ion channels and mediate fast excitatory transmission in the central and peripheral nervous systems. Among the different existing receptor subtypes, the homomeric α7 nAChR has attracted considerable attention because of its possible implication in several neurological and psychiatric disorders, including cognitive decline associated with Alzheimer's disease or schizophrenia. Allosteric modulators of ligand-gated ion channels are of particular interest as therapeutic agents, as they modulate receptor activity without affecting normal fluctuations of synaptic neurotransmitter release. Here, we used X-ray crystallography and surface plasmon resonance spectroscopy of α7-acetylcholine-binding protein (AChBP), a humanized chimera of a snail AChBP, which has 71% sequence similarity with the extracellular ligand-binding domain of the human α7 nAChR, to investigate the structural determinants of allosteric modulation. We extended previous observations that an allosteric site located in the vestibule of the receptor offers an attractive target for receptor modulation. We introduced seven additional humanizing mutations in the vestibule-located binding site of AChBP to improve its suitability as a model for studying allosteric binding. Using a fragment-based screening approach, we uncovered an allosteric binding site located near the β8-β9 loop, which critically contributes to coupling ligand binding to channel opening in human α7 nAChR. This work expands our understanding of the topology of allosteric binding sites in AChBP and, by extrapolation, in the human α7 nAChR as determined by electrophysiology measurements. Our insights pave the way for drug design strategies targeting nAChRs involved in ion channel-mediated disorders.

CheckMyMetal: a macromolecular metal-binding validation tool

Metals are essential in many biological processes, and metal ions are modeled in roughly 40% of the macromolecular structures in the Protein Data Bank (PDB). However, a significant fraction of these structures contain poorly modeled metal-binding sites. CheckMyMetal (CMM) is an easy-to-use metal-binding site validation server for macromolecules that is freely available at http://csgid.org/csgid/metal_sites. The CMM server can detect incorrect metal assignments as well as geometrical and other irregularities in the metal-binding sites. Guidelines for metal-site modeling and validation in macromolecules are illustrated by several practical examples grouped by the type of metal. These examples show CMM users (and crystallographers in general) problems they may encounter during the modeling of a specific metal ion.

AceDRG: a stereochemical description generator for ligands

The program AceDRG is designed for the derivation of stereochemical information about small molecules. It uses local chemical and topological environment-based atom typing to derive and organize bond lengths and angles from a small-molecule database: the Crystallography Open Database (COD). Information about the hybridization states of atoms, whether they belong to small rings (up to seven-membered rings), ring aromaticity and nearest-neighbour information is encoded in the atom types. All atoms from the COD have been classified according to the generated atom types. All bonds and angles have also been classified according to the atom types and, in a certain sense, bond types. Derived data are tabulated in a machine-readable form that is freely available from CCP4. AceDRG can also generate stereochemical information, provided that the basic bonding pattern of a ligand is known. The basic bonding pattern is perceived from one of the computational chemistry file formats, including SMILES, mmCIF, SDF MOL and SYBYL MOL2 files. Using the bonding chemistry, atom types, and bond and angle tables generated from the COD, AceDRG derives the `ideal' bond lengths, angles, plane groups, aromatic rings and chirality information, and writes them to an mmCIF file that can be used by the refinement program REFMAC5 and the model-building program Coot. Other refinement and model-building programs such as PHENIX and BUSTER can also use these files. AceDRG also generates one or more coordinate sets corresponding to the most favourable conformation(s) of a given ligand. AceDRG employs RDKit for chemistry perception and for initial conformation generation, as well as for the interpretation of SMILES strings, SDF MOL and SYBYL MOL2 files.