An omics perspective to the molecular mechanisms of anticancer metallo-drugs in the computational microscope era

Rationalizing the structural basis of organic-platinum hybrid complexes binding towards quadruplex-duplex hybrids through all-atom simulations

Guanine-rich sequences containing complementary base pairs can fold into non-canonical quadruplex-duplex hybrid (QDH) conformations. These structures possess unique structural features, leading to the presence of a peculiar binding pocket that can be distinguished from a canonical double helix or a G-quadruplex (G4) structure. Recently, two organic-metal hybrid platinum complexes, able to selectively and strongly recognize a particular type of QDH with a lateral duplex stem-loop, were reported in the literature. However, solution structures are not available for all the investigated compounds, leaving unanswered questions on the structural traits underlying the different binding affinity of these complexes. In this work, we address this gap using all-atom simulations to unravel the key features driving the high selectivity of these organic‑platinum hybrid complexes at an atomistic level. In particular, their binding affinity depends on a delicate balance between the extended π-π stacking interactions performed in the G4-duplex binding pocket and the capacity to form stable hydrogen bonds with the surrounding nucleobases. Thus, our findings provide essential insights to guide the rational design of novel compounds that selectively target QDH structures.

Unraveling the Molecular Basis for G-Quadruplex-Binders to ALS/FTD-Associated G4C2 Repeats of the C9orf72 Gene

The most recurrent familial cause of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) is the presence of an abnormal number of intronic GGGGCC (G4C2) repetitions in the C9orf72 gene, which has been proposed to drive ALS/FTD pathogenesis. Recently, it has been shown that such G4C2 repetitions can fold into G-quadruplex (G4) secondary structures. These G4s have been selectively stabilized by small-molecule binders, furnishing proof-of-principle that targeting these non-canonical nucleic acid sequences represents a novel and effective therapeutic strategy to tackle neurodegenerative disorders. However, precise information on the mechanism of action of these compounds is still lacking. Here, by performing in silico investigations, we unraveled the molecular basis for the selectivity of a series of known structurally related C9orf72 G4-binders. Moreover, we investigated the binding properties of a strong and selective metal-based G4 stabilizer, the AuI bis-N-heterocyclic carbene (NHC) complex - Au(TMX)2 - showing that it moderately stabilizes G4C2 G4 RNA by Förster resonance energy transfer (FRET) DNA melting assays. Using metadynamics (metaD) simulations, the Au(TMX)2 binding mode and the associated free-energy landscape were also evaluated. This information paves the way for developing improved compounds to tackle ALS/FTD neurodegenerative disorders.

Frontiers of metal-coordinating drug design

Introduction: The occurrence of metal ions in biomolecules is required to exert vital cellular functions. Metal-containing biomolecules can be modulated by small-molecule inhibitors targeting their metal-moiety. As well, the discovery of cisplatin ushered the rational discovery of metal-containing-drugs. The use of both drug types exploiting metal-ligand interactions is well established to treat distinct pathologies. Therefore, characterizing and leveraging metal-coordinating drugs is a pivotal, yet challenging, part of medicinal chemistry.Area covered: Atomic-level simulations are increasingly employed to overcome the challenges met by traditional drug-discovery approaches and to complement wet-lab experiments in elucidating the mechanisms of drugs' action. Multiscale simulations, allow deciphering the mechanism of metal-binding inhibitors and metallo-containing-drugs, enabling a reliable description of metal-complexes in their biological environment. In this compendium, the authors review selected applications exploiting the metal-ligand interactions by focusing on understanding the mechanism and design of (i) inhibitors targeting iron and zinc-enzymes, and (ii) ruthenium and gold-based anticancer agents targeting the nucleosome and aquaporin protein, respectively.Expert opinion: The showcased applications exemplify the current role and the potential of atomic-level simulations and reveal how their synergic use with experiments can contribute to uncover fundamental mechanistic facets and exploit metal-ligand interactions in medicinal chemistry.

Targeting Orthosteric and Allosteric Pockets of Aromatase via Dual-Mode Novel Azole Inhibitors

Breast cancer (BC) is the most diffused cancer type in women and the second leading cause of death among the female population. Effective strategies to fight estrogen responsive (ER+) BC, which represents 70% of all BC cases, rely on estrogen deprivation, via the inhibition of the aromatase enzyme, or the modulation of its cognate estrogen receptor. Current clinical therapies significantly increased patient survival time. Nevertheless, the onset of resistance in metastatic BC patients undergoing prolonged treatments is becoming a current clinical challenge, urgently demanding to devise innovative strategies. In this context, here we designed, synthesized, and performed in vitro inhibitory tests on the aromatase enzyme and distinct ER+/ER- BC cell line types of novel azole bridged xanthones. These compounds are active in the low μM range and behave as dual-mode inhibitors, targeting both the orthosteric and the allosteric sites of the enzyme placed along one access channel. Classical and quantum-classical molecular dynamics simulations of the new compounds, as compared with selected steroidal and nonsteroidal inhibitors, provide a rationale to the observed inhibitory potency and supply the guidelines to boost the activity of inhibitors able to exploit coordination to iron and occupation of the access channel to modulate estrogen production.

Dynamics Insights Into the Gain of Flexibility by Helix-12 in ESR1 as a Mechanism of Resistance to Drugs in Breast Cancer Cell Lines

Incidents of breast cancer (BC) are on the rise on a daily basis and have proven to be the most prevelant cause of death for women in both developed and developing countries. Among total BC cases diagnosed after menopause, 70% of cases are Estrogen Receptor (ER) positive (ER-positive or ER+). Mutations in the LBD (ligand-binding domain) of the ER have recently been reported to be the major cause of resistance to potent antagonists. In this study, the experimentally reported mutations K303R, E380Q, V392I, S463P, V524E, P535H, P536H, Y537C, Y537N, Y537S, and D538G were analyzed, and the most significant mutations were shortlisted based on multiple analyses. Initial analyses, such as mCSM stability, occluded depth analysis, mCSM-binding affinity, and FoldX energy changes shortlisted only six mutations as being highly resistant. Finally, simulations of force field-based molecular dynamics (MD on wild type (WT) ERα) on six mERα variants (E380Q, S463P, Y537S, Y537C, Y537N, and D538G) were carried out to justify mechanism of the resistance. It was observed that these mutations increased the flexibility of the H12. A bonding analysis suggested that previously reported important residue His524 lost bonding upon mutation. Other parameters, such as PCA (principal component analysis), DCCM (dynamics cross-correlation), and FEL (free energy landscape), verified that the shortlisted mutations affect the H12 helix, which opens up the co-activator binding conformation. These results provide deep insight into the mechanism of relative resistance posed to fulvestrant due to mutations in breast cancer. This study will facilitate further understanding of the important aspects of designing specific and more effective drugs.

Unraveling the Impact of Cysteine-to-Serine Mutations on the Structural and Functional Properties of Cu(I)-Binding Proteins

Appropriate maintenance of Cu(I) homeostasis is an essential requirement for proper cell function because its misregulation induces the onset of major human diseases and mortality. For this reason, several research efforts have been devoted to dissecting the inner working mechanism of Cu(I)-binding proteins and transporters. A commonly adopted strategy relies on mutations of cysteine residues, for which Cu(I) has an exquisite complementarity, to serines. Nevertheless, in spite of the similarity between these two amino acids, the structural and functional impact of serine mutations on Cu(I)-binding biomolecules remains unclear. Here, we applied various biochemical and biophysical methods, together with all-atom simulations, to investigate the effect of these mutations on the stability, structure, and aggregation propensity of Cu(I)-binding proteins, as well as their interaction with specific partner proteins. Among Cu(I)-binding biomolecules, we focused on the eukaryotic Atox1-ATP7B system, and the prokaryotic CueR metalloregulator. Our results reveal that proteins containing cysteine-to-serine mutations can still bind Cu(I) ions; however, this alters their stability and aggregation propensity. These results contribute to deciphering the critical biological principles underlying the regulatory mechanism of the in-cell Cu(I) concentration, and provide a basis for interpreting future studies that will take advantage of cysteine-to-serine mutations in Cu(I)-binding systems.

Post-Translational Regulation of CYP450s Metabolism As Revealed by All-Atoms Simulations of the Aromatase Enzyme

Phosphorylation by kinases enzymes is a widespread regulatory mechanism able of rapidly altering the function of target proteins. Among these are cytochrome P450s (CYP450), a superfamily of enzymes performing the oxidation of endogenous and exogenous substrates thanks to the electron supply of a redox partner. In spite of its pivotal role, the molecular mechanism by which phosphorylation modulates CYP450s metabolism remains elusive. Here by performing microsecond-long all-atom molecular dynamics simulations, we disclose how phosphorylation regulates estrogen biosynthesis, catalyzed by the Human Aromatase (HA) enzyme. Namely, we unprecedentedly propose that HA phosphorylation at Y361 markedly stabilizes its adduct with the flavin mononucleotide domain of CYP450s reductase (CPR), the redox partner of microsomal CYP450s, and a variety of other proteins. With CPR present at physiological conditions in a limiting ratio with respect to its multiple oxidative partners, the enhanced stability of the CPR/HA adduct may favor HA in the competition with the other proteins requiring CPR's electron supply, ultimately facilitating the electron transfer and estrogen biosynthesis. As a result, our work elucidates at atomic-level the post-translational regulation of CYP450s catalysis. Given the potential for rational clinical management of diseases associated with steroid metabolism disorders, unraveling this mechanism is of utmost importance, and raises the intriguing perspective of exploiting this knowledge to devise novel therapies.

The Catalytic Mechanism of Steroidogenic Cytochromes P450 from All-Atom Simulations: Entwinement with Membrane Environment, Redox Partners, and Post-Transcriptional Regulation

Cytochromes P450 (CYP450s) promote the biosynthesis of steroid hormones with major impact on the onset of diseases such as breast and prostate cancers. By merging distinct functions into the same catalytic scaffold, steroidogenic CYP450s enhance complex chemical transformations with extreme efficiency and selectivity. Mammalian CYP450s and their redox partners are membrane-anchored proteins, dynamically associating to form functional machineries. Mounting evidence signifies that environmental factors are strictly intertwined with CYP450s catalysis. Atomic-level simulations have the potential to provide insights into the catalytic mechanism of steroidogenic CYP450s and on its regulation by environmental factors, furnishing information often inaccessible to experimental means. In this review, after an introduction of computational methods commonly employed to tackle these systems, we report the current knowledge on three steroidogenic CYP450s—CYP11A1, CYP17A1, and CYP19A1—endowed with multiple catalytic functions and critically involved in cancer onset. In particular, besides discussing their catalytic mechanisms, we highlight how the membrane environment contributes to (i) regulate ligand channeling through these enzymes, (ii) modulate their interactions with specific protein partners, (iii) mediate post-transcriptional regulation induced by phosphorylation. The results presented set the basis for developing novel therapeutic strategies aimed at fighting diseases originating from steroid metabolism dysfunction.

A Computational Assay of Estrogen Receptor α Antagonists Reveals the Key Common Structural Traits of Drugs Effectively Fighting Refractory Breast Cancers

Somatic mutations of the Estrogen Receptor α (ERα) occur with an up to 40% incidence in ER sensitive breast cancer (BC) patients undergoing prolonged endocrine treatments. These polymorphisms are implicated in acquired resistance, disease relapse, and increased mortality rates, hence representing a current major clinical challenge. Here, multi-microseconds (12.5 µs) molecular dynamics simulations revealed that recurrent ERα polymorphisms (i. e. L536Q, Y537S, Y537N, D538G) (mERα) are constitutively active in their apo form and that they prompt the selection of an agonist (active)-like conformation even upon antagonists binding. Interestingly, our simulations rationalize, for the first time, the efficacy profile of (pre)clinically used Selective Estrogen Receptor Modulators/Downregulators (SERMs/SERDs) against these variants, enlightening, at atomistic level of detail, the key common structural traits needed by drugs able to effectively fight refractory BC types. This knowledge represents a key advancement for mechanism-based therapeutics targeting resistant ERα isoforms, potentially allowing the community to move a step closer to ‘precision medicine’ calibrated on patients’ genetic profiles and disease progression.