Development and Application of a Nonbonded Cu(2+) Model That Includes the Jahn-Teller Effect

Engineering Synthetic Electron Transfer Chains from Metallopeptide Membranes

The energetic and geometric features enabling redox chemistry across the copper cupredoxin fold contain key components of electron transfer chains (ETC), which have been extended here by templating the cross-β bilayer assembly of a synthetic nonapeptide, HHQALVFFA-NH2 (K16A), with copper ions. Similar to ETC cupredoxin plastocyanin, these assemblies contain copper sites with blue-shifted (λmax 573 nm) electronic transitions and strongly oxidizing reduction potentials. Electron spin echo envelope modulation and X-ray absorption spectroscopies define square planar Cu(II) sites containing a single His ligand. Restrained molecular dynamics of the cross-β peptide bilayer architecture support metal ion coordination stabilizing the leaflet interface and indicate that the relatively high reduction potential is not simply the result of distorted coordination geometry (entasis). Cyclic voltammetry (CV) supports a charge-hopping mechanism across multiple copper centers placed 10-12 Å apart within the assembled peptide leaflet interface. This metal-templated scaffold accordingly captures the electron shuttle and cupredoxin functionality in a peptide membrane-localized electron transport chain.

Molecular Insights into Cu/Zn Metal Response to the Amyloid β-Peptide (1-42)

Aβ1-40 peptide and Aβ1-42 peptide are the building units of beta-amyloid plaques present in Alzheimer's disease (AD)-affected brain. The binding affinity of various divalent metal ions such as Cu and Zn present in AD-affected brain with different amino acids available in Aβ-peptide became the focus to explore their role in soluble neurotoxic oligomer formation. Cu2+ metal ions are known to enhance the neurotoxicity of the Aβ1-42 peptide by catalyzing the formation of soluble neurotoxic oligomers. The competitive preference of both Cu2+ and Zn2+ simultaneously to interact with the Aβ-peptide is unknown. The divalent Cu and Zn ions were inserted in explicit aqueous Aβ1-42 peptide configurations to get insights into the binding competence of these metal ions with peptides using classical molecular dynamics (MD) simulations. The metal-ion interactions reveal that competitive binding preferences of various peptide sites become metal-ion-specific and differ significantly. For Cu2+, interactions are found to be more significant with respect to those of Asp-7, His-6, Glu-11, and His-14. Asp-1, Glu-3, Asp-7, His-6, Glu-11, and His-13 amino acid residues show higher affinity toward Zn2+ ions. MD simulations show notable variation in the solvent-accessible surface area in the hydrophobic region of the peptide. Infinitesimal mobility was obtained for Zn2+ compared to Cu2+ in an aqueous solution and Cu2+ diffusivity deviated significantly at different time scales, proving its labile features in aqueous Aβ1-42 peptides.

Prion Protein Octarepeat Domain Forms Transient β-Sheet Structures upon Residue-Specific Binding to Cu(II) and Zn(II) Ions

Misfolding of the cellular prion protein (PrP^C) is associated with the development of fatal neurodegenerative diseases called transmissible spongiform encephalopathies (TSEs). Metal ions appear to play a crucial role in PrP^C misfolding. PrP^C is a combined Cu(II) and Zn(II) metal-binding protein, where the main metal-binding site is located in the octarepeat (OR) region. Thus, the biological function of PrP^C may involve the transport of divalent metal ions across membranes or buffering concentrations of divalent metal ions in the synaptic cleft. Recent studies have shown that an excess of Cu(II) ions can result in PrP^C instability, oligomerization, and/or neuroinflammation. Here, we have used biophysical methods to characterize Cu(II) and Zn(II) binding to the isolated OR region of PrP^C. Circular dichroism (CD) spectroscopy data suggest that the OR domain binds up to four Cu(II) ions or two Zn(II) ions. Binding of the first metal ion results in a structural transition from the polyproline II helix to the β-turn structure, while the binding of additional metal ions induces the formation of β-sheet structures. Fluorescence spectroscopy data indicate that the OR region can bind both Cu(II) and Zn(II) ions at neutral pH, but under acidic conditions, it binds only Cu(II) ions. Molecular dynamics simulations suggest that binding of either metal ion to the OR region results in the formation of β-hairpin structures. As the formation of β-sheet structures can be a first step toward amyloid formation, we propose that high concentrations of either Cu(II) or Zn(II) ions may have a pro-amyloid effect in TSE diseases.

Computational Analysis of Triazole-Based Kojic Acid Analogs as Tyrosinase Inhibitors by Molecular Dynamics and Free Energy Calculations

Molecular docking, molecular dynamics (MD) simulations and the linear interaction energy (LIE) method were used here to predict binding modes and free energy for a set of 1,2,3-triazole-based KA analogs as potent inhibitors of Tyrosinase (TYR), a key metalloenzyme of the melanogenesis process. Initially, molecular docking calculations satisfactorily predicted the binding mode of evaluated KA analogs, where the KA part overlays the crystal conformation of the KA inhibitor into the catalytic site of TYR. The MD simulations were followed by the LIE method, which reproduced the experimental binding free energies for KA analogs with an r² equal to 0.97, suggesting the robustness of our theoretical model. Moreover, the van der Waals contributions performed by some residues such as Phe197, Pro201, Arg209, Met215 and Val218 are responsible for the binding recognition of 1,2,3-triazole-based KA analogs in TYR catalytic site. Finally, our calculations provide suitable validation of the combination of molecular docking, MD, and LIE approaches as a powerful tool in the structure-based drug design of new and potent TYR inhibitors.

Refining details of the structural and electronic properties of the CuB site in pMMO enzyme through sequential molecular dynamics/CPKS-EPR calculations

This work investigated the structural and electronic properties of the copper mononuclear site of the PmoB part of the pMMO enzyme at the molecular level. We propose that the Cu_B catalytic site in the soluble portion of pMMO at room temperature and under physiological conditions is a mononuclear copper complex in a distorted octahedral arrangement with the residues His³³, His¹³⁷, and His¹³⁹ on the equatorial base and two water molecules on the axial axis. Our view was based on the molecular dynamics results and DFT calculations of the electronic paramagnetic resonance parameters and comparisons with experimental EPR data. This new proposed model for the Cu_B site brings additional support concerning the recent experimental evidence, which pointed out that a saturated coordination sphere of the copper ion in the Cu_B center is an essential factor that makes it less efficient than the Cu_C site in the methane oxidation. Therefore, according to the Cu_B site model proposed here, an additional step involving a displacement of at least one water molecule of the copper coordination sphere by the O₂ molecule prior to its activation must be necessary. This scenario is less likely to occur in the Cu_C center once this one is buried in the alpha-helices, which are part of the pMMO structure bound to the membrane wall, and consequently located in a less solvent-exposed region. In addition, we also present a simple and efficient sequential S-MD/CPKS protocol to compute EPR parameters that can, in principle, be expanded for the study of other copper-containing proteins.

A Polarizable Cationic Dummy Metal Ion Model

A novel locally polarizable multisite model based on the original cation dummy atom (CDA) model is described for molecular dynamics simulations of ions in condensed phases. Polarization effects are introduced by the electronegativity equalization model (EEM) method where charges on the metal ion and its dummy atoms can fluctuate to respond to the environment. This model includes explicit polarization and ion-induced interactions and can be coupled with nonpolarizable or polarizable water models, making it more transferable to simpler force fields. This approach allows us to enhance the original fixed charge CDA model where the charge distribution cannot adapt to the local solvent structure. To illustrate the new CDApol model, we examined properties of the Zn2+, Al3+, and Zr4+ ions in aqueous solution. The polarizable model and Lennard-Jones parameters were refined for octahedrally coordinated Zn2+, Al3+, and Zr4+ CDAs to reproduce thermodynamic and geometrical properties. Using this locally polarizable model, we were able to obtain the experimental hydration free energy, ion-oxygen distance, and coordination number coupled with the standard 12-6 Lennard-Jones model. This model can be used in myriad additional applications where local polarization and charge transfer is important.

Regulating the Spin State of FeIII Enhances the Magnetic Effect of the Molecular Catalysis Mechanism

Aqueous-phase oxygen evolution reaction (OER) is the bottleneck of water splitting. The formation of the O-O bond involves the generation of paramagnetic oxygen molecules from the diamagnetic hydroxides. The spin configurations might play an important role in aqueous-phase molecular electrocatalysis. However, spintronic electrocatalysis is almost an uncultivated land for the exploration of the oxygen molecular catalysis process. Herein, we present a novel magnetic Fe^III site spin-splitting strategy, wherein the electronic structure and spin states of the Fe^III sites are effectively induced and optimized by the Jahn-Teller effect of Cu²⁺. The theoretical calculations and operando attenuated total reflectance-infrared Fourier transform infrared (ATR FT-IR) reveal the facilitation for the O-O bond formation, which accelerates the production of O₂ from OH^- and improves the OER activity. The Cu₁-Ni₆Fe₂-LDH catalyst exhibits a low overpotential of 210 mV at 10 mA cm^-2 and a low Tafel slope (33.7 mV dec^-1), better than those of the initial Cu₀-Ni₆Fe₂-LDHs (278 mV, 101.6 mV dec^-1). With the Cu²⁺ regulation, we have realized the transformation of NiFe-LDHs from ferrimagnets to ferromagnets and showcase that the OER performance of Cu-NiFe-LDHs significantly increases compared with that of NiFe-LDHs under the effect of a magnetic field for the first time. The magnetic-field-assisted Cu₁-Ni₆Fe₂-LDHs provide an ultralow overpotential of 180 mV at 10 mA cm^-2, which is currently one of the best OER performances. The combination of the magnetic field and spin configuration provides new principles for the development of high-performance catalysts and understandings of the catalytic mechanism from the spintronic level.

Copper(II) import and reduction are dependent on His-Met clusters in the extracellular amino terminus of human copper transporter-1

Copper(I) is an essential metal for all life forms. Though Cu(II) is the most abundant and stable state, its reduction to Cu(I) via an unclear mechanism is prerequisite for its bioutilization. In eukaryotes, the copper transporter-1 (CTR1) is the primary high-affinity copper importer, although its mechanism and role in Cu(II) reduction remain uncharacterized. Here we show that extracellular amino-terminus of human CTR1 contains two methionine-histidine clusters and neighboring aspartates that distinctly bind Cu(I) and Cu(II) preceding its import. We determined that hCTR1 localizes at the basolateral membrane of polarized MDCK-II cells and that its endocytosis to Common-Recycling-Endosomes is regulated by reduction of Cu(II) to Cu(I) and subsequent Cu(I) coordination by the methionine cluster. We demonstrate the transient binding of both Cu(II) and Cu(I) during the reduction process is facilitated by aspartates that also act as another crucial determinant of hCTR1 endocytosis. Mutating the first Methionine cluster (⁷Met-Gly-Met⁹) and Asp¹³ abrogated copper uptake and endocytosis upon copper treatment. This phenotype could be reverted by treating the cells with reduced and nonreoxidizable Cu(I). We show that histidine clusters, on other hand, bind Cu(II) and are crucial for hCTR1 functioning at limiting copper. Finally, we show that two N-terminal His-Met-Asp clusters exhibit functional complementarity, as the second cluster is sufficient to preserve copper-induced CTR1 endocytosis upon complete deletion of the first cluster. We propose a novel and detailed mechanism by which the two His-Met-Asp residues of hCTR1 amino-terminus not only bind copper, but also maintain its reduced state, crucial for intracellular uptake.

Rational Design of Nonbonded Point Charge Models for Highly Charged Metal Cations with Lennard-Jones 12-6 Potential

Here, we developed nonbonded point charge models using a simple Lennard-Jones (LJ) 12-6 potential for highly charged metal cations (18 trivalent and 6 tetravalent ions) for use with 11 water models of TIP3P, OPC3, SPC/E, SPC/E_b, TIP3P-FB, a99SB-disp, TIP4P-Ew, OPC, TIP4P/2005, TIP4P-D, and TIP4P-FB. The designed models simultaneously reproduce the hydration free energy (HFE) and ion-oxygen distance (IOD) in the first hydration shell with an error within 1 kcal/mol and 0.01 Å on average, respectively, and yield reasonable coordination numbers for most cations. Such performance is equivalent to the previously reported point charge models using a more complex 12-6-4 LJ-type potential, while the LJ R parameters of our models are much close to Shannon's revised effective ion radii than that of the 12-6-4 models. Our designed models overestimate the diffusion constants of several trivalent ions by 5-68%. The performance in predicting osmotic coefficients of trivalent chlorides in aqueous solution depends on the salt type. A calibration of cation-anion interacting LJ parameters reproduces the experimental osmotic coefficients of an AlCl₃ solution at 0.2-3.0 mol/L. The effectiveness of our new models is further demonstrated by simulating a metalloprotein system with four force field/water combinations. This work facilitates accurate modeling of metal-containing systems by a variety of force fields and water models in aqueous solution.

Developing and Assessing Nonbonded Dummy Models of Magnesium Ion with Different Hydration Free Energy References

A large diversity in the targeted hydration free energies (HFEs) during model parameterization of metal ions was reported in the literature with a difference by dozens of kcal/mol. Here, we developed a series of nonbonded dummy models of the Mg²⁺ ion targeting different HFE references in TIP3P water, followed by assessments of the designed models in the simulations of MgCl₂ solution and biological systems. Together with the comparison of existing models, we conclude that the difference in the targeted HFEs has a limited influence on the model performance, while the usability of these models differs from case to case. The feasibility of reproducing more properties of Mg²⁺ such as diffusion constants and water exchange rates using a nonbonded dummy model is demonstrated. Underestimated activity derivative and osmotic coefficient of MgCl₂ solutions in high concentration reveal a necessity for further optimization of ion-pair interactions. The developed dummy models are applicable to metal coordination with Asp, Glu, and His residues in metalloenzymes, and the performance in predicting monodentate or bidentate binding modes of Asp/Glu residues depends on the complexity of metal centers and the choice of protein force fields. When both the binding modes coexist, the nonbonded dummy models outperform point charge models, probably in need of considering polarization of metal-binding residues by, for instance, charge calibration in classical force fields. This work is valuable for the use and further development of magnesium ion models for simulations of metal-containing systems with good accuracy.

Gonçalves Vianez JL, Santana da Costa K, de Molfetta FA, Nahum Alves C. Analysis of Kojic Acid Derivatives as Competitive Inhibitors of Tyrosinase: A Molecular Modeling Approach

Tyrosinases belong to the functional copper-containing proteins family, and their structure contains two copper atoms, in the active site, which are coordinated by three histidine residues. The biosynthesis of melanin in melanocytes has two stages depending on the actions of the natural substrates L-DOPA and L-tyrosine. The dysregulation of tyrosinase is involved in skin cancer initiation. In the present study, using molecular modeling tools, we analyzed the inhibition activity of tyrosinase activity using kojic acid (KA) derivatives designed from aromatic aldehydes and malononitrile. All derivatives showed conformational affinity to the enzyme active site, and a favorable distance to chelate the copper ion, which is essential for enzyme function. Molecular dynamics simulations revealed that the derivatives formed promising complexes, presenting stable conformations with deviations between 0.2 and 0.35 Å. In addition, the investigated KA derivatives showed favorable binding free energies. The most stable KA derivatives showed the following binding free energies: -17.65 kcal mol-1 (D6), -18.07 kcal mol-1 (D2), -18.13 (D5) kcal mol-1, and -10.31 kcal mol-1 (D4). Our results suggest that these derivatives could be potent competitive inhibitors of the natural substrates of L-DOPA (-12.84 kcal mol-1) and L-tyrosine (-9.04 kcal mol-1) in melanogenesis.

Characterization of putative transcriptional regulator (PH0140) and its distal homologue

In this study, a phylogenetic tree was constructed using 1854 sequences of various Lrp/AnsC (FFRPs) and ArsR proteins from pathogenic and non-pathogenic organisms. Despite having sequence similarities, FFRPs and ArsR proteins functioning differently as a transcriptional regulator and de-repressor in the presence of exogenous amino acids and metal ions, respectively. To understand these functional differences, the structures of various FFRPs and ArsR proteins (134 sequences) were modeled. Several ArsR proteins exhibited high similarity to the FFRPs while in few proteins, unusual structural folds were observed. However, the Helix-turn-Helix (HTH) domains are common among them and the ligand-binding domains are structurally dissimilar suggest the differences in their binding preferences. Despite low sequence conservation, most of these proteins revealed negatively charged surfaces in the active site pockets. Representative structures (PH0140 and TtArsR protein) from FFRPs and ArsR protein families were considered and evaluated for their functional differences using molecular modeling studies. Our earlier study has explained the binding preference of exogenous Tryptophan and the related transcriptional regulatory mechanism of PH0140 protein. In this study, a Cu²⁺ ion-induced de-repression mechanism of the TtArsR-DNA complex was characterized through docking and molecular dynamics. Further, the proteins were purified and their efficiency for sensing Tryptophan and Cu²⁺ ions were analyzed using cyclic voltammetry. Overall, the study explores the structural evolution and functional difference of FFRPs and ArsR proteins that present the possibilities of PH0140 and TtArsR as potential bio-sensory molecules.

A Rare Natural Benzo[k,l]xanthene as a Turn-Off Fluorescent Sensor for Cu2+ Ion

Rapid and efficient analyses of copper ions are crucial to providing key information for Cu²⁺ in living cells because of their biological importance. In this study, we reported one new turn-off fluorescent sensor for Cu²⁺ with a benzo[k,l]xanthene core, which served as an efficient cation sensor for copper ion over a wide range of other cations (Na⁺, K⁺, Ag⁺, Hg²⁺, Cd²⁺, Co²⁺, Ni²⁺, Zn²⁺, Mg²⁺, and Fe³⁺) owing to the catechol group in the aromatic core. The sensor showed selectivity for Cu²⁺ over other ions; the logK_β for Cu²⁺ binding to compound 1 had a value of 13.265. In the presence of Cu²⁺, sensor 1 provided significant fluorescence decrement; Co²⁺, and Ni²⁺ caused a fluorescence decrement when employed at a higher concentration than Cu²⁺, while Na⁺, K⁺, Hg²⁺, Cd²⁺, Zn²⁺, and Mg²⁺ metal ions produced only minor changes in fluorescence intensity. Fluorescence experiments demonstrate that compound 1 may have an application as a fluorescent probe for detecting Cu²⁺ with a limit of detection of 0.574 µM.

Evaluating the Performance of a Non-Bonded Cu2+ Model Including Jahn-Teller Effect into the Binding of Tyrosinase Inhibitors

Tyrosinase (TYR) is a metalloenzyme classified as a type-3 copper protein, which is involved in the synthesis of melanin through a catalytic process beginning with the conversion of the amino acid l-Tyrosine (l-Tyr) to l-3,4-dihydroxyphenylalanine (l-DOPA). It plays an important role in the mechanism of melanogenesis in various organisms including mammals, plants, and fungi. Herein, we used a combination of computational molecular modeling techniques including molecular dynamic (MD) simulations and the linear interaction energy (LIE) model to evaluate the binding free energy of a set of analogs of kojic acid (KA) in complex with TYR. For the MD simulations, we used a dummy model including the description of the Jahn–Teller effect for Cu²⁺ ions in the active site of this enzyme. Our results show that the LIE model predicts the TYR binding affinities of the inhibitor in close agreement to experimental results. Overall, we demonstrate that the classical model provides a suitable description of the main interactions between analogs of KA and Cu²⁺ ions in the active site of TYR.

Systematic Parametrization of Divalent Metal Ions for the OPC3, OPC, TIP3P-FB, and TIP4P-FB Water Models

Divalent metal ions play important roles in biological and materials systems. Molecular dynamics simulation is an efficient tool to investigate these systems at the microscopic level. Recently, four new water models (OPC3, OPC, TIP3P-FB, and TIP4P-FB) have been developed and better represent the physical properties of water than previous models. Metal ion parameters are dependent on the water model employed, making it necessary to develop metal ion parameters for select new water models. In the present work, we performed parameter scanning for the 12-6 Lennard-Jones nonbonded model of divalent metal ions in conjunction with the four new water models as well as four previous water models (TIP3P, SPC/E, TIP4P, and TIP4P-Ew). We found that these new three-point and four-point water models provide comparable or significantly improved performance for the simulation of divalent metal ions when compared to previous water models in the same category. Among all eight water models, the OPC3 water model yields the best performance for the simulation of divalent metal ions in the aqueous phase when using the 12-6 model. On the basis of the scanning results, we independently parametrized the 12-6 model for 24 divalent metal ions with each of the four new water models. As noted previously, the 12-6 model still fails to simultaneously reproduce the experimental hydration free energy (HFE) and ion-oxygen distance (IOD) values even with these new water models. To solve this problem, we parametrized the 12-6-4 model for the 16 divalent metal ions for which we have both experimental HFE and IOD values for each of the four new water models. The final parameters are able to reproduce both the experimental HFE and IOD values accurately. To validate the transferability of our parameters, we carried out benchmark calculations to predict the energies and geometries of ion-water clusters as well as the ion diffusivity coefficient of Mg²⁺. By comparison to quantum chemical calculations and experimental data, these results show that our parameters are well designed and have excellent transferability. The metal ion parameters for the 12-6 and 12-6-4 models reported herein can be employed in simulations of various biological and materials systems when using the OPC3, OPC, TIP3P-FB, or TIP4P-FB water model.

Site-Specific Interactions with Copper Promote Amyloid Fibril Formation for λ6aJL2-R24G

Light-chain amyloidosis (AL) is one of the most common systemic amyloidoses, and it is characterized by the deposition of immunoglobulin light chain (LC) variable domains as insoluble amyloid fibers in vital organs and tissues. The recombinant protein 6aJL2-R24G contains λ6a and JL2 germline genes and also contains the Arg24 by Gly substitution. This mutation is present in 25% of all amyloid-associated λ6 LC cases, reduces protein stability, and increases the propensity to form amyloid fibers. In this study, it was found that the interaction of 6aJL2-R24G with Cu(II) decreases the thermal stability of the protein and accelerates the amyloid fibril formation, as observed by fluorescence spectroscopy. Isothermal calorimetry titration showed that Cu(II) binds to the protein with micromolar affinity. His99 may be one of the main Cu(II) interaction sites, as observed by nuclear magnetic resonance spectroscopy. The binding of Cu(II) to His99 induces larger fluctuations of the CDR1 and loop C″, as shown by molecular dynamics simulations. Thus, Cu(II) binding may be inducing the loss of interactions between CDR3 and CDR1, making the protein less stable and more prone to form amyloid fibers. This study provides insights into the mechanism of metal-induced aggregation of the 6aJL2-R24G protein and sheds light on the bio-inorganic understanding of AL disease.

Comparing Nonbonded Metal Ion Models in the Divalent Cation Binding Protein PsaA

Divalent metal cations are essential for many biological processes; however, accurately modeling divalent metal ions has proved a significant challenge for molecular dynamics force fields. Here we show that the choice of ion model influences the observed dynamics in PsaA, a metal binding protein from Streptococcus pneumoniae. We conduct extensive unbiased simulations and free energy calculations of PsaA bound to its cognate ligand Mn²⁺ and inhibitory ligand Zn²⁺ using three nonbonded ion models: a 12-6 model, a 12-6-4 model, and a multisite model. The observed coordination geometries and metal binding dynamics are sensitive to the choice of ion model, with the most dramatic differences observed in free energy calculations of ion release. We show that the conformational ensemble of Mn-bound PsaA is more similar to the crystallographic metal bound open state. This work extends the current model of PsaA metal binding and provides a framework for the rationalization of experimentally determined metal binding behavior. Our findings support the use of the 12-6-4 ion model for further simulations of divalent cation binding proteins.

Structures of the solvated copper(ii) ion in ammonia at various temperatures

We investigated theoretically the structures and relative stabilities of the solvated copper(ii) ion in ammonia, Cu²⁺(NH₃)_n, n = 1–10.

Understanding Carbamate Formation Reaction Thermochemistry of Amino Acids as Solvents for Postcombustion CO2 Capture

The carbamate stability constant for a data set of 10 amino acids, having potential for being postcombustion CO₂ capture (PCC) solvents, has been calculated using various implicit and explicit solvation shell models. This work also includes an extensive study of gas-phase free energy and enthalpy for the amino acid carbamate formation reaction with the Hartree Fock method, density functional methods [B3LYP/6-311++G(d,p)], and composite methods (G3MP2B3, G3MP2, CBS-QB3, and G4MP2). Ideal PCC solvent properties require finding a profitable tradeoff between various thermodynamic and system optimization parameters. Benchmark gaseous-phase and solution-phase thermodynamic properties given in this work can help in making informed decisions when choosing promising PCC solvents. The temperature dependency of the carbamate stability constant of amino acids is predicted using PCM and SM8T implicit solvation models. PCC is a temperature swing absorption-desorption process, and the high-temperature sensitivity of the ln K_c^AmCOO- value is of vital importance in attaining cost-efficient processes.

Intrinsically disordered proteins in various hypotheses on the pathogenesis of Alzheimer's and Parkinson's diseases

Amyloid-β (Aβ) and α-synuclein (αS) are two intrinsically disordered proteins (IDPs) at the centers of the pathogenesis of Alzheimer's and Parkinson's diseases, respectively. Different hypotheses have been proposed for explanation of the molecular mechanisms of the pathogenesis of these two diseases, with these two IDPs being involved in many of these hypotheses. Currently, we do not know, which of these hypothesis is more accurate. Experiments face challenges due to the rapid conformational changes, fast aggregation processes, solvent and paramagnetic effects in studying these two IDPs in detail. Furthermore, pathological modifications impact their structures and energetics. Theoretical studies using computational chemistry and computational biology have been utilized to understand the structures and energetics of Aβ and αS. In this chapter, we introduce Aβ and αS in light of various hypotheses, and discuss different experimental and theoretical techniques that are used to study these two proteins along with their weaknesses and strengths. We suggest that a promising solution for studying Aβ and αS at the center of varying hypotheses could be provided by developing new techniques that link quantum mechanics, statistical mechanics, thermodynamics, bioinformatics to machine learning. Such new developments could also lead to development in experimental techniques.

Simulations of CYP51A from Aspergillus fumigatus in a model bilayer provide insights into triazole drug resistance

Azole antifungal drugs target CYP51A in Aspergillus fumigatus by binding with the active site of the protein, blocking ergosterol biosynthesis. Resistance to azole antifungal drugs is now common, with a leucine to histidine amino acid substitution at position 98 the most frequent, predominantly conferring resistance to itraconazole, although cross-resistance has been reported in conjunction with other mutations. In this study, we create a homology model of CYP51A using a recently published crystal structure of the paralog protein CYP51B. The derived structures, wild type, and L98H mutant are positioned within a lipid membrane bilayer and subjected to molecular dynamics simulations in order improve the accuracy of both models. The structural analysis from our simulations suggests a decrease in active site surface from the formation of hydrogen bonds between the histidine substitution and neighboring polar side chains, potentially preventing the binding of azole drugs. This study yields a biologically relevant structure and set of dynamics of the A. fumigatus Lanosterol 14 alpha-demethylase enzyme and provides further insight into azole antifungal drug resistance.

A unified picture of the covalent bond within quantum-accurate force fields: From organic molecules to metallic complexes' reactivity

Computational studies of chemical processes taking place over extended size and time scales are inaccessible by electronic structure theories and can be tackled only by atomistic models such as force fields. These have evolved over the years to describe the most diverse systems. However, as we improve the performance of a force field for a particular physical/chemical situation, we are also moving away from a unified description. Here, we demonstrate that a unified picture of the covalent bond is achievable within the framework of machine learning-based force fields. Ridge regression, together with a representation of the atomic environment in terms of bispectrum components, can be used to map a general potential energy surface for molecular systems at chemical accuracy. This protocol sets the ground for the generation of an accurate and universal class of potentials for both organic and organometallic compounds with no specific assumptions on the chemistry involved.

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.

Oxidase or peptidase? A computational insight into a putative aflatoxin oxidase from Armillariella tabescens

Aflatoxin oxidase (AFO), an enzyme isolated from Armillariella tabescens, has been reported to degrade aflatoxin B1 (AFB1). However, recent studies reported sequence and structure similarities with the dipeptidyl peptidase III (DPP III) family of enzymes and confirmed peptidase activity toward DPP III substrates. In light of these investigations, an extensive computational study was performed in order to improve understanding of the AFO functions. Steered MD simulations revealed long-range domain motions described as protein opening, characteristic for DPPs III and necessary for substrate binding. Newly identified open and partially open forms of the enzyme closely resemble those of the human DPP III orthologue. Docking of a synthetic DPP III substrate Arg₂ -2-naphthylamide revealed a binding mode similar to the one found in crystal structures of human DPP III complexes with peptides with the S1 and S2 subsites' amino acid residues conserved. On the other hand, no energetically favorable AFB1 binding mode was detected, suggesting that aflatoxins are not good substrates of AFO. High plasticity of the zinc ion coordination sphere within the active site, consistent with that of up to date studied DPPs III, was observed as well. A detailed electrostatic analysis of the active site revealed a predominance of negatively charged regions, unsuitable for the binding of the neutral AFB1. The present study is in line with the most recent experimental study on this enzyme, both suggesting that AFO is a typical member of the DPP III family.

Hydrolytic zinc metallopeptides using a computational multi-state design approach

Combination of multi-state design and long-timescale conformational dynamics as a powerful strategy to obtain metalloenzymes.

Assessing the Performance of the Nonbonded Mg2+ Models in a Two-Metal-Dependent Ribonuclease

Magnesium ions (Mg²⁺), abundant in living cells, are essential for biomolecular structure, dynamics, and function. The biological importance of Mg²⁺ has motivated continuous development and improvement of various Mg²⁺ models for molecular dynamics (MD) simulations during the last decades. There are four types of nonbonded Mg²⁺ models: the point charge models based on a 12-6 or 12-6-4 type Lennard-Jones (LJ) potential, and the multisite models based on a 12-6 or 12-6-4 LJ potential. Here, we systematically assessed the performance of these four types of nonbonded Mg²⁺ models (21 models in total) in terms of maintaining a challenging intermediate state configuration captured in the structure of a prototypical two-metal-ion RNase H complex with an RNA/DNA hybrid. Our data demonstrate that the 12-6-4 multisite models, which account for charge-induced dipole interactions, perform the best in reproducing all the unique coordination modes in this intermediate state and maintaining the correct carboxylate denticity. Our benchmark work provides a useful guideline for MD simulations and structural refinement of Mg²⁺-containing biomolecular systems.

Single-Ion Thermodynamics from First Principles: Calculation of the Absolute Hydration Free Energy and Single-Electrode Potential of Aqueous Li+ Using ab Initio Quantum Mechanical/Molecular Mechanical Molecular Dynamics Simulations

A recently proposed thermodynamic integration (TI) approach formulated in the framework of quantum mechanical/molecular mechanical molecular dynamics (QM/MM MD) simulations is applied to study the structure, dynamics, and absolute intrinsic hydration free energy Δ_s G_M⁺,wat^◦ of the Li⁺ ion at a correlated ab initio level of theory. Based on the results, standard values (298.15 K, ideal gas at 1 bar, ideal solute at 1 molal) for the absolute intrinsic hydration free energy [Formula: see text] of the proton, the surface electric potential jump χ_wat^◦ upon entering bulk water, and the absolute single-electrode potential [Formula: see text] of the reference hydrogen electrode are calculated to be -1099.9 ± 4.2 kJ·mol^-1, 0.13 ± 0.08 V, and 4.28 ± 0.04 V, respectively, in excellent agreement with the standard values recommended by Hünenberger and Reif on the basis of an extensive evaluation of the available experimental data (-1100 ± 5 kJ·mol^-1, 0.13 ± 0.10 V, and 4.28 ± 0.13 V). The simulation results for Li⁺ are also compared to those for Na⁺ and K⁺ from a previous study in terms of relative hydration free energies ΔΔ_s G_M⁺,wat^◦ and relative electrode potentials [Formula: see text]. The calculated values are found to agree extremely well with the experimental differences in standard conventional hydration free energies ΔΔ_s G_M⁺,wat^• and redox potentials [Formula: see text]. The level of agreement between simulation and experiment, which is quantitative within error bars, underlines the substantial accuracy improvement achieved by applying a highly demanding QM/MM approach at the resolution-of-identity second-order Møller-Plesset perturbation (RIMP2) level over calculations relying on purely molecular mechanical or density functional theory (DFT) descriptions. A detailed analysis of the structural and dynamical properties of the Li⁺ hydrate indicates that a correct description of the solvation structure and dynamics is achieved as well at this level of theory. Consideration of the QM/MM potential-energy components also shows that the partitioning into QM and MM zones does not induce any significant energetic artifact for the system considered.

Empirical Valence Bond Simulations Suggest a Direct Hydride Transfer Mechanism for Human Diamine Oxidase

Diamine oxidase (DAO) is an enzyme involved in the regulation of cell proliferation and the immune response. This enzyme performs oxidative deamination in the catabolism of biogenic amines, including, among others, histamine, putrescine, spermidine, and spermine. The mechanistic details underlying the reductive half-reaction of the DAO-catalyzed oxidative deamination which leads to the reduced enzyme cofactor and the aldehyde product are, however, still under debate. The catalytic mechanism was proposed to involve a prototropic shift from the substrate-Schiff base to the product-Schiff base, which includes the rate-limiting cleavage of the Cα-H bond by the conserved catalytic aspartate. Our detailed mechanistic study, performed using a combined quantum chemical cluster approach with empirical valence bond simulations, suggests that the rate-limiting cleavage of the Cα-H bond involves direct hydride transfer to the topaquinone cofactor-a mechanism that does not involve the formation of a Schiff base. Additional investigation of the D373E and D373N variants supported the hypothesis that the conserved catalytic aspartate is indeed essential for the reaction; however, it does not appear to serve as the catalytic base, as previously suggested. Rather, the electrostatic contributions of the most significant residues (including D373), together with the proximity of the Cu2+ cation to the reaction site, lower the activation barrier to drive the chemical reaction.

Challenges and advances in the computational modeling of biological phosphate hydrolysis

Phosphate ester hydrolysis is fundamental to many life processes, and has been the topic of substantial experimental and computational research effort. However, even the simplest of phosphate esters can be hydrolyzed through multiple possible pathways that can be difficult to distinguish between, either experimentally, or computationally. Therefore, the mechanisms of both the enzymatic and non-enzymatic reactions have been historically controversial. In the present contribution, we highlight a number of technical issues involved in reliably modeling these computationally challenging reactions, as well as proposing potential solutions. We also showcase examples of our own work in this area, discussing both the non-enzymatic reaction in aqueous solution, as well as insights obtained from the computational modeling of organophosphate hydrolysis and catalytic promiscuity amongst enzymes that catalyze phosphoryl transfer.

Computer simulations of the catalytic mechanism of wild-type and mutant β-phosphoglucomutase

β-Phosphoglucomutase (β-PGM) has served as an important model system for understanding biological phosphoryl transfer.

Extending the Nonbonded Cationic Dummy Model to Account for Ion-Induced Dipole Interactions

Modeling metalloproteins often requires classical molecular dynamics (MD) simulations in order to capture their relevant motions, which in turn necessitates reliable descriptions of the metal centers involved. One of the most successful approaches to date is provided by the "cationic dummy model", where the positive charge of the metal ion is transferred toward dummy particles that are bonded to the central metal ion in a predefined coordination geometry. While this approach allows for ligand exchange, and captures the correct electrostatics as demonstrated for different divalent metal ions, current dummy models neglect ion-induced dipole interactions. In the present work, we resolve this weakness by taking advantage of the recently introduced 12-6-4 type Lennard-Jones potential to include ion-induced dipole interactions. We revise our previous dummy model for Mg²⁺ and demonstrate that the resulting model can simultaneously reproduce the experimental solvation free energy and metal-ligand distances without the need for artificial restraints or bonds. As ion-induced dipole interactions become particularly important for highly charged metal ions, we develop dummy models for the biologically relevant ions Al³⁺, Fe³⁺, and Cr³⁺. Finally, the effectiveness of our new models is demonstrated in MD simulations of several diverse (and highly challenging to simulate) metalloproteins.

Free-Energy Calculations of Ionic Hydration Consistent with the Experimental Hydration Free Energy of the Proton

Computational free-energy correction strategies and the choice of experimental proton hydration free energy, ΔG_s^*(H⁺), are analyzed to investigate the apparent controversy in experimental thermodynamics of ionic hydration. Without corrections, the hydration free-energy (ΔG_hyd) calculations match experiments with ΔG_s^*(H⁺) = -1064 kJ/mol as reference. Using the Galvani surface potential the resulting (real) ΔG_hyd are consistent with ΔG_s^*(H⁺) = -1098 kJ/mol. When applying, in an ad hoc manner, the discrete solvent correction, ΔG_hyd matching the "consensus" ΔG_s^*(H⁺) of -1112 kJ/mol are obtained. This analysis rationalizes reports on ΔG_hyd calculations for ions using different experimental references. For neutral amino acid side chains ΔG_hyd are independent of the water model, whereas there are large differences in ΔG_hyd due to the water model for charged species, suggesting that long-range ordering of water around ions yields an important contribution to the ΔG_hyd. These differences are reduced significantly when applying consistent corrections, but to obtain the most accurate results it is recommended to use the water model belonging to the force field.

Development of a multisite model for Ni(II) ion in solution from thermodynamic and kinetic data

Force-field parameters are developed for a multisite model of Ni(II) ions to be used in molecular dynamics simulations combined to enhanced sampling methods. The performances of two charge-partitioning schemes are validated by taking into account structural, thermodynamic, and kinetic observables. One of the two models, featuring partial charges on the dummy atoms only, matches both Ni(II) free energy of solvation and water exchange rates. Such model is particularly suited to study complexation events at a fully dynamic description. © 2017 Wiley Periodicals, Inc.

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.

Excitation Dependent Phosphorous Property and New Model of the Structured Green Luminescence in ZnO

The copper induced green luminescence (GL) with two sets of fine structures in ZnO crystal has been found for several decades (i.e., R. Dingle, Phys. Rev. Lett. 23, 579 (1969)), but the physical origin of the doublet still remains as an open question up to now. In this paper, we provide new insight into the mechanism of the structured GL band in terms of new experimental findings and theoretical calculations. It is found, for the first time, that the GL signal exhibits persistent afterglow for tens of minutes after the switch-off of below-band-gap excitation light but it cannot occur under above-band-gap excitation. Such a phosphorous property may be interpreted as de-trapping and feeding of electrons from a shallow trapping level via the conduction band to the Cu-related luminescence centers where the Cu³⁺ ion is proposed to work as the final state of the GL emission. From first-principles calculation, such a Cu³⁺ ion in wurtzite ZnO prefers a high spin 3d⁸ state with two non-degenerated half-filled orbitals due to the Jahn-Teller effect, probably leading to the double structures in photoluminescence spectrum. Therefore, this model gives a comprehensively new understanding on the mechanism of the structured GL band in ZnO.

Characterization of Mn(II) ion binding to the amyloid-β peptide in Alzheimer's disease

Growing evidence links neurodegenerative diseases to metal exposure. Aberrant metal ion concentrations have been noted in Alzheimer's disease (AD) brains, yet the role of metals in AD pathogenesis remains unresolved. A major factor in AD pathogenesis is considered to be aggregation of and amyloid formation by amyloid-β (Aβ) peptides. Previous studies have shown that Aβ displays specific binding to Cu(II) and Zn(II) ions, and such binding has been shown to modulate Aβ aggregation. Here, we use nuclear magnetic resonance (NMR) spectroscopy to show that Mn(II) ions also bind to the N-terminal part of the Aβ(1-40) peptide, with a weak binding affinity in the milli- to micromolar range. Circular dichroism (CD) spectroscopy, solid state atomic force microscopy (AFM), fluorescence spectroscopy, and molecular modeling suggest that the weak binding of Mn(II) to Aβ may not have a large effect on the peptide's aggregation into amyloid fibrils. However, identification of an additional metal ion displaying Aβ binding reveals more complex AD metal chemistry than has been previously considered in the literature.

Probing the mechanisms for the selectivity and promiscuity of methyl parathion hydrolase

Diverse organophosphate hydrolases have convergently evolved the ability to hydrolyse man-made organophosphates. Thus, these enzymes are attractive model systems for studying the factors shaping enzyme functional evolution. Methyl parathion hydrolase (MPH) is an enzyme from the metallo-β-lactamase superfamily, which hydrolyses a wide range of organophosphate, aryl ester and lactone substrates. In addition, MPH demonstrates metal-ion-dependent selectivity patterns. The origins of this remain unclear, but are linked to open questions about the more general role of metal ions in functional evolution and divergence within enzyme superfamilies. Here, we present detailed mechanistic studies of the paraoxonase and arylesterase activities of MPH complexed with five different transition metal ions, and demonstrate that the hydrolysis reactions proceed via similar pathways and transition states. However, while it is possible to discern a clear structural origin for the selectivity between different substrates, the selectivity between different metal ions appears to lie instead in the distinct electrostatic properties of the metal ions themselves, which causes subtle changes in transition state geometries and metal-metal distances at the transition state rather than significant structural changes in the active site. While subtle, these differences can be significant for shaping the metal-ion-dependent activity patterns observed for this enzyme.This article is part of the themed issue 'Multiscale modelling at the physics-chemistry-biology interface'.