Atomistic molecular simulations of Aß-Zn conformational ensembles.

The amyloid-forming Aß peptide is able to interact with metal cations to form very stable complexes that influence fibril formation and contribute to the onset of Alzheimer's disease. Multiple structures of peptides derived from Aß in complex with different metals have been resolved experimentally to provide an atomic-level description of the metal-protein interactions. However, Aß is intrinsically disordered, and hence more amenable to an ensemble description. Molecular dynamics simulations can now reach the timescales needed to generate ensembles for these type of complexes. However, this requires accurate force fields both for the protein and the protein-metal interactions. Here we use state-of-the-art methods to generate force field parameters for the Zn(II) cations in a set of complexes formed by different Aß variants and combine them with the Amber99SB*-ILDN optimized force field. Upon comparison of NMR experiments with the simulation results, further optimized with a Bayesian/Maximum entropy approach, we provide an accurate description of the molecular ensembles for most Aß-metal complexes. We find that the resulting conformational ensembles are more heterogeneous than the NMR models deposited in the Protein Data Bank.

All-Purpose Measure of Electron Correlation for Multireference Diagnostics.

We present an analytical relationship between two natural orbital occupancy-based indices, IND¯ and INDmax, and two established electron correlation metrics: the leading term of a configuration interaction expansion, c0, and the D2 diagnostic. Numerical validation revealed that IND¯ and INDmax can effectively substitute for c0 and D2, respectively. These indices offer three distinct advantages: (i) they are universally applicable across all electronic structure methods, (ii) their interpretation is more intuitive, and (iii) they can be readily incorporated into the development of hybrid electronic structure methods. Additionally, we draw a distinction between correlation measures and correlation diagnostics, establishing MP2 and CCSD numerical thresholds for INDmax, which are to be used as a multireference diagnostic. Our findings further demonstrate that establishing thresholds for other electronic structure methods can be easily accomplished using small data sets.

Rules governing metal coordination in Aß-Zn(II) complex models from quantum mechanical calculations.

Transition metals directly contribute to the neurotoxicity of the aggregates of the amyloid-forming Aß peptide. The understanding and rationalization of the coordination modes of metals to Aß amyloid is, therefore, of paramount importance to understand the capacity of a given metal to promote peptide aggregation. Experimentally, multiple Aß-metal structures have been resolved, which exhibit different modes of coordination in both the monomeric and oligomeric forms of Aß. Although Zn(II) metalloproteins are very abundant and often involve cysteine residues in the first coordination shell, in the case of Aß-Zn(II), though, Zn(II) is coordinated by glutamic/aspartic acid and/or histidine residues exclusively, making for an interesting case study. Here we present a systematic analysis of the underlying chemistry on Aß-Zn(II) coordination, where relative stabilities of different coordination arrangements indicate that a mixture of Glu/Asp and His residues is favored. A detailed comparison between different coordination shell geometries shows that tetrahedral coordination is generally favored in the aqueous phase. Our calculations show an interplay between dative covalent interactions and electrostatics which explains the observed trends. Multiple structures deposited in the Protein Data Bank support our findings, suggesting that the trends found in our work may be transferable to other Zn(II) metalloproteins with this type of coordination.

Influence of metal binding on the conformational landscape of neurofilament peptides.

In order to understand the preferred modes of chelation in metal-binding peptides, quantum mechanical calculations can be used to compute energies, resulting in a hierarchy of binding affinities. These calculations often produce increasing stabilization energies the higher the coordination of the complex. However, as the coordination of a metal increases, the conformational freedom of the polypeptide chain is inevitably reduced, resulting in an entropic penalty. Estimating the magnitude of this penalty from the many different degrees of freedom of biomolecular systems is very challenging, and as a result this contribution to the free energy is often ignored. Here we explore this problem focusing on a family of phosphorylated neuropeptides that bind to aluminum. We find that there is a general negative correlation between both stabilization energy and entropy. Our results suggest that a subtle interplay between enthalpic and entropic forces will determine the population of the most favourable species. Additionally, we discuss the requirements for a possible "Metal Ion Hypothesis" based on our findings.

One- and Two-Electron Reductions in MiniSOG and their Implication in Catalysis.

The unconventional bioorthogonal catalytic activation of anticancer metal complexes by flavin and flavoproteins photocatalysis has been reported recently. The reactivity is based on a two-electron redox reaction of the photoactivated flavin. Furthermore, when it comes to flavoproteins, we recently reported that site mutagenesis can modulate and improve this catalytic activity in the mini Singlet Oxygen Generator protein (SOG). In this paper, we analyze the reductive half-reaction in different miniSOG environments by means of density functional theory. We report that the redox properties of flavin and the resulting reactivity of miniSOG is modulated by specific mutations, which is in line with the experimental results in the literature. This modulation can be attributed to the fundamental physicochemical properties of the system, specifically (i) the competition of single and double reduction of the flavin and (ii) the probability of electron transfer from the protein to the flavin. These factors are ultimately linked to the stability of flavin's electron-accepting orbitals in different coordination modes.

The Wigner localization of interacting electrons in a one-dimensional harmonic potential.

In this work, we study the Wigner localization of interacting electrons that are confined to a quasi-one-dimensional harmonic potential using accurate quantum chemistry approaches. We demonstrate that the Wigner regime can be reached using small values of the confinement parameter. To obtain physical insight in our results, we analyze them with a semi-analytical model for two electrons. Thanks to electronic-structure properties such as the one-body density and the particle-hole entropy, we are able to define a path that connects the Wigner regime to the Fermi-gas regime by varying the confinement parameter. In particular, we show that the particle-hole entropy, as a function of the confinement parameter, smoothly connects the two regimes. Moreover, it exhibits a maximum that could be interpreted as the transition point between the localized and delocalized regimes.

Unravelling the binding affinity and selectivity of molybdenum(II) phenanthroline complexes with DNA G-quadruplexes by using linear-scaling DFT studies. The important role of ancillary ligands.

We have used near linear-scaling density functional theory (LS-DFT) methods including dispersion, for the first time, to study the interaction of two isomers, equatorial (Eq) and axial (Ax), of the [Mo(η3-C3H5)Br(CO)2(phen)] metal complex with the DNA G-quadruplexes (GQ) to gain insight into its cytotoxicity. The LMKLL/DZDP level of calculation, which includes van der Waals contributions, with the SIESTA software was used to treat by means of first-principles computations the whole biological studied model system with ∼1000 atoms. Computed formation energies point to systems containing the Ax isomer as the most stable although the nearest system in energy containing the Eq isomer is only 7.5 kcal mol-1 above. On the other hand, the energy decomposition analysis (EDA) favours interaction energies for the systems containing the Eq isomer. However, when solvent effects are taken into account the systems containing the Ax isomer are again the most stable. This Ax isomer was found interacting by means of end-stacking with the GQ and surprisingly totally inside the non-canonical secondary structure, where all the ligands of the metal complex produce several weak interactions with the DNA structure. On the other hand, the Eq isomer prefers to interact from outside by means of intercalation in which the ancillary ligands also have some role in the interaction. Such features and comparison with the results regarding the interaction of the [Mo(η3-C3H5)Br(CO)2(phen)] metal complex with duplex DNA suggest that the [Mo(η3-C3H5)Br(CO)2(phen)] would have a higher affinity and eventual selectivity for non-canonical DNA GQ structures.

Hydrogen Tunneling in Catalytic Hydrolysis and Alcoholysis of Silanes.

An unprecedented quantum tunneling effect has been observed in catalytic Si-H bond activations at room temperature. The cationic hydrido-silyl-iridium(III) complex, {Ir[SiMe(o-C6 H4 SMe)2 ](H)(PPh3 )(THF)}[BArF 4 ], has proven to be a highly efficient catalyst for the hydrolysis and the alcoholysis of organosilanes. When triethylsilane was used as a substrate, the system revealed the largest kinetic isotopic effect (KIESi-H/Si-D =346±4) ever reported for this type of reaction. This unexpectedly high KIE, measured at room temperature, together with the calculated Arrhenius preexponential factor ratio (AH /AD =0.0004) and difference in the observed activation energy [(E a D -E a H )=34.07 kJ mol-1 ] are consistent with the participation of quantum tunneling in the catalytic process. DFT calculations have been used to unravel the reaction pathway and identify the rate-determining step. Aditionally, isotopic effects were considered by different methods, and tunneling effects have been calculated to be crucial in the process.

Study of the DNA binding mechanism and in vitro activity against cancer cells of iron(III) and aluminium(III) kojic acid derivative complexes.

Metal ions have unique electrochemical and spectroscopical properties that cannot be attained by purely organic compounds. Most of the metal ions are toxic to humans, but paradoxically, metallodrugs are used in medicine as therapeutics and theranostics. Metallodrugs are eliminated in urine and faeces, and therefore release toxic metals and ligands into aquatic ecosystems, thereby raising concerns regarding environmental risks. The use of metallodrugs based on essential metal ions (i.e., iron, copper and zinc), instead of toxic ions, is a new alternative with minor hazards. Kojic acid is an Asperigillus oryzae metabolite of low toxicity used in the food and cosmetics industries. Its derivatives form stable complexes with iron(III) ions, which bind effectively to DNA and inhibit DNA polymerization. The iron(III)/S2 ligand complexes reduce in vitro colon carcinoma (Caco2) cell viability and significantly decrease the cell number. The kojic acid derivative complexes with iron(III) presented here are an alternative to the currently used platinum complexes in cancer therapy.

Semi-empirical and linear-scaling DFT methods to characterize duplex DNA and G-quadruplexes in the presence of interacting small molecules.

The computational study of DNA and its interaction with ligands is a highly relevant area of research, with significant consequences for developing new therapeutic strategies. However, the computational description of such large and complex systems requires considering interactions of different types simultaneously in a balanced way, such as non-covalent weak interactions (namely hydrogen bonds and stacking), metal-ligand interactions, polarisation and charge transfer effects. All these considerations imply a real challenge for computational chemistry. The possibility of studying large biological systems using quantum methods for the entire system requires significant computational resources, with improvements in parallelisation and optimisation of theoretical strategies. Computational methods, such as Linear-Scaling Density Functional Theory (LS-DFT) and DLPNO-CCSD(T), may allow performing ab initio quantum mechanics calculations, including the electronic structure for large biological systems, in a reasonable computing time. In this work, we study the interaction of small molecules and cations with DNA (both duplex DNA and G-quadruplexes), comparing different computational methods: a LS-DFT method at the LMKLL/DZDP level of theory, semi-empirical methods (PM6-DH2 and PM7), mixed QM/MM, and DLPNO-CCSD(T). Our goal is to demonstrate the adequacy of LS-DFT to treat the different types of interactions present in DNA-dependent systems. We show that LMKLL/DZDP using SIESTA can yield very accurate geometries and energetics in all the different systems considered in this work: duplex DNA (dDNA), phenanthroline intercalating dDNA, G-quadruplexes, and metal-G-tetrads considering alkaline metals of different sizes. As far as we know, this is the first time that full G-quadruplex geometry optimisations have been carried out using a DFT method thanks to its linear-scaling capabilities. Moreover, we show that LS-DFT provides high-quality structures, and some semi-empirical Hamiltonians can also yield suitable geometries. However, DLPNO-CCSD(T) and LS-DFT are the only methods that accurately describe interaction energies for all the systems considered in our study.

Influence of the Nonprotein Amino Acid Mimosine in Peptide Conformational Propensities from Novel Amber Force Field Parameters.

Mimosine is a nonprotein amino acid derived from plants known for its ability to bind to divalent and trivalent metal cations such as Zn2+, Ni2+, Fe2+, or Al3+. This results in interesting antimicrobial and anticancer properties, which make mimosine a promising candidate for therapeutic applications. One possibility is to incorporate mimosine into synthetic short peptide drugs. However, how this amino acid affects the peptide structure is not well understood, reducing our ability to design effective therapeutic compounds. In this work, we used computer simulations to understand this question. We first built parameters for the mimosine residue to be used in combination with two classical force fields of the Amber family. Then, we used atomistic molecular dynamics simulations with the resulting parameter sets to evaluate the influence of mimosine in the structural propensities for this amino acid. We compared the results of these simulations with homologous peptides, where mimosine is replaced by either phenylalanine or tyrosine. We found that the strong dipole in mimosine induces a preference for conformations where the amino acid rings are stacked over more extended conformations. We validated our results using quantum mechanical calculations, which provide a robust foundation for the outcome of our classical simulations.

Flavin-mediated photoactivation of Pt(IV) anticancer complexes: computational insights on the catalytic mechanism.

The mechanism for the photocatalytic activation of Pt(IV) anticancer prodrugs by riboflavin in the presence of NADH has been investigated by DFT. In the first step of the reaction, the oxidation kinetics of NADH to afford the catalytically active riboflavin hydroquinone is dramatically favoured by generation of the flavin triplet excited state. In the triplet, formation of a π-π stacked adduct promotes the hydride transfer from NADH to riboflavin with an almost barrierless pathway (2.7 kcal mol-1). In the singlet channel, conversely, the process is endergonic and requires overcoming a higher activation energy (19.2 kcal mol-1). In the second half of the reaction, the reduction of the studied Pt(IV) complexes by riboflavin hydroquinone occurs via an inner sphere mechanism, displaying free energy barriers smaller than 10 kcal mol-1. Pt reduction by bioreductants such as NADH and ascorbate involve instead less stabilized transition states (22.2-38.3 kcal mol-1), suggesting that riboflavin hydroquinone is an efficient reducing agent for Pt(IV) derivatives in biological settings.

Thymosin ß4 Is an Endogenous Iron Chelator and Molecular Switcher of Ferroptosis.

Thymosin ß4 (Tß4) was extracted forty years agofrom calf thymus. Since then, it has been identified as a G-actin binding protein involved in blood clotting, tissue regeneration, angiogenesis, and anti-inflammatory processes. Tß4 has also been implicated in tumor metastasis and neurodegeneration. However, the precise roles and mechanism(s) of action of Tß4 in these processes remain largely unknown, with the binding of the G-actin protein being insufficient to explain these multi-actions. Here we identify for the first time the important role of Tß4 mechanism in ferroptosis, an iron-dependent form of cell death, which leads to neurodegeneration and somehow protects cancer cells against cell death. Specifically, we demonstrate four iron2+ and iron3+ binding regions along the peptide and show that the presence of Tß4 in cell growing medium inhibits erastin and glutamate-induced ferroptosis in the macrophage cell line. Moreover, Tß4 increases the expression of oxidative stress-related genes, namely BAX, hem oxygenase-1, heat shock protein 70 and thioredoxin reductase 1, which are downregulated during ferroptosis. We state the hypothesis that Tß4 is an endogenous iron chelator and take part in iron homeostasis in the ferroptosis process. We discuss the literature data of parallel involvement of Tß4 and ferroptosis in different human pathologies, mainly cancer and neurodegeneration. Our findings confronted with literature data show that controlled Tß4 release could command on/off switching of ferroptosis and may provide novel therapeutic opportunities in cancer and tissue degeneration pathologies.

Can aluminum, a non-redox metal, alter the thermodynamics of key biological redox processes? The DPPH-QH2 radical scavenging reaction as a test case.

The increased bioavailability of aluminum has led to a concern about its toxicity on living systems. Among the most important toxic effects, it has been proven that aluminum increases oxidative stress in biological systems, a controversial fact, however, due to its non-redox nature. In the present work, we characterize in detail how aluminum can alter redox equilibriums by analyzing its effects on the thermodynamics of the redox scavenging reaction between DPPH., a radical compound often used as a reactive oxygen species model, and hydroquinones, a potent natural antioxidant. For the first time, theoretical and experimental redox potentials within aluminum biochemistry are directly compared. Our results fully agree with experimental reduction and oxidation potentials, unequivocally revealing how aluminum alters the spontaneity of the reaction by stabilizing the reduction of DPPH⋅ to DPPH- and promoting a proton transfer to the diazine moiety, leading to the production of a DPPH-H species. The capability of aluminum to modify redox potentials shown here confirms previous experimental findings on the role of aluminum to interfere with free radical scavenging reactions, affecting the natural redox processes of living organisms.

Markov state models from hierarchical density-based assignment.

Markov state models (MSMs) have become one of the preferred methods for the analysis and interpretation of molecular dynamics (MD) simulations of conformational transitions in biopolymers. While there is great variation in terms of implementation, a well-defined workflow involving multiple steps is often adopted. Typically, molecular coordinates are first subjected to dimensionality reduction and then clustered into small "microstates," which are subsequently lumped into "macrostates" using the information from the slowest eigenmodes. However, the microstate dynamics is often non-Markovian, and long lag times are required to converge the relevant slow dynamics in the MSM. Here, we propose a variation on this typical workflow, taking advantage of hierarchical density-based clustering. When applied to simulation data, this type of clustering separates high population regions of conformational space from others that are rarely visited. In this way, density-based clustering naturally implements assignment of the data based on transitions between metastable states, resulting in a core-set MSM. As a result, the state definition becomes more consistent with the assumption of Markovianity, and the timescales of the slow dynamics of the system are recovered more effectively. We present results of this simplified workflow for a model potential and MD simulations of the alanine dipeptide and the FiP35 WW domain.

New Insights on the Interaction of Phenanthroline Based Ligands and Metal Complexes and Polyoxometalates with Duplex DNA and G-Quadruplexes.

This work provides new insights from our team regarding advances in targeting canonical and non-canonical nucleic acid structures. This modality of medical treatment is used as a form of molecular medicine specifically against the growth of cancer cells. Nevertheless, because of increasing concerns about bacterial antibiotic resistance, this medical strategy is also being explored in this field. Up to three strategies for the use of DNA as target have been studied in our research lines during the last few years: (1) the intercalation of phenanthroline derivatives with duplex DNA; (2) the interaction of metal complexes containing phenanthroline with G-quadruplexes; and (3) the activity of Mo polyoxometalates and other Mo-oxo species as artificial phosphoesterases to catalyze the hydrolysis of phosphoester bonds in DNA. We demonstrate some promising computational results concerning the favorable interaction of these small molecules with DNA that could correspond to cytotoxic effects against tumoral cells and microorganisms. Therefore, our results open the door for the pharmaceutical and medical applications of the compounds we propose.

Mechanistic Insights into Promoted Hydrolysis of Phosphoester Bonds by MoO2Cl2(DMF)2.

A phosphoester bond is a crucial structural block in biological systems, whose occurrence is regulated by phosphatases. Molybdenum compounds have been reported to be active in phosphate ester hydrolysis of model phosphates. Specifically, MoO2Cl2(DMF)2 is active in the hydrolysis of para-nitrophenyl phosphate (pNPP), leading to heteropolyoxometalate structures. We use density functional theory (DFT) to clarify the mechanism by which these species promote the hydrolysis of the phosphoester bond. The present calculations give insight into several key aspects of this reaction: (i) the speciation of this complex prior to interaction with the phosphate (DMF release, Mo-Cl hydrolysis, and pH influence on the speciation), (ii) the competition between phosphate addition and the molybdate nucleation process, (iii) and the mechanisms by which some plausible active species promote this hydrolysis in different conditions. We described thoroughly two different pathways depending on the nucleation possibilities of the molybdenum complex: one mononuclear mechanism, which is preferred in conditions in which very low complex concentrations are used, and another dinuclear mechanism, which is preferred at higher concentrations.

Kojic acid derivatives as double face ligands for metal and phosphate ions.

A family of combined Kojic acid and polyamine derivatives has been synthesized as phosphate anion and metal ion ligands. The stoichiometry, stability and structure of the ion/ligand adducts were determined by 1H NMR spectroscopy, potentiometry, EXAFS and DFT calculations. The presented dual ligands bind effectively both phosphate anions and metal ions and could be used as effective ion receptors in challenging water conditions in the broad pH range. A careful analysis of the heatmaps of the stability constants allows to choose the most appropriate ligand for the ion for qualitative and/or quantitative analysis in water, without analyte pre-treatment. Extremely high-water solubility (>0.6 M) and ion(s)/ligand stability of the adducts in the pH 3-11 are the greatest advantages of the presented here molecules over other known ion sensors. The presented here molecules represent an innovative class of dual metal/anion ligands, with perspective of medical and environmental use.

Enhancing the Photocatalytic Conversion of Pt(IV) Substrates by Flavoprotein Engineering.

Our recent work demonstrates that certain flavoproteins can catalyze the redox activation of Pt(IV) prodrug complexes under light irradiation. Herein, we used site-directed mutagenesis on the mini singlet oxygen generator (mSOG) to modulate the photocatalytic activity of this flavoprotein toward two model Pt(IV) substrates. Among the prepared mutants, Q103V mSOG displayed enhanced catalytic efficiency as a result of its longer triplet excited-state lifetime. This study shows, for the first time, that protein engineering can improve the catalytic capacity of a protein toward metal-containing substrates.

Probing the Catalytically Active Species in POM-Catalysed DNA-Model Hydrolysis*.

Phosphoester hydrolysis is an important chemical step in DNA repair. One archetypal molecular model of phosphoesters is para-nitrophenylphosphate (pNPP). It has been shown previously that the presence of molecular metal oxide [Mo7 O24 ]6- may catalyse the hydrolysis of pNPP through the partial decomposition of polyoxomolybdate framework resulting in a [(PO4 )2 Mo5 O15 ]6- product. Real-time monitoring of the catalytic system using electrospray ionisation mass spectrometry (ESI-MS) provided a glance into the species present in the reaction mixture and identification of potential catalytic candidates. Following up on the obtained spectrometric data, Density Functional Theory (DFT) calculations were carried out to characterise the hypothetical intermediate [Mo5 O15 (pNPP)2 (H2 O)6 ]6- that would be required to form under the hypothesised transformation. Surprisingly, our results point to the dimeric [Mo2 O8 ]4- anion resulting from the decomposition of [Mo7 O24 ]6- as the active catalytic species involved in the hydrolysis of pNPP rather than the originally assumed {Mo5 O15 } species. A similar study was carried out involving the same species but substituting Mo by W. The mechanism involving W species showed a higher barrier and less stable products in agreement with the non-catalytic effect found in experimental results.