The estimation of the redox potentials of biologically relevant systems by means of theoretical-computational approaches still represents a challenge. In fact, the size of these systems typically does not allow a full quantum-mechanical treatment needed to describe electron loss/gain in such a complex environment, where the redox process takes place. Therefore, a number of different theoretical strategies have been developed so far to make the calculation of the redox free energy feasible with current computational resources. In this review, we provide a survey of such theoretical-computational approaches used in this context, highlighting their physical principles and discussing their advantages and limitations. Several examples of these approaches applied to the estimation of the redox potentials of both proteins and nucleic acids are described and critically discussed. Finally, general considerations on the most promising strategies are reported.

First ML/MM-based nonadiabatic dynamics study modelling interactions of environment and electronic states of ML-region with electrostatic embedding.

Substrate binding to the cytochrome P450 OleP is coupled to a large open-to-closed transition that remodels the active site, minimizing its exposure to the external solvent. When the aglycone substrate binds, a small empty cavity is formed between the I and G helices, the BC loop, and the substrate itself, where solvent molecules accumulate mediating substrate-enzyme interactions. Herein, we analyzed the role of this cavity in substrate binding to OleP by producing three mutants (E89Y, G92W, and S240Y) to decrease its volume. The crystal structures of the OleP mutants in the closed state bound to the aglycone 6DEB showed that G92W and S240Y occupied the cavity, providing additional contact points with the substrate. Conversely, mutation E89Y induces a flipped-out conformation of this amino acid side chain, that points towards the bulk, increasing the empty volume. Equilibrium titrations and molecular dynamic simulations indicate that the presence of a bulky residue within the cavity impacts the binding properties of the enzyme, perturbing the conformational space explored by the complexes. Our data highlight the relevance of this region in OleP substrate binding and suggest that it represents a key substrate-protein contact site to consider in the perspective of redirecting its activity towards alternative compounds.

Anthracycline doxorubicin hydrochloride (DX) is a positively charged fluorescent drug, which in water self-associates into non-fluorescent antiparallel dimers upon increasing concentration and/or ionic strength. The positive charge of DX allows for complexation with negatively charged polymers and drug carriers. The fluorescence of DX following complexation with polyanion polystyrene sulfonate (PSS) is studied here. The fluorescence emission of DX decreases in the presence of PSS, being almost completely quenched when the ratio (R) of PSS monomers-to-DX molecules is larger than 10. Increasing R values over 30 results in a progressive recovery of fluorescence. The circular dichroism of PSS–DX complexes shows inverted characteristic bands of DX dimers suggesting the presence of parallel dimers at a concentration of DX below dimerization in water. Molecular dynamics studies corroborate a preferential orientation of DX into parallel dimers when interacting with PSS and show that DX molecules interact with a binding pocket of PSS monomers rather than with one single monomer. Increasing the ionic strength results in a recovery of fluorescence without an apparent release of DX from the PSS–DX complex as shown by DOSY NMR. PSS acts as a template for concentrating DX, triggering dimerisation and orienting DX molecules with their charged groups facing the negatively charged PSS monomers.