Thr729 in human topoisomerase I modulates anti-cancer drug resistance by altering protein domain communications as suggested by molecular dynamics simulations

Uncovering the selectivity mechanism of phosphodiesterase 7A/8A inhibitors through computational studies

The expression of phosphodiesterase 7A (PDE7A) and phosphodiesterase 8A (PDE8) genes is integral to human signaling pathways, and the inhibition of PDE7A has been associated with the onset of various diseases, including effects on the immune system and nervous system. The development of PDE7 selective inhibitors can promote research on immune and nervous system diseases, such as multiple sclerosis, chronic inflammation, and autoimmune responses. PDE8A is expressed alongside PDE8B, and its inhibitory mechanism is still unclear. Studying the mechanisms of selective inhibitors against different PDE subtypes is crucial to prevent potential side effects, such as nausea and cardiac toxicity, and the sequence similarity of the two protein subtypes was 55.9%. Therefore, it is necessary to investigate the differences of both subtypes' ligand binding sites. Selective inhibitors of two proteins were chosen to summarize the reason for their selectivity through molecular docking, molecular dynamics simulation, alanine scanning mutagenesis, and MM-GBSA calculation. We found that Phe384PDE7A, Leu401PDE7A, Gln413PDE7A, Tyr419PDE7A, and Phe416PDE7A in the active site positively contribute to the selectivity towards PDE7A. Additionally, Asn729PDE8A, Phe767PDE8A, Gln778PDE8A, and Phe781PDE8A positively contribute to the selectivity towards PDE8A.

In silico study to identify novel NEK7 inhibitors from natural sources by a combination strategy

Cancer poses a significant global health challenge and significantly contributes to mortality. NEK7, related to the NIMA protein kinase family, plays a crucial role in spindle assembly and cell division. The dysregulation of NEK7 is closely linked to the onset and progression of various cancers, especially colon and breast cancer, making it a promising target for cancer therapy. Nevertheless, the shortage of high-quality NEK7 inhibitors highlights the need for new therapeutic strategies. In this study, we utilized a multidisciplinary approach, including virtual screening, molecular docking, pharmacokinetics, molecular dynamics simulations (MDs), and MM/PBSA calculations, to evaluate natural compounds as NEK7 inhibitors comprehensively. Through various docking strategies, we identified three natural compounds: (-)-balanol, digallic acid, and scutellarin. Molecular docking revealed significant interactions at residues such as GLU112 and ALA114, with docking scores of -15.054, -13.059, and -11.547 kcal/mol, respectively, highlighting their potential as NEK7 inhibitors. MDs confirmed the stability of these compounds at the NEK7-binding site. Hydrogen bond analysis during simulations revealed consistent interactions, supporting their strong binding capacity. MM/PBSA analysis identified other crucial amino acids contributing to binding affinity, including ILE20, VAL28, ILE75, LEU93, ALA94, LYS143, PHE148, LEU160, and THR161, crucial for stabilizing the complex. This research demonstrated that these compounds exceeded dabrafenib in binding energy, according to MM/PBSA calculations, underscoring their effectiveness as NEK7 inhibitors. ADME/T predictions showed lower oral toxicity for these compounds, suggesting their potential for further development. This study highlights the promise of these natural compounds as bases for creating more potent derivatives with significant biological activities, paving the way for future experimental validation.

Structural insight to human Retinoid X receptor alpha-Thyroid hormone receptor beta heterodimer by molecular modelling and MD-simulation studies: role of conserved water molecules

The Retinoid X receptor alpha-Thyroid hormone receptor beta (RXRα-THRβ) heterodimer plays an important role in physiological function of humans specially in the growth and development. Extensive MD-simulation studies on the aquated complexes of modelled RXRα-THRβ heterodimer with DNA-duplex have indicated the role of some conserved/semiconserved water molecules in the complexation process in presence or absence of Triiodothyronine (T3) and 9-cis retinoic acid (9CR) in the respective Ligand Binding Domain (LBD) domain. Among the seventeen conserved/semi-conserved water molecules, the W1-W4 water centers have been observed to mediate the interaction between the residues of A-chain (DBD of RXR) to consensus sequence (C-chain) of DNA. The W5-W8 water centers involve in recognition of the residues of B-chain (DBD of THR) to C-chain of DNA. The W9-W13 centers have connected the different residues of B-chain (THR) to D-chain of DNA through H-bonds, whereas W14-W17 water molecules were involved in the interaction of A-chain's (RXR) residues to D-chain of DNA. In our previous study with homodimeric THRβ from Rattus norvegicus we have identified fifteen conserved water molecules at the DNA-DBD interface. Moreover, the conformational flexibility of Met313 (in the LBD of THR) from open to close form in presence or absence of T3 molecule in the holo and Apo-protein may provide a plausible rational on the possible role of that residue to acts as gate which could restrict the solvent molecules to enter into the hydrophobic T3-binding pocket of LBD during the absence of ligand molecule and thus could help the stabilization of that domain in THRβ structure.Communicated by Ramaswamy H. Sarma.

Dynamozones are the most obvious sign of the evolution of conformational dynamics in HIV-1 protease

Proteins are not static but are flexible molecules that can adopt many different conformations. The HIV-1 protease is an important target for the development of therapies to treat AIDS, due to its critical role in the viral life cycle. We investigated several dynamics studies on the HIV-1 protease families to illustrate the significance of examining the dynamic behaviors and molecular motions for an entire understanding of their dynamics-structure-function relationships. Using computer simulations and principal component analysis approaches, the dynamics data obtained revealed that: (i) The flap regions are the most obvious sign of the evolution of conformational dynamics in HIV-1 protease; (ii) There are dynamic structural regions in some proteins that contribute to the biological function and allostery of proteins via appropriate flexibility. These regions are a clear sign of the evolution of conformational dynamics of proteins, which we call dynamozones. The flap regions are one of the most important dynamozones members that are critical for HIV-1 protease function. Due to the existence of other members of dynamozones in different proteins, we propose to consider dynamozones as a footprint of the evolution of the conformational dynamics of proteins.

Remdesivir and Ledipasvir among the FDA-Approved Antiviral Drugs Have Potential to Inhibit SARS-CoV-2 Replication

The rapid spread of the virus, the surge in the number of deaths, and the unavailability of specific SARS-CoV-2 drugs thus far necessitate the identification of drugs with anti-COVID-19 activity. SARS-CoV-2 enters the host cell and assembles a multisubunit RNA-dependent RNA polymerase (RdRp) complex of viral nonstructural proteins that plays a substantial role in the transcription and replication of the viral genome. Therefore, RdRp is among the most suitable targets in RNA viruses. Our aim was to investigate the FDA approved antiviral drugs having potential to inhibit the viral replication. The methodology adopted was virtual screening and docking of FDA-approved antiviral drugs into the RdRp protein. Top hits were selected and subjected to molecular dynamics simulations to understand the dynamics of RdRp in complex with these drugs. The antiviral activity of the drugs against SARS-CoV-2 was assessed in Vero E6 cells. Notably, both remdesivir (half-maximal effective concentration (EC50) 6.6 μM, 50% cytotoxicity concentration (CC50) > 100 µM, selectivity index (SI) = 15) and ledipasvir (EC50 34.6 μM, CC50 > 100 µM, SI > 2.9) exerted antiviral action. This study highlights the use of direct-acting antiviral drugs, alone or in combination, for better treatments of COVID-19.

Molecular modeling and molecular dynamics simulation studies on thyroid hormone receptor from Rattus norvegicus: role of conserved water molecules

Thyroid hormone receptor (THR) belongs to the nuclear receptor (NR) superfamily that is activated by binding of appropriate ligand molecules (thyroid hormones). These receptors directly bind to specific DNA sequences for gene expression, which is essential for metabolism, homeostasis, and the development of organisms, making it an important drug target. Extensive MD-simulation studies of triiodothyronine (T3) docked modeled rnTHRβ1 structures have indicated the presence of twelve conserved water molecules at the DNA-DBD (DNA binding domain) interface. The W1-W5 water centers have been involved in the recognition between the A-chain of DBD to C-chain of DNA, W6 and W7 mediated the interaction between A-chain of DBD and D-chain of DNA, W8 and W9 recognized the B-chain of DBD and C-chain of DNA, and W9-W12 centers conjugated the residues of B-chain of DBD to D-chain of DNA through hydrogen bonds. The conformation flexibility of Phe272 and Met313 residues in the absence of T3 at the LBD (ligand-binding domain) region have been observed and reported.

The benzylisoquinoline alkaloids, berberine and coptisine, act against camptothecin-resistant topoisomerase I mutants

DNA replication inhibitors are utilized extensively in studies of molecular biology and as chemotherapy agents in clinical settings. The inhibition of DNA replication often triggers double-stranded DNA breaks (DSBs) at stalled DNA replication sites, resulting in cytotoxicity. In East Asia, some traditional medicines are administered as anticancer drugs, although the mechanisms underlying their pharmacological effects are not entirely understood. In this study, we screened Japanese herbal medicines and identified two benzylisoquinoline alkaloids (BIAs), berberine and coptisine. These alkaloids mildly induced DSBs, and this effect was dependent on the function of topoisomerase I (Topo I) and MUS81-EME1 structure-specific endonuclease. Biochemical analysis revealed that the action of BIAs involves inhibiting the catalytic activity of Topo I rather than inducing the accumulation of the Topo I-DNA complex, which is different from the action of camptothecin (CPT). Furthermore, the results showed that BIAs can act as inhibitors of Topo I, even against CPT-resistant mutants, and that the action of these BIAs was independent of CPT. These results suggest that using a combination of BIAs and CPT might increase their efficiency in eliminating cancer cells.

Exploring the dynamic mechanism of allosteric drug SHP099 inhibiting SHP2E69K

The E69K mutation is one of the most frequent protein tyrosine phosphatase-2 (SHP2) mutations in leukemia, and it can cause the increase in the protein activity. Recent studies have shown that the E69K mutation was fairly sensitive to the allosteric inhibitor of SHP2 (SHP099). However, the molecular mechanism of the allosteric drug SHP099 inhibiting SHP2^E69K remains unclear. Thus, the molecular dynamic simulations and the post-dynamics analyses (RMSF, PCA, DCCM, RIN and the binding free energies) for SHP2^WT, SHP2^WT-SHP099, SHP2^E69K and SHP2^E69K-SHP099 were carried out, respectively. Owing to the strong binding affinity of SHP099 to residues Thr219 and Arg220, the flexibility of linker region (residues Val209-Arg231) was reduced. Moreover, the presence of SHP099 kept the autoinhibition state of the SHP2 protein through enhancing the interactions between the linker region and Q loop in PTP domain, such as Thr219/Val490, Thr219/Asn491, Arg220/Ile488 and Leu254/Asn491. In addition, it was found that the residues (Thr219, Arg220, Leu254 and Asn491) might be the key residues responsible for the conformational changes of protein. Overall, this study may provide an important basis for understanding how the SHP099 effectively inhibited the SHP2^E69K activity at the molecular level.

Identification of the potential dual inhibitor of protein tyrosine phosphatase sigma and leukocyte common antigen-related phosphatase by virtual screen, molecular dynamic simulations and post-analysis

Owing to their inhibitory role in regulating oligodendrocyte differentiation and apoptosis, protein tyrosine phosphatase sigma (PTPσ) and leukocyte common antigen-related phosphatase (LAR) play a crucial potential role in treating spinal cord injury (SCI) disease. In this research, the computer aided drug design (CADD) methods were applied to discover the potential dual-target drug involving virtual screen, molecular docking and molecular dynamic simulation. Initially, the top 20 compounds with higher docking score than the positive controls (ZINC13749892, ZINC14516161) were virtually screened out from NCI and ZINC databases, and then were submitted in ADMET to predict their drug properties. Among these potential compounds, ZINC72417086 showed a higher docking score and satisfied Lipinski's rule of five. In addition, the post-analysis demonstrated that when ZINC72417086 bound to PTPσ and LAR, it could stable proteins conformations and destroy the residues interactions between P-loop and other loop regions in active pocket. Meanwhile, residue ARG1595 and ARG1528 could play a crucial role in in the inhibition of PTPσ and LAR, respectively. This research offered a novel approach for rapid discovery of dual-target leads compounds to treat SCI.Communicated by Ramaswamy H. Sarma.

Structural and dynamic analysis of G558R mutation in chicken TSHR gene shows altered signal transduction and corroborates its role as a domestication gene

The thyroid-stimulating hormone receptor (TSHR) has been indicated as a putative domestication gene in chicken. Comparison of WGS identified a variant in residue 558 of the transmembrane domain (TM) of TSHR, where the domestic chicken (GGD) presents an arginine, whereas the red jungle fowl (RJF) shares a conserved glycine with other vertebrates. This variant has been demonstrated to be associated with phenotypes that are important for domestication and related to thyroid regulation, such as less fearful behavior, reduced aggressive behavior and reduced dependence on seasonal reproduction in GGD as compared with RJF. By means of molecular dynamics simulations, we highlighted the structural and dynamic differences of variant Gly558Arg in the TSHR TM domain. Alterations in TM helix flexibility, structure and protein overall motion are described. The so-called 'arginine snorkeling' of residue 568 in GGD is observed and we hypothesize it as the originating force that produces the observed whole-protein perturbation in the helix bundle dynamics, capable of altering the TSHR signal transduction. The results are discussed in the context of their implications for a better understanding of biological mechanisms in chicken under control of the thyroid, such as body metabolism, as well as for their usefulness in biomedical research.

Mathematical modeling and computational prediction of cancer drug resistance

Diverse forms of resistance to anticancer drugs can lead to the failure of chemotherapy. Drug resistance is one of the most intractable issues for successfully treating cancer in current clinical practice. Effective clinical approaches that could counter drug resistance by restoring the sensitivity of tumors to the targeted agents are urgently needed. As numerous experimental results on resistance mechanisms have been obtained and a mass of high-throughput data has been accumulated, mathematical modeling and computational predictions using systematic and quantitative approaches have become increasingly important, as they can potentially provide deeper insights into resistance mechanisms, generate novel hypotheses or suggest promising treatment strategies for future testing. In this review, we first briefly summarize the current progress of experimentally revealed resistance mechanisms of targeted therapy, including genetic mechanisms, epigenetic mechanisms, posttranslational mechanisms, cellular mechanisms, microenvironmental mechanisms and pharmacokinetic mechanisms. Subsequently, we list several currently available databases and Web-based tools related to drug sensitivity and resistance. Then, we focus primarily on introducing some state-of-the-art computational methods used in drug resistance studies, including mechanism-based mathematical modeling approaches (e.g. molecular dynamics simulation, kinetic model of molecular networks, ordinary differential equation model of cellular dynamics, stochastic model, partial differential equation model, agent-based model, pharmacokinetic–pharmacodynamic model, etc.) and data-driven prediction methods (e.g. omics data-based conventional screening approach for node biomarkers, static network approach for edge biomarkers and module biomarkers, dynamic network approach for dynamic network biomarkers and dynamic module network biomarkers, etc.). Finally, we discuss several further questions and future directions for the use of computational methods for studying drug resistance, including inferring drug-induced signaling networks, multiscale modeling, drug combinations and precision medicine.

Exploring the effect of N308D mutation on protein tyrosine phosphatase-2 cause gain-of-function activity by a molecular dynamics study

One of the most common protein tyrosine phosphatase-2 (SHP2) mutations in Noonan syndrome is the N308D mutation, and it increases the activity of the protein. However, the molecular basis of the activation of N308D mutation on SHP2 conformations is poorly understood. Here, molecular dynamic simulations were performed on SHP2 and SHP2-N308D to explore the effect of N308D mutation on SHP2 cause gain of function activity, respectively. The principal component analysis, dynamic cross-correlation map, secondary structure analysis, residue interaction networks, and solvent accessible surface area analysis suggested that the N308D mutation distorted the residues interactions network between the allosteric site (residue Gly244-Gly246) and C-SH2 domain, including the hydrogen bond formation and the binding energy. Meanwhile, the activity of catalytic site (residue Gly503-Val505) located in the Q-loop in mutant increased due to this region's high fluctuations. Therefore, the substrate had more chances to access to the catalytic activity site of the precision time protocol domain of SHP2-N308D, which was easy to be exposed. In addition, we had speculated that the Lys244 located in the allosteric site was the key residue which lead to the protein conformation changes. Consequently, overall calculations presented in this study ultimately provide a useful understanding of the increased activity of SHP2 caused by the N308D mutation.

Structural insights into the Middle East respiratory syndrome coronavirus 4a protein and its dsRNA binding mechanism

Middle East respiratory syndrome coronavirus (MERS-CoV) has evolved to navigate through the sophisticated network of a host's immune system. The immune evasion mechanism including type 1 interferon and protein kinase R-mediated antiviral stress responses has been recently attributed to the involvement of MERS-CoV protein 4a (p4a) that masks the viral dsRNA. However, the structural mechanism of how p4a recognizes and establishes contacts with dsRNA is not well explained. In this study, we report a dynamic mechanism deployed by p4a to engage the viral dsRNA and make it unavailable to the host immune system. Multiple variants of p4a-dsRNA were created and investigated through extensive molecular dynamics procedures to highlight crucial interfacial residues that may be used as potential pharmacophores for future drug development. The structural analysis revealed that p4a exhibits a typical αβββα fold structure, as found in other dsRNA-binding proteins. The α1 helix and the β1-β2 loop play a crucial role in recognizing and establishing contacts with the minor grooves of dsRNA. Further, mutational and binding free energy analyses suggested that in addition to K63 and K67, two other residues, K27 and W45, might also be crucial for p4a-dsRNA stability.

Advantages of an optical nanosensor system for the mechanistic analysis of a novel topoisomerase I targeting drug: a case study

The continuous need for the development of new small molecule anti-cancer drugs calls for easily accessible sensor systems for measuring the effect of vast numbers of new drugs on their potential cellular targets. Here we demonstrate the use of an optical DNA biosensor to unravel the inhibitory mechanism of a member of a new family of small molecule human topoisomerase I inhibitors, the so-called indeno-1,5-naphthyridines. By analysing human topoisomerase I catalysis on the biosensor in the absence or presence of added drug complemented with a few traditional assays, we demonstrate that the investigated member of the indeno-1,5-naphthyridine family inhibited human topoisomerase I activity by blocking enzyme-DNA dissociation. To our knowledge, this represents the first characterized example of a small molecule drug that inhibits a post-ligation step of catalysis. The elucidation of a completely new and rather surprising drug mechanism-of-action using an optical real time sensor highlights the value of this assay system in the search for new topoisomerase I targeting small molecule drugs.

Inhibition of human DNA topoisomerase IB by nonmutagenic ruthenium(II)-based compounds with antitumoral activity

Herein we synthesized two new ruthenium(II) compounds [Ru(pySH)(bipy)(dppb)]PF6 (1) and [Ru(HSpym)(bipy)(dppb)]PF6 (2) that are analogs to an antitumor agent recently described, [Ru(SpymMe2)(bipy)(dppb)]PF6 (3), where [(Spy) = 2-mercaptopyridine anion; (Spym) = 2-mercaptopyrimidine anion and (SpymMe2) = 4,6-dimethyl-2-mercaptopyrimidine anion]. In vitro cell culture experiments revealed significant anti-proliferative activity for 1-3 against HepG2 and MDA-MB-231 tumor cells, higher than the standard anti-cancer drugs doxorubicin and cisplatin. No mutagenicity is detected when compounds are evaluated by cytokinesis-blocked micronucleus cytome and Ames test in the presence and absence of S9 metabolic activation from rat liver. Interaction studies show that compounds 1-3 can bind to DNA through electrostatic interactions and to albumin through hydrophobic interactions. The three compounds are able to inhibit the DNA supercoiled relaxation mediated by human topoisomerase IB (Top1). Compound 3 is the most efficient Top1 inhibitor and the inhibitory effect is enhanced upon pre-incubation with the enzyme. Analysis of different steps of Top1 catalytic cycle indicates that 3 inhibits the cleavage reaction impeding the binding of the enzyme to DNA and slows down the religation reaction. Molecular docking shows that 3 preferentially binds closer to the residues of the active site when Top1 is free and lies on the DNA groove downstream of the cleavage site in the Top1-DNA complex. Thus, 3 can be considered in further studies for a possible use as an anticancer agent.

Human topoisomerase IB is a target of a thiosemicarbazone copper(II) complex

The human topoisomerase IB inhibition and the antiproliferative activity of 3-(4-bromophenyl)-1-pyridin-2-ylprop-2-en-1-one thiosemicarbazone HPyCT4BrPh alone and its copper(II) complex [Cu(PyCT4BrPh)Cl] was investigated. [Cu(PyCT4BrPh)Cl] inhibits both the DNA cleavage and religation step of the enzyme, whilst the ligand alone does not display any effect. In addition we show that coordination to copper(II) improves the cytotoxicity of HPyCT4BrPh against THP-1 leukemia and MCF-7 breast cancer cells. The data indicate that the copper(II) thiosemicarbazone complex may hit human topoisomerase IB and that metal coordination can be useful to improve cytotoxicity of this versatile class of compounds.

Geometrical constraints limiting the poly(ADP-ribose) conformation investigated by molecular dynamics simulation

Poly(ADP-ribosylation) is a post-transductional modification that regulates protein's function. Most of the proteins subjected to this control mechanism belong to machineries involved in DNA damage repair, or DNA interacting proteins. Poly(ADP-ribose) polymers are long chains of even 100 monomer length that can be branched at several positions but, not withstanding its importance, nothing is known concerning its structure. To understand, which are the geometrical parameters that confer to the polymer the structural constraints that determine its interaction with the target proteins, we have performed molecular dynamics of three chains of different length, made by 5, 25, and 30 units, the last one being branched. Analysis of the simulations allowed us to identify the main intra- and inter-monomer dihedral angles that govern the structure of the polymer that however, does not reach a unique definite conformation.

Smyd3 open & closed lock mechanism for substrate recruitment: The hinge motion of C-terminal domain inferred from μ-second molecular dynamics simulations

BACKGROUND

The human lysine methyltransferase Smyd3, a member of the SET and MYND domain containing protein family, harbors methylation activity on both histone and non-histone targets in a tightly regulated manner. The mechanism of how Smyd3 dynamically regulates substrate recognition is still not fully unveiled.

METHODS

Here, we employed molecular dynamics simulations on full length human Smyd3, performed to a total of 1.2 μ-second, in the presence (holo) and absence (apo) of the S-Adenosyl methionine (AdoMet) cofactor. The dynamical features of Smyd3 in apo and holo states have been examined and compared via examining geometrical and electrostatic properties.

RESULTS

The results show a distinct dynamics of the C-terminal domain (CTD) in the two states. In the apo state, the CTD undergoes a large hinge like motion and samples more opened configurations, thus acting like a loosened clamp and resulting in expanded substrate binding crevice. In the holo state, the CTD exhibits a restricted motion while the overall structure remains compact, mimicking a closed clamp. This leads to a localized increase in the negative potential at the substrate binding cleft. Further, solvent accessibility of critical residues at the target lysine access channel, important for methylation activity, is increased.

CONCLUSIONS

We postulate that AdoMet cofactor acts like a key and locks Smyd3 in a closed conformation. In effect, the cofactor binding restricts the elasticity of the CTD, presenting a compact substrate binding cleft with high negative potential, which may have implications on substrate recruitment via long range electrostatics.

GENERAL SIGNIFICANCE

The deletion of the CTD from Smyd3 has been shown to abolish the basal histone methylation activity. Our study highlights the importance of the CTD elasticity in shaping the substrate binding site for recognition and supports the previously proposed role of the CTD in stabilizing the active site for methylation activity.

Structural and Dynamic Characterization of the C313Y Mutation in Myostatin Dimeric Protein, Responsible for the "Double Muscle" Phenotype in Piedmontese Cattle

The knowledge of the molecular effects of the C313Y mutation, responsible for the “double muscle” phenotype in Piedmontese cattle, can help understanding the actual mechanism of phenotype determination and paves the route for a better modulation of the positive effects of this economic important phenotype in the beef industry, while minimizing the negative side effects, now inevitably intersected. The structure and dynamic behavior of the active dimeric form of Myostatin in cattle was analyzed by means of three state-of-the-art Molecular Dynamics simulations, 200-ns long, of wild-type and C313Y mutants. Our results highlight a role for the conserved Arg333 in establishing a network of short and long range interactions between the two monomers in the wild-type protein that is destroyed upon the C313Y mutation even in a single monomer. Furthermore, the native protein shows an asymmetry in residue fluctuation that is absent in the double monomer mutant. Time window analysis on further 200-ns of simulation demonstrates that this is a characteristic behavior of the protein, likely dependent on long range communications between monomers. The same behavior, in fact, has already been observed in other mutated dimers. Finally, the mutation does not produce alterations in the secondary structure elements that compose the characteristic TGF-β cystine-knot motif.

Structural Mechanism behind Distinct Efficiency of Oct4/Sox2 Proteins in Differentially Spaced DNA Complexes

The octamer-binding transcription factor 4 (Oct4) and sex-determining region Y (SRY)-box 2 (Sox2) proteins induce various transcriptional regulators to maintain cellular pluripotency. Most Oct4/Sox2 complexes have either 0 base pairs (Oct4/Sox2^0bp) or 3 base pairs (Oct4/Sox2^3bp) separation between their DNA-binding sites. Results from previous biochemical studies have shown that the complexes separated by 0 base pairs are associated with a higher pluripotency rate than those separated by 3 base pairs. Here, we performed molecular dynamics (MD) simulations and calculations to determine the binding free energy and per-residue free energy for the Oct4/Sox2^0bp and Oct4/Sox2^3bp complexes to identify structural differences that contribute to differences in induction rate. Our MD simulation results showed substantial differences in Oct4/Sox2 domain movements, as well as secondary-structure changes in the Oct4 linker region, suggesting a potential reason underlying the distinct efficiencies of these complexes during reprogramming. Moreover, we identified key residues and hydrogen bonds that potentially facilitate protein-protein and protein-DNA interactions, in agreement with previous experimental findings. Consequently, our results confess that differential spacing of the Oct4/Sox2 DNA binding sites can determine the magnitude of transcription of the targeted genes during reprogramming.

Recovery of the poisoned topoisomerase II for DNA religation: coordinated motion of the cleavage core revealed with the microsecond atomistic simulation

Type II topoisomerases resolve topological problems of DNA double helices by passing one duplex through the reversible double-stranded break they generated on another duplex. Despite the wealth of information in the cleaving operation, molecular understanding of the enzymatic DNA ligation remains elusive. Topoisomerase poisons are widely used in anti-cancer and anti-bacterial therapy and have been employed to entrap the intermediates of topoisomerase IIβ with religatable DNA substrate. We removed drug molecules from the structure and conducted molecular dynamics simulations to investigate the enzyme-mediated DNA religation. The drug-unbound intermediate displayed transitions toward the resealing-compliant configuration: closing distance between the cleaved DNA termini, B-to-A transformation of the double helix, and restoration of the metal-binding motif. By mapping the contact configurations and the correlated motions between enzyme and DNA, we identified the indispensable role of the linker preceding winged helix domain (WHD) in coordinating the movements of TOPRIM, the nucleotide-binding motifs, and the bound DNA substrate during gate closure. We observed a nearly vectorial transition in the recovery of the enzyme and identified the previously uncharacterized roles of Asn508 and Arg677 in DNA rejoining. Our findings delineate the dynamic mechanism of the DNA religation conducted by type II topoisomerases.

Structural insight of DNA topoisomerases I from camptothecin-producing plants revealed by molecular dynamics simulations

DNA topoisomerase I (Top1) catalyzes changes in DNA topology by cleaving and rejoining one strand of the double stranded (ds)DNA. Eukaryotic Top1s are the cellular target of the plant-derived anticancer indole alkaloid camptothecin (CPT), which reversibly stabilizes the Top1-dsDNA complex. However, CPT-producing plants, including Camptotheca acuminata, Ophiorrhiza pumila and Ophiorrhiza liukiuensis, are highly resistant to CPT because they possess point-mutated Top1. Here, the adaptive convergent evolution is reported between CPT production ability and mutations in their Top1, as a universal resistance mechanism found in all tested CPT-producing plants. This includes Nothapodytes nimmoniana, one of the major sources of CPT. To obtain a structural insight of the resistance mechanism, molecular dynamics simulations of CPT- resistant and -sensitive plant Top1s complexed with dsDNA and topotecan (a CPT derivative) were performed, these being compared to that for the CPT-sensitive human Top1. As a result, two mutations, Val617Gly and Asp710Gly, were identified in O. pumila Top1 and C. acuminata Top1, respectively. The substitutions at these two positions, surprisingly, are the same as those found in a CPT derivative-resistant human colon adenocarcinoma cell line. The results also demonstrated an increased linker flexibility of the CPT-resistant Top1, providing an additional explanation for the resistance mechanism found in CPT-producing plants. These mutations could reflect the long evolutionary adaptation of CPT-producing plant Top1s to confer a higher degree of resistance.

Mutation of Gly717Phe in human topoisomerase 1B has an effect on enzymatic function, reactivity to the camptothecin anticancer drug and on the linker domain orientation

Human topoisomerase 1B controls the topological state of supercoiled DNA allowing the progression of fundamental cellular processes. The enzyme, which is the unique molecular target of the natural anticancer compound camptothecin, acts by cleaving one DNA strand and forming a transient protein-DNA covalent adduct. In this work the role of the Gly717 residue, located in a α-helix structure bridging the active site and the linker domain, has been investigated mutating it in Phe. The mutation gives rise to drug resistance in vivo as observed through a viability assay of yeast cells. In vitro activity assays show that the mutant is characterized by a fast religation rate, only partially reduced by the presence of the drug. Comparative molecular dynamics simulations of the native and mutant proteins indicate that the mutation of Gly717 affects the motion orientation of the linker domain, changing its interaction with the DNA substrate, likely affecting the strand rotation and religation rate. The mutation also causes a slight rearrangement of the active site and of the drug binding site, providing an additional explanation for the lowered effect of camptothecin toward the mutant.

Simulations of DNA topoisomerase 1B bound to supercoiled DNA reveal changes in the flexibility pattern of the enzyme and a secondary protein-DNA binding site

Human topoisomerase 1B has been simulated covalently bound to a negatively supercoiled DNA minicircle, and its behavior compared to the enzyme bound to a simple linear DNA duplex. The presence of the more realistic supercoiled substrate facilitates the formation of larger number of protein–DNA interactions when compared to a simple linear duplex fragment. The number of protein–DNA hydrogen bonds doubles in proximity to the active site, affecting all of the residues in the catalytic pentad. The clamp over the DNA, characterized by the salt bridge between Lys369 and Glu497, undergoes reduced fluctuations when bound to the supercoiled minicircle. The linker domain of the enzyme, which is implicated in the controlled relaxation of superhelical stress, also displays an increased number of contacts with the minicircle compared to linear DNA. Finally, the more complex topology of the supercoiled DNA minicircle gives rise to a secondary DNA binding site involving four residues located on subdomain III. The simulation trajectories reveal significant changes in the interactions between the enzyme and the DNA for the more complex DNA topology, which are consistent with the experimental observation that the protein has a preference for binding to supercoiled DNA.

Solvent dependency of the UV-Vis spectrum of indenoisoquinolines: role of keto-oxygens as polarity interaction probes

Indenoisoquinolines are the most promising non-campthotecins topoisomerase IB inhibitors. We present an integrated experimental/computational investigation of the UV-Vis spectra of the IQNs parental compound (NSC314622) and two of its derivatives (NSC724998 and NSC725776) currently undergoing Phase I clinical trials. In all the three compounds a similar dependence of the relative absorption intensities at 270 nm and 290 nm on solvent polarity is found. The keto-oxygens in positions 5 and 11 of the molecular scaffold of the molecule are the principal chromophores involved in this dependence. Protic interactions on these sites are also found to give rise to absorptions at wavelength <250 nm observed in water solution, due to the stabilization of highly polarized tautomers of the molecule. These results suggest that the keto-oxygens are important polarizable groups that can act as useful interactors with the molecular receptor, providing at the same time an useful fingerprint for the monitoring of the drug binding to topoisomerase IB.

Sequence selectivity of the cleavage sites induced by topoisomerase I inhibitors: a molecular dynamics study

Topoisomerase IB (Top1) inhibitors, such as camptothecin (CPT), stabilize the Top1-DNA cleavage complex in a DNA sequence-dependent manner. The sequence selectivity of Top1 inhibitors is important for targeting specific genomic sequences of therapeutic value. However, the molecular mechanisms underlying this selectivity remain largely unknown. We performed molecular dynamics simulations to delineate structural, dynamic and energetic features that contribute to the differential sequence selectivity of the Top1 inhibitors. We found the sequence selectivity of CPT to be highly correlated with the drug binding energies, dynamic and structural properties of the linker domain. Chemical insights, gained by per-residue binding energy analysis revealed that the non-polar interaction between CPT and nucleotide at the +1 position of the cleavage site was the major (favorable) contributor to the total binding energy. Mechanistic insights gained by a potential of mean force analysis implicated that the drug dissociation step was associated with the sequence selectivity. Pharmaceutical insights gained by our molecular dynamics analyses explained why LMP-776, an indenoisoquinoline derivative under clinical development at the National Institutes of Health, displays different sequence selectivity when compared with camptothecin and its clinical derivatives.

Molecular mechanism of the camptothecin resistance of Glu710Gly topoisomerase IB mutant analyzed in vitro and in silico

BACKGROUND

DNA topoisomerases are key enzymes that modulate the topological state of DNA through the breaking and rejoining of DNA strands. Human topoisomerase IB can be inhibited by several compounds that act through different mechanisms, including clinically used drugs, such as the derivatives of the natural compound camptothecin that reversibly bind the covalent topoisomerase-DNA complex, slowing down the religation of the cleaved DNA strand, thus inducing cell death. Three enzyme mutations, which confer resistance to irinotecan in an adenocarcinoma cell line, were recently identified but the molecular mechanism of resistance was unclear.

METHODS

The three resistant mutants have been investigated in S. cerevisiae model system following their viability in presence of increasing amounts of camptothecin. A systematical analysis of the different catalytic steps has been made for one of these mutants (Glu710Gly) and has been correlated with its structural-dynamical properties studied by classical molecular dynamics simulation.

RESULTS

The three mutants display a different degree of camptothecin resistance in a yeast cell viability assay. Characterization of the different steps of the catalytic cycle of the Glu710Gly mutant indicated that its resistance is related to a high religation rate that is hardly affected by the presence of the drug. Analysis of the dynamic properties through simulation indicate that the mutant displays a much lower degree of correlation in the motion between the different protein domains and that the linker almost completely loses its correlation with the C-terminal domain, containing the active site tyrosine.

CONCLUSIONS

These results indicate that a fully functional linker is required to confer camptothecin sensitivity to topoisomerase I since the destabilization of its structural-dynamical properties is correlated to an increase of religation rate and drug resistance.

Structural and functional interactions of the prostate cancer suppressor protein NKX3.1 with topoisomerase I

NKX3.1 (NK3 homeobox 1) is a prostate tumour suppressor protein with a number of activities that are critical for its role in tumour suppression. NKX3.1 mediates the cellular response to DNA damage by interacting with ATM (ataxia telangiectasia mutated) and by activation of topoisomerase I. In the present study we characterized the interaction between NKX3.1 and topoisomerase I. The NKX3.1 homeodomain binds to a region of topoisomerase I spanning the junction between the core and linker domains. Loss of the topoisomerase I N-terminal domain, a region for frequent protein interactions, did not affect binding to NKX3.1 as was shown by the activation of Topo70 (N-terminal truncated topoisomerase I) in vitro. In contrast, NKX3.1 interacts with the enzyme reconstituted from peptide fragments of the core and linker active site domains, but inhibits the DNA-resolving activity of the reconstituted enzyme in vitro. The effect of NKX3.1 on both Topo70 and the reconstituted enzyme was seen in the presence and absence of camptothecin. Neither NKX3.1 nor CPT (camptothecin) had an effect on the interaction of the other with topoisomerase I. Therefore the interactions of NKX3.1 and CPT with the linker domain of topoisomerase I are mutually exclusive. However, in cells the effect of NKX3.1 on topoisomerase binding to DNA sensitized the cells to cellular toxicity and the induction of apoptosis by low doses of CPT. Lastly, topoisomerase I is important for the effect of NKX3.1 on cell survival after DNA damage as topoisomerase knockdown blocked the effect of NKX3.1 on clonogenicity after DNA damage. Therefore NKX3.1 and topoisomerase I interact in vitro and in cells to affect the CPT sensitivity and DNA-repair functions of NKX3.1.

Replacement of the human topoisomerase linker domain with the plasmodial counterpart renders the enzyme camptothecin resistant

A human/plasmodial hybrid enzyme, generated by swapping the human topoisomerase IB linker domain with the corresponding domain of the Plasmodium falciparum enzyme, has been produced and characterized. The hybrid enzyme displays a relaxation activity comparable to the human enzyme, but it is characterized by a much faster religation rate. The hybrid enzyme is also camptothecin resistant. A 3D structure of the hybrid enzyme has been built and its structural-dynamical properties have been analyzed by molecular dynamics simulation. The analysis indicates that the swapped plasmodial linker samples a conformational space much larger than the corresponding domain in the human enzyme. The large linker conformational variability is then linked to important functional properties such as an increased religation rate and a low drug reactivity, demonstrating that the linker domain has a crucial role in the modulation of the topoisomerase IB activity.

Molecular dynamics simulations show how the FMRP Ile304Asn mutation destabilizes the KH2 domain structure and affects its function

Mutations or deletions of FMRP, involved in the regulation of mRNA metabolism in brain, lead to the Fragile X syndrome (FXS), the most frequent form of inherited intellectual disability. A severe manifestation of the disease has been associated with the Ile304Asn mutation, located on the KH2 domain of the protein. Several hypotheses have been proposed to explain the possible molecular mechanism responsible for the drastic effect of this mutation in humans. Here, we performed a molecular dynamics simulation and show that the Ile304Asn mutation destabilizes the hydrophobic core producing a partial unfolding of two α-helices and a displacement of a third one. The affected regions show increased residue flexibility and motion. Molecular docking analysis revealed strongly reduced binding to a model single-stranded nucleic acid in agreement with known data that the two partially unfolded helices form the RNA-binding surface. The third helix, which we show here to be also affected, is involved in the PAK1 protein interaction. These two functional binding sites on the KH2 domain do not overlap spatially, and therefore, they can simultaneously bind their targets. Since the Ile304Asn mutation affects both binding sites, this may justify the severe clinical manifestation observed in the patient in which both mRNA metabolism activity and cytoskeleton remodeling would be affected.

Characterization of DNA topoisomerase-1 in Spodoptera exigua for toxicity evaluation of camptothecin and hydoxy-camptothecin

Camptothecin (CPT), a plant alkaloid originally isolated from the native Chinese tree, Camptotheca acuminate, exerts the toxic effect by targeting eukaryotic DNA topoisomerase 1 (DNA Topo1). Besides as potent anti-cancer agents, CPT and its derivatives are now being explored as potential pesticides for insect control. In this study, we assessed their toxicity to an insect homolog, the Topo1 protein from beet armyworms (Spodoptera exigua Hübner), a worldwide pest of many important crops. The S. exigua Topo1 gene contains an ORF of 2790 base pairs that is predicted to encode a polypeptide of 930 amino acids. The deduced polypeptide exhibits polymorphism at residue sites V420, L530, A653 and T729 (numbered according to human Topo1) among insect species, which are predicted to confer sensitivity to CPT. The DNA relaxation activity of this protein was subsequently examined using a truncated form that contained the residues 337–930 and was expressed in bacteria BL21 cells. The purified protein retained the ability to relax double-stranded DNA and was susceptible to CPT and its derivative hydroxy-camptothecin (HCPT) in a dose-dependent manner. The same inhibitory effect was also found on the native Topo1 extracted from IOZCAS-Spex-II cells, a cell line established from beet armyworms. Additionally, CPT and HCPT treatment reduced the steady accumulation of Topo1 protein despite the increased mRNA expression in response to the treatment. Our studies provide information of the S. exigua Topo1 gene and its amino acid polymorphism in insects and uncover some clues about potential mechanisms of CPT toxicity against insect pests. These results also are useful for development of more effective Topo1-targeted CPT insecticides in the future.

A derivative of the natural compound kakuol affects DNA relaxation of topoisomerase IB inhibiting the cleavage reaction

Topoisomerases IB are anticancer and antimicrobial targets whose inhibition by several natural and synthetic compounds has been documented over the last three decades. Here we show that kakuol, a natural compound isolated from the rhizomes of Asarum sieboldii, and a derivative analogue are able to inhibit the DNA relaxation mediated by the human enzyme. The analogue is the most efficient one and the inhibitory effect is enhanced upon pre-incubation with the enzyme. Analysis of the different steps of the catalytic cycle indicates that the inhibition occurs at the cleavage level and does not prevent DNA binding. Molecular docking shows that the compound preferentially binds near the active site at the bottom of the catalytic residue Tyr723, providing an atomistic explanation for its inhibitory activity.

Structural and dynamical effects induced by the anticancer drug topotecan on the human topoisomerase I - DNA complex

BACKGROUND

Human topoisomerase I catalyzes the relaxation of DNA supercoils in fundamental cell processes like transcription, replication and chromosomal segregation. It is the only target of the camptothecin family of anticancer drugs. Among these, topotecan has been used to treat lung and ovarian carcinoma for several years. Camptothecins reversibly binds to the covalent intermediate DNA-enzyme, stabilizing the cleavable complex and reducing the religation rate. The stalled complex then collides with the progression of the replication fork, producing lethal double strand DNA breaks and eventually cell death.

METHODOLOGY/PRINCIPAL FINDINGS

Long lasting molecular dynamics simulations of the DNA-topoisomerase I binary complex and of the DNA-topoisomerase-topotecan ternary complex have been performed and compared. The conformational space sampled by the binary complex is reduced by the presence of the drug, as observed by principal component and cluster analyses. This conformational restraint is mainly due to the reduced flexibility of residues 633-643 (the region connecting the linker to the core domain) that causes an overall mobility loss in the ternary complex linker domain. During the simulation, DNA/drug stacking interactions are fully maintained, and hydrogen bonds are maintained with the enzyme. Topotecan keeps the catalytic residue Lys532 far from the DNA, making it unable to participate to the religation reaction. Arg364 is observed to interact with both the B and E rings of topotecan with two stable direct hydrogen bonds. An interesting constrain exerted by the protein on the geometrical arrangement of topotecan is also observed.

CONCLUSIONS/SIGNIFICANCE

Atomistic-scale understanding of topotecan interactions with the DNA-enzyme complex is fundamental to the explaining of its poisonous effect and of the drug resistance observed in several single residue topoisomerase mutants. We observed significant alterations due to topotecan in both short-range interactions and long-range protein domain communications.

Erybraedin C, a natural compound from the plant Bituminaria bituminosa, inhibits both the cleavage and religation activities of human topoisomerase I

The interaction of human topoisomerase I and erybraedin C, a pterocarpan purified from the plant Bituminaria bituminosa, that was shown to have an antitumour activity, was investigated through enzymatic activity assays and molecular docking procedures. Erybraedin C is able to inhibit both the cleavage and the religation steps of the enzyme reaction. In both cases, pre-incubation of the drug with the enzyme is required to produce a complete inhibition. Molecular docking simulations indicate that, when interacting with the enzyme alone, the preferential drug-binding site is localized in proximity to the active Tyr723 residue, with one of the two prenilic groups close to the active-site residues Arg488 and His632, essential for the catalytic reaction. When interacting with the cleavable complex, erybraedin C interacts with both the enzyme and DNA in a way similar to that found for topotecan. This is the first example of a natural compound able to act on both the cleavage and religation reaction of human topoisomerase I.

Evidence of the crucial role of the linker domain on the catalytic activity of human topoisomerase I by experimental and simulative characterization of the Lys681Ala mutant

The functional and structural-dynamical properties of the Lys681Ala mutation in the human topoisomerase IB linker domain have been investigated by catalytic assays and molecular dynamics simulation. The mutant is characterized by a comparable cleavage and a strongly reduced religation rate when compared to the wild type protein. The mutant also displays perturbed linker dynamics, as shown by analysis of the principal components of the motion, and a reduced electrostatic interaction with DNA. Inspection of the inter atomic distances in proximity of the active site shows that in the mutant the distance between the amino group of Lys532 side chain and the 5′ OH of the scissile phosphate is longer than the wild type enzyme, providing an atomic explanation for the reduced religation rate of the mutant. Taken together these results indicate the existence of a long range communication between the linker domain and the active site region and points out the crucial role of the linker in the modulation of the catalytic activity.

Conjugated eicosapentaenoic acid inhibits human topoisomerase IB with a mechanism different from camptothecin

Conjugated eicosapentaenoic acid (cEPA) has been found to have antitumor effects which has been ascribed to their ability to inhibit DNA topoisomerases and DNA polymerases. We here show that cEPA inhibits the catalytic activity of human topoisomerase I, but unlike camptothecin it does not stabilize the cleavable complex, indicating a different mechanism of action. cEPA inhibits topoisomerase by impeding the catalytic cleavage of the DNA substrate as demonstrated using specific oligonucleotide substrates, and prevents the stabilization of the cleavable complex by camptothecin. Preincubation of the inhibitor with the enzyme is required to obtain complete inhibition. Molecular docking simulations indicate that the preferred cEPA binding site is proximal to the active site with the carboxylic group strongly interacting with the positively charged K443 and K587. Taken together the results indicate that cEPA inhibitor does not prevent DNA binding but inhibits DNA cleavage, binding in a region close to the topoisomerase active site.