Insights into Cu(I) exchange in HAH1 using quantum mechanical and molecular simulations

Unraveling Copper Exchange in the Atox1-Cu(I)-Mnk1 Heterodimer: A Simulation Approach

Copper, an essential metal for various cellular processes, requires tight regulation to prevent cytotoxicity. Intracellular pathways crucial for maintaining optimal copper levels involve soluble and membrane transporters, namely, metallochaperones and P-type ATPases, respectively. In this study, we used a simulation workflow based on free-energy perturbation (FEP) theory and parallel bias metadynamics (PBMetaD) to predict the Cu(I) exchange mechanism between the human Cu(I) chaperone, Atox1, and one of its two physiological partners, ATP7A. ATP7A, also known as the Menkes disease protein, is a transmembrane protein and one of the main copper-transporting ATPases. It pumps copper into the trans-Golgi network for the maturation of cuproenzymes and is also essential for the efflux of excess copper across the plasma membrane. In this analysis, we utilized the nuclear magnetic resonance (NMR) structure of the Cu(I)-mediated complex between Atox1 and the first soluble domain of the Menkes protein (Mnk1) as a starting point. Independent free-energy simulations were conducted to investigate the dissociation of both Atox1 and Mnk1. The calculations revealed that the two dissociations require free energy values of 6.3 and 6.2 kcal/mol, respectively, following a stepwise dissociation mechanism.

Role of Explicit Solvation in Computational Modeling of Chemical Reactions: Mechanism of Cu(I) Transfer Between Thiolate-Based Chelators in Water

Solvation plays important roles in controlling the thermodynamic and kinetic aspects of chemical reactions. The conventional approaches to treat solvation via electronic structure methods are likely to become inadequate, when the reacting solutes have strong electrostatic and hydrogen bonding interactions with the solvent and undergo significant structural changes during the course of the reaction. In this article, we present evidence of such solvent and structural effects in the computational study of the Cu(I) transfer reaction between thiolate-based chelators dithiobutylamine (DTBA) and dithiotheritol (DTT) in water, inspired from biological copper trafficking phenomena. We propose a general solution to the problem by combining classical molecular dynamics (MD) simulations of the bulk system and static quantum chemistry calculations. The fluctuating solvation shell was estimated from MD, and energetics was assessed by averaging QM energies of a series of molecular clusters constructed from the MD snapshots. Applying this approach, we propose a reaction pathway with estimates of relative intermediate stabilities and barriers, which suggest the overall reaction to be reversible in nature and likely to go through both two and three coordinated intermediates, confirming previous studies of similar protein analogues. An interesting fact that emerged from our study was the strong indication that the rate-determining step is the deprotonation of initial thiol bound Cu(I) complex, without involving any Cu(I)-S bonds. The proposed method will lead to a better treatment of solvations, and these mechanistic insights will aid our understanding of biological copper(I) trafficking.

Disrupting Cu trafficking as a potential therapy for cancer

Copper ions play a crucial role in various cellular biological processes. However, these copper ions can also lead to toxicity when their concentration is not controlled by a sophisticated copper-trafficking system. Copper dys-homeostasis has been linked to a variety of diseases, including neurodegeneration and cancer. Therefore, manipulating Cu-trafficking to trigger selective cancer cell death may be a viable strategy with therapeutic benefit. By exploiting combined in silico and experimental strategies, we identified small peptides able to bind Atox1 and metal-binding domains 3-4 of ATP7B proteins. We found that these peptides reduced the proliferation of cancer cells owing to increased cellular copper ions concentration. These outcomes support the idea of harming copper trafficking as an opportunity for devising novel anti-cancer therapies.

Copper binding leads to increased dynamics in the regulatory N-terminal domain of full-length human copper transporter ATP7B

ATP7B is a human copper-transporting P_1B-type ATPase that is involved in copper homeostasis and resistance to platinum drugs in cancer cells. ATP7B consists of a copper-transporting core and a regulatory N-terminal tail that contains six metal-binding domains (MBD1-6) connected by linker regions. The MBDs can bind copper, which changes the dynamics of the regulatory domain and activates the protein, but the underlying mechanism remains unknown. To identify possible copper-specific structural dynamics involved in transport regulation, we constructed a model of ATP7B spanning the N-terminal tail and core catalytic domains and performed molecular dynamics (MD) simulations with (holo) and without (apo) copper ions bound to the MBDs. In the holo protein, MBD2, MBD3 and MBD5 showed enhanced mobilities, which resulted in a more extended N-terminal regulatory region. The observed separation of MBD2 and MBD3 from the core protein supports a mechanism where copper binding activates the ATP7B protein by reducing interactions among MBD1-3 and between MBD1-3 and the core protein. We also observed an increased interaction between MBD5 and the core protein that brought the copper-binding site of MBD5 closer to the high-affinity internal copper-binding site in the core protein. The simulation results assign specific, mechanistic roles to the metal-binding domains involved in ATP7B regulation that are testable in experimental settings.

Living organisms depend upon active transport against gradients across biological membranes for survival. Such transport can be accomplished by ATP-dependent membrane protein transporters for which the activity must be regulated to maintain optimal concentrations in the cellular compartments. The regulatory mechanisms often involve structural responses inherent to the protein structure, which because of their dynamic nature can be hard to assess experimentally. A prime example is regulation of cellular copper levels by a copper-binding tail in the human copper transporter ATP7B. Dysregulation can cause severe diseases, for example the copper metabolism disorder Wilson’s disease, which is caused by mutations in ATP7B regulation machinery. Due to the practical difficulties in working with membrane proteins, most studies of ATP7B have been conducted in the absence of the membrane-bound protein core. Here, we used computer simulations of full-length ATP7B to study how structural dynamics in the regulatory tail differ between copper-bound and copper-free states. Copper induced increased dynamics in the tail, resulting in an overall movement towards the ion-binding site in the protein core. The simulations identified several, hitherto not reported, interactions between the regulatory tail and the protein core that can be targeted experimentally to enhance our understanding of this medically relevant regulatory mechanism.

The conformational plasticity of the selectivity filter methionines controls the in-cell Cu(I) uptake through the CTR1 transporter

Copper is a trace element vital to many cellular functions. Yet its abnormal levels are toxic to cells, provoking a variety of severe diseases. The high affinity copper transporter 1 (CTR1), being the main in-cell copper [Cu(I)] entry route, tightly regulates its cellular uptake via a still elusive mechanism. Here, all-atoms simulations unlock the molecular terms of Cu(I) transport in eukaryotes disclosing that the two methionine (Met) triads, forming the selectivity filter, play an unprecedented dual role both enabling selective Cu(I) transport and regulating its uptake rate thanks to an intimate coupling between the conformational plasticity of their bulky side chains and the number of bound Cu(I) ions. Namely, the Met residues act as a gate reducing the Cu(I) import rate when two ions simultaneously bind to CTR1. This may represent an elegant autoregulatory mechanism through which CTR1 protects the cells from excessively high, and hence toxic, in-cell Cu(I) levels. Overall, our outcomes resolve fundamental questions in CTR1 biology and open new windows of opportunity to tackle diseases associated with an imbalanced copper uptake.

The Advantages of EPR Spectroscopy in Exploring Diamagnetic Metal Ion Binding and Transfer Mechanisms in Biological Systems

Electron paramagnetic resonance (EPR) spectroscopy has emerged as an ideal biophysical tool to study complex biological processes. EPR spectroscopy can follow minor conformational changes in various proteins as a function of ligand or protein binding or interactions with high resolution and sensitivity. Resolving cellular mechanisms, involving small ligand binding or metal ion transfer, is not trivial and cannot be studied using conventional biophysical tools. In recent years, our group has been using EPR spectroscopy to study the mechanism underlying copper ion transfer in eukaryotic and prokaryotic systems. This mini-review focuses on our achievements following copper metal coordination in the diamagnetic oxidation state, Cu(I), between biomolecules. We discuss the conformational changes induced in proteins upon Cu(I) binding, as well as the conformational changes induced in two proteins involved in Cu(I) transfer. We also consider how EPR spectroscopy, together with other biophysical and computational tools, can identify the Cu(I)-binding sites. This work describes the advantages of EPR spectroscopy for studying biological processes that involve small ligand binding and transfer between intracellular proteins.

Computing Metal-Binding Proteins for Therapeutic Benefit

Over one third of biomolecules rely on metal ions to exert their cellular functions. Metal ions can play a structural role by stabilizing the structure of biomolecules, a functional role by promoting a wide variety of biochemical reactions, and a regulatory role by acting as messengers upon binding to proteins regulating cellular metal-homeostasis. These diverse roles in biology ascribe critical implications to metal-binding proteins in the onset of many diseases. Hence, it is of utmost importance to exhaustively unlock the different mechanistic facets of metal-binding proteins and to harness this knowledge to rationally devise novel therapeutic strategies to prevent or cure pathological states associated with metal-dependent cellular dysfunctions. In this compendium, we illustrate how the use of a computational arsenal based on docking, classical, and quantum-classical molecular dynamics simulations can contribute to extricate the minutiae of the catalytic, transport, and inhibition mechanisms of metal-binding proteins at the atomic level. This knowledge represents a fertile ground and an essential prerequisite for selectively targeting metal-binding proteins with small-molecule inhibitors aiming to (i) abrogate deregulated metal-dependent (mis)functions or (ii) leverage metal-dyshomeostasis to selectively trigger harmful cells death.

Cu(I) Controls Conformational States in Human Atox1 Metallochaperone: An EPR and Multiscale Simulation Study

Atox1 is a human copper metallochaperone that is responsible for transferring copper ions from the main human copper transporter, hCtr1, to ATP7A/B in the Golgi apparatus. Atox1 interacts with the Ctr1 C-terminal domain as a dimer, although it transfers the copper ions to ATP7A/B in a monomeric form. The copper binding site in the Atox1 dimer involves Cys12 and Cys15, while Lys60 was also suggested to play a role in the copper binding. We recently showed that Atox1 can adopt various conformational states, depending on the interacting protein. In the current study, we apply EPR experiments together with hybrid quantum mechanics–molecular mechanics molecular dynamics simulations using a recently developed semiempirical density functional theory approach, to better understand the effect of Atox1’s conformational states on copper coordination. We propose that the flexibility of Atox1 occurs owing to protonation of one or more of the cysteine residues, and that Cys15 is an important residue for Atox1 dimerization, while Cys12 is a critical residue for Cu(I) binding. We also show that Lys60 electrostatically stabilizes the Cu(I)–Atox1 dimer.

Molecular Dynamics Simulations Based on Newly Developed Force Field Parameters for Cu2+ Spin Labels Provide Insights into Double-Histidine-Based Double Electron-Electron Resonance

Electron paramagnetic resonance (EPR) in combination with the recently developed double-histidine (dHis)-based Cu²⁺ spin labeling has provided valuable insights into protein structure and conformational dynamics. To relate sparse distance constraints measured by EPR to protein fluctuations in solution, modeling techniques are needed. In this work, we have developed force field parameters for Cu²⁺-nitrilotriacetic and Cu²⁺-iminodiacetic acid spin labels. We employed molecular dynamics (MD) simulations to capture the atomic-level details of dHis-labeled protein fluctuations. The interspin distances extracted from 200 ns MD trajectories show good agreement with the experimental results. The MD simulations also illustrate the dramatic rigidity of the Cu²⁺ labels compared to the standard nitroxide spin label. Further, the relative orientations between spin-labeled sites were measured to provide insight into the use of double electron-electron resonance (DEER) methods for such labels. The relative mean angles, as well as the standard deviations of the relative angles, agree well in general with the spectral simulations published previously. The fluctuations of relative orientations help rationalize why orientation selectivity effects are minimal at X-band frequencies, but observable at the Q-band for such labels. In summary, the results show that by combining the experimental results with MD simulations precise information about protein conformations as well as flexibility can be obtained.

Unraveling the Impact of Cysteine-to-Serine Mutations on the Structural and Functional Properties of Cu(I)-Binding Proteins

Appropriate maintenance of Cu(I) homeostasis is an essential requirement for proper cell function because its misregulation induces the onset of major human diseases and mortality. For this reason, several research efforts have been devoted to dissecting the inner working mechanism of Cu(I)-binding proteins and transporters. A commonly adopted strategy relies on mutations of cysteine residues, for which Cu(I) has an exquisite complementarity, to serines. Nevertheless, in spite of the similarity between these two amino acids, the structural and functional impact of serine mutations on Cu(I)-binding biomolecules remains unclear. Here, we applied various biochemical and biophysical methods, together with all-atom simulations, to investigate the effect of these mutations on the stability, structure, and aggregation propensity of Cu(I)-binding proteins, as well as their interaction with specific partner proteins. Among Cu(I)-binding biomolecules, we focused on the eukaryotic Atox1-ATP7B system, and the prokaryotic CueR metalloregulator. Our results reveal that proteins containing cysteine-to-serine mutations can still bind Cu(I) ions; however, this alters their stability and aggregation propensity. These results contribute to deciphering the critical biological principles underlying the regulatory mechanism of the in-cell Cu(I) concentration, and provide a basis for interpreting future studies that will take advantage of cysteine-to-serine mutations in Cu(I)-binding systems.

Transition Metal Ion Interactions with Disordered Amyloid-β Peptides in the Pathogenesis of Alzheimer's Disease: Insights from Computational Chemistry Studies

Monomers and oligomers of the amyloid-β peptide aggregate to form the fibrils found in the brains of Alzheimer's disease patients. These monomers and oligomers are largely disordered and can interact with transition metal ions, affecting the mechanism and kinetics of amyloid-β aggregation. Due to the disordered nature of amyloid-β, its rapid aggregation, as well as solvent and paramagnetic effects, experimental studies face challenges in the characterization of transition metal ions bound to amyloid-β monomers and oligomers. The details of the coordination chemistry between transition metals and amyloid-β obtained from experiments remain debated. Furthermore, the impact of transition metal ion binding on the monomeric or oligomeric amyloid-β structures and dynamics are still poorly understood. Computational chemistry studies can serve as an important complement to experimental studies and can provide additional knowledge on the binding between amyloid-β and transition metal ions. Many research groups conducted first-principles calculations, ab initio molecular dynamics simulations, quantum mechanics/classical mechanics simulations, and classical molecular dynamics simulations for studying the interplay between transition metal ions and amyloid-β monomers and oligomers. This review summarizes the current understanding of transition metal interactions with amyloid-β obtained from computational chemistry studies. We also emphasize the current view of the coordination chemistry between transition metal ions and amyloid-β. This information represents an important foundation for future metal ion chelator and drug design studies aiming to combat Alzheimer's disease.

The pivotal role of MBD4–ATP7B in the human Cu(i) excretion path as revealed by EPR experiments and all-atom simulations

Atox1–MBD4 interaction mediates the in-cell Cu(i) concentration.

Trapping intermediates in metal transfer reactions of the CusCBAF export pump of Escherichia coli

Escherichia coli CusCBAF represents an important class of bacterial efflux pump exhibiting selectivity towards Cu(I) and Ag(I). The complex is comprised of three proteins: the CusA transmembrane pump, the CusB soluble adaptor protein, and the CusC outer-membrane pore, and additionally requires the periplasmic metallochaperone CusF. Here we used spectroscopic and kinetic tools to probe the mechanism of copper transfer between CusF and CusB using selenomethionine labeling of the metal-binding Met residues coupled to RFQ-XAS at the Se and Cu edges. The results indicate fast formation of a protein-protein complex followed by slower intra-complex metal transfer. An intermediate coordinated by ligands from each protein forms in 100 ms. Stopped-flow fluorescence of the capping CusF-W44 tryptophan that is quenched by metal transfer also supports this mechanism. The rate constants validate a process in which shared-ligand complex formation assists protein association, providing a driving force that raises the rate into the diffusion-limited regime.

Metal Ion Modeling Using Classical Mechanics

Metal ions play significant roles in numerous fields including chemistry, geochemistry, biochemistry, and materials science. With computational tools increasingly becoming important in chemical research, methods have emerged to effectively face the challenge of modeling metal ions in the gas, aqueous, and solid phases. Herein, we review both quantum and classical modeling strategies for metal ion-containing systems that have been developed over the past few decades. This Review focuses on classical metal ion modeling based on unpolarized models (including the nonbonded, bonded, cationic dummy atom, and combined models), polarizable models (e.g., the fluctuating charge, Drude oscillator, and the induced dipole models), the angular overlap model, and valence bond-based models. Quantum mechanical studies of metal ion-containing systems at the semiempirical, ab initio, and density functional levels of theory are reviewed as well with a particular focus on how these methods inform classical modeling efforts. Finally, conclusions and future prospects and directions are offered that will further enhance the classical modeling of metal ion-containing systems.

Structural and Functional Characterization of a Lytic Polysaccharide Monooxygenase with Broad Substrate Specificity

The recently discovered lytic polysaccharide monooxygenases (LPMOs) carry out oxidative cleavage of polysaccharides and are of major importance for efficient processing of biomass. NcLPMO9C from Neurospora crassa acts both on cellulose and on non-cellulose β-glucans, including cellodextrins and xyloglucan. The crystal structure of the catalytic domain of NcLPMO9C revealed an extended, highly polar substrate-binding surface well suited to interact with a variety of sugar substrates. The ability of NcLPMO9C to act on soluble substrates was exploited to study enzyme-substrate interactions. EPR studies demonstrated that the Cu(2+) center environment is altered upon substrate binding, whereas isothermal titration calorimetry studies revealed binding affinities in the low micromolar range for polymeric substrates that are due in part to the presence of a carbohydrate-binding module (CBM1). Importantly, the novel structure of NcLPMO9C enabled a comparative study, revealing that the oxidative regioselectivity of LPMO9s (C1, C4, or both) correlates with distinct structural features of the copper coordination sphere. In strictly C1-oxidizing LPMO9s, access to the solvent-facing axial coordination position is restricted by a conserved tyrosine residue, whereas access to this same position seems unrestricted in C4-oxidizing LPMO9s. LPMO9s known to produce a mixture of C1- and C4-oxidized products show an intermediate situation.

Predicting the coordination number within copper chaperones: Atox1 as case study

The concentration of copper ions in biological systems is tightly regulated by metallochaperone proteins which are responsible for Cu(I) delivery to designated locations in the cell. These proteins contain a unique motif (MXCXXC) that binds Cu(I) very tightly and specifically but at the same time allows efficient metal transfer to target proteins that often contain a similar copper binding motif. It was found that binding to Cu(I) is achieved through the two cysteine residues in a low coordination number of 2-3 due to possible binding of a third external ligand. Understanding copper transport requires better understanding of copper coordination. Here we therefore focused on establishing a computational method that can predict the coordination number of copper in copper chaperones. The method is shown to be successful in predicting the coordination of Cu(I) within the human copper chaperone (Atox1). Based on the results, a possible rationale for this unique Cu(I) dicoordination in Atox1 is suggested.

The fictile coordination chemistry of cuprous-thiolate sites in copper chaperones

Copper plays vital roles in the active sites of cytochrome oxidase and in several other enzymes essential for human health. Copper is also highly toxic when dysregulated; because of this an elaborate array of accessory proteins have evolved which act as intracellular carriers or chaperones for the copper ions. In most cases chaperones transport cuprous copper. This review discusses some of the chemistry of these copper sites, with a view to some of the structural factors in copper coordination which are important in the biological function of these chaperones. The coordination chemistry and accessible geometries of the cuprous oxidation state are remarkably plastic and we discuss how this may relate to biological function. This article is part of a Special Issue entitled: Biogenesis/Assembly of Respiratory Enzyme Complexes.

Copper-transfer mechanism from the human chaperone Atox1 to a metal-binding domain of Wilson disease protein

The molecular details of how copper (Cu) is transferred from the human Cu chaperone Atox1 to metal-binding domains (MBDs) of P(1B)-type ATPases are still unclear. Here, we use a computational approach, employing quantum mechanics/molecular mechanics (QM/MM) methods, to shed light on the reaction mechanism [probable intermediates, Cu(I) coordination geometries, activation barriers, and energetics] of Cu(I) transfer from Atox1 to the fourth MBD of Wilson disease protein (WD4). Both Atox1 and WD4 have solvent-exposed metal-binding motifs with two Cys residues that coordinate Cu(I). After assessing the existence of all possible 2-, 3- and 4-coordinate Cu-intermediate species, one dominant reaction path emerged. First, without activation barrier, WD4's Cys1 binds Cu(I) in Atox1 to form a 3-coordinated intermediate. Next, with an activation barrier of about 9.5 kcal/mol, a second 3-coordinated intermediate forms that involves both of the Cys residues in WD4 and Cys1 of Atox1. This species can then form the product by decoordination of Atox1's Cys1 (barrier of about 8 kcal/mol). Overall, the Cu-transfer reaction from Atox1 to WD4 appears to be kinetically accessible but less energetically favorable (DeltaE = 7.7 kcal/mol). Our results provide unique insights into the molecular mechanism of protein-mediated Cu(I) transfer in the secretory pathway and are in agreement with existing experimental data.

Comparative analysis of the performance of commonly available density functionals in the determination of geometrical parameters for zinc complexes

A set of 44 Zinc-ligand bond-lengths and of 60 ligand-metal-ligand bond angles from 10 diverse transition-metal complexes, representative of the coordination spheres of typical biological Zn systems, were used to evaluate the performance of a total of 18 commonly available density functionals in geometry determination. Five different basis sets were considered for each density functional, namely two all-electron basis sets (a double-zeta and triple-zeta formulation) and three basis sets including popular types of effective-core potentials: Los Alamos, Steven-Basch-Krauss, and Stuttgart-Dresden. The results show that there are presently several better alternatives to the popular B3LYP density functional for the determination of Zn-ligand bond-lengths and angles. BB1K, MPWB1K, MPW1K, B97-2 and TPSS are suggested as the strongest alternatives for this effect presently available in most computational chemistry software packages. In addition, the results show that the use of effective-core potentials (in particular Stuttgart-Dresden) has a very limited impact, in terms of accuracy, in the determination of metal-ligand bond-lengths and angles in Zinc-complexes, and is a good and safe alternative to the use of an all-electron basis set such as 6-31G(d) or 6-311G(d,p).

Structural organization of human Cu-transporting ATPases: learning from building blocks

Copper-transporting ATPases (Cu-ATPases) ATP7A and ATP7B play an essential role in human physiological function. Their primary function is to deliver copper to the secretory pathway and export excess copper from the cell for removal or further utilization. Cells employ Cu-ATPases in numerous physiological processes that include the biosynthesis of copper-dependent enzymes, lactation, and response to hypoxia. Biochemical studies of human Cu-ATPases and their orthologs have demonstrated that Cu-ATPases share many common structural and mechanistic characteristics with other members of the P-type ATPase family. Nevertheless, the Cu-ATPases have a unique coordinate environment for their ligands, copper and ATP, and additional domains that are required for sophisticated regulation of their intracellular localization and activity. Here, we review recent progress that has been made in understanding the structure of Cu-ATPases from the analysis of their individual domains and orthologs from microorganisms, and speculate about the implications of these findings for the function and regulation of human copper pumps.

Structure and dynamics of Cu(I) binding in copper chaperones Atox1 and CopZ: a computer simulation study

Copper chaperones deliver reduced copper (i.e., Cu(I)) to metal-binding domains of P-type ATPases in the cytoplasm of a range of organisms. Both chaperones and target domains have a ferredoxin-like fold and metal-binding motifs involving two Cys residues. Here, we investigated the Cu-binding geometry and structural dynamics of two homologous Cu(I) chaperones, Homo sapiens Atox1 and Bacillus subtilis CopZ, using a combination of quantum mechanical-molecular mechanics (QM-MM) and classical molecular dynamics (MD) methods. Our QM-MM optimized geometries for the holo- proteins suggested that Cu(I) in Atox1 favors a linear Cys(S)-Cu-Cys(S) arrangement but that this angle is close to 150 degrees in CopZ. Classical MD simulations suggest that both Atox1 and CopZ apo- forms have an increased conformational flexibility as compared to the respective holo- forms. This difference is most pronounced in CopZ and correlates with a lower in vitro thermal stability. Both average fluctuation (i.e., rmsd) and radius of gyration data demonstrate that the effects of Cu(I) coordination extend throughout the proteins. Distinct deviations between the two homologues were found in protein-solvent interactions, entropy of Cu(I) binding, and apo-protein Cys-Cys distance distributions. Our in silico results provide new insights into copper chaperone behavior with direct implications for copper transport mechanisms in vivo.