Predicting Conformational Ensembles of Intrinsically Disordered Proteins: From Molecular Dynamics to Machine Learning.

Intrinsically disordered proteins and regions (IDP/IDRs) are ubiquitous across all domains of life. Characterized by a lack of a stable tertiary structure, IDP/IDRs populate a diverse set of transiently formed structural states that can promiscuously adapt upon binding with specific interaction partners and/or certain alterations in environmental conditions. This malleability is foundational for their role as tunable interaction hubs in core cellular processes such as signaling, transcription, and translation. Tracing the conformational ensemble of an IDP/IDR and its perturbation in response to regulatory cues is thus paramount for illuminating its function. However, the conformational heterogeneity of IDP/IDRs poses several challenges. Here, we review experimental and computational methods devised to disentangle the conformational landscape of IDP/IDRs, highlighting recent computational advances that permit proteome-wide scans of IDP/IDRs conformations. We briefly evaluate selected computational methods using the disordered N-terminal of the human copper transporter 1 as a test case and outline further challenges in IDP/IDRs ensemble prediction.

Molecular dynamics simulations for the structure-based drug design: targeting small-GTPases proteins.

INTRODUCTION: Molecular Dynamics (MD) simulations can support mechanism-based drug design. Indeed, MD simulations by capturing biomolecule motions at finite temperatures can reveal hidden binding sites, accurately predict drug-binding poses, and estimate the thermodynamics and kinetics, crucial information for drug discovery campaigns. Small-Guanosine Triphosphate Phosphohydrolases (GTPases) regulate a cascade of signaling events, that affect most cellular processes. Their deregulation is linked to several diseases, making them appealing drug targets. The broad roles of small-GTPases in cellular processes and the recent approval of a covalent KRas inhibitor as an anticancer agent renewed the interest in targeting small-GTPase with small molecules. AREA COVERED: This review emphasizes the role of MD simulations in elucidating small-GTPase mechanisms, assessing the impact of cancer-related variants, and discovering novel inhibitors. EXPERT OPINION: The application of MD simulations to small-GTPases exemplifies the role of MD simulations in the structure-based drug design process for challenging biomolecular targets. Furthermore, AI and machine learning-enhanced MD simulations, coupled with the upcoming power of quantum computing, are promising instruments to target elusive small-GTPases mutations and splice variants. This powerful synergy will aid in developing innovative therapeutic strategies associated to small-GTPases deregulation, which could potentially be used for personalized therapies and in a tissue-agnostic manner to treat tumors with mutations in small-GTPases.

Third Metal Ion Dictates the Catalytic Activity of the Two-Metal-Ion Pre-Ribosomal RNA-Processing Machinery.

Nucleic acid processing enzymes use a two-Mg2+-ion motif to promote the formation and cleavage of phosphodiester bonds. Yet, recent evidence demonstrates the presence of spatially conserved second-shell cations surrounding the catalytic architecture of proteinaceous and RNA-dependent enzymes. The RNase mitochondrial RNA processing (MRP) complex, which cleaves the ribosomal RNA (rRNA) precursor at the A3 cleavage site to yield mature 5'-end of 5.8S rRNA, hosts in the catalytic core one atypically-located Mg2+ ion, in addition to the ions forming the canonical catalytic motif. Here, we employ biased quantum classical molecular dynamics simulations of RNase MRP to discover that the third Mg2+ ion inhibits the catalytic process. Instead, its displacement in favour of a second-shell monovalent K+ ion propels phosphodiester bond cleavage by enabling the formation of a specific hydrogen bonding network that mediates the essential proton transfer step. This study points to a direct involvement of a transient K+ ion in the catalytic cleavage of the phosphodiester bond and implicates cation trafficking as a general mechanism in nucleic acid processing enzymes and ribozymes.

Neutral rhodol-based dyes expressing localization in mitochondria.

Neutral rhodol-based red emitters are shown to efficiently localize in mitochondria, as demonstrated by confocal microscopy and co-localization studies. A simple model is proposed to explain the localization mechanism of neutral molecules. The model takes into account the strong coupling between the molecular dipole moment and the electric field of the inner mitochondrial membrane.

Assessing the mechanism of fast-cycling cancer-associated mutations of Rac1 small Rho GTPase.

Rho-GTPases proteins function as molecular switches alternating from an active to an inactive state upon Guanosine triphosphate (GTP) binding and hydrolysis to Guanosine diphosphate (GDP). Among them, Rac subfamily regulates cell dynamics, being overexpressed in distinct cancer types. Notably, these proteins are object of frequent cancer-associated mutations at Pro29 (P29S, P29L, and P29Q). To assess the impact of these mutations on Rac1 structure and function, we performed extensive all-atom molecular dynamics simulations on wild-type (wt) and oncogenic isoforms of this protein in GDP- and GTP-bound states. Our results unprecedentedly elucidate that P29Q/S-induced structural and dynamical perturbations of Rac1 core domain weaken the binding of the catalytic site Mg2+ ion, and reduce the GDP residence time within protein, enhancing the GDP/GTP exchange rate and Rac1 activity. This broadens our knowledge of the role of cancer-associated mutations on small GTPases mechanism supplying valuable information for future drug discovery efforts targeting specific Rac1 isoforms.

Tracing Allostery in the Spliceosome Ski2-like RNA Helicase Brr2.

RNA ATPases/helicases remodel substrate RNA-protein complexes in distinct ways. The different RNA ATPases/helicases, taking part in the spliceosome complex, reshape the RNA/RNA-protein contacts to enable premature-mRNA splicing. Among them, the bad response to refrigeration 2 (Brr2) helicase promotes U4/U6 small nuclear (sn)RNA unwinding via ATP-driven translocation of the U4 snRNA strand, thus playing a pivotal role during the activation, catalytic, and disassembly phases of splicing. The plastic Brr2 architecture consists of an enzymatically active N-terminal cassette (N-cassette) and a structurally similar but inactive C-terminal cassette (C-cassette). The C-cassette, along with other allosteric effectors and regulators, tightly and timely controls Brr2's function via an elusive mechanism. Here, microsecond-long molecular dynamics simulations, dynamical network theory, and community network analysis are combined to elucidate how allosteric effectors/regulators modulate the Brr2 function. We unexpectedly reveal that U4 snRNA itself acts as an allosteric regulator, amplifying the cross-talk of distal Brr2 domains and triggering a conformational reorganization of the protein. Our findings offer fundamental understanding into Brr2's mechanism of action and broaden our knowledge on the sophisticated regulatory mechanisms by which spliceosome ATPases/helicases control gene expression. This includes their allosteric regulation exerted by client RNA strands, a mechanism that may be broadly applicable to other RNA-dependent ATPases/helicases.

Frontiers and Challenges of Computing ncRNAs Biogenesis, Function and Modulation.

Non-coding RNAs (ncRNAs), generated from nonprotein coding DNA sequences, constitute 98-99% of the human genome. Non-coding RNAs encompass diverse functional classes, including microRNAs, small interfering RNAs, PIWI-interacting RNAs, small nuclear RNAs, small nucleolar RNAs, and long non-coding RNAs. With critical involvement in gene expression and regulation across various biological and physiopathological contexts, such as neuronal disorders, immune responses, cardiovascular diseases, and cancer, non-coding RNAs are emerging as disease biomarkers and therapeutic targets. In this review, after providing an overview of non-coding RNAs' role in cell homeostasis, we illustrate the potential and the challenges of state-of-the-art computational methods exploited to study non-coding RNAs biogenesis, function, and modulation. This can be done by directly targeting them with small molecules or by altering their expression by targeting the cellular engines underlying their biosynthesis. Drawing from applications, also taken from our work, we showcase the significance and role of computer simulations in uncovering fundamental facets of ncRNA mechanisms and modulation. This information may set the basis to advance gene modulation tools and therapeutic strategies to address unmet medical needs.

Monovalent metal ion binding promotes the first transesterification reaction in the spliceosome.

Cleavage and formation of phosphodiester bonds in nucleic acids is accomplished by large cellular machineries composed of both protein and RNA. Long thought to rely on a two-metal-ion mechanism for catalysis, structure comparisons revealed many contain highly spatially conserved second-shell monovalent cations, whose precise function remains elusive. A recent high-resolution structure of the spliceosome, essential for pre-mRNA splicing in eukaryotes, revealed a potassium ion in the active site. Here, we employ biased quantum mechanics/ molecular mechanics molecular dynamics to elucidate the function of this monovalent ion in splicing. We discover that the K+ ion regulates the kinetics and thermodynamics of the first splicing step by rigidifying the active site and stabilizing the substrate in the pre- and post-catalytic state via formation of key hydrogen bonds. Our work supports a direct role for the K+ ion during catalysis and provides a mechanistic hypothesis likely shared by other nucleic acid processing enzymes.

Assessing the Binding Mode of a Splicing Modulator Stimulating Pre-mRNA Binding to the Plastic U2AF2 Splicing Factor.

RNA recognition motifs (RRMs) play a pivotal role in RNA metabolism and the regulation of gene expression. Owing to their plasticity and fuzziness, targeting RRM/RNA interfaces with small molecules is a daunting challenge for drug discovery campaigns. The U2AF2 splicing factor, which recognizes the polypyrimidine (polyPy) sequence of premature messenger (pre-m)RNA, exhibits a dynamic architecture consisting of two RRMs joined by a disordered linker. An inhibitor, NSC-194308, was shown to enhance the binding of pre-mRNA to U2AF2, selectively triggering cell death in leukemia cell lines containing spliceosome mutations. The NSC-194308 binding mode remains elusive; yet, unraveling its knowledge may offer intriguing insights for effectively targeting U2AF2 and other flexible protein/protein/RNA interfaces with small molecules. To infer plausible NSC-194308 binding poses to U2AF2, here, we applied and benchmarked the performance of static and dynamic docking approaches, elucidating the molecular basis of the NSC-194308-induced pre-mRNA stabilization on U2AF2. We demonstrate that introducing dynamic effects is mandatory to assess the binding mode of the inhibitors when they target plastic and modular architectures, such as those formed by interacting RRMs. The latter are widespread across RNA binding proteins; therefore, this mechanism may be broadly applicable to discover new therapeutics aimed at selectively modulating the RNA function by targeting protein/protein/RNA interfaces.

Exploring the Anticancer Activity of Tamoxifen-Based Metal Complexes Targeting Mitochondria.

Two new 'hybrid' metallodrugs of Au(III) (AuTAML) and Cu(II) (CuTAML) were designed featuring a tamoxifen-derived pharmacophore to ideally synergize the anticancer activity of both the metal center and the organic ligand. The compounds have antiproliferative effects against human MCF-7 and MDA-MB 231 breast cancer cells. Molecular dynamics studies suggest that the compounds retain the binding activity to estrogen receptor (ERα). In vitro and in silico studies showed that the Au(III) derivative is an inhibitor of the seleno-enzyme thioredoxin reductase, while the Cu(II) complex may act as an oxidant of different intracellular thiols. In breast cancer cells treated with the compounds, a redox imbalance characterized by a decrease in total thiols and increased reactive oxygen species production was detected. Despite their different reactivities and cytotoxic potencies, a great capacity of the metal complexes to induce mitochondrial damage was observed as shown by their effects on mitochondrial respiration, membrane potential, and morphology.

Cancer-Related Mutations Alter RNA-Driven Functional Cross-Talk Underlying Premature-Messenger RNA Recognition by Splicing Factor SF3b.

The pillar of faithful premature-messenger (pre-mRNA) splicing is the precise recognition of key intronic sequences by specific splicing factors. The heptameric splicing factor 3b (SF3b) recognizes the branch point sequence (BPS), a key part of the 3' splice site. SF3b contains SF3B1, a protein holding recurrent cancer-associated mutations. Among these, K700E, the most-frequent SF3B1 mutation, triggers aberrant splicing, being primarily implicated in hematologic malignancies. Yet, K700E and the BPS recognition site are 60 Å apart, suggesting the existence of an allosteric cross-talk between the two distal spots. Here, we couple molecular dynamics simulations and dynamical network theory analysis to unlock the molecular terms underpinning the impact of SF3b splicing factor mutations on pre-mRNA selection. We establish that by weakening and remodeling interactions of pre-mRNA with SF3b, K700E scrambles RNA-mediated allosteric cross-talk between the BPS and the mutation site. We propose that the altered allostery contributes to cancer-associated missplicing by mutated SF3B1. This finding broadens our comprehension of the elaborate mechanisms underlying pre-mRNA metabolism in eukaryotes.

Molecular Dynamics Simulations Elucidate the Molecular Basis of Pre-mRNA Translocation by the Prp2 Spliceosomal Helicase.

The spliceosome machinery catalyzes precursor-messenger RNA (pre-mRNA) splicing by undergoing at each splicing cycle assembly, activation, catalysis, and disassembly processes, thanks to the concerted action of specific RNA-dependent ATPases/helicases. Prp2, a member of the DExH-box ATPase/helicase family, harnesses the energy of ATP hydrolysis to translocate a single pre-mRNA strand in the 5' to 3' direction, thus promoting spliceosome remodeling to its catalytic-competent state. Here, we established the functional coupling between ATPase and helicase activities of Prp2. Namely, extensive multi-µs molecular dynamics simulations allowed us to unlock how, after pre-mRNA selection, ATP binding, hydrolysis, and dissociation induce a functional typewriter-like rotation of the Prp2 C-terminal domain. This movement, endorsed by an iterative swing of interactions established between specific Prp2 residues with the nucleobases at 5'- and 3'-ends of pre-mRNA, promotes pre-mRNA translocation. Notably, some of these Prp2 residues are conserved in the DExH-box family, suggesting that the translocation mechanism elucidated here may be applicable to all DExH-box helicases.

Monovalent Ionic Atmosphere Modulates the Selection of Suboptimal RNA Sequences by Splicing Factors' RNA Recognition Motifs.

The U2AF2 splicing factor is involved in the RNA recognition of the pre-mRNA poly-pyrimidine signaling sequence. This protein contains two RRM domains connected by a flexible linker, which ensure the preferential selection of a poly-uridine sequence over a poly-cytosine one. In this work, all-atom simulations provide insights into the U2AF2 recognition mechanism and on the features underlying its selectivity. Our outcomes show that U2AF2's RNA recognition is driven by cooperative events modulated by RNA-protein and RNA-ion interactions. Stunningly, monovalent ions contribute to mediating the binding of the weakly binding polyC strand, thus contributing to the selection of suboptimal poly-pyrimidine tracts. This finding broadens our understanding of the diverse traits tuning splicing factors' selectivity and adaptability to precisely handle and process diverse pre-mRNA sequences.

Switching from Aromatase Inhibitors to Dual Targeting Flavonoid-Based Compounds for Breast Cancer Treatment.

Despite the significant outcomes attained by scientific research, breast cancer (BC) still represents the second leading cause of death in women. Estrogen receptor-positive (ER+) BC accounts for the majority of diagnosed BCs, highlighting the disruption of estrogenic signalling as target for first-line treatment. This goal is presently pursued by inhibiting aromatase (AR) enzyme or by modulating Estrogen Receptor (ER) α. An appealing strategy for fighting BC and reducing side effects and resistance issues may lie in the design of multifunctional compounds able to simultaneously target AR and ER. In this paper, previously reported flavonoid-related potent AR inhibitors were suitably modified with the aim of also targeting ERα. As a result, homoisoflavone derivatives 3b and 4a emerged as well-balanced submicromolar dual acting compounds. An extensive computational study was then performed to gain insights into the interactions the best compounds established with the two targets. This study highlighted the feasibility of switching from single-target compounds to balanced dual-acting agents, confirming that a multi-target approach may represent a valid therapeutic option to counteract ER+ BC. The homoisoflavone core emerged as a valuable natural-inspired scaffold for the design of multifunctional compounds.

Molecular Basis of RNA-Driven ATP Hydrolysis in DExH-Box Helicases.

The spliceosome machinery catalyzes precursor messenger (pre-m)RNA splicing. In each cycle, the spliceosome experiences massive compositional and conformational remodeling fueled by the concerted action of specific RNA-dependent ATPases/helicases. Intriguingly, these enzymes are allosterically activated to perform ATP hydrolysis and trigger helicase activity only upon pre-mRNA binding. Yet, the molecular mechanism underlying the RNA-driven regulation of their ATPase function remains elusive. Here, we focus on the Prp2 ATPase/helicase which contributes to reshaping the spliceosome into its catalytic competent state. By performing classical and quantum-classical molecular dynamics simulations, we unprecedentedly unlock the molecular terms governing the Prp2 ATPase/helicase function. Namely, we dissect the molecular mechanism of ATP hydrolysis, and we disclose that RNA binding allosterically triggers the formation of a set of interactions linking the RNA binding tunnel to the catalytic site. This activates the Prp2's ATPase function by optimally placing the nucleophilic water and the general base of the enzymatic process to perform ATP hydrolysis. The key structural motifs, mechanically coupling RNA gripping and the ATPase/helicase functions, are conserved across all DExH-box helicases. This mechanism could thus be broadly applicable to all DExH-box helicase family.

Role of Monovalent Ions in the NKCC1 Inhibition Mechanism Revealed through Molecular Simulations.

The secondary active Na-K-Cl cotransporter 1 (NKCC1) promotes electroneutral uptake of two chloride ions, one sodium ion and one potassium ion. NKCC1 regulates Cl- homeostasis, thus being implicated in transepithelial water transport and in neuronal excitability. Aberrant NKCC1 transport is linked to a variety of human diseases. The loop diuretic drugs bumetanide, furosemide, azosemide and ethacrynic acid target NKCC1, but are characterized by poor selectivity leading to severe side effects. Despite its therapeutic importance, the molecular details of the NKCC1 inhibition mechanism remain unclear. Using all-atom simulations, we predict a putative binding mode of these drugs to the zebrafish (z) and human (h) NKCC1 orthologs. Although differing in their specific interactions with NKCC1 and/or monovalent ions, all drugs can fit within the same cavity and engage in hydrophobic interactions with M304/M382 in z/hNKCC1, a proposed ion gating residue demonstrated to be key for bumetanide binding. Consistent with experimental evidence, all drugs take advantage of the K+/Na+ ions, which plastically respond to their binding. This study not only provides atomic-level insights useful for drug discovery campaigns of more selective/potent NKCC1 inhibitors aimed to tackle diseases related to deregulated Cl- homeostasis, but it also supplies a paradigmatic example of the key importance of dynamical effects when drug binding is mediated by monovalent ions.

Intrinsically disordered ectodomain modulates ion permeation through a metal transporter.

The function of many channels and transporters is enriched by the conformational plasticity of intrinsically disordered regions (IDRs). Copper transporter 1 (Ctr1) is the main entry point for Cu(I) ions in eukaryotes and contains IDRs both at its N-terminal (Nterm) and C-terminal ends. The former delivers copper ions from the extracellular matrix to the selectivity filter in the Ctr1 lumen. However, the molecular mechanism of this process remains elusive due to Nterm's disordered nature. Here, we combine advanced molecular dynamics simulations and circular dichroism experiments to show that Cu(I) ions and a lipidic environment drive the insertion of the Nterm into the Ctr1 selectivity filter, causing its opening. Through a lipid-aided conformational switch of one of the transmembrane helices, the conformational change of the selectivity filter propagates down to the cytosolic gate of Ctr1. Taken together, our results elucidate how conformational variability of IDRs modulates ion transport.

Single-digit nanomolar inhibitors lock the aromatase active site via a dualsteric targeting strategy.

The most frequently diagnosed breast cancer (BC) type in women expresses estrogen receptor (ER) and depends on estrogens for its growth, being classified as ER positive (ER+). The gold standard therapy for the treatment of this tumor relies on the inhibition of the aromatase enzyme, which catalyzes estrogen biosynthesis. Despite the clinical success of current aromatase inhibitors (AIs), after prolonged therapeutic regimens, BC ER + patients experience acquired resistance and disease relapse. This points up the urgent need for a newer generation of AIs able to overcome resistance issues, while mitigating toxicity and side effects of current therapies. Here we performed the synthesis, biological evaluation, and extensive structural characterization by advanced molecular simulation methods of a new generation of dualsteric non-steroidal AIs, which simultaneously target the enzyme's active and allosteric sites. Notably, 3d, the most active AI of the series, exhibits a single-digit nM potency (IC50 2 nM). A detailed inspection of its binding mode reveals that the ancillary alkoxy chain predatorily takes advantage of the small hydrophobic cavities lining the allosteric site, triggering a remodeling of its residues and completely sealing the active site access-channel. As a result, the inhibitor is effectively locked in. This study sets a conceptual basis to develop a new generation of AIs exploiting a dualsteric targeting strategy.

Disrupting Cu trafficking as a potential therapy for cancer.

Copper ions play a crucial role in various cellular biological processes. However, these copper ions can also lead to toxicity when their concentration is not controlled by a sophisticated copper-trafficking system. Copper dys-homeostasis has been linked to a variety of diseases, including neurodegeneration and cancer. Therefore, manipulating Cu-trafficking to trigger selective cancer cell death may be a viable strategy with therapeutic benefit. By exploiting combined in silico and experimental strategies, we identified small peptides able to bind Atox1 and metal-binding domains 3-4 of ATP7B proteins. We found that these peptides reduced the proliferation of cancer cells owing to increased cellular copper ions concentration. These outcomes support the idea of harming copper trafficking as an opportunity for devising novel anti-cancer therapies.

Role of computational and structural biology in the development of small-molecule modulators of the spliceosome.

INTRODUCTION: RNA splicing is a pivotal step of eukaryotic gene expression during which the introns are excised from the precursor (pre-)RNA and the exons are joined together to form mature RNA products (i.e a protein-coding mRNA or long non-coding (lnc)RNAs). The spliceosome, a complex ribonucleoprotein machine, performs pre-RNA splicing with extreme precision. Deregulated splicing is linked to cancer, genetic, and neurodegenerative diseases. Hence, the discovery of small-molecules targeting core spliceosome components represents an appealing therapeutic opportunity. AREA COVERED: Several atomic-level structures of the spliceosome and distinct splicing-modulators bound to its protein/RNA components have been solved. Here, we review recent advances in the discovery of small-molecule splicing-modulators, discuss opportunities and challenges for their therapeutic applicability, and showcase how structural data and/or all-atom simulations can illuminate key facets of their mechanism, thus contributing to future drug-discovery campaigns. EXPERT OPINION: This review highlights the potential of modulating pre-RNA splicing with small-molecules, and anticipates how the synergy of computer and wet-lab experiments will enrich our understanding of splicing regulation/deregulation mechanisms. This information will aid future structure-based drug-discovery efforts aimed to expand the currently limited portfolio of selective splicing-modulators.