The role of electrostatics in TrxR electron transfer mechanism: A computational approach

Fighting Multidrug Resistance with Ruthenium-Cyclopentadienyl Compounds: Unveiling the Mechanism of P-gp Inhibition

The search for more effective and selective drugs to overcome cancer multidrug resistance is urgent. As such, a new series of ruthenium-cyclopentadienyl ("RuCp") compounds with the general formula [Ru(η5-C5H4R)(4,4'-R'-2,2'-bipy)(PPh3)] were prepared and fully characterized. All compounds were evaluated toward non-small cell lung cancer cells with different degrees of cisplatin sensitivity (A549, NCI-H2228, Calu-3, and NCI-H1975), showing better cytotoxicity than the first-line chemotherapeutic drug cisplatin. Compounds 2 and 3 (R' = -OCH3; R = CHO (2) or CH2OH (3)) further inhibited the activity of P-gp and MRP1 efflux pumps by impairing their catalytic activity. Molecular docking calculations identified the R-site P-gp pocket as the preferred one, which was further validated using site-directed mutagenesis experiments in P-gp. Altogether, our results unveil the first direct evidence of the interaction between P-gp and "RuCp" compounds in the modulation of P-gp activity and establish them as valuable candidates to circumvent cancer MDR.

Analysis of Pleurotin binding to human thioredoxin reductase using docking and molecular dynamics simulation

Thioredoxin reductase (TrxR) has been considered a potential target for cancer chemotherapy. It acts by controlling the redox homeostasis of human cells and, therefore, interfering in its function may trigger apoptosis, which is a crucial tumor suppression mechanism. Despite the great effort in the search for TrxR inhibitors, none was approved for human therapy. In the present study a virtual screening for natural organic compounds is discussed for a set of 72 compounds with known IC-50 for TrxR inhibition. The results suggest the Pleurotin, a naphthoquinone obtained from Hohenbuehelia grisea fungus, as a potential TrxR inhibitor, which acts by binding to the active site of the enzyme, between the N- and C-terminal domains. The presence of the ligand blocks the approximation of the C-terminal arm to the N-terminal, which is an essential step of the enzyme function. Besides, the two equivalent binding sites of TrxR were explored, by docking two ligands simultaneously. The results indicate that both sites have an allosteric correlation and, the presence of the ligand in one site may interfere, or even prevent, the binding of the second ligand at the other site. All these findings are quantitatively discussed based on the analysis of long molecular dynamics trajectories, which provides a full description of the ligand-receptor binding modes, average binding energies and conformational changes.Communicated by Ramaswamy H. Sarma.

Novel US-CpHMD Protocol to Study the Protonation-Dependent Mechanism of the ATP/ADP Carrier

We have designed a protocol combining constant-pH molecular dynamics (CpHMD) simulations with an umbrella sampling (US) scheme (US-CpHMD) to study the mechanism of ADP/ATP transport (import and export) by their inner mitochondrial membrane carrier protein [ADP/ATP carrier (AAC)]. The US scheme helped overcome the limitations of sampling the slow kinetics involved in these substrates' transport, while CpHMD simulations provided an unprecedented realism by correctly capturing the associated protonation changes. The import of anionic substrates along the mitochondrial membrane has a strong energetic disadvantage due to a smaller substrate concentration and an unfavorable membrane potential. These limitations may have created an evolutionary pressure on AAC to develop specific features benefiting the import of ADP. In our work, the potential of mean force profiles showed a clear selectivity in the import of ADP compared to ATP, while in the export, no selectivity was observed. We also observed that AAC sequestered both substrates at longer distances in the import compared to the export process. Furthermore, only in the import process do we observe transient protonation of both substrates when going through the AAC cavity, which is an important advantage to counteract the unfavorable mitochondrial membrane potential. Finally, we observed a substrate-induced disruption of the matrix salt-bridge network, which can promote the conformational transition (from the C- to M-state) required to complete the import process. This work unraveled several important structural features where the complex electrostatic interactions were pivotal to interpreting the protein function and illustrated the potential of applying the US-CpHMD protocol to other transport processes involving membrane proteins.

pK a Calculations in Membrane Proteins from Molecular Dynamics Simulations

The conformational changes of membrane proteins are crucial to their function and usually lead to fluctuations in the electrostatic environment of the protein surface. A very effective way to quantify these changes is by calculating the pK _a values of the protein's titratable residues, which can be regarded as electrostatic probes. To achieve this, we need to take advantage of the fast and reliable pK _a calculators developed for globular proteins and adapt them to include the explicit effects of membranes. Here, we provide a detailed linear response approximation protocol that uses our own software (PypKa) to calculate reliable pK _a values from short MD simulations of membrane proteins.

Computational Analysis of the Interactions between the S100B Extracellular Chaperone and Its Amyloid β Peptide Client

S100B is an astrocytic extracellular Ca²⁺-binding protein implicated in Alzheimer’s disease, whose role as a holdase-type chaperone delaying Aβ₄₂ aggregation and toxicity was recently uncovered. Here, we employ computational biology approaches to dissect the structural details and dynamics of the interaction between S100B and Aβ₄₂. Driven by previous structural data, we used the Aβ_25–35 segment, which recapitulates key aspects of S100B activity, as a starting guide for the analysis. We used Haddock to establish a preferred binding mode, which was studied with the full length Aβ using long (1 μs) molecular dynamics (MD) simulations to investigate the structural dynamics and obtain representative interaction complexes. From the analysis, Aβ-Lys28 emerged as a key candidate for stabilizing interactions with the S100B binding cleft, in particular involving a triad composed of Met79, Thr82 and Glu86. Binding constant calculations concluded that coulombic interactions, presumably implicating the Lys28(Aβ)/Glu86(S100B) pair, are very relevant for the holdase-type chaperone activity. To confirm this experimentally, we examined the inhibitory effect of S100B over Aβ aggregation at high ionic strength. In agreement with the computational predictions, we observed that electrostatic perturbation of the Aβ-S100B interaction decreases anti-aggregation activity. Altogether, these findings unveil features relevant in the definition of selectivity of the S100B chaperone, with implications in Alzheimer’s disease.

Role of the Enzymatic Environment in the Reactivity of the AuIII-C^N^C Anticancer Complexes

The action mechanism of anticancer gold(III) complexes is a multi-step process and depends on their redox stability. First, the gold(III) complex undergoes a ligand exchange reaction in the presence of cellular thiols, such as those available in the active site of the enzyme TrxR, and then, the Au^III → Au^I reduction occurs. Most experimental and theoretical studies describe these processes under chemical conditions without considering the enzyme structure effect. In the present study, molecular models are proposed for the [Au^III(C^N^C)(SHCys-R)]⁺ adduct, with the [Au^III(C^N^C)]⁺ moiety bonded to the Cys498 residue in the C-terminal arm of the TrxR. This one represents the product of the first ligand exchange reaction. Overall, our results suggest that the exchange of the auxiliary ligand (for instance, Cl^- to S-R) plays a primary role in increasing the reduction potential, with the enzyme structure having a small effect. The parent compound [Au^III(C^N^C)Cl] has E° = -1.20 V, which enlarges to -0.72 V for [Au^III(C^N^C)CH₃SH]⁺ and to -0.65 V for the largest model studied, Au-trx. In addition to the effect of the enzyme structure on the redox stability, we also analyze the Au transfer to the enzyme using a small peptide model (a tetramer). This reaction is dependent on the Cys497 protonation state. Thermodynamics and kinetic analysis suggests that the C^N^C ligand substitution by Cys497 is an exergonic process, with an energy barrier estimated at 20.2 kcal mol^-1. The complete transfer of the Au ion to the enzyme's active site would lead to a total loss of enzyme activity, generating oxidative damage and, consequently, cancer cell death.

Targeting the Thioredoxin System as a Strategy for Cancer Therapy

Thioredoxin reductase (TrxR) participates in the regulation of redox reactions in organisms. It works mainly via its substrate molecule, thioredoxin, to maintain the redox balance and regulate signal transduction, which controls cell proliferation, differentiation, death, and other important physiological processes. In recent years, increasing evidence has shown that the overactivation of TrxR is related to the development of tumors. The exploration of TrxR-targeted antitumor drugs has attracted wide attention and is expected to provide new therapies for cancer treatment. In this perspective, we highlight the specific relationship between TrxR and apoptotic signaling pathways. The cytoplasm and mitochondria both contain TrxR, resulting in the activation of apoptosis. TrxR activity influences reactive oxygen species (ROS) and further regulates the inflammatory signaling pathway. In addition, we discuss representative TrxR inhibitors with anticancer activity and analyze the challenges in developing TrxR inhibitors as anticancer drugs.