Drug metabolism of thiophene containing substrates by cytochrome P450s (CYP450) leads to toxic side effects, for example, nephrotoxicity (suprofen, ticlopidine), hepatotoxicity (tienilic acid), thrombotic thrombocytopenic purpura (clopidogrel), and aplastic anemia (ticlopidine). The origin of toxicity in these cases has been attributed to two different CYP450 mediated metabolic reactions: S-oxidation and epoxidation. In this work, the molecular level details of the bioinorganic chemistry associated with the generation of these competitive reactions are reported. Density functional theory was utilized (i) to explore the molecular mechanism for S-oxidation and epoxidation using the radical cationic center Cpd I [(iron(IV)-oxo-heme porphine system with SH(-) as the axial ligand, to mimic CYP450s] as the model oxidant, (ii) to establish the 3D structures of the reactants, transition states, and products on both the metabolic pathways, and (iii) to examine the potential energy (PE) profile for both the pathways to determine the energetically preferred toxic metabolite formation. The energy barrier required for S-oxidation was observed to be 14.75 kcal/mol as compared to that of the epoxidation reaction (13.23 kcal/mol) on the doublet PE surface of Cpd I. The formation of the epoxide metabolite was found to be highly exothermic (-23.24 kcal/mol), as compared to S-oxidation (-8.08 kcal/mol). Hence, on a relative scale the epoxidation process was observed to be thermodynamically and kinetically more favorable. The energy profiles associated with the reactions of the S-oxide and epoxide toxic metabolites were also explored. This study helps in understanding the CYP450-catalyzed toxic reactions of drugs containing the thiophene ring at the atomic level.

No abstract

Severing proteins are nanomachines from the AAA+ (ATPases associated with various cellular activities) superfamily whose function is to remodel the largest cellular filaments, microtubules. The standard AAA+ machines adopt hexameric ring structures for functional reasons, while being primarily monomeric in the absence of the nucleotide. Both major severing proteins, katanin and spastin, are believed to follow this trend. However, studies proposed that they populate lower-order oligomers in the presence of cofactors, which are functionally relevant. Our simulations show that the preferred oligomeric assembly is dependent on the binding partners and on the type of severing protein. Essential dynamics analysis predicts that the stability of an oligomer is dependent on the strength of the interface between the helical bundle domain (HBD) of a monomer and the convex face of the nucleotide binding domain (NBD) of a neighboring monomer. Hot spots analysis found that the region consisting of the HBD tip and the C-terminal (CT) helix is the only common element between the allosteric networks responding to nucleotide, substrate, and intermonomer binding. Clustering analysis indicates the existence of multiple pathways for the transition between the secondary structure of the HBD tip in monomers and the structure(s) it adopts in oligomers.

This is the author manuscript accepted for publication and has undergone full peer review but has not been through the copyediting, typesetting, pagination and proofreading process, which may lead to differences between this version and the Version of Record. Please cite this article as

Essential cellular processes of microtubule disassembly and protein degradation, which span lengths from tens of μm to nm, are mediated by specialized molecular machines with similar hexameric structure and function. Our molecular simulations at atomistic and coarse-grained scales show that both the microtubule-severing protein spastin and the caseinolytic protease ClpY, accomplish spectacular unfolding of their diverse substrates, a microtubule lattice and dihydrofolate reductase (DHFR), by taking advantage of mechanical anisotropy in these proteins. Unfolding of wild-type DHFR requires disruption of mechanically strong β-sheet interfaces near each terminal, which yields branched pathways associated with unzipping along soft directions and shearing along strong directions. By contrast, unfolding of circular permutant DHFR variants involves single pathways due to softer mechanical interfaces near terminals, but translocation hindrance can arise from mechanical resistance of partially unfolded intermediates stabilized by β-sheets. For spastin, optimal severing action initiated by pulling on a tubulin subunit is achieved through specific orientation of the machine versus the substrate (microtubule lattice). Moreover, changes in the strength of the interactions between spastin and a microtubule filament, which can be driven by the tubulin code, lead to drastically different outcomes for the integrity of the hexameric structure of the machine.