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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Toxic and Physiological Metal Uptake and Release by Human Serum Transferrin.

An atomistic understanding of metal transport in the human body is critical to anticipate the side effects of metal-based therapeutics and holds promise for new drugs and drug delivery designs. Human serum transferrin (hTF) is a central part of the transport processes because of its ubiquitous ferrying of physiological Fe(III) and other transition metals to tightly controlled parts of the body. There is an atomistic mechanism for the uptake process with Fe(III), but not for the release process, or for other metals. This study provides initial insight into these processes for a range of transition metals-Ti(IV), Co(III), Fe(III), Ga(III), Cr(III), Fe(II), Zn(II)-through fully atomistic, extensive quantum mechanical/discrete molecular dynamics sampling and provides, to our knowledge, a new technique we developed to calculate relative binding affinities between metal cations and the protein. It identifies protonation of Tyr188 as a trigger for metal release rather than protonation of Lys206 or Lys296. The study identifies the difficulty of metal release from hTF as potentially related to cytotoxicity. Simulations identify a few critical interactions that stabilize the metal binding site in a flexible, nuanced manner.

Computational Studies on the Binding Preferences of Molybdenum(II) Phenanthroline Complexes with Duplex DNA. The Important Role of the Ancillary Ligands.

The interaction of two isomers, equatorial (Eq) and axial (Ax), of the [Mo(η3-C3H5)Br(CO)2(phen)] metal complex with DNA was studied by using large-scaling density functional theory methods including dispersion for the whole system, represented as a d(AGACGTCT)2 DNA octamer, to gain insight into its experimentally found cytotoxicity. Three different modes of interaction were considered: (1) minor groove (mg) binding, (2) intercalation through the major groove (MG), and (3) the apparently unexpected intercalation via the mg. Computed formation energies, energy decomposition analysis, solvation energies, and noncovalent interaction analysis explain the preference for Eq and Ax isomers of the complex for intercalation via the mg. π-π interactions of the phenanthroline (phen) flat ligand that appear in the intercalation mode and do not exist for the mg binding mode suggest the preference of [Mo(η3-C3H5)Br(CO)2(phen)] for intercalation. On the other hand, the role of the ancillary ligands is crucial for better interaction of the metal complex including phen than when the phen ligand alone is considered because of their additional interactions with base pairs (bps). The role of the ancillary ligands is enhanced when intercalation takes place through the mg because such ligands are able to interact not only with bps but also with the sugar and phosphate backbone, whereas for intercalation through the MG, the interaction of these ligands is only with bps. This feature explains the preference of [Mo(η3-C3H5)Br(CO)2(phen)] for intercalation via the mg in crystal structures. Finally, the solvation penalty is more important for intercalation through the mg than via the MG, which suggests a subtle mechanism involving weak interactions with solvent molecules to explain the selectivity for intercalation in solution to answer the MG versus mg question.

Aluminum and Fenton reaction: how can the reaction be modulated by speciation? A computational study using citrate as a test case.

The pro-oxidant ability of aluminum is behind many of the potential toxic effects of this exogenous element in the human organism. Although the overall process is still far from being understood at the molecular level, the well known ability of aluminum to promote the Fenton reaction is mediated through the formation of stable aluminum-superoxide radical complexes. However, the properties of metal complexes are highly influenced by the speciation of the metal. In this paper, we investigate the effect that speciation could have on the pro-oxidant activity of aluminum. We choose citrate as a test case, because it is the main low-molecular-mass chelator of aluminum in blood serum, forming very stable aluminum-citrate complexes. The influence of citrate in the interaction of aluminum with the superoxide radical is investigated, determining how the formation of aluminum-citrate complexes affects the promotion of the Fenton reaction. The results indicate that citrate increases the stability of the aluminum-superoxide complexes through the formation of ternary compounds, and that the Fenton reaction is even more favorable when aluminum is chelated to citrate. Nevertheless, our results demonstrate that overall, citrate may prevent the pro-oxidant activity of aluminum: on one hand, in an excess of citrate, the formation of 1 : 2 aluminum-citrate complexes is expected. On the other hand, the chelation of iron by citrate makes the reduction of iron thermodynamically unfavorable. In summary, the results suggest that citrate can have both a promotion and protective role, depending on subtle factors, such as initial concentration, non-equilibrium behavior and the exchange rate of ligands in the first shell of the metals.

Photosensitization mechanism of Cu(ii) porphyrins.

This work presents the mechanism of the photoinduced generation of reactive oxygen species (ROS) by paramagnetic copper porphyrins in aqueous solution. Electronic structure calculations within the framework of the (time-dependent) density functional theory, (TD)DFT, reveal the details regarding the development of the atomistic and electronic structures of the copper porphyrin in solution along the set of chemical reactions accessible upon photoactivation. This study identifies the key parameters controlling the feasibility of the various reaction pathways that drive the formation of specific reactive oxygen species, ROS, i.e. superoxide, peroxyl and hydroxyl radicals. An important outcome of our results is the rationalization of how the water solvent molecules play a crucial role in most steps of the overall reaction. The present study is illustrated by focusing on one specific copper porphyrin for which precise experimental data have recently been measured, and can readily be generalized to the whole family of paramagnetic porphyrins. The conclusions of this work shed light on the rational design of metalloporphyrins as photosensitizers for photodynamic therapy.

The Bond Order of C2 from a Strictly N-Representable Natural Orbital Energy Functional Perspective.

The bond order of the ground electronic state of the carbon dimer has been analyzed in the light of natural orbital functional theory calculations carried out with an approximate, albeit strictly N-representable, energy functional. Three distinct solutions have been found from the Euler equations of the minimization of the energy functional with respect to the natural orbitals and their occupation numbers, which expand upon increasing values of the internuclear coordinate. In the close vicinity of the minimum energy region, two of the solutions compete around a discontinuity point. The former, corresponding to the absolute minimum energy, features two valence natural orbitals of each of the following symmetries, σ, σ*, π and π*, and has three bonding interactions and one antibonding interaction, which is very suggestive of a bond order large than two but smaller than three. The latter, features one σ-σ* linked pair of natural orbitals and three degenerate pseudo-bonding like orbitals, paired each with one triply degenerate pseudo-antibonding orbital, which points to a bond order larger than three. When correlation effects, other than Hartree-Fock for example, between the paired natural orbitals are accounted for, this second solution vanishes yielding a smooth continuous dissociation curve. Comparison of the vibrational energies and electron ionization energies, calculated on this curve, with their corresponding experimental marks, lend further support to a bond order for C2 intermediate between acetylene and ethylene.

Cholesterol-Ceramide Interactions in Phospholipid and Sphingolipid Bilayers As Observed by Positron Annihilation Lifetime Spectroscopy and Molecular Dynamics Simulations.

Free volume voids in lipid bilayers can be measured by positron annihilation lifetime spectroscopy (PALS). This technique has been applied, together with differential scanning calorimetry and molecular dynamics (MD) simulations, to study the effects of cholesterol (Chol) and ceramide (Cer) on free volume voids in sphingomyelin (SM) or dipalmitoylphosphatidylcholine (DPPC) bilayers. Binary lipid samples with Chol were studied (DPPC:Chol 60:40, SM:Chol 60:40 mol ratio), and no phase transition was detected in the 20-60 °C range, in agreement with calorimetric data. Chol-driven liquid-ordered phase showed an intermediate free volume void size as compared to gel and fluid phases. For SM and SM:Cer (85:15 mol:mol) model membranes measured in the 20-60 °C range the gel-to-fluid phase transition could be observed with a related increase in free volume, which was more pronounced for the SM:Cer sample. MD simulations suggest a hitherto unsuspected lipid tilting in SM:Cer bilayers but not in pure SM. Ternary samples of DPPC:Cer:Chol (54:23:23) and SM:Cer:Chol (54:23:23) were measured, and a clear pattern of free volume increase was observed in the 20-60 °C because of the gel-to-fluid transition. Interestingly, MD simulations showed a tendency of Cer to change its distribution along the membrane to make room for Chol in ternary mixtures. The results suggest that the gel phase formed in these ternary mixtures is stabilized by Chol-Cer interactions.