Corrigendum to "Recognition motifs for importin 4 [(L)PPRS(G/P)P] and importin 5 [KP(K/Y)LV] binding, identified by bio-informatic simulation and experimental in vitro validation" [Comput Struct Biotechnol J 20 (2022) 5952-5961].

[This corrects the article DOI: 10.1016/j.csbj.2022.10.015.].

Recognition motifs for importin 4 [(L)PPRS(G/P)P] and importin 5 [KP(K/Y)LV] binding, identified by bio-informatic simulation and experimental in vitro validation.

Nuclear translocation of large proteins is mediated through karyopherins, carrier proteins recognizing specific motifs of cargo proteins, known as nuclear localization signals (NLS). However, only few NLS signals have been reported until now. In the present work, NLS signals for Importins 4 and 5 were identified through an unsupervised in silico approach, followed by experimental in vitro validation. The sequences LPPRS(G/P)P and KP(K/Y)LV were identified and are proposed as recognition motifs for Importins 4 and 5 binding, respectively. They are involved in the trafficking of important proteins into the nucleus. These sequences were validated in the breast cancer cell line T47D, which expresses both Importins 4 and 5. Elucidating the complex relationships of the nuclear transporters and their cargo proteins is very important in better understanding the mechanism of nuclear transport of proteins and laying the foundation for the development of novel therapeutics, targeting specific importins.

FEPrepare: A Web-Based Tool for Automating the Setup of Relative Binding Free Energy Calculations.

Relative binding free energy calculations in drug design are becoming a useful tool in facilitating lead binding affinity optimization in a cost- and time-efficient manner. However, they have been limited by technical challenges such as the manual creation of large numbers of input files to set up, run, and analyze free energy simulations. In this Application Note, we describe FEPrepare, a novel web-based tool, which automates the setup procedure for relative binding FEP calculations for the dual-topology scheme of NAMD, one of the major MD engines, using OPLS-AA force field topology and parameter files. FEPrepare provides the user with all necessary files needed to run a FEP/MD simulation with NAMD. FEPrepare can be accessed and used at https://feprepare.vi-seem.eu/.

The sequence [EKRKI(E/R)(K/L/R/S/T)] is a nuclear localization signal for importin 7 binding (NLS7).

BACKGROUND: Nuclear translocation of large proteins is mediated through specific protein carriers, collectively named karyopherins (importins, exportins and adaptor proteins). Cargo proteins are recognized by importins through specific motifs, known as nuclear localization signals (NLS). However, only the NLS recognized by importin α and transportin (M9 NLS) have been identified so far METHODS: An unsupervised in silico approach was used, followed by experimental validation. RESULTS: We identified the sequence EKRKI(E/R)(K/L/R/S/T) as an NLS signal for importin 7 recognition. This sequence was validated in the breast cancer cell line T47D, which expresses importin 7. Finally, we verified that importin 7-mediated nuclear protein transport is affected by cargo protein phosphorylation. CONCLUSIONS: The NLS sequence for importin 7 was identified and we propose this approach as an identification method of novel specific NLS sequences for ß-karyopherin family members. GENERAL SIGNIFICANCE: Elucidating the complex relationships of the nuclear transporters and their cargo proteins may help in laying the foundation for the development of novel therapeutics, targeting specific importins, with an immediate translational impact.

ChemBioServer 2.0: an advanced web server for filtering, clustering and networking of chemical compounds facilitating both drug discovery and repurposing.

SUMMARY: ChemBioServer 2.0 is the advanced sequel of a web server for filtering, clustering and networking of chemical compound libraries facilitating both drug discovery and repurposing. It provides researchers the ability to (i) browse and visualize compounds along with their physicochemical and toxicity properties, (ii) perform property-based filtering of compounds, (iii) explore compound libraries for lead optimization based on perfect match substructure search, (iv) re-rank virtual screening results to achieve selectivity for a protein of interest against different protein members of the same family, selecting only those compounds that score high for the protein of interest, (v) perform clustering among the compounds based on their physicochemical properties providing representative compounds for each cluster, (vi) construct and visualize a structural similarity network of compounds providing a set of network analysis metrics, (vii) combine a given set of compounds with a reference set of compounds into a single structural similarity network providing the opportunity to infer drug repurposing due to transitivity, (viii) remove compounds from a network based on their similarity with unwanted substances (e.g. failed drugs) and (ix) build custom compound mining pipelines. AVAILABILITY AND IMPLEMENTATION: http://chembioserver.vi-seem.eu.

NanoCrystal: A Web-Based Crystallographic Tool for the Construction of Nanoparticles Based on Their Crystal Habit.

Modeling of nanoparticles is an essential first step to assess their capacities for different uses such as energy storage and drug delivery. However, creating an initial starting conformation for modeling and simulation is tedious because every crystalline material grows with a different crystal habit. In this application note, we describe NanoCrystal, a novel web-based crystallographic tool that creates nanoparticle models from any crystal structure guided by their preferred equilibrium shape under standard conditions according to the Wulff morphology (crystal habit). Users can upload a cif file, define the Miller indices and their corresponding minimum surface energies according to the Wulff construction of a particular crystal, and specify the size of the nanocrystal. As a result, the nanoparticle is constructed and visualized, and the coordinates of the atoms are output to the user. NanoCrystal can be accessed at http://nanocrystal.vi-seem.edu/ .

Using physics-based pose predictions and free energy perturbation calculations to predict binding poses and relative binding affinities for FXR ligands in the D3R Grand Challenge 2.

Computer-aided drug design has become an integral part of drug discovery and development in the pharmaceutical and biotechnology industry, and is nowadays extensively used in the lead identification and lead optimization phases. The drug design data resource (D3R) organizes challenges against blinded experimental data to prospectively test computational methodologies as an opportunity for improved methods and algorithms to emerge. We participated in Grand Challenge 2 to predict the crystallographic poses of 36 Farnesoid X Receptor (FXR)-bound ligands and the relative binding affinities for two designated subsets of 18 and 15 FXR-bound ligands. Here, we present our methodology for pose and affinity predictions and its evaluation after the release of the experimental data. For predicting the crystallographic poses, we used docking and physics-based pose prediction methods guided by the binding poses of native ligands. For FXR ligands with known chemotypes in the PDB, we accurately predicted their binding modes, while for those with unknown chemotypes the predictions were more challenging. Our group ranked #1st (based on the median RMSD) out of 46 groups, which submitted complete entries for the binding pose prediction challenge. For the relative binding affinity prediction challenge, we performed free energy perturbation (FEP) calculations coupled with molecular dynamics (MD) simulations. FEP/MD calculations displayed a high success rate in identifying compounds with better or worse binding affinity than the reference (parent) compound. Our studies suggest that when ligands with chemical precedent are available in the literature, binding pose predictions using docking and physics-based methods are reliable; however, predictions are challenging for ligands with completely unknown chemotypes. We also show that FEP/MD calculations hold predictive value and can nowadays be used in a high throughput mode in a lead optimization project provided that crystal structures of sufficiently high quality are available.

DIANA-miRGen v3.0: accurate characterization of microRNA promoters and their regulators.

microRNAs (miRNAs) are small non-coding RNAs that actively fine-tune gene expression. The accurate characterization of the mechanisms underlying miRNA transcription regulation will further expand our knowledge regarding their implication in homeostatic and pathobiological networks. Aim of DIANA-miRGen v3.0 (http://www.microrna.gr/mirgen) is to provide for the first time accurate cell-line-specific miRNA gene transcription start sites (TSSs), coupled with genome-wide maps of transcription factor (TF) binding sites in order to unveil the mechanisms of miRNA transcription regulation. To this end, more than 7.3 billion RNA-, ChIP- and DNase-Seq next generation sequencing reads were analyzed/assembled and combined with state-of-the-art miRNA TSS prediction and TF binding site identification algorithms. The new database schema and web interface facilitates user interaction, provides advanced queries and innate connection with other DIANA resources for miRNA target identification and pathway analysis. The database currently supports 276 miRNA TSSs that correspond to 428 precursors and >19M binding sites of 202 TFs on a genome-wide scale in nine cell-lines and six tissues of Homo sapiens and Mus musculus.

Biological and computational evaluation of resveratrol inhibitors against Alzheimer's disease.

It has been reported that beta amyloid induces production of radical oxygen species and oxidative stress in neuronal cells, which in turn upregulates ß-secretase (BACE-1) expression and beta amyloid levels, thereby propagating oxidative stress and increasing neuronal injury. A series of resveratrol derivatives, known to be inhibitors of oxidative stress-induced neuronal cell death (oxytosis) were biologically evaluated against BACE-1 using homogeneous time-resolved fluorescence (TRF) assay. Correlation between oxytosis inhibitory and BACE-1 inhibitory activity of resveratrol derivatives was statistically significant, supporting the notion that BACE-1 may act as pivotal mediator of neuronal cell oxytosis. Four of the biologically evaluated resveratrol analogs demonstrated considerably higher activity than resveratrol in either assay. The discovery of some "hits" led us to initiate detailed docking studies associated with Molecular Dynamics in order to provide a plausible explanation for the experimental results and understand their molecular basis of action.

Hydrogen bonding in 1-butyl- and 1-ethyl-3-methylimidazolium chloride ionic liquids.

A detailed investigation of hydrogen bonding in the pure ionic liquids [C4C1im]Cl and [C2C1im]Cl has been carried out using primarily molecular dynamics techniques. Analyses of the individual atom-atom pair radial distribution functions, and in particular those for C···Cl(-), have revealed that hydrogen bonding to the first methylene or methyl units of the substituent groups is important. Multiple geometric criteria for defining a hydrogen bond have been applied, and in particular the choice of the cutoff angle has been carefully examined. The interpretation of hydrogen bonding within these ionic liquids is highly angle dependent, and justification is provided for why it may be appropriate to employ a wider angle criteria than the 30° used for water or alcohol systems. The different types of hydrogen bond formed are characterized, and "top" conformations where the Cl anion resides above (or below) the imidazolium ring are investigated. The number of hydrogen bonds undertaken by each hydrogen atom (and the chloride anion) is quantified, and the propensity to form zero, one, or two hydrogen bonds is established. The effects of an increase in temperature on the static hydrogen bonding are also briefly examined.

Solvation structure and dynamics of cis- and trans-1,2 dichloroethene isomers in supercritical carbon dioxide. A molecular dynamics simulation study.

Molecular dynamics simulation techniques have been employed to investigate the solvation structure and dynamics in dilute mixtures of cis- and trans-1,2-dichloroethene in supercritical carbon dioxide. The calculations were performed for state points along a near-critical isotherm (1.02 T(c)) over a wide range of densities, using new developed optimized potential models for both isomers. The similarities and differences in the solvation structures around each isomer have been presented and discussed. The local density augmentation and enhancement factors of CO(2) around the isomers have been found significantly larger than the corresponding values for pure supercritical CO(2). The dynamic local density reorganization has been investigated and related to previously proposed relaxation mechanisms. The density dependence of the calculated self-diffusion coefficients has revealed the existence of a plateau in the region of 0.7-1.1 ρ(c), where the local density augmentation exhibits the maximum value. The reorientational dynamics of the C═C bond vector have been also studied, exhibiting significant differences between the two isomers in the case of the second-order Legendre time correlation functions.

Conformational and solvation studies via computer simulation of the novel large scale diastereoselectively synthesized phosphinic MMP inhibitor RXP03 diluted in selected solvents.

Structure-activity relationship studies, regarding the influence of side chains of phosphinic pseudotripeptidic inhibitors on matrix metalloproteinases (MMPs), provided potent and selective inhibitors for this family of structurally and functionally related proteases. Among them, phosphinic pseudopeptide CbzPhepsi[P(O)(OH)CH(2)] phenylpropyl TrpNH(2), known as RXP03, has been extensively used for in vivo and in vitro studies so far. The large quantities of RXP03 required for in vivo studies, as well as the necessity for diastereoisomeric purity, motivated us to further explore and develop an efficient synthetic methodology, which allows separation of the four diastereoisomers of RXP03 based on the astonishing observed differences in solubility of the four isomers in various solvents. This fact prompted us to examine theoretically the conformational differences of these four isomers via computer simulations in the solvents used experimentally. Given the fact that the four examined diastereoisomeric forms of the phosphinic peptides exhibit different behavior in terms of potency and selectivity profiles toward zinc-metalloproteases, this theoretical study provides valuable information on the conformation of phosphinic inhibitors and therefore improves the design and synthesis of active structures. The differences in solubility of RXP03 diastereoisomers in the used solvents were examined in terms of intra- and intermolecular structure. It is found that the different solubility of the RRS and RSS diastereoisomers in EtOH is a result of the different number of hydrogen bonds formed by each isomer with EtOH molecules. In the case of SRS and SSS in Et(2)O, their different solubility might be attributed to the different intramolecular hydrogen bonds formed on these diastereoisomers.

The effect of intermolecular interactions on local density inhomogeneities and related dynamics in pure supercritical fluids. A comparative molecular dynamics simulation study.

The effect of intermolecular interactions of different strength on the local density inhomogeneities in pure supercritical fluids (scfs), with different intramolecular structure, was investigated by employing molecular dynamics (MD) simulation techniques. The simulations were performed at state points along an isotherm close to the critical temperature of each system (T(r) = T/T(c) = 1.03). The molecular fluids under study have been chosen on the basis of the electrostatic character of their intermolecular interactions as follows: monatomic, dipolar and hydrogen bonding (HB), quadrupolar, and octupolar. In the case of dipolar scfs, their HB nature when present was systematically explored and related to the behavior of the created local density inhomogeneities at all densities. The results obtained reveal strong influence of the dipolar and HB interactions of the investigated systems upon the local density augmentation. We found that this effect is fairly larger in the case of the dipolar and HB fluids (H2O, CH3OH, and NH3) compared to those for the non-dipolar ones (Xe, CH4, CO2, and N2). In the case of sc CO2, the dependence of the local density augmentation on the bulk density is in agreement with available experimental data as also reported previously. The estimated average number of hydrogen bonds per molecule (nHB) in these HB fluids shows an analogue nonlinear trend compared to the behavior of the average coordination numbers Nco(rho) of a particle with bulk density. The local density dynamics of the first and second solvation shell of each fluid were further analyzed and related to our previously proposed [Skarmoutsos, I.; Samios, J. J. Chem. Phys. 2007, 126, 044503] different time-scale relaxation mechanisms. Finally, the effect of the different strength of the molecular interactions corresponding to these fluids upon the local density dynamics has also been revealed in the behavior of the predicted appropriate time correlation functions and their corresponding correlation times.

Investigation of the local composition enhancement and related dynamics in supercritical CO2-cosolvent mixtures via computer simulation: the case of ethanol in CO2.

The supercritical mixture ethanol-carbon dioxide (EtOH-CO2) with mole fraction of ethanol X(EtOH) congruent with 0.1 was investigated at 348 K, by employing the molecular dynamics simulation technique in the canonical ensemble. The local intermolecular structure of the fluid was studied in terms of the calculated appropriate pair radial distribution functions. The estimated average local coordination numbers and mole fractions around the species in the mixture reveal the existence of local composition enhancement of ethanol around the ethanol molecules. This finding indicates the nonideal mixing behavior of the mixture due to the existence of aggregation between the ethanol molecules. Furthermore, the local environment redistribution dynamics have been explored by analyzing the time correlation functions (TCFs) of the total local coordination number (solvent, cosolvent) around the cosolvent molecules in appropriate parts. The analysis of these total TCFs in the auto-(solvent-solvent, cosolvent-cosolvent) and cross-(solvent-cosolvent, cosolvent-solvent) TCFs has shown that the time dependent redistribution process of the first solvation shell of ethanol is mainly determined by the redistribution of the CO2 solvent molecules. These results might be explained on the basis of the CO2-CO2 and EtOH-CO2 intermolecular forces, which are sufficiently weaker in comparison to the EtOH-EtOH hydrogen bonding interactions, creating in this way a significantly faster redistribution of the CO2 molecules in comparison with EtOH. Finally, the self-diffusion coefficients and the single reorientational dynamics of both the cosolvent and solvent species in the mixture have been predicted and discussed in relationship with the local environment around the species, which in the case of the EtOH molecules seem to be strongly affected.

Pressure and temperature dependence of the hydrogen bonding in supercritical ethanol: a computer simulation study.

As a step toward deeper insight on the "hydrogen bonding" in supercritical ethanol (scEtOH), we carried out NVT molecular dynamics simulations of the fluid over a wide range of temperatures and pressures. The fluid was studied at SC conditions for which thermodynamic and spectroscopic (NMR, infrared, Raman, dielectric) data are available. The various site-site pair distribution functions (pdf's) were calculated, and their temperature and pressure dependence was obtained. It was found that over the thermodynamic conditions investigated here, scEtOH remains highly structured. Moreover, the characteristic behavior of the first peaks in H-H, O-O, and H-O pdf's reveals that hydrogen bonds still exist in scEtOH. The analysis focuses also on the reorientational dynamics of the bond unit vectors O-H, C-O, and of the permanent dipole moment of the molecules as well as the total dipole moment of the sample. The corresponding Legendre time correlation functions were discussed in connection to the "hydrogen bonding" in the fluid and in the context of experimental results. Specifically, the behavior of the O-H dynamics exhibits the well-known associative nature of the molecules in the system. A further analysis of the hydrogen bonds was carried out, and the degree of aggregation (average number of H-bonds per molecule) was obtained and compared with results from NMR chemical shift studies. Also the estimated monomer and free O-H groups in the fluid were compared with results from IR and Raman vibrational spectroscopy. The percentage analysis fi of the liquid and scEtOH molecules, with i = 0, 1, 2, 3, ... hydrogen bonds per molecule, has been obtained. The results show the existence of small, linear-chain oligomers formed mainly by two molecules, whereas the number of the three body oligomers, and specifically that of four body oligomers in the sample, is relatively small.