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.

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.

The experimental absorption measurements and an advanced treatment of quantum-mechanical and molecular dynamics data are here combined to provide a complete picture of the absorption behavior doxorubicin hydrochloride (DX) in different environments.

TipA, a MerR family transcription factor from Streptomyces lividans , promotes antibiotic resistance by sequestering broad-spectrum thiopeptide-based antibiotics, thus counteracting their inhibitory effect on ribosomes. TipAS, a minimal binding motif which is expressed as an isoform of TipA, harbors a partially disordered N-terminal subdomain that folds upon binding multiple antibiotics. The extent and nature of the underlying molecular heterogeneity in TipAS that shapes its promiscuous folding–function landscape is an open question and is critical for understanding antibiotic-sequestration mechanisms. Here, combining equilibrium and time-resolved experiments, statistical modeling, and simulations, we show that the TipAS native ensemble exhibits a pre-equilibrium between binding-incompetent and binding-competent substates, with the fully folded state appearing only as an excited state under physiological conditions. The binding-competent state characterized by a partially structured N-terminal subdomain loses structure progressively in the physiological range of temperatures, swells on temperature increase, and displays slow conformational exchange across multiple conformations. Binding to the bactericidal antibiotic thiostrepton follows a combination of induced-fit and conformational-selection-like mechanisms, via partial binding and concomitant stabilization of the binding-competent substate. These ensemble features are evolutionarily conserved across orthologs from select bacteria that infect humans, underscoring the functional role of partial disorder in the native ensemble of antibiotic-sequestering proteins belonging to the MerR family.

(1) Background: the theoretical modelling of reactions occurring in liquid phase is a research line of primary importance both in theoretical–computational chemistry and in the context of organic and biological chemistry. Here we present the modelling of the kinetics of the hydroxide-promoted hydrolysis of phosphoric diesters. (2) Method: the theoretical–computational procedure involves a hybrid quantum/classical approach based on the perturbed matrix method (PMM) in conjunction with molecular mechanics. (3) Results: the presented study reproduces the experimental data both in the rate constants and in the mechanistic aspects (C–O bond vs. O–P bond reactivity). The study suggests that the basic hydrolysis of phosphodiesters occurs through a concerted $A_N{D}_N$ mechanism, with no formation of penta-coordinated species as reaction intermediates. (4) Conclusions: the presented approach, despite the approximations, is potentially applicable to a large number of bimolecular transformations in solution and therefore leads the way to a fast and general method to predict the rate constants and reactivities/selectivities in complex environments.