Factors Governing Metal−Ligand Distances and Coordination Geometries of Metal Complexes

The Role of Axial Ligand in Determining the Mg2+/TM2+ (TM = Fe, Mn, Cu, Zn) Selectivity in Chlorophyll

Life on Earth is sustained due to plants, algae, and cyanobacteria's exquisite ability to store energy from light and convert it into a form that can be utilized by living organisms. In this process, the compound of immense significance, chlorophyll, actively participates. It consists of a chlorin unit with different substituents, depending on the type of chlorophyll molecule. The main element, however, is the magnesium cation bound to the center of the tetrapyrrole macrocycle. One additional axial ligand of diverse nature is often attached to the cation to complement the structure of the complex. The present study aims at elucidating the role of the additional ligand for determining the metal selectivity in chlorophyll molecules by implementing a widely applied methodology based on Density Functional Theory. The obtained results provide information about the thermodynamic outcome of the competition between Mg2+ and a series of transition biogenic metal cations such as Fe2+, Mn2+, Cu2+, and Zn2+ in model chlorophyll metal binding sites in two environments of different polarity, delineating the key factors that contribute to the process in the greatest extent. The results of calculations are discussed in light of known experimental data and hence shed light on the metal selectivity process in this system.

Flexibility-Induced Collective Behavior Drives Symmetry Breaking in Discrimination of Undesired Ions

Structure flexibility is essential for the biological function of proteins. At the same time, many proteins need to discriminate ligands with subtle differences, with one example being ion selectivity. Investigating the mechanisms by which flexible proteins achieve such precise discrimination is crucial for advancing our understanding of their functions. In this work, we study transporter KCC4, which undergoes continuous conformation changes during ion transport and can realize K+ over Na+ selectivity. Our findings reveal that the center of the binding site no longer represents a stable equilibrium for the undesired Na+, and its binding mode exhibits bifurcation. Interestingly, protein conformation fluctuation can induce collective behavior throughout the entire binding region, which contributes to this bifurcation. Thus, the symmetry of the binding mode decreases from the inherent T d symmetry to a C2v symmetry, and the binding stability of Na+ is largely reduced. A similar phenomenon is observed in a GPCR, β2-AR, where a less favored ligand forms a biased binding mode with reduced stability. The mechanism underlying the selectivity in such flexible regions could be interpreted as spontaneous symmetry breaking, which may represent a general mechanism by which flexible proteins achieve efficient ligand discrimination.

Direct Identification of α-Bisabolol Enantiomers in an Essential Oil Using a Combined Ion Mobility-Mass Spectrometry/Quantum Chemistry Approach

Enantiomer-specific identification of chiral molecules in natural extracts is a challenging task, as many routine analytical techniques fail to provide selectivity in multicomponent mixtures. Here we describe an alternative approach, based on the combination of ion mobility-mass spectrometry (IM-MS) and quantum chemistry (QM), for the direct enantiomers differentiation in crude essential oils. The identification of α-bisabolol enantiomers contained in the raw essential oil (EO) from the Corsican Xanthium italicum fruits is reported as a proof-of-concept. Accordingly, IM-MS experiments performed in Ag+-doped methanol revealed the presence of both (+)- and (-)-α-bisabolol in the EO, while molecular simulations provided the structures of the two α-bisabolol enantiomer silver(I) adducts.

Discovery, Characterization, and Structure of a Cell Active PAD2 Inhibitor Acting through a Novel Allosteric Mechanism

Peptidyl arginine deiminases (PADs) are important enzymes in many diseases, especially those involving inflammation and autoimmunity. Despite many years of effort, developing isoform-specific inhibitors has been a challenge. We describe herein the discovery of a potent, noncovalent PAD2 inhibitor, with selectivity over PAD3 and PAD4, from a DNA-encoded library. The biochemical and biophysical characterization of this inhibitor and two noninhibitory binders indicated a novel, Ca2+ competitive mechanism of inhibition. This was confirmed via X-ray crystallographic analysis. Finally, we demonstrate that this inhibitor selectively inhibits PAD2 in a cellular context.

Disorder-to-order active site capping regulates the rate-limiting step of the inositol pathway

Myo-inositol-1-phosphate synthase (MIPS) catalyzes the NAD+-dependent isomerization of glucose-6-phosphate (G6P) into inositol-1-phosphate (IMP), controlling the rate-limiting step of the inositol pathway. Previous structural studies focused on the detailed molecular mechanism, neglecting large-scale conformational changes that drive the function of this 240 kDa homotetrameric complex. In this study, we identified the active, endogenous MIPS in cell extracts from the thermophilic fungus Thermochaetoides thermophila. By resolving the native structure at 2.48 Å (FSC = 0.143), we revealed a fully populated active site. Utilizing 3D variability analysis, we uncovered conformational states of MIPS, enabling us to directly visualize an order-to-disorder transition at its catalytic center. An acyclic intermediate of G6P occupied the active site in two out of the three conformational states, indicating a catalytic mechanism where electrostatic stabilization of high-energy intermediates plays a crucial role. Examination of all isomerases with known structures revealed similar fluctuations in secondary structure within their active sites. Based on these findings, we established a conformational selection model that governs substrate binding and eventually inositol availability. In particular, the ground state of MIPS demonstrates structural configurations regardless of substrate binding, a pattern observed across various isomerases. These findings contribute to the understanding of MIPS structure-based function, serving as a template for future studies targeting regulation and potential therapeutic applications.

Spontaneous ligand loss by soft landed [Ni(bpy)3]2+ ions on perfluorinated self-assembled monolayer surfaces

Transition metal (TM) complexes are widely used in catalysis, photochemical energy conversion, and sensing. Understanding factors that affect ligand loss from TM complexes at interfaces is important both for generating catalytically-active undercoordinated TM complexes and for controlling the degradation pathways of photosensitizers and photoredox catalysts. Herein, we demonstrate that well-defined TM complexes prepared on surfaces using ion soft landing undergo substantial structural rearrangements resulting in ligand loss and formation of both stable and reactive undercoordinated species. We employ nickel bipyridine (Ni-bpy) cations as a model system and explore their structural reorganization on surfaces using a combination of experimental and computational approaches. The controlled preparation of surface layers by mass-selected deposition of [Ni(bpy)3]2+ cations provides insights into the chemical reactivity of these species on surfaces. Both surface characterization using mass spectrometry and electronic structure calculations using density functional theory (DFT) indicate that [Ni(bpy)3]2+ undergoes a substantial geometry distortion on surfaces in comparison with its gas-phase structure. This distortion reduces the ligand binding energy and facilitates the formation of the undercoordinated [Ni(bpy)2]2+. Additionally, charge reduction by the soft landed [Ni(bpy)3]2+ facilitates ligand loss. We observe that ligand loss is inhibited by co-depositing [Ni(bpy)3]2+ with a stable anion such as closo-dodecaborate dianion, [B12F12]2-. The strong electrostatic interaction between [Ni(bpy)3]2+ and [B12F12]2- diminishes the distortion of the cation due to interactions with the surface. This interaction stabilizes the soft landed cation by reducing the extent of charge reduction and its structural reorganization. Overall, this study shows the intricate interplay of charge state, ion surface interactions, and stabilization by counterions on the structure and reactivity of metal complexes on surfaces. The combined experimental and computational approach used in this study offers detailed insights into factors that affect the integrity and stability of active species relevant to energy production and catalysis.

Diverse crystalline protein scaffolds through metal-dependent polymorphism

As protein crystals are increasingly finding diverse applications as scaffolds, controlled crystal polymorphism presents a facile strategy to form crystalline assemblies with controllable porosity with minimal to no protein engineering. Polymorphs of consensus tetratricopeptide repeat proteins with varying porosity were obtained through co-crystallization with metal salts, exploiting the innate metal ion geometric requirements. A single structurally exposed negative amino acid cluster was responsible for metal coordination, despite the abundance of negatively charged residues. Density functional theory calculations showed that while most of the crystals were the most thermodynamically stable assemblies, some were kinetically trapped states. Thus, crystalline porosity diversity is achieved and controlled with metal coordination, opening a new scope in the application of proteins as biocompatible protein-metal-organic frameworks (POFs). In addition, metal-dependent polymorphic crystals allow direct comparison of metal coordination preferences.

Utilising a Proton-Responsive 1,8-Naphthyridine Ligand for the Synthesis of Bimetallic Palladium and Platinum Complexes

We present four proton-responsive palladium and platinum complexes, [MCl2 (R PONNHO)] (M=Pd, Pt; R=i Pr, t Bu) synthesised by complexation of PdCl2 or PtCl2 (COD) with the 1,8-naphthyridine ligand R PONNHO. Deprotonation of [MCl2 (tBu PONNHO)] switches ligand coordination from mono- to dinucleating, offering a synthetic pathway to bimetallic PdII and PtII complexes [M2 Cl2 (tBu PONNO)2 ]. Two-electron reduction gives planar MI -MI complexes [M2 (tBu PONNO)2 ] (M=Pd, Pt) containing a metal-metal bond. In contrast to the related nickel system that forms a metallophosphorane [Ni2 (tBu PONNOPONNO)], an unusual phosphinite binding mode is observed in [M2 (tBu PONNO)2 ] containing close phosphinite-naphthyridinone P⋅⋅⋅O interactions, which is investigated spectroscopically, crystallographically and computationally. The presented proton-responsive and structurally-responsive R PONNHO and bimetallic R PONNO complexes offer a novel platform for future explorations of metal-ligand and metal-metal cooperativity with palladium and platinum.

Reanalysis of a μ opioid receptor crystal structure reveals a covalent adduct with BU72

BACKGROUND

The first crystal structure of the active μ opioid receptor (μOR) exhibited several unexplained features. The ligand BU72 exhibited many extreme deviations from ideal geometry, along with unexplained electron density. I previously showed that inverting the benzylic configuration resolved these problems, establishing revised stereochemistry of BU72 and its analog BU74. However, another problem remains unresolved: additional unexplained electron density contacts both BU72 and a histidine residue in the N-terminus, revealing the presence of an as-yet unidentified atom.

RESULTS

These short contacts and uninterrupted density are inconsistent with non-covalent interactions. Therefore, BU72 and μOR form a covalent adduct, rather than representing two separate entities as in the original model. A subsequently proposed magnesium complex is inconsistent with multiple lines of evidence. However, oxygen fits the unexplained density well. While the structure I propose is tentative, similar adducts have been reported previously in the presence of reactive oxygen species. Moreover, known sources of reactive oxygen species were present: HEPES buffer, nickel ions, and a sequence motif that forms redox-active nickel complexes. This motif contacts the unexplained density. The adduct exhibits severe strain, and the tethered N-terminus forms contacts with adjacent residues. These forces, along with the nanobody used as a G protein substitute, would be expected to influence the receptor conformation. Consistent with this, the intracellular end of the structure differs markedly from subsequent structures of active μOR bound to Gi protein.

CONCLUSIONS

Later Gi-bound structures are likely to be more accurate templates for ligand docking and modelling of active G protein-bound μOR. The possibility of reactions like this should be considered in the choice of protein truncation sites and purification conditions, and in the interpretation of excess or unexplained density.

High-resolution crystal structure of the Mu8.1 conotoxin from Conus mucronatus

Marine cone snails produce a wealth of peptide toxins (conotoxins) that bind their molecular targets with high selectivity and potency. Therefore, conotoxins constitute valuable biomolecular tools with a variety of biomedical purposes. The Mu8.1 conotoxin from Conus mucronatus is the founding member of the newly identified saposin-like conotoxin class of conotoxins and has been shown to target Cav2.3, a voltage-gated calcium channel. Two crystal structures have recently been determined of Mu8.1 at 2.3 and 2.1 Å resolution. Here, a high-resolution crystal structure of Mu8.1 was determined at 1.67 Å resolution in the high-symmetry space group I4122. The asymmetric unit contained one molecule, with a symmetry-related molecule generating a dimer equivalent to that observed in the two previously determined structures. The high resolution allows a detailed atomic analysis of a water-filled cavity buried at the dimer interface, revealing a tightly coordinated network of waters that shield a lysine residue (Lys55) with a predicted unusually low side-chain pKa value. These findings are discussed in terms of a potential functional role of Lys55 in target interaction.

Mechanisms of inward transmembrane proton translocation

Proton transport is indispensable for cell life. It is believed that molecular mechanisms of proton movement through different types of proton-conducting molecules have general universal features. However, elucidation of such mechanisms is a challenge. It requires true-atomic-resolution structures of all key proton-conducting states. Here we present a comprehensive function-structure study of a light-driven bacterial inward proton pump, xenorhodopsin, from Bacillus coahuilensis in all major proton-conducting states. The structures reveal that proton translocation is based on proton wires regulated by internal gates. The wires serve as both selectivity filters and translocation pathways for protons. The cumulative results suggest a general concept of proton translocation. We demonstrate the use of serial time-resolved crystallography at a synchrotron source with sub-millisecond resolution for rhodopsin studies, opening the door for principally new applications. The results might also be of interest for optogenetics since xenorhodopsins are the only alternative tools to fire neurons.

Benchmarking biomolecular force field-based Zn2+ for mono- and bimetallic ligand binding sites

Zn²⁺ is one of the most versatile biologically available metal ions, but accurate modeling of Zn²⁺ -containing metalloproteins at the biomolecular force field level can be challenging. Since most Zn²⁺ models are parameterized in bulk solvent, in-depth knowledge about their performance in a protein environment is limited. Thus, we systematically investigate here the behavior of non-polarizable Zn²⁺ models for their ability to reproduce experimentally determined metal coordination and ligand binding in metalloproteins. The benchmarking is performed in challenging environments, including mono- (carbonic anhydrase II) and bimetallic (metallo-β-lactamase VIM-2) ligand binding sites. We identify key differences in the performance between the Zn²⁺ models with regard to the preferred ligating atoms (charged/non-charged), attraction of water molecules, and the preferred coordination geometry. Based on these results, we suggest suitable simulation conditions for varying Zn²⁺ site geometries that could guide the further development of biomolecular Zn²⁺ models.

Lanthanides as Calcium Mimetic Species in Calcium-Signaling/Buffering Proteins: The Effect of Lanthanide Type on the Ca2+/Ln3+ Competition

Lanthanides, the 14 4f-block elements plus Lanthanum, have been extensively used to study the structure and biochemical properties of metalloproteins. The characteristics of lanthanides within the lanthanide series are similar, but not identical. The present research offers a systematic investigation of the ability of the entire Ln3+ series to substitute for Ca2+ in biological systems. A well-calibrated DFT/PCM protocol is employed in studying the factors that control the metal selectivity in biological systems by modeling typical calcium signaling/buffering binding sites and elucidating the thermodynamic outcome of the competition between the “alien” La3+/Ln3+ and “native” Ca2+, and La3+ − Ln3+ within the lanthanide series. The calculations performed reveal that the major determinant of the Ca2+/Ln3+ selectivity in calcium proteins is the net charge of the calcium binding pocket; the more negative the charge, the higher the competitiveness of the trivalent Ln3+ with respect to its Ca2+ contender. Solvent exposure of the binding site also influences the process; buried active centers with net charge of −4 or −3 are characterized by higher Ln3+ over Ca2+ selectivity, whereas it is the opposite for sites with overall charge of −1. Within the series, the competition between La3+ and its fellow lanthanides is determined by the balance between two competing effects: electronic (favoring heavier lanthanides) and solvation (generally favoring the lighter lanthanides).

CMM-An enhanced platform for interactive validation of metal binding sites

Metal ions bound to macromolecules play an integral role in many cellular processes. They can directly participate in catalytic mechanisms or be essential for the structural integrity of proteins and nucleic acids. However, their unique nature in macromolecules can make them difficult to model and refine, and a substantial portion of metal ions in the PDB are misidentified or poorly refined. CheckMyMetal (CMM) is a validation tool that has gained widespread acceptance as an essential tool for researchers working on metal-macromolecule complexes. CMM can be used during structure determination or to validate metal binding sites in structural models within the PDB. The functionalities of CMM have recently been greatly enhanced and provide researchers with additional information that can guide modeling decisions. The new version of CMM shows metals in the context of electron density maps and allows for on-the-fly refinement of metal binding sites. The improvements should increase the reproducibility of biomedical research. The web server is available at https://cmm.minorlab.org.

Probing and Manipulating Noncovalent Interactions in Functional Polymeric Systems

Noncovalent interactions, which usually feature tunable strength, reversibility, and environmental adaptability, have been recognized as driving forces in a variety of biological and chemical processes, contributing to the recognition between molecules, the formation of molecule clusters, and the establishment of complex structures of macromolecules. The marriage of noncovalent interactions and conventional covalent polymers offers the systems novel mechanical, physicochemical, and biological properties, which are highly dependent on the binding mechanisms of the noncovalent interactions that can be illuminated via quantification. This review systematically discusses the nanomechanical characterization of typical noncovalent interactions in polymeric systems, mainly through direct force measurements at microscopic, nanoscopic, and molecular levels, which provide quantitative information (e.g., ranges, strengths, and dynamics) on the binding behaviors. The fundamental understandings of intermolecular and interfacial interactions are then correlated to the macroscopic performances of a series of noncovalently bonded polymers, whose functions (e.g., stimuli-responsiveness, self-healing capacity, universal adhesiveness) can be customized through the manipulation of the noncovalent interactions, providing insights into the rational design of advanced materials with applications in biomedical, energy, environmental, and other engineering fields.

Stereocontrolled Synthesis of C20S-C26 and C20R-C26 Fragments of Amphidinolide L

Amphidinolide L is a cytotoxic macrolide isolated from marine symbiotic dinoflagellates of the genus Amphidinium. While its planar structure and the absolute stereochemistry of the C21-C26 part have been determined, six stereocenters have remained unassigned. Aiming at structure determination, we have developed a synthetic route to the C20S-C26 and C20R-C26 fragments via the Li-mediated stereocontrolled aldol reaction. Two aldehydes, 16 with the C22-hydroxy group and 19 with the C22-TES ether, were synthesized from lactone 4. The aldol reactions using the Li-enolate of 4-methyl-2-pentanone in THF provided the C20S-C26 fragment 20 from 16 and a 1:3.5 mixture of the C20-C26 fragment 22 favoring the C20R-isomer. Mechanistic studies based on an extensive search of transition states in explicit solvents indicated that the C20S-isomer would be generated via a tri-solvated transition state, while the C20R-isomer would be formed via a di-solvated transition state. The calculation emphasizes the importance of the coordination network as a higher-order complex composed of solvent molecules, aldehyde, enolate, and Li atoms in the reaction of 16 to minimize steric interactions but maximize the stabilizing effect by the coordination of solvents. The presence of the rotationally free aldehyde in the reaction of 19 results in moderate diastereoselectivity.

Molecular and in vivo studies of a glutamate-class prolyl-endopeptidase for coeliac disease therapy

The digestion of gluten generates toxic peptides, among which a highly immunogenic proline-rich 33-mer from wheat α-gliadin, that trigger coeliac disease. Neprosin from the pitcher plant is a reported prolyl endopeptidase. Here, we produce recombinant neprosin and its mutants, and find that full-length neprosin is a zymogen, which is self-activated at gastric pH by the release of an all-β pro-domain via a pH-switch mechanism featuring a lysine plug. The catalytic domain is an atypical 7+8-stranded β-sandwich with an extended active-site cleft containing an unprecedented pair of catalytic glutamates. Neprosin efficiently degrades both gliadin and the 33-mer in vitro under gastric conditions and is reversibly inactivated at pH > 5. Moreover, co-administration of gliadin and the neprosin zymogen at the ratio 500:1 reduces the abundance of the 33-mer in the small intestine of mice by up to 90%. Neprosin therefore founds a family of eukaryotic glutamate endopeptidases that fulfils requisites for a therapeutic glutenase.

Celiac disease is characterized by intolerance to gluten, a cereal protein. Here, the authors show that neprosin, a glutamate peptidase from the pitcher plant, efficiently cleaves gluten components under physiological conditions in vitro and in the gut of mice.

Factors allowing small monovalent Li+ to displace Ca2+ in proteins

Because Li⁺ and Ca²⁺ differ in both charge and size, the possibility that monovalent Li⁺ could dislodge the bulkier, divalent Ca²⁺ in Ca²⁺ proteins had not been considered. However, our recent density functional theory/continuum dielectric calculations predicted that Li⁺ could displace the native Ca²⁺ from the C2 domain of cytosolic PKCα/γ. This would reduce electrostatic interactions between the Li⁺-bound C2 domain and the membrane, consistent with experimental studies showing that Li⁺ can inhibit the translocation of cytoplasmic PKC to membranes. Besides the trinuclear Ca²⁺-site in the PKCα/γ C2 domain, it is not known whether other Ca²⁺-sites in human proteins may be susceptible to Li⁺ substitution. Furthermore, it is unclear what factors determine the outcome of the competition between divalent Ca²⁺ and monovalent Li⁺. Here we show that the net charge of residues in the first and second coordination shell is a key determinant of the selectivity for divalent Ca²⁺ over monovalent Li⁺ in proteins: neutral/anionic Ca²⁺-carboxylate sites are protected against Li⁺ attack. They are further protected by outer-shell Asp^-/Glu^- and the protein matrix rigidifying the Ca²⁺-site or limiting water entry. In contrast, buried, cationic Ca²⁺-sites surrounded by Arg⁺/Lys⁺, which are found in the C2 domains of PKCα/γ, as well as certain synaptotagmins, are prone to Li⁺ attack.

Metal content and kinetic properties of yeast RNA lariat debranching enzyme Dbr1

In eukaryotic cells, intron lariats produced by the spliceosome contain a 2'5' phosphodiester linkage. The RNA lariat debranching enzyme, Dbr1, is the only enzyme known to hydrolyze this bond. Dbr1 is a member of the metallophosphoesterase (MPE) family of enzymes, and recent X-ray crystal structures and biochemistry data demonstrate that Dbr1 from Entamoeba histolytica uses combinations of Mn²⁺, Zn²⁺, and Fe²⁺ as enzymatic cofactors. Here, we examine the kinetic properties and metal dependence of the Dbr1 homolog from Saccharomyces cerevisiae (yDbr1). Elemental analysis measured stoichiometric quantities of Fe and Zn in yDbr1 purified following heterologous expression E. coli We analyzed the ability of Fe²⁺, Zn²⁺, and Mn²⁺ to reconstitute activity in metal-free apoenzyme. Purified yDbr1 was highly active, turning over substrate at 5.6 sec^-1, and apo-yDbr1 reconstituted with Fe²⁺ was the most active species, turning over at 9.2 sec^-1 We treated human lymphoblastoid cells with the iron-chelator deferoxamine and measured a twofold increase in cellular lariats. These data suggest that Fe is an important biological cofactor for Dbr1 enzymes.

A Nd6 molecular butterfly: a unique all-in-one material for SMM, MCE and maiden photosensitized opto-electronic device fabrication

Besides iron, ironically neodymium (Nd) is the most ubiquitously used metal for magnetic purposes, even among the lanthanides, when it comes to the field of molecular magnetism, yet it ranks among the least studied metals. However, strong apathy towards this magnetic lanthanide means that vital information will be missed, which is required for the advancement of the subject. Herein, we have successfully demonstrated the usefulness of a hexanuclear neodymium complex as a magnetic material, and also in electronic device fabrication. A {NdIII6} cage with an aesthetically pleasing butterfly topology was synthesized using a rather non-conventional N-rich pyridyl-pyrazolyl based ligand. The cage shows single molecule magnet (SMM) properties, with an effective energy barrier, U_eff, value of 3.4 K and relaxation time, τ₀, of 3.1 × 10^-4 s, originating from an unusual occurrence of metal centres with different coordination environments. Furthermore, magnetic studies reveal significant cyrogenic magnetic cooling, with a magnetic entropy change of 8.28 J kg^-1 K^-1 at 5 T and 3 K. To the best of our knowledge, the titular compound is the only example of a Nd-complex that exhibits concomitant magnetocaloric effect (MCE) and SMM properties. Complete active space self-consistent field (CASSCF) calculations were carried out to shed light on the origin of the magnetic anisotropy and magnetic relaxation of the compound. The same uniqueness is also true for the first electronic investigation carried out on the Nd complex. The maiden electronic device fabricated using the Nd complex shows an interesting intertwining of electronic and optical features, which contribute towards its improved photosensitized optoelectronic data.

Divalent-Metal-Ion Selectivity of the CRISPR-Cas System-Associated Cas1 Protein: Insights from Classical Molecular Dynamics Simulations and Electronic Structure Calculations

CRISPR-associated protein 1 (Cas1) is a universally conserved essential metalloenzyme of the clustered regularly interspaced short palindromic repeat (CRISPR) immune system of prokaryotes (bacteria, archaea) that can cut and integrate a part of viral DNA to its host genome with the help of other proteins. The integrated DNA acts as a memory of viral infection, which can be transcribed to RNA and stop future infection by recognition (based on the RNA/DNA complementarity principle) followed by protein-mediated degradation of the viral DNA. It has been proposed that the presence of a single manganese (Mn²⁺) ion in a conserved divalent-metal-ion binding pocket (key residues: E190, H254, D265, D268) of Cas1 is crucial for its function. Cas1-mediated DNA degradation was proposed to be hindered by metal substitution, metal chelation, or mutation of the binding pocket residues. Cas1 is active toward dsDNA degradation with both Mn²⁺ and Mg²⁺. X-ray structures of Cas1 revealed an intricate atomic interaction network of the divalent-metal-ion binding pocket and opened up the possibility of modeling related metal ions (viz., Mg²⁺, Ca²⁺) in the binding pocket of wild-type (WT) and mutated Cas1 proteins for computational analysis, which includes (1) quantitative estimation of the energetics of the divalent-metal-ion preference and (2) exploring the structural and dynamical aspects of the protein in response to divalent-metal-ion substitution or amino acid mutation. Using the X-ray structure of the Cas1 protein from Pseudomonas aeruginosa as a template (PDB 3GOD), we performed (∼2.23 μs) classical molecular dynamics (MD) simulations to compare structural and dynamical differences between Mg²⁺- and Ca²⁺-bound binding pockets of wild-type (WT) and mutant (E190A, H254A, D265A, D268A) Cas1. Furthermore, reduced binding pocket models were generated from X-ray and molecular dynamics (MD) trajectories, and the resulting structures were subjected to quantum chemical calculations. Results suggest that Cas1 prefers Mg²⁺ binding relative to Ca²⁺ and the preference is the strongest for WT and the weakest for the D268A mutant. Quantum chemical calculations indicate that Mn²⁺ is the most preferred relative to both Mg²⁺ and Ca²⁺ in the wild-type and mutant Cas1. Substitution of Mg²⁺ by Ca²⁺ does not alter the interaction network between Cas1 and the divalent metal ion but increases the wetness of the binding pocket by introducing a single water molecule in the first coordination shell of the latter. The strength of metal-ion preference (Mg²⁺ versus Ca²⁺) seems to be dependent on the solvent accessibility of the divalent-metal-ion binding pocket, strongest for wild-type Cas1 (in which the metal-ion binding pocket is dry, which includes two water molecules) and the weakest for the D268A mutant (in which the metal-ion binding pocket is wet, which includes four water molecules).

Spectroscopic and computational investigations of organometallic complexation of group 12 transition metals by methanobactins from Methylocystis sp. SB2

Methanotrophic bacteria catalyze the aerobic oxidation of methane to methanol using Cu-containing enzymes, thereby exerting a modulating influence on the global methane cycle. To facilitate the acquisition of Cu ions, some methanotrophic bacteria secrete small modified peptides known as "methanobactins," which strongly bind Cu and function as an extracellular Cu recruitment relay, analogous to siderophores and Fe. In addition to Cu, methanobactins form complexes with other late transition metals, including the Group 12 transition metals Zn, Cd, and Hg, although the interplay among solution-phase configurations, metal interactions, and the spectroscopic signatures of methanobactin-metal complexes remains ambiguous. In this study, the complexation of Zn, Cd, and Hg by methanobactin from Methylocystis sp. strain SB2 was studied using a combination of absorbance, fluorescence, extended x-ray absorption fine structure (EXAFS) spectroscopy, and time-dependent density functional theory (TD-DFT) calculations. We report changes in sample absorbance and fluorescence spectral dynamics, which occur on a wide range of experimental timescales and characterize a clear stoichiometric complexation dependence. Mercury L3-edge EXAFS and TD-DFT calculations suggest a linear model for HgS coordination, and TD-DFT suggests a tetrahedral model for Zn2+ and Cd2+. We observed an enhancement in the fluorescence of methanobactin upon interaction with transition metals and propose a mechanism of complexation-hindered isomerization drawing inspiration from the wild-type Green Fluorescent Protein active site. Collectively, our results represent the first combined computational and experimental spectroscopy study of methanobactins and shed new light on molecular interactions and dynamics that characterize complexes of methanobactins with Group 12 transition metals.

Calcium in Signaling: Its Specificity and Vulnerabilities toward Biogenic and Abiogenic Metal Ions

Divalent calcium ion (Ca²⁺) plays an indispensable role as a second messenger in a myriad of signal transduction processes. Of utmost importance for the faultless functioning of calcium-modulated signaling proteins is their binding selectivity of the native metal cation over rival biogenic/abiogenic metal ion contenders in the intra/extracellular fluids. In this Perspective, we summarize recent findings on the competition between the cognate Ca²⁺ and other biogenic or abiogenic divalent cations for binding to Ca²⁺-signaling proteins or organic cofactors. We describe the competition between the two most abundant intracellular biogenic metal ions (Mg²⁺ and Ca²⁺) for Ca²⁺-binding sites in signaling proteins, followed by the rivalry between native Ca²⁺ and "therapeutic" Li⁺ as well as "toxic" Pb²⁺. We delineate the key factors governing the rivalry between the native and non-native cations in proteins and highlight key implications for the biological performance of the respective proteins/organic cofactors.

Computational Discovery of Transition-metal Complexes: From High-throughput Screening to Machine Learning

Transition-metal complexes are attractive targets for the design of catalysts and functional materials. The behavior of the metal-organic bond, while very tunable for achieving target properties, is challenging to predict and necessitates searching a wide and complex space to identify needles in haystacks for target applications. This review will focus on the techniques that make high-throughput search of transition-metal chemical space feasible for the discovery of complexes with desirable properties. The review will cover the development, promise, and limitations of "traditional" computational chemistry (i.e., force field, semiempirical, and density functional theory methods) as it pertains to data generation for inorganic molecular discovery. The review will also discuss the opportunities and limitations in leveraging experimental data sources. We will focus on how advances in statistical modeling, artificial intelligence, multiobjective optimization, and automation accelerate discovery of lead compounds and design rules. The overall objective of this review is to showcase how bringing together advances from diverse areas of computational chemistry and computer science have enabled the rapid uncovering of structure-property relationships in transition-metal chemistry. We aim to highlight how unique considerations in motifs of metal-organic bonding (e.g., variable spin and oxidation state, and bonding strength/nature) set them and their discovery apart from more commonly considered organic molecules. We will also highlight how uncertainty and relative data scarcity in transition-metal chemistry motivate specific developments in machine learning representations, model training, and in computational chemistry. Finally, we will conclude with an outlook of areas of opportunity for the accelerated discovery of transition-metal complexes.

Trinuclear Calcium Site in the C2 Domain of PKCα/γ Is Prone to Lithium Attack

Lithium (Li⁺) is the first-line therapy for bipolar disorder and a candidate drug for various diseases such as amyotrophic lateral sclerosis, multiple sclerosis, and stroke. Despite being the captivating subject of many studies, the mechanism of lithium's therapeutic action remains unclear. To date, it has been shown that Li⁺ competes with Mg²⁺ and Na⁺ to normalize the activity of inositol and neurotransmitter-related signaling proteins, respectively. Furthermore, Li⁺ may co-bind with Mg²⁺-loaded adenosine or guanosine triphosphate to alter the complex's susceptibility to hydrolysis and mediate cellular signaling. Bipolar disorder patients exhibit abnormally high cytosolic Ca²⁺ levels and protein kinase C (PKC) hyperactivity that can be downregulated by long-term Li⁺ treatment. However, the possibility that monovalent Li⁺ could displace the bulkier divalent Ca²⁺ and inhibit PKC activity has not been considered. Here, using density functional theory calculations combined with continuum dielectric methods, we show that Li⁺ may displace the native dication from the positively charged trinuclear site in the C2 domain of cytosolic PKCα/γ. This would affect the membrane-docking ability of cytosolic PKCα/γ and reduce the abnormally high membrane-associated active PKCα/γ levels, thus downregulating the PKC hyperactivity found in bipolar patients.

Evaluating the Multifunctionality of a New Modulator of Zinc-Induced Aβ Aggregation Using a Novel Computational Approach

The high concentration of zinc metal ions in Aβ aggregations is one of the most cited hallmarks of Alzheimer's disease (AD), and several substantial pieces of evidence emphasize the key role of zinc metal ions in the pathogenesis of AD. In this study, while designing a multifunctional peptide for simultaneous targeting Aβ aggregation and chelating the zinc metal ion, a novel and comprehensive approach is introduced for evaluating the multifunctionality of a multifunctional drugs based on computational methods. The multifunctional peptide consists of inhibitor and chelator domains, which are included in the C-terminal hydrophobic region of Aβ, and the first four amino acids of human albumin. The ability of the multifunctional peptide in zinc ion chelation has been investigated using molecular dynamics (MD) simulations of the peptide-zinc interaction for 300 ns, and Bennett's acceptance ratio (BAR) method has been used to accurately calculate the chelation free energy. Data analysis demonstrates that the peptide chelating domain can be stably linked to the zinc ion. Besides, the introduced method used for evaluating chelation and calculating the free energy of peptide binding to zinc ions was successfully validated by comparison with previous experimental and theoretical published data. The results indicate that the multifunctional peptide, coordinating with the zinc metal ion, can be effective in Aβ inhibition by preserving the native helical structure of the Aβ₄₂ monomer as well as disrupting the β-sheet structure of Aβ₄₂ aggregates. Detailed assessments of the Aβ₄₂-peptide interactions elucidate that the inhibition of Aβ is achieved by considerable hydrophobic interactions and hydrogen bonding between the multifunctional peptide and the hydrophobic Aβ regions, along with interfering in stable bridges formed inside the Aβ aggregate.

DeepCys: Structure-based multiple cysteine function prediction method trained on deep neural network: Case study on domains of unknown functions belonging to COX2 domains

Cysteine (Cys) is the most reactive amino acid participating in a wide range of biological functions. In-silico predictions complement the experiments to meet the need of functional characterization. Multiple Cys function prediction algorithm is scarce, in contrast to specific function prediction algorithms. Here we present a deep neural network-based multiple Cys function prediction, available on web-server (DeepCys) (https://deepcys.herokuapp.com/). DeepCys model was trained and tested on two independent datasets curated from protein crystal structures. This prediction method requires three inputs, namely, PDB identifier (ID), chain ID and residue ID for a given Cys and outputs the probabilities of four cysteine functions, namely, disulphide, metal-binding, thioether and sulphenylation and predicts the most probable Cys function. The algorithm exploits the local and global protein properties, like, sequence and secondary structure motifs, buried fractions, microenvironments and protein/enzyme class. DeepCys outperformed most of the multiple and specific Cys function algorithms. This method can predict maximum number of cysteine functions. Moreover, for the first time, explicitly predicts thioether function. This tool was used to elucidate the cysteine functions on domains of unknown functions belonging to cytochrome C oxidase subunit-II like transmembrane domains. Apart from the web-server, a standalone program is also available on GitHub (https://github.com/vam-sin/deepcys).

Second Coordination Sphere Effects in an Evolved Ru Complex Based on Highly Adaptable Ligand Results in Rapid Water Oxidation Catalysis

A new Ru complex containing the deprotonated 2,2':6',2''-terpyridine-6,6''-diphosphonic acid (H₄tPa) and pyridine (py) of general formula [Ru^II(H₃tPa-κ-N³O)(py)₂]⁺, 2⁺, has been prepared and thoroughly characterized by means of spectroscopic and electrochemical techniques, X-ray diffraction analysis, and density functional theory (DFT) calculations. Complex 2⁺ presents a dynamic behavior in the solution that involves the synchronous coordination and the decoordination of the dangling phosphonic groups of the tPa^4- ligand. However, at oxidation state IV, complex 2⁺ becomes seven coordinated with the two phosphonic groups now bonded to the metal center. Further, at this oxidation state at neutral and basic pH, the Ru complex undergoes the coordination of an exogenous OH^- group from the solvent that leads to an intramolecular aromatic O atom insertion into the CH bond of one of the pyridyl groups, forming the corresponding phenoxo-phosphonate Ru complex [Ru^III(tPaO-κ-N²O_PO_C)(py)₂]^2-, 4^2-, where tPaO^5- is the 3-(hydroxo-[2,2':6',2''-terpyridine]-6,6''-diyl)bis(phosphonate) ligand. This new in situ generated Ru complex, 4^2-, has been isolated and spectroscopically and electrochemically characterized. In addition, a crystal structure has been also obtained using single-crystal X-ray diffraction techniques. Complex 4^2- turns out to be an exceptional water oxidation catalyst achieving record maximum turnover frequencies (TOF_max) on the order of 16 000 s^-1. A mechanistic analysis complemented with DFT calculations has also been carried out, showing the critical role of intramolecular second coordination sphere effects exerted by the phosphonate groups in lowering the activation energy at the rate-determining step.

A Computational and Structural Database Study of the Metal-Carbene Bond in Groups IA, IIA, and IIIA Imidazol-2-Ylidene Complexes

Imidazol-2-ylidenes are important N-heterocyclic carbenes which have become universal ligands in organometallic and coordination chemistry. Generally classified as σ-donor ligands, these compounds have been used to stabilize various metal complexes which hitherto were less stable in their catalytic processes. Herein, the number and distribution of group IA, group IIA, and group IIIA metal-imidazol-2-ylidene complexes retrieved from the Cambridge Structural Database (CSD) are assessed. The data showed that the mean M-Ccarbene bond length increases with increasing ionic size but is similar across each diagonal. Dominant factors such as Lewis acidity and electrostatic attractions were found to control the bonding modes of the respective ions. Generally, the metal ions show preference for tetrahedral coordination with larger cations forming complexes with higher coordination numbers. For their high number of entries (101), tetrahedrally coordinated boron complexes with various electron withdrawing and electron donating groups were studied computationally at the DFT/B3LYP level of theory. The strength of the B-Ccarbene bond was found to depend on steric interactions between bulky groups on the borenium atom and substituents on the N-positions of the imidazol-2-ylidene ligand. This observation was further confirmed by estimation of the binding energy, natural charge, and the electron distribution in the B-Ccarbene bond.

Why Cellular Di/Triphosphates Preferably Bind Mg2+ and Not Ca2

Di/triphosphates perform a multitude of essential tasks, being important components of many vital organic cofactors such as adenosine/guanosine di/triphosphate (ADP/GDP, ATP/GTP), flavin adenine dinucleotide, and nicotinamide adenine dinucleotide and its phosphate derivative. They are generally bound to cations inside cells, in particular Mg²⁺ in the case of ATP/GTP. Yet how their metal-binding modes depend on the number, charge, and solvent exposure of the polyphosphate group and how Mg²⁺and Ca²⁺ dications that coexist in cellular fluids compete for di/triphosphates in biological systems remain elusive. Using density functional theory calculations combined with a polarizable continuum model, we have determined the relative free energies and stabilities of the different binding modes of di- and triphosphate groups to Mg²⁺ and Ca²⁺. We show that the thermodynamic outcome of the competition between Mg²⁺ and Ca²⁺ for cellular di/triphosphates depends mainly on the oligomericity/charge and metal-binding mode of the phosphate ligand as well as the solvent exposure of the binding site. Increasing the charge and thus denticity of the phosphate ligand from bi- to tridentate in a buried binding pocket enhances the affinity of the host system for the stronger charge acceptor, Mg²⁺. The cellular di/triphosphates's intrinsic properties and the protein matrix allowing them to bind a dication bi/tridentately, along with the higher cytosolic concentration of Mg²⁺ compared to Ca²⁺, enables Mg²⁺ to outcompete Ca²⁺ in binding to these highly charged anions. This suggests an explanation for why nature has chosen Mg²⁺ but not Ca²⁺ to perform most of the essential tasks associated with biological triphosphates.

Identification of a Novel Multifunctional Ligand for Simultaneous Inhibition of Amyloid-Beta (Aβ42) and Chelation of Zinc Metal Ion

Zinc binding to β-amyloid structure could promote amyloid-β aggregation, as well as reactive oxygen species (ROS) production, as suggested in many experimental and theoretical studies. Therefore, the introduction of multifunctional drugs capable of chelating zinc metal ion and inhibiting Aβ aggregation is a promising strategy in the development of AD treatment. The present study has evaluated the efficacy of a new bifunctional peptide drug using molecular docking and molecular dynamics (MD) simulations. This drug comprises two different domains, an inhibitor domain, obtained from the C-terminal hydrophobic region of Aβ, and a Zn²⁺ chelating domain, derived from rapeseed meal, merge with a linker. The multifunctionality of the ligand was evaluated using a comprehensive set of MD simulations spanning up to 3.2 μs including Aβ relaxation, ligand-Zn²⁺ bilateral interaction, and, more importantly, ligand-Zn²⁺-Aβ₄₂ trilateral interactions. Analysis of the results strongly indicated that the bifunctional ligand can chelate zinc metal ion and avoid Aβ aggregation simultaneously. The present study illustrated that the proposed ligand has considerable hydrophobic interactions and hydrogen bonding with monomeric Aβ in the presence of zinc metal ion. Therefore, in light of these considerable interactions and contacts, the α-helical structure of Aβ has been enhanced, while the β-sheet formation is prevented and the α-helix native structure is protected. Furthermore, the analysis of interactions between Aβ and ligand-zinc complex revealed that the zinc metal ion is coordinated to Met13, the ending residue of the ligand and merely one residue in Aβ. The results have proven the previous experimental and theoretical findings in the literature about Aβ interactions with zinc metal ion and also Aβ interactions with the first domain of the proposed ligand. Moreover, the current research has evaluated the chelation using MD simulation and linear interaction energy (LIE) methods, and the result has been satisfactorily verified with previous experimental and theoretical (DFT) studies.

Free and Bound Therapeutic Lithium in Brain Signaling

Lithium, a first-line therapy for bipolar disorder, is effective in preventing suicide and new depressive/manic episodes. Yet, how this beguilingly simple monocation with only two electrons could yield such profound therapeutic effects remains unclear. An in-depth understanding of lithium's mechanisms of actions would help one to develop better treatments limiting its adverse side effects and repurpose lithium for treating traumatic brain injury and chronic neurodegenerative diseases. In this Account, we begin with a comparison of the physicochemical properties of Li⁺ and its key native rivals, Na⁺ and Mg²⁺, to provide physical grounds for their competition in protein binding sites. Next, we review the abnormal signaling pathways and proteins found in bipolar patients, who generally have abnormally high intracellular Na⁺ and Ca²⁺ concentrations, high G-protein levels, and hyperactive phosphatidylinositol signaling and glycogen synthase kinase-3β (GSK3β) activity. We briefly summarize experimental findings on how lithium, at therapeutic doses, modulates these abnormal signaling pathways and proteins. Following this survey, we address the following aspects of lithium's therapeutic actions: (1) Can Li⁺ displace Na⁺ from the allosteric Na⁺-binding sites in neurotransmitter transporters and G-protein coupled receptors (GPCRs); if so, how would this affect the host protein's function? (2) Why are certain Mg²⁺-dependent enzymes targeted by Li⁺? (3) How does Li⁺ binding to Mg²⁺-bound ATP/GTP (denoted as NTP) in solution affect the cofactor's conformation and subsequent recognition by the host protein? (4) How do NTP-Mg-Li complexes modulate the properties of the respective cellular receptors and signal-transducing proteins? We show that Li⁺ may displace Na⁺ from allosteric Na⁺-binding sites in certain GPCRs and stabilize inactive conformations, preventing these receptors from relaying signal to the respective G-proteins. It may also displace Mg²⁺ in enzymes containing highly cationic Mg²⁺-binding sites such as GSK3β, but not in enzymes containing Mg²⁺-binding sites with low or zero charge. We further show that Li⁺ binding to Mg²⁺-NTP in water does not alter the NTP conformation, which is locked by all three phosphates binding to Mg²⁺. However, bound lithium in the form of [NTP-Mg-Li]^2- dianions can activate or inhibit the host protein depending on the NTP-binding pocket's shape, which determines the metal-binding mode: The ATP-binding pocket's shape in the P2X receptor is complementary to the native ATP-Mg solution conformation and nicely fits [ATP-Mg-Li]^2-. However, since the ATP βγ phosphates bind Li⁺, bimetallic [ATP-Mg-Li]^2- may be more resistant to hydrolysis than the native cofactor, enabling ATP to reside longer in the binding site and elicit a prolonged P2X response. In contrast, the elongated GTP-binding pockets in G-proteins allow only two GTP phosphates to bind Mg²⁺, so the GTP conformation is no longer "triply-locked". Consequently, Li⁺ binding to GTP-Mg can significantly alter the native cofactor's structure, lowering the activated G-protein level, thus attenuating hyperactive G-protein-mediated signaling in bipolar patients. In summary, we have presented a larger "connected" picture of lithium's diverse effects based on its competition as a free monocation with native cations or as a phosphate-bound polyanionic complex modulating the host protein function.

Harnessing Ion-Binding Sites for GPCR Pharmacology

Endogenous ions play important roles in the function and pharmacology of G-protein coupled receptors (GPCRs). Historically the evidence for ionic modulation of GPCR function dates to 1973 with studies of opioid receptors, where it was demonstrated that physiologic concentrations of sodium allosterically attenuated agonist binding. This Na+-selective effect was distinct from effects of other monovalent and divalent cations, with the latter usually counteracting sodium's negative allosteric modulation of binding. Since then, numerous studies documenting the effects of mono- and divalent ions on GPCR function have been published. While ions can act selectively and nonselectively at many sites in different receptors, the discovery of the conserved sodium ion site in class A GPCR structures in 2012 revealed the unique nature of Na+ site, which has emerged as a near-universal site for allosteric modulation of class A GPCR structure and function. In this review, we synthesize and highlight recent advances in the functional, biophysical, and structural characterization of ions bound to GPCRs. Taken together, these findings provide a molecular understanding of the unique roles of Na+ and other ions as GPCR allosteric modulators. We will also discuss how this knowledge can be applied to the redesign of receptors and ligand probes for desired functional and pharmacological profiles. SIGNIFICANCE STATEMENT: The function and pharmacology of GPCRs strongly depend on the presence of mono and divalent ions in experimental assays and in living organisms. Recent insights into the molecular mechanism of this ion-dependent allosterism from structural, biophysical, biochemical, and computational studies provide quantitative understandings of the pharmacological effects of drugs in vitro and in vivo and open new avenues for the rational design of chemical probes and drug candidates with improved properties.

Imidazol-2-ylidene stabilized tetrahedral cobalt carbonyl complexes: A computational and structural database study

A combined crystallographic database and computational chemistry analyses of tetrahedral cobalt (Co) carbonyl (CO) complexes bearing substituted imidazol-2-ylidene (Im) ligands were studied. These group of compounds are potential catalysts in the hydroformylation of alkenes as well as effective precursors for both chemical vapour deposition and atomic layer deposition in integrated circuit manufacture. A substructure search revealed seventeen (17) crystal structures in the Cambridge Structural Database (CSD). Generally, there was a positive correlation between the Im-Co bond length and the Co-CO bond length implying that a strong Im-CO bonding interaction could give rise to a more stable carbonyl complex. The gas phase structures of these complexes were further studied computationally at the BP86 level of the density functional theory. The def2-TZVPP basis set was used for cobalt with 6-31G (d,p) for all other atoms. It was observed that the Co-CO bond length increased with increase in bulkiness of alkyl substituents at the ortho positions of the Im ligands. Subsequent natural bond orbital (NBO) and chemical reactivity parameter analyses revealed that the chemical softness, electronic chemical potential and the ease of ionization also increased with increase in bulkiness of the alkyl substituents. This property has useful applications in catalytic systems and for the design of precursors in integrated circuit manufacture.

A Million Crystal Structures: The Whole Is Greater than the Sum of Its Parts

The founding in 1965 of what is now called the Cambridge Structural Database (CSD) has reaped dividends in numerous and diverse areas of chemical research. Each of the million or so crystal structures in the database was solved for its own particular reason, but collected together, the structures can be reused to address a multitude of new problems. In this Review, which is focused mainly on the last 10 years, we chronicle the contribution of the CSD to research into molecular geometries, molecular interactions, and molecular assemblies and demonstrate its value in the design of biologically active molecules and the solid forms in which they are delivered. Its potential in other commercially relevant areas is described, including gas storage and delivery, thin films, and (opto)electronics. The CSD also aids the solution of new crystal structures. Because no scientific instrument is without shortcomings, the limitations of CSD research are assessed. We emphasize the importance of maintaining database quality: notwithstanding the arrival of big data and machine learning, it remains perilous to ignore the principle of garbage in, garbage out. Finally, we explain why the CSD must evolve with the world around it to ensure it remains fit for purpose in the years ahead.

A Fragment Quantum Mechanical Method for Metalloproteins

An accurate energy calculation of metalloprotein is of crucial importance and also a theoretical challenge. In this work, a metal molecular fractionation with conjugate caps (metal-MFCC) approach is developed for efficient linear-scaling quantum calculation of potential energy and atomic forces of metalloprotein. In this approach, the potential energy of a given protein is calculated by a linear combination of potential energies of the neighboring residues, two-body interaction energy between non-neighboring residues that are spatially in close contact and the potential energy of the metal binding group. The calculation of each fragment is embedded in a field of point charges representing the remaining protein environment. Numerical studies were carried out to check the performance of this method, and the calculated potential energies and atomic forces all show excellent agreement with the full system calculations at the M06-2X/6-31G(d) level. By combining the energy calculation with molecular dynamic simulation, we performed an ab initio structural optimization for a zinc finger protein with high efficiency. The present metal-MFCC approach is linear-scaling with a low prefactor and trivially parallelizable. The individual fragment typically contains about 50 atoms, and it is thus possible to be calculated at higher levels of the quantum chemistry method. This fragment method can be routinely applied to perform structural optimization and ab initio molecular dynamic simulation for metalloproteins of any size.

Activation Stoichiometry and Pore Architecture of TRPA1 Probed with Channel Concatemers

The nociceptor ion channel TRPA1 detects a wide range of hazardous chemicals, including reactive electrophiles such as cinnamaldehyde, which gate the channel allowing Na⁺ and Ca²⁺ entry. TRPA1 assembles as a tetramer, with a central pore within which an aspartate residue (D918) determines Ca²⁺ permeability. Here, we report that introduction of histidine at this position, D918H, makes TRPA1 channels sensitive to block by nanomolar concentration of Zn²⁺ and can be used to functionally tag subunits in concatemers. Concatemers with increasing numbers of D918H subunits display increasing sensitivity to Zn²⁺ inhibition, indicating that the four side chains at position 918 of the tetramer directly coordinate Zn²⁺ and other permeating divalent cations. In the published structure of TRPA1, this requires a rearrangement of the pore region which may represent the true open state of the channel. Concatemeric channels containing subunits mutated to be insensitive to reactive electrophiles (C622S) could be activated by cinnamaldehyde when as few as two subunits contained intact ligand binding sites. Activation upon liganding of just two of the four possible subunits may represent an optimal strategy to rapidly and reliably detect noxious chemicals.

Data on PAGE analysis and MD simulation for the interaction of endonuclease Apn1 from Saccharomyces cerevisiae with DNA substrates containing 5,6-dihydrouracyl and 2-aminopurine

This article presents new data on nucleotide incision repair (NIR) activity of apurinic/apyrimidinic endonuclease Apn1 of Saccharomyces cerevisiae, which is known as a key player of the base excision DNA repair (BER) pathway, see “Yeast structural gene (APN1) for the major apurinic endonuclease: homology to Escherichia coli endonuclease IV” [1], “Abasic sites in DNA: repair and biological consequences in Saccharomyces cerevisiae” [2] and “Characterisation of new substrate specificities of Escherichia coli and Saccharomyces cerevisiae AP endonucleases” [3]. The characterization of NIR activity of wild type Apn1 and mutant form Ape1 H83A were made by denaturing PAGE analysis, and MD simulations of Apn1 complexed with DNA containing 5,6-dihydro-2′-deoxyuridine (DHU) and 2-aminopurine (2-aPu) residues. This data article is associated to the manuscript titled “Apurinic/apyrimidinic endonuclease Apn1 from Saccharomyces cerevisiae is recruited to the nucleotide incision repair pathway: kinetic and structural features” [4].

Published by Elsevier Inc. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).

Divalent metal ions control activity and inhibition of protein kinases

Protein kinases are key enzymes in the regulation of eukaryotic signal transduction. As metalloenzymes they employ divalent cations for catalysis and regulation. We used the catalytic (C) subunit of cAMP-dependent protein kinase (PKA) as a model protein to investigate the role of a variety of physiologically or pathophysiologically relevant divalent metal ions in distinct steps within the catalytic cycle. It is established that divalent metal ions play a crucial role in co-binding of nucleotides and also assist in catalysis. Our studies reveal that besides the physiologically highly relevant magnesium, metals like zinc and manganese can assist in phosphoryl transfer, however, only a few support efficient substrate turnover (turnover catalysis). Those trace metals allow for substrate binding and phosphotransfer but hamper product release. We further established the unique role of magnesium as the common biologically relevant divalent metal ideally enabling (co-) substrate binding and orientation. Magnesium allows stable substrate binding and, on the other hand accelerates product release, thus being extremely efficient in turnover catalysis. We extended our studies to non-catalytic functions of protein kinases looking at pseudokinases, a subfamily of protein kinases inherently lacking critical residues for catalysis. Recently, pseudokinases have been linked to human diseases. Some pseudokinases are still capable of binding metal ions, yet have lost the ability to transfer the phosphoryl group from ATP to a given substrate. Here metal ions stabilize an active like, though catalytically unproductive conformation and are therefore crucial to maintain non-catalytic function. Finally, we demonstrate for the canonical kinase PKA that the trace metal manganese alone can stabilize protein kinases in an active like conformation allowing them to bind substrates even in the absence of nucleotides.