Granger causality based on vector time series and quaternion algebra with possible applications to molecular dynamics data analysis.

Causal analysis plays a significant role in physics, chemistry, and biology. Dynamics of complex (bio)molecular and nanosystems, from the microscopic to the macroscopic scale, are characterized by time-dependent vectors such as positions, forces, momenta, angular momenta, or torques. Identification and analysis of causal relationships between these time-dependent signals is an important problem in the multidimensional time-series analysis and is of great practical importance in describing the properties of such dynamical systems, and to understanding their functionality. For linear stochastic systems characterized by multidimensional scalar signals, Granger proposed a simple procedure to detect causal relationships, called Granger causality. In this study we extended this formalism to vector signals representing physical vector quantities. For this purpose, we used quaternion algebra, where vector signals are treated as time-dependent quaternions. The developed analytical model is based on the autoregressive formalism. This formalism (Q-MVAR) and its numerical implementation were validated using two simple dynamic models: a rigid body model represented by a benzenelike molecular fragment, interacting with a short-range harmonic potential with a wall, as well as a system of three model atomic balls moving inside a soft spherical surface and interacting with long range electrostatic forces. Although the motivation to these studies was the analysis of classical motions in complex (bio)molecular systems, described with a mechanical model and based on molecular dynamics (MD) simulations, in particular coarse-grained ones, it should be noted that the developed extended formalism can be applied to any system composed of many rigid elements that interact with arbitrary potentials and are characterized by complex internal motions. A description of the detailed procedure for calculating causality measures is provided in the Appendices of the Supplemental Material. This formalism and the prototype of its numerical implementation can be further developed and applied in many different fields of physical, natural, and engineering sciences.

DAMA: a method for computing multiple alignments of protein structures using local structure descriptors.

MOTIVATION: The well-known fact that protein structures are more conserved than their sequences forms the basis of several areas of computational structural biology. Methods based on the structure analysis provide more complete information on residue conservation in evolutionary processes. This is crucial for the determination of evolutionary relationships between proteins and for the identification of recurrent structural patterns present in biomolecules involved in similar functions. However, algorithmic structural alignment is much more difficult than multiple sequence alignment. This study is devoted to the development and applications of DAMA-a novel effective environment capable to compute and analyze multiple structure alignments. RESULTS: DAMA is based on local structural similarities, using local 3D structure descriptors and thus accounts for nearest-neighbor molecular environments of aligned residues. It is constrained neither by protein topology nor by its global structure. DAMA is an extension of our previous study (DEDAL) which demonstrated the applicability of local descriptors to pairwise alignment problems. Since the multiple alignment problem is NP-complete, an effective heuristic approach has been developed without imposing any artificial constraints. The alignment algorithm searches for the largest, consistent ensemble of similar descriptors. The new method is capable to capture most of the biologically significant similarities present in canonical test sets and is discriminatory enough to prevent the emergence of larger, but meaningless, solutions. Tests performed on the test sets, including protein kinases, demonstrate DAMA's capability of identifying equivalent residues, which should be very useful in discovering the biological nature of proteins similarity. Performance profiles show the advantage of DAMA over other methods, in particular when using a strict similarity measure QC, which is the ratio of correctly aligned columns, and when applying the methods to more difficult cases. AVAILABILITY AND IMPLEMENTATION: DAMA is available online at http://dworkowa.imdik.pan.pl/EP/DAMA. Linux binaries of the software are available upon request. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Exploring Covalent Docking Mechanisms of Boron-Based Inhibitors to Class A, C and D ß-Lactamases Using Time-dependent Hybrid QM/MM Simulations.

Recently, molecular covalent docking has been extensively developed to design new classes of inhibitors that form chemical bonds with their biological targets. This strategy for the design of such inhibitors, in particular boron-based inhibitors, holds great promise for the vast family of ß-lactamases produced, inter alia, by Gram-negative antibiotic-resistant bacteria. However, the description of covalent docking processes requires a quantum-mechanical approach, and so far, only a few studies of this type have been presented. This study accurately describes the covalent docking process between two model inhibitors - representing two large families of inhibitors based on boronic-acid and bicyclic boronate scaffolds, and three ß-lactamases which belong to the A, C, and D classes. Molecular fragments containing boron can be converted from a neutral, trigonal, planar state with sp2 hybridization to the anionic, tetrahedral sp3 state in a process sometimes referred to as morphing. This study applies multi-scale modeling methods, in particular, the hybrid QM/MM approach which has predictive power reaching well beyond conventional molecular modeling. Time-dependent QM/MM simulations indicated several structural changes and geometric preferences, ultimately leading to covalent docking processes. With current computing technologies, this approach is not computationally expensive, can be used in standard molecular modeling and molecular design works, and can effectively support experimental research which should allow for a detailed understanding of complex processes important to molecular medicine. In particular, it can support the rational design of covalent boron-based inhibitors for ß-lactamases as well as for many other enzyme systems of clinical relevance, including SARS-CoV-2 proteins.

Structural characterisation of inhibitory and non-inhibitory MMP-9-TIMP-1 complexes and implications for regulatory mechanisms of MMP-9.

MMP-9 plays a number of important physiological functions but is also responsible for many pathological processes, including cancer invasion, metastasis, and angiogenesis. It is, therefore, crucial to understand its enzymatic activity, including activation and inhibition mechanisms. This enzyme may also be partially involved in the "cytokine storm" that is characteristic of COVID-19 disease (SARS-CoV-2), as well as in the molecular mechanisms responsible for lung fibrosis. Due to the variety of processing pathways involving MMP-9 in biological systems and its uniqueness due to the O-glycosylated domain (OGD) and fibronectin-like (FBN) domain, specific interactions with its natural TIMP-1 inhibitor should be carefully studied, because they differ significantly from other homologous systems. In particular, earlier experimental studies have indicated that the newly characterised circular form of a proMMP-9 homotrimer exhibits stronger binding properties to TIMP-1 compared to its monomeric form. However, molecular structures of the complexes and the binding mechanisms remain unknown. The purpose of this study is to fill in the gaps in knowledge. Molecular modelling methods are applied to build the inhibitory and non-inhibitory MMP-9-TIMP-1 complexes, which allows for a detailed description of these structures and should allow for a better understanding of the regulatory processes in which MMP-9 is involved.

A novel formulation of an approximate valence bond model (AVB2) and its application to the tautomeric forms of porphyrin and porphycene.

This study presents a novel formulation of the approximate valence bond method, which can be applied as a very fast generator of the molecular potential energy function. The AVB2 model was formulated and parameterized for porphyrin and porphycene using results of quantum mechanical computations at the B3LYP/6-31G (d, p) level. The DFT potential energy, its gradients, and the Hessian-matrix elements, as well as effective atomic charges at local energy minima and transition states, were used for the parameterization of the AVB2 Hamiltonian matrix. The AVB2 method, and in particular its anharmonic version, very well reproduce the potential energy maps for all representative geometries of the studied systems, including harmonic frequencies, and possible proton translocations. For validation of the method, we performed molecular dynamics simulations for isolated molecules accounting for internal double proton transfer processes, which are strongly correlated with changes of the electronic charge density. The simulated power spectra were compared with the experimental infrared spectra. More precise simulations of IR spectra at the classical and quantum dynamics levels, as well as extensions of the AVB2 parameterization to electronic excited states, are the subject of further research.

Surface markers of cancer stem-like cells of ovarian cancer and their clinical relevance.

Cancer stem-like cells (CSLCs) are defined as cancer cells with stem cell characteristics. Although CSLCs constitute no more than a few percent of the tumor mass, they play important roles in cancer chemo-resistance, metastasis and disease recurrence. Ovarian cancer (OC) is considered the most aggressive gynecological malignancy in which the role of CSLCs is of major significance, although it remains to be specified. The studies describing ovarian CSLC phenotype vary in the definition of the molecular pattern of expression of the main markers such as CD133, CD44, CD117, and CD24. Stem-like features of OC have been shown to correlate with the clinical course of the disease and permit diagnosis, prognosis and treatment outcome to be improved. Identification of CSLC markers could provide hallmarks which, related to the chemo-resistance of the disease, will facilitate treatment selection. This review describes recent advances in research on stem-like cell status in OC, mainly focusing on surface markers of CSLCs and their clinical relevance.

Synthesis, structural and antimicrobial studies of type II topoisomerase-targeted copper(II) complexes of 1,3-disubstituted thiourea ligands.

A series of Cu(II) complexes of 3-(trifluoromethyl)phenylthiourea derivatives was synthesized. Their structural properties were investigated by spectroscopic techniques (infrared and electron paramagnetic resonance), as well as molecular modeling. All studied coordination compounds are mononuclear complexes containing two chelating ligands bonded to the metal cation via S and deprotonated N atoms. The new chelates were evaluated for their antimicrobial potency. The complex of 1-(3,4-dichlorophenyl)-3-[3-(trifluoromethyl)phenyl]thiourea (3) presented the highest activity against Gram-positive pathogens, even stronger than the activity of its non-complexed counterpart and the reference drug. The compound also prevented the biofilm formation of methicillin-resistant and standard strains of staphylococcal cocci. The title derivatives were found to be effective inhibitors of DNA gyrase and topoisomerase IV isolated from Staphylococcus aureus. The binding modes of the ligand L3 with DNA gyrase and topoisomerase IV were presented.

Estimation of Posturographic Trajectory Using k-Nearest Neighbors Classifier in Patients with Rheumatoid Arthritis and Osteoarthritis.

Rheumatoid arthritis (RA) and osteoarthritis (OA) are common rheumatic diseases and account for a significant percentage of disability. Posturography is a method that assesses postural stability and quantitatively evaluates postural sways. The objective of this study was to estimate posturographic trajectories applying pattern recognition algorithms. To this end, k-nearest neighbors (k-NN) classifier was used to differentiate between healthy subjects and patients with OA and RA. The following parameters of trajectories were computed: radius of sways, developed area, total length, and two directional components of sways: length of left-right and forward-backward motions. Posturographic tests were applied with eyes open and closed, and with biofeedback control. We found that in RA, the radius of sways, the trajectory area, and the biofeedback coordination were related to the patients' condition. The trajectory dynamics in OA patients were smaller compared to those in RA patients. The smallest misclassification errors were observed after feature selection in the biofeedback test compared with the eyes open and closed tests. We conclude that the estimation of posturographic trajectory with k-NN classifier could be helpful in monitoring the condition of RA patients.

Inferring Molecular Processes Heterogeneity from Transcriptional Data.

RNA microarrays and RNA-seq are nowadays standard technologies to study the transcriptional activity of cells. Most studies focus on tracking transcriptional changes caused by specific experimental conditions. Information referring to genes up- and downregulation is evaluated analyzing the behaviour of relatively large population of cells by averaging its properties. However, even assuming perfect sample homogeneity, different subpopulations of cells can exhibit diverse transcriptomic profiles, as they may follow different regulatory/signaling pathways. The purpose of this study is to provide a novel methodological scheme to account for possible internal, functional heterogeneity in homogeneous cell lines, including cancer ones. We propose a novel computational method to infer the proportion between subpopulations of cells that manifest various functional behaviour in a given sample. Our method was validated using two datasets from RNA microarray experiments. Both experiments aimed to examine cell viability in specific experimental conditions. The presented methodology can be easily extended to RNA-seq data as well as other molecular processes. Moreover, it complements standard tools to indicate most important networks from transcriptomic data and in particular could be useful in the analysis of cancer cell lines affected by biologically active compounds or drugs.

Cinnamic Acid Derivatives as Inhibitors of Oncogenic Protein Kinases--Structure, Mechanisms and Biomedical Effects.

Cinnamic acid belongs to phenolic-acid class of polyphenols, one of the most abundant plant secondary metabolites. These substances are widely studied because of plethora of their biological activities. In particular, their inhibition of protein kinases contributes to the pleiotropic effects in the cell. Protein kinases are essential in controlling cell signaling networks. Selective targeting of oncogenic protein kinases increases clinical anticancer efficacy. Cinnamic acid and related compounds have inspired researchers in the design of numerous synthetic and semisynthetic inhibitors of oncogenic protein kinases for the past three decades. Interest in cinnamoyl-scaffold-containing compounds revived in recent years, which was stimulated by modern drug design and discovery methodologies such as in vitro and in silico HTS. This review presents cinnamic acid derivatives and analogs for which direct inhibition of protein kinases was identified. We also summarize significance of the above protein kinase families - validated or promising targets for anticancer therapies. The inhibition mode may vary from ATP-competitive, through bisubstrate-competitive and mixedcompetitive, to non-competitive one. Kinase selectivity is often correlated with subtle chemical modifications, and may also be steered by an additional non-cinnamoyl fragment of the inhibitor. Specific cinnamic acid congeners may synergize their effects in the cell by a wider range of activities, like suppression of additional enzymes, e.g. deubiquitinases, influencing the same signaling pathways (e.g. JAK2/STAT). Cinnamic acid, due to its biological and physicochemical properties, provides nature-inspired ideas leading to novel inhibitors of oncogenic protein kinases and related enzymes, capable to target a variety of cancer cells.

Toward the identification of molecular cogs.

Computer simulations of molecular systems allow determination of microscopic interactions between individual atoms or groups of atoms, as well as studies of intramolecular motions. Nevertheless, description of structural transformations at the mezoscopic level and identification of causal relations associated with these transformations is very difficult. Structural and functional properties are related to free energy changes. Therefore, to better understand structural and functional properties of molecular systems, it is required to deepen our knowledge of free energy contributions arising from molecular subsystems in the course of structural transformations. The method presented in this work quantifies the energetic contribution of each pair of atoms to the total free energy change along a given collective variable. Next, with the help of a genetic clustering algorithm, the method proposes a division of the system into two groups of atoms referred to as molecular cogs. Atoms which cooperate to push the system forward along a collective variable are referred to as forward cogs, and those which work in the opposite direction as reverse cogs. The procedure was tested on several small molecules for which the genetic clustering algorithm successfully found optimal partitionings into molecular cogs. The primary result of the method is a plot depicting the energetic contributions of the identified molecular cogs to the total Potential of Mean Force (PMF) change. Case-studies presented in this work should help better understand the implications of our approach, and were intended to pave the way to a future, publicly available implementation.

ResiCon: a method for the identification of dynamic domains, hinges and interfacial regions in proteins.

MOTIVATION: Structure of most proteins is flexible. Identification and analysis of intramolecular motions is a complex problem. Breaking a structure into relatively rigid parts, the so-called dynamic domains, may help comprehend the complexity of protein's mobility. We propose a new approach called ResiCon (Residue Contacts analysis), which performs this task by applying a data-mining analysis of an ensemble of protein configurations and recognizes dynamic domains, hinges and interfacial regions, by considering contacts between residues. RESULTS: Dynamic domains found by ResiCon are more compact than those identified by two other popular methods: PiSQRD and GeoStaS. The current analysis was carried out using a known reference set of 30 NMR protein structures, as well as molecular dynamics simulation data of flap opening events in HIV-1 protease. The more detailed analysis of HIV-1 protease dataset shows that ResiCon identified dynamic domains involved in structural changes of functional importance. AVAILABILITY AND IMPLEMENTATION: The ResiCon server is available at URL: http://dworkowa.imdik.pan.pl/EP/ResiCon. CONTACT: pawel@bioexploratorium.pl SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

WeBIAS: a web server for publishing bioinformatics applications.

BACKGROUND: One of the requirements for a successful scientific tool is its availability. Developing a functional web service, however, is usually considered a mundane and ungratifying task, and quite often neglected. When publishing bioinformatic applications, such attitude puts additional burden on the reviewers who have to cope with poorly designed interfaces in order to assess quality of presented methods, as well as impairs actual usefulness to the scientific community at large. RESULTS: In this note we present WeBIAS-a simple, self-contained solution to make command-line programs accessible through web forms. It comprises a web portal capable of serving several applications and backend schedulers which carry out computations. The server handles user registration and authentication, stores queries and results, and provides a convenient administrator interface. WeBIAS is implemented in Python and available under GNU Affero General Public License. It has been developed and tested on GNU/Linux compatible platforms covering a vast majority of operational WWW servers. Since it is written in pure Python, it should be easy to deploy also on all other platforms supporting Python (e.g. Windows, Mac OS X). Documentation and source code, as well as a demonstration site are available at http://bioinfo.imdik.pan.pl/webias . CONCLUSIONS: WeBIAS has been designed specifically with ease of installation and deployment of services in mind. Setting up a simple application requires minimal effort, yet it is possible to create visually appealing, feature-rich interfaces for query submission and presentation of results.

Analogs of cinnamic acid benzyl amide as nonclassical inhibitors of activated JAK2 kinase.

Scaffold-based analogs of cinnamic acid benzyl amide (CABA) exhibit pleiotropic effects in cancer cells, and their exact molecular mechanism of action is under investigation. The present study is part of our systemic analysis of interactions of CABA analogs with their molecular targets. These compounds were shown to inhibit Janus kinase 2 (JAK2)/signal transducer and activator of transcription 3 (STAT3) and JAK2/signal transducer and activator of transcription 5 (STAT5) signaling and thus are attractive scaffolds for anticancer drug design. To identify the potential mechanisms of action of this class of compounds, direct interactions of the selected CABA analogs with JAK2 kinase were examined. Inhibition of JAK2 enzymatic activity was assessed, and molecular modeling studies of selected compounds-(E)-2-cyano-N-[(S)-1-phenylethyl]-3-(pyridin-2-yl)acrylamide (WP1065), (E)-2-cyano-N-[(S)-1-phenylbutyl]- 3-(3-bromopyridin-2-yl)acrylamide (WP1130), and (E)-2-cyano-N-[(S)-1,4-diphenylbutyl]-3-(3-bromopyridin-2-yl)acrylamide (WP1702)-in the JAK2 kinase domain were used to support interpretation of the experimental data. Our results indicated that the tested CABA analogs are nonclassical inhibitors of activated (phosphorylated) JAK2, although markedly weaker than clinically tested ATP-competitive JAK2 inhibitors. Relatively small structural changes in the studied compounds affected interactions with JAK2, and their mode of action ranged from allosteric-noncompetitive to bisubstratecompetitive. These results demonstrated that direct inhibition of JAK2 enzymatic activity by the WP1065 (half-maximal inhibitory concentration [IC50] = 14.8 µM), WP1130 (IC50 = 3.8 µM), and WP1702 (IC50 = 2.9 µM) potentially contributes, albeit minimally, to suppression of the JAK2/STAT signaling pathways in cancer cells and that additional specific structural modifications may amplify JAK2-inhibitory effects.

Development of novel molecular probes of the Rio1 atypical protein kinase.

The RIO kinases are essential protein factors required for the synthesis of new ribosomes in eukaryotes. Conserved in archaeal organisms as well, RIO kinases are among the most ancient of protein kinases. Their exact molecular mechanisms are under investigation and progress of this research would be significantly improved with the availability of suitable molecular probes that selectively block RIO kinases. RIO kinases contain a canonical eukaryotic protein kinase fold, but also display several unusual structural features that potentially create opportunity for the design of selective inhibitors. In an attempt to identify structural leads to target the RIO kinases, a series of pyridine caffeic acid benzyl amides (CABA) were tested for their ability to inhibit the autophosphorylation activity of Archeaoglobus fulgidus Rio1 (AfRio1). Screening of a small library of CABA molecules resulted in the identification of four compounds that measurably inhibited AfRio1 activity. Additional biochemical characterization of binding and inhibition activity of these compounds demonstrated an ATP competitive inhibition mode, and allowed identification of the functional groups that result in the highest binding affinity. In addition, docking of the compound to the structure of Rio1 and determination of the X-ray crystal structure of a model compound (WP1086) containing the desired functional groups allowed detailed analysis of the interactions between these compounds and the enzyme. Furthermore, the X-ray crystal structure demonstrated that these compounds stabilize an inactive form of the enzyme. Taken together, these results provide an important step in identification of a scaffold for the design of selective molecular probes to study molecular mechanisms of Rio1 kinases in vitro and in vivo. In addition, it provides a rationale for the future design of potent inhibitors with drug-like properties targeting an inactive form of the enzyme. This article is part of a Special Issue entitled: Inhibitors of Protein Kinases (2012).

Thermal fluctuations and infrared spectra of the formamide-formamidine complex.

Combined effects of hydrogen bonding and thermal fluctuations on the structure and infrared spectra of the formamide-formamidine dimer, FM···FI, are studied using ab initio molecular dynamics simulations. The equilibrium structure of the dimer is stabilized by two hydrogen bonds that form a pattern reminiscent of that found in the adenine-thymine base pair. The structure of the hydrogen bonds at 300 K is subject to large fluctuations, with the hydrogen atom being tightly bound to the donor in the covalent bonding scenario. The hydrogen bond acceptor has a tendency to detach farther away from the D-H pair, approaching the dimer dissociation limit. Moreover, the N-H···O hydrogen bond breaks occasionally, thus giving rise to an "open" structure of the dimer, while the N-H···N bond stays largely intact at this temperature. Thermal fluctuations result in the minor red shifts of the monomer vibrational frequencies indicative of the anharmonicity of the potential energy surface. In contrast, the IR frequencies of the two symmetric NH(2) vibrational modes of the FM and FI monomers are shifted substantially toward the red upon hydrogen bond formation in the FM···FI dimer. Dynamical effects studied here are relevant, in particular, to the hydrogen bonding of nucleic acids at finite temperatures.

Novel small molecular inhibitors disrupt the JAK/STAT3 and FAK signaling pathways and exhibit a potent antitumor activity in glioma cells.

JAK (Janus kinase)/STAT (signal transducers and activators of transcription) signaling is involved in the regulation of cell growth, differentiation and apoptosis. Constitutive activation of STATs, in particular STAT3, is observed in a large number of human tumors, including gliomas and may contribute to oncogenesis by stimulating cell proliferation and preventing apoptosis, thus it emerges as a promising target for anti-cancer therapy. To investigate the therapeutic potential of blocking STAT3 in glioma cells a set of small synthetic molecules - caffeic acid derivatives, structurally related to AG490 was screened for its ability to inhibit STAT3. Inhibitor 2 (E)-2-cyano-N-[(S)-1-phenylethyl]-3-(pyridin-2-yl)acrylamide was the most effective in inhibition of JAK/STAT3 signaling and at doses ≥ 25 µM significantly reduced the level of phosphorylated JAK1, JAK2 and STAT3 (at Tyr705) and downregulated the expression of known STAT3 targets. In treated cells we observed rapid detachment and rounding of cells associated with reduction of focal adhesion kinase phosphorylation and activity, followed by upregulation of phosphorylated p38, JNK and ERK1/2 levels. Accumulation of cells with fragmented DNA, increases of the cleaved caspase 3 and fragmented PARP levels were detected 24 h after the treatment suggesting ongoing apoptotic cell death. Three human malignant glioblastoma cell lines defective in tumor suppressors TP53 and/or PTEN were susceptible to inhibitor 2 that induced the programmed cell death. Global gene expression profiling revealed modulation of numerous genes in cells treated with inhibitor 2 revealing novel, potential JAK/STAT targets. Our study demonstrates that suitably modified caffeic acid molecules exhibit significant cytotoxic potential toward glioma cells.

Vibrations and hydrogen bonding in porphycene.

Combined use of IR, Raman, neutron scattering and fluorescence measurements for porphycene isolated in helium nanodroplets, supersonic jet and cryogenic matrices, as well as for solid and liquid solutions, resulted in the assignments of almost all of 108 fundamental vibrations. The puzzling feature of porphycene is the apparent lack of the N-H stretching band in the IR spectrum, predicted to be the strongest of all bands by standard harmonic calculations. Theoretical modeling of the IR spectra, based on ab initio molecular dynamics simulations, reveals that the N-H stretching mode should appear as an extremely broad band in the 2250-3000 cm(-1) region. Coupling of the N-H stretching vibration to other modes is discussed in the context of multidimensional character of intramolecular double hydrogen transfer in porphycene. The analysis can be generalized to other strongly hydrogen-bonded systems.

Activities of topoisomerase I in its complex with SRSF1.

Human DNA topoisomerase I (topo I) catalyzes DNA relaxation and phosphorylates SRSF1. Whereas the structure of topo I complexed with DNA has been resolved, the structure of topo I in the complex with SRSF1 and structural determinants of topo I activities in this complex are not known. The main obstacle to resolving the structure is a contribution of unfolded domains of topo I and SRSF1 in formation of the complex. To overcome this difficulty, we employed a three-step strategy: identifying the interaction regions, modeling the complex, and validating the model with biochemical methods. The binding sites in both topo I and SRSF1 are localized in the structured regions as well as in the unfolded domains. One observes cooperation between the binding sites in topo I but not in SRSF1. Our results indicate two features of the unfolded RS domain of SRSF1 containing phosphorylated residues that are critical for the kinase activity of topo I: its spatial arrangement relative to topo I and the organization of its sequence. The efficiency of phosphorylation of SRSF1 depends on the length and flexibility of the spacer between the two RRM domains that uniquely determine an arrangement of the RS domain relative to topo I. The spacer also influences inhibition of DNA nicking, a prerequisite for DNA relaxation. To be phosphorylated, the RS domain has to include a short sequence recognized by topo I. A lack of this sequence in the mutants of SRSF1 or its spatial inaccessibility in SRSF9 makes them inadequate as topo I/kinase substrates.