Metal-binding sites are designed to achieve optimal mechanical and signaling properties

Formation of Supplementary Metal-Binding Centers in Proteins under Stress Conditions

In many proteins, supplementary metal-binding centers appear under stress conditions. They are known as aberrant or atypical sites. Physico-chemical properties of proteins are significantly changed after such metal binding, and very stable protein aggregates are formed, in which metals act as "cross-linking" agents. Supplementary metal-binding centers in proteins often arise as a result of posttranslational modifications caused by reactive oxygen and nitrogen species and reactive carbonyl compounds. New chemical groups formed as a result of these modifications can act as ligands for binding metal ions. Special attention is paid to the role of cysteine SH-groups in the formation of supplementary metal-binding centers, since these groups are the main target for the action of reactive species. Supplementary metal binding centers may also appear due to unmasking of amino acid residues when protein conformation changing. Appearance of such centers is usually considered as a pathological process. Such unilateral approach does not allow to obtain an integral view of the phenomenon, ignoring cases when formation of metal complexes with altered proteins is a way to adjust protein properties, activity, and stability under the changed redox conditions. The role of metals in protein aggregation is being studied actively, since it leads to formation of non-membranous organelles, liquid condensates, and solid conglomerates. Some proteins found in such aggregates are typical for various diseases, such as Alzheimer's and Huntington's diseases, amyotrophic lateral sclerosis, and some types of cancer.

Regulation of zinc-dependent enzymes by metal carrier proteins

Zn²⁺ ions are essential in many physiological processes, including enzyme catalysis, protein structural stabilization, and the regulation of many proteins. The affinities of proteins for Zn²⁺ ions span several orders of magnitude, with catalytic Zn²⁺ ions generally held more tightly than structural or regulatory ones. Metal carrier proteins, most of which are not specific for Zn²⁺, bind these ions with a broad range of affinities that overlap those of catalytic, structural, and regulatory Zn²⁺ ions and are thought to be responsible for distributing the metal through most cells, tissues, and fluid compartments. While little is known about how many proteins obtain or release these ions, there is now considerable experimental evidence suggesting that metal carrier proteins may be responsible for transferring metals to and from some Zn²⁺-dependent proteins, thus serving as a major regulatory factor for them. In this review, the biological roles of Zn²⁺ and structures of Zn²⁺ binding sites are examined, and experimental evidence demonstrating the direct participation of metal carrier proteins in enzyme regulation is discussed. Mechanisms of metal ion transfer are also offered, and the potential physiological significance of this phenomenon is explored.

Epoxyqueuosine Reductase QueH in the Biosynthetic Pathway to tRNA Queuosine Is a Unique Metalloenzyme

Queuosine is a structurally unique and functionally important tRNA modification, widely distributed in eukaryotes and bacteria. The final step of queuosine biosynthesis is the reduction/deoxygenation of epoxyqueuosine to form the cyclopentene motif of the nucleobase. The chemistry is performed by the structurally and functionally characterized cobalamin-dependent QueG. However, the queG gene is absent from several bacteria that otherwise retain queuosine biosynthesis machinery. Members of the IPR003828 family (previously known as DUF208) have been recently identified as nonorthologous replacements of QueG, and this family was renamed QueH. Here, we present the structural characterization of QueH from Thermotoga maritima. The structure reveals an unusual active site architecture with a [4Fe-4S] metallocluster along with an adjacent coordinated iron metal. The juxtaposition of the cofactor and coordinated metal ion predicts a unique mechanism for a two-electron reduction/deoxygenation of epoxyqueuosine. To support the structural characterization, in vitro biochemical and genomic analyses are presented. Overall, this work reveals new diversity in the chemistry of iron/sulfur-dependent enzymes and novel insight into the last step of this widely conserved tRNA modification.

Cardiovascular Risk Stratification Based on Oxidative Stress for Early Detection of Pathology

Aims: Current cardiovascular (CV) risk prediction algorithms are able to quantify the individual risk of CV disease. However, CV risk in young adults is underestimated due to the high dependency of age in biomarker-based algorithms. Because oxidative stress is associated with CV disease, we sought to examine CV risk stratification in young adults based on oxidative stress to approach the discovery of new markers for early detection of pathology. Results: Young adults were stratified into (i) healthy controls, (ii) subjects with CV risk factors, and (iii) patients with a reported CV event. Plasma samples were analyzed using FASILOX, a novel approach to interrogate the dynamic thiol redox proteome. We also analyzed irreversible oxidation by targeted searches using the Uniprot database. Irreversible oxidation of cysteine (Cys) residues was greater in patients with reported CV events than in healthy subjects. These results also indicate that oxidation is progressive. Moreover, we found that glutathione reductase and glutaredoxin 1 proteins are differentially expressed between groups and are proteins involved in antioxidant response, which is in line with the impaired redox homeostasis in CV disease. Innovation: This study, for the first time, describes the oxidative stress (reversible and irreversible Cys oxidation) implication in human plasma according to CV risk stratification. Conclusion: The identification of redox targets and the quantification of protein and oxidative changes might help to better understand the role of oxidative stress in CV disease, and aid stratification for CV events beyond traditional prognostic and diagnostic markers. Antioxid. Redox Signal. 35, 602-617.

Identification of Allosteric Effects in Proteins by Elastic Network Models

Allostery is a fundamental regulatory mechanism in the majority of biological processes of molecular machines. Allostery is well-known as a dynamic-driven process, and thus, the molecular mechanism of allosteric signal transmission needs to be established. Elastic network models (ENMs) provide efficient methods for investigating the intrinsic dynamics and allosteric communication pathways in proteins. In this chapter, two ENM methods including Gaussian network model (GNM) coupled with Markovian stochastic model, as well as the anisotropic network model (ANM), were introduced to identify allosteric effects in hemoglobins. Techniques on model parameters, scripting and calculation, analysis, and visualization are shown step by step.

Transcriptome Analysis of Genes Involved in Cold Hardiness of Peach Tree (Prunus persica) Shoots during Cold Acclimation and Deacclimation

We analyzed the transcriptomes in the shoots of five-year-old 'Soomee' peach trees (Prunus persica) during cold acclimation (CA), from early CA (end of October) to late CA (middle of January), and deacclimation (DA), from late CA to late DA (middle of March), to identify the genes involved in cold hardiness. Cold hardiness of the shoots increased from early to late CA, but decreased from late CA to late DA, as indicated by decreased and increased the median lethal temperature (LT₅₀), respectively. Transcriptome analysis identified 17,208 assembled transcripts during all three stages. In total, 1891 and 3008 transcripts were differentially expressed with a |fold change| > 2 (p < 0.05) between early and late CA, and between late CA and late DA, respectively. Among them, 1522 and 2830, respectively, were functionally annotated with gene ontology (GO) terms having a greater proportion of differentially expressed genes (DEGs) associated with molecular function than biological process or cellular component categories. The biochemical pathways best represented both periods from early to late CA and from late CA to late DA were 'metabolic pathway' and 'biosynthesis of secondary metabolites'. We validated these transcriptomic results by performing reverse transcription quantitative polymerase chain reaction on the selected DEGs showing significant fold changes. The relative expressions of the selected DEGs were closely related to the LT₅₀ values of the peach tree shoots: 'Soomee' shoots exhibited higher relative expressions of the selected DEGs than shoots of the less cold-hardy 'Odoroki' peach trees. Irrespective of the cultivar, the relative expressions of the DEGs that were up- and then down-regulated during CA, from early to late CA, and DA, from late CA to late DA, were more closely correlated with cold hardiness than those of the DEGs that were down- and then up-regulated. Therefore, our results suggest that the significantly up- and then down-regulated DEGs are associated with cold hardiness in peach tree shoots. These DEGs, including early light-induced protein 1, chloroplastic, 14-kDa proline-rich protein DC2.15, glutamate dehydrogenase 2, and triacylglycerol lipase 2, could be candidate genes to determine cold hardiness.

Structural Model of the ETR1 Ethylene Receptor Transmembrane Sensor Domain

The structure, mechanism of action and copper stoichiometry of the transmembrane sensor domain of the plant ethylene receptor ETR1 and homologs have remained elusive, hampering the understanding on how the perception of the plant hormone ethylene is transformed into a downstream signal. We generated the first structural model of the transmembrane sensor domain of ETR1 by integrating ab initio structure prediction and coevolutionary information. To refine and independently validate the model, we determined protein-related copper stoichiometries on purified receptor preparations and explored the helix arrangement by tryptophan scanning mutagenesis. All-atom molecular dynamics simulations of the dimeric model reveal how ethylene can bind proximal to the copper ions in the receptor, illustrating the initial stages of the ethylene perception process.

Interpreting the Dynamics of Binding Interactions of snRNA and U1A Using a Coarse-Grained Model

The binding interactions of small nuclear RNAs (snRNA) and the associated protein factors are critical to the function of spliceosomes in alternatively splicing primary RNA transcripts. Although molecular dynamics simulations are a powerful tool to interpret the mechanism of biological processes, the atomic-level simulations are, however, too expensive and with limited accuracy for the large-size systems, such as snRNA-protein complexes. We extend the coarse-grained Gaussian network model, which models the RNA-protein complexes as a harmonic chain of Cα, P, and O4' atoms, to investigating the impact of the snRNA-binding interaction on the dynamic stability of the human U1A protein, which is a major component of the spliceosomal U1 small nuclear ribonucleoprotein particle. The results reveal that the first and third loops and the C-terminal helix regions of the U1A domain undergo a significant loss of flexibility upon the RNA binding due to the forming of mostly electrostatic and hydrogen bond interactions with RNA 5' stem and loop. By examining the residues whose mutations significantly change the binding free energy between U1A and snRNA, the Gaussian network model-based calculations show that not only the residues at the binding sites that are traditionally considered to play a major role in U1A-RNA association but also those residues that are far away from the RNA-binding interface can participate in the long-range allosteric signal transmission; these calculations are quantitatively consistent with the data observed in the recent snRNA binding experiments. The study demonstrates a useful avenue to utilize the simplified elastic network model to investigate the dynamics characteristics of the biologically important macromolecular interactions.

Comparative Study of Elastic Network Model and Protein Contact Network for Protein Complexes: The Hemoglobin Case

The overall topology and interfacial interactions play key roles in understanding structural and functional principles of protein complexes. Elastic Network Model (ENM) and Protein Contact Network (PCN) are two widely used methods for high throughput investigation of structures and interactions within protein complexes. In this work, the comparative analysis of ENM and PCN relative to hemoglobin (Hb) was taken as case study. We examine four types of structural and dynamical paradigms, namely, conformational change between different states of Hbs, modular analysis, allosteric mechanisms studies, and interface characterization of an Hb. The comparative study shows that ENM has an advantage in studying dynamical properties and protein-protein interfaces, while PCN is better for describing protein structures quantitatively both from local and from global levels. We suggest that the integration of ENM and PCN would give a potential but powerful tool in structural systems biology.

Comparison of the Internal Dynamics of Metalloproteases Provides New Insights on Their Function and Evolution

Metalloproteases have evolved in a vast number of biological systems, being one of the most diverse types of proteases and presenting a wide range of folds and catalytic metal ions. Given the increasing understanding of protein internal dynamics and its role in enzyme function, we are interested in assessing how the structural heterogeneity of metalloproteases translates into their dynamics. Therefore, the dynamical profile of the clan MA type protein thermolysin, derived from an Elastic Network Model of protein structure, was evaluated against those obtained from a set of experimental structures and molecular dynamics simulation trajectories. A close correspondence was obtained between modes derived from the coarse-grained model and the subspace of functionally-relevant motions observed experimentally, the later being shown to be encoded in the internal dynamics of the protein. This prompted the use of dynamics-based comparison methods that employ such coarse-grained models in a representative set of clan members, allowing for its quantitative description in terms of structural and dynamical variability. Although members show structural similarity, they nonetheless present distinct dynamical profiles, with no apparent correlation between structural and dynamical relatedness. However, previously unnoticed dynamical similarity was found between the relevant members Carboxypeptidase Pfu, Leishmanolysin, and Botulinum Neurotoxin Type A, despite sharing no structural similarity. Inspection of the respective alignments shows that dynamical similarity has a functional basis, namely the need for maintaining proper intermolecular interactions with the respective substrates. These results suggest that distinct selective pressure mechanisms act on metalloproteases at structural and dynamical levels through the course of their evolution. This work shows how new insights on metalloprotease function and evolution can be assessed with comparison schemes that incorporate information on protein dynamics. The integration of these newly developed tools, if applied to other protein families, can lead to more accurate and descriptive protein classification systems.

Study, by use of coarse-grained models, of the functionally crucial residues and allosteric pathway of anesthetic regulation of the Gloeobacter violaceus ligand-gated ion channel

Although pentameric ligand-gated ion channels (pLGICs) have been found to be the targets of general anesthetics, the mechanism of the effects of anesthetics on pLGICs remains elusive. pLGICs from Gloeobacter violaceus (GLIC) can be inhibited by the anesthetic ketamine. X-ray crystallography has shown that the ketamine binding site is distant from the channel gate of the GLIC. It is still not clear how ketamine controls the function of the GLIC by long-range allosteric regulation. In this work, the functionally crucial residues and allosteric pathway of anesthetic regulation of the GLIC were identified by use of a coarse-grained thermodynamic method developed by our group. In our method, the functionally crucial sites were identified as the residues thermodynamically coupled with binding of ketamine. The results from calculation were highly consistent with experimental data. Our study aids understanding of the mechanism of the anesthetic action of ketamine on the GLIC by long-range allosteric modulation.

Prediction of allosteric sites on protein surfaces with an elastic-network-model-based thermodynamic method

Allostery is a rapid and efficient way in many biological processes to regulate protein functions, where binding of an effector at the allosteric site alters the activity and function at a distant active site. Allosteric regulation of protein biological functions provides a promising strategy for novel drug design. However, how to effectively identify the allosteric sites remains one of the major challenges for allosteric drug design. In the present work, a thermodynamic method based on the elastic network model was proposed to predict the allosteric sites on the protein surface. In our method, the thermodynamic coupling between the allosteric and active sites was considered, and then the allosteric sites were identified as those where the binding of an effector molecule induces a large change in the binding free energy of the protein with its ligand. Using the proposed method, two proteins, i.e., the 70 kD heat shock protein (Hsp70) and GluA2 alpha-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid (AMPA) receptor, were studied and the allosteric sites on the protein surface were successfully identified. The predicted results are consistent with the available experimental data, which indicates that our method is a simple yet effective approach for the identification of allosteric sites on proteins.

Identification of functionally key residues in AMPA receptor with a thermodynamic method

AMPA receptor mediates the fast excitatory synaptic transmission in the central nervous system, and it is activated by the binding of glutamate that results in the opening of the transmembrane ion channel. In the present work, the thermodynamic method developed by our group was improved and then applied to identify the functionally key residues that regulate the glutamate-binding affinity of AMPA receptor. In our method, the key residues are identified as those whose perturbation largely changes the ligand binding free energy of the protein. It is found that besides the ligand binding sites, other residues distant from the binding cleft can also influence the glutamate binding affinity through a long-range allosteric regulation. These allosteric sites include the hinge region of the ligand binding cleft, the dimer interface of the ligand binding domain, the linkers between the ligand binding domain and the transmembrane domain, and the interface between the N-terminal domain and the ligand binding domain. Our calculation results are consistent with the available experimental data. The results are helpful for our understanding of the mechanism of long-range allosteric communication in the AMPA receptor and the mechanism of channel opening triggered by glutamate binding.

Molecular basis for manganese sequestration by calprotectin and roles in the innate immune response to invading bacterial pathogens

The S100A8/S100A9 heterodimer calprotectin (CP) functions in the host response to pathogens through a mechanism termed "nutritional immunity." CP binds Mn(2+) and Zn(2+) with high affinity and starves bacteria of these essential nutrients. Combining biophysical, structural, and microbiological analysis, we identified the molecular basis of Mn(2+) sequestration. The asymmetry of the CP heterodimer creates a single Mn(2+)-binding site from six histidine residues, which distinguishes CP from all other Mn(2+)-binding proteins. Analysis of CP mutants with altered metal-binding properties revealed that, despite both Mn(2+) and Zn(2+) being essential metals, maximal growth inhibition of multiple bacterial pathogens requires Mn(2+) sequestration. These data establish the importance of Mn(2+) sequestration in defense against infection, explain the broad-spectrum antimicrobial activity of CP relative to other S100 proteins, and clarify the impact of metal depletion on the innate immune response to infection.

Sequence evolution correlates with structural dynamics

Biochemical activity and core stability are essential properties of proteins, maintained usually by conserved amino acids. Structural dynamics emerged in recent years as another essential aspect of protein functionality. Structural dynamics enable the adaptation of the protein to binding substrates and to undergo allosteric transitions, while maintaining the native fold. Key residues that mediate structural dynamics would thus be expected to be conserved or exhibit coevolutionary patterns at least. Yet, the correlation between sequence evolution and structural dynamics is yet to be established. With recent advances in efficient characterization of structural dynamics, we are now in a position to perform a systematic analysis. In the present study, a set of 34 enzymes representing various folds and functional classes is analyzed using information theory and elastic network models. Our analysis shows that the structural regions distinguished by their coevolution propensity as well as high mobility are predisposed to serve as substrate recognition sites, whereas residues acting as global hinges during collective dynamics are often supported by conserved residues. We propose a mobility scale for different types of amino acids, which tends to vary inversely with amino acid conservation. Our findings suggest the balance between physical adaptability (enabled by structure-encoded motions) and chemical specificity (conferred by correlated amino acid substitutions) underlies the selection of a relatively small set of versatile folds by proteins.

New structural and functional contexts of the Dx[DN]xDG linear motif: insights into evolution of calcium-binding proteins

Binding of calcium ions (Ca²⁺) to proteins can have profound effects on their structure and function. Common roles of calcium binding include structure stabilization and regulation of activity. It is known that diverse families – EF-hands being one of at least twelve – use a Dx[DN]xDG linear motif to bind calcium in near-identical fashion. Here, four novel structural contexts for the motif are described. Existing experimental data for one of them, a thermophilic archaeal subtilisin, demonstrate for the first time a role for Dx[DN]xDG-bound calcium in protein folding. An integrin-like embedding of the motif in the blade of a β-propeller fold – here named the calcium blade – is discovered in structures of bacterial and fungal proteins. Furthermore, sensitive database searches suggest a common origin for the calcium blade in β-propeller structures of different sizes and a pan-kingdom distribution of these proteins. Factors favouring the multiple convergent evolution of the motif appear to include its general Asp-richness, the regular spacing of the Asp residues and the fact that change of Asp into Gly and vice versa can occur though a single nucleotide change. Among the known structural contexts for the Dx[DN]xDG motif, only the calcium blade and the EF-hand are currently found intracellularly in large numbers, perhaps because the higher extracellular concentration of Ca²⁺ allows for easier fixing of newly evolved motifs that have acquired useful functions. The analysis presented here will inform ongoing efforts toward prediction of similar calcium-binding motifs from sequence information alone.

Cooperative dynamics of proteins unraveled by network models

Recent years have seen a significant increase in the number of computational studies that adopted network models for investigating biomolecular systems dynamics and interactions. In particular, elastic network models have proven useful in elucidating the dynamics and allosteric signaling mechanisms of proteins and their complexes. Here we present an overview of two most widely used elastic network models, the Gaussian Network Model (GNM) and Anisotropic Network Model (ANM). We illustrate their use in (i) explaining the anisotropic response of proteins observed in external pulling experiments, (ii) identifying residues that possess high allosteric potentials, and demonstrating in this context the propensity of catalytic sites and metal-binding sites for enabling efficient signal transduction, and (iii) assisting in structure refinement, molecular replacement and comparative modeling of ligand-bound forms via efficient sampling of energetically favored conformers.

Structural basis for metal sensing by CnrX

CnrX is the metal sensor and signal modulator of the three-protein transmembrane signal transduction complex CnrYXH of Cupriavidus metallidurans CH34 that is involved in the setup of cobalt and nickel resistance. We have determined the atomic structure of the soluble domain of CnrX in its Ni-bound, Co-bound, or Zn-bound form. Ni and Co ions elicit a biological response, while the Zn-bound form is inactive. The structures presented here reveal the topology of intraprotomer and interprotomer interactions and the ability of metal-binding sites to fine-tune the packing of CnrX dimer as a function of the bound metal. These data suggest an allosteric mechanism to explain how the complex is switched on and how the signal is modulated by Ni or Co binding. These results provide clues to propose a model for signal propagation through the membrane in the complex.