Background MicroRNAs (miRNAs) are noncoding RNA molecules heavily involved in human tumors, in which few of them circulating the human body. Finding a tumor-associated signature of miRNA, that is, the minimum miRNA entities to be measured for discriminating both different types of cancer and normal tissues, is of utmost importance. Feature selection techniques applied in machine learning can help however they often provide naive or biased results. Results An ensemble feature selection strategy for miRNA signatures is proposed. miRNAs are chosen based on consensus on feature relevance from high-accuracy classifiers of different typologies. This methodology aims to identify signatures that are considerably more robust and reliable when used in clinically relevant prediction tasks. Using the proposed method, a 100-miRNA signature is identified in a dataset of 8023 samples, extracted from TCGA. When running eight-state-of-the-art classifiers along with the 100-miRNA signature against the original 1046 features, it could be detected that global accuracy differs only by 1.4%. Importantly, this 100-miRNA signature is sufficient to distinguish between tumor and normal tissues. The approach is then compared against other feature selection methods, such as UFS, RFE, EN, LASSO, Genetic Algorithms, and EFS-CLA. The proposed approach provides better accuracy when tested on a 10-fold cross-validation with different classifiers and it is applied to several GEO datasets across different platforms with some classifiers showing more than 90% classification accuracy, which proves its cross-platform applicability. Conclusions The 100-miRNA signature is sufficiently stable to provide almost the same classification accuracy as the complete TCGA dataset, and it is further validated on several GEO datasets, across different types of cancer and platforms. Furthermore, a bibliographic analysis confirms that 77 out of the 100 miRNAs in the signature appear in lists of circulating miRNAs used in cancer studies, in stem-loop or mature-sequence form. The remaining 23 miRNAs offer potentially promising avenues for future research.

The G-protein coupled estrogen receptor 1 GPER/GPR30 is a transmembrane seven-helix (7TM) receptor involved in the growth and proliferation of breast cancer. Due to the absence of a crystal structure of GPER/GPR30, in this work, molecular modeling studies have been carried out to build a three-dimensional structure, which was subsequently refined by molecular dynamics (MD) simulations (up to 120 ns). Furthermore, we explored GPER/GPR30's molecular recognition properties by using reported agonist ligands (G1, estradiol (E2), tamoxifen, and fulvestrant) and the antagonist ligands (G15 and G36) in subsequent docking studies. Our results identified the E2 binding site on GPER/GPR30, as well as other receptor cavities for accepting large volume ligands, through GPER/GPR30 π-π, hydrophobic, and hydrogen bond interactions. Snapshots of the MD trajectory at 14 and 70 ns showed almost identical binding motifs for G1 and G15. It was also observed that C107 interacts with the acetyl oxygen of G1 (at 14 ns) and that at 70 ns the residue E275 interacts with the acetyl group and with the oxygen from the other agonist whereas the isopropyl group of G36 is oriented toward Met141, suggesting that both C107 and E275 could be involved in the protein activation. This contribution suggest that GPER1 has great structural changes which explain its great capacity to accept diverse ligands, and also, the same ligand could be recognized in different binding pose according to GPER structural conformations.

Circulating microRNAs (miRNA) are small noncoding RNA molecules that can be detected in bodily fluids without the need for major invasive procedures on patients. miRNAs have shown great promise as biomarkers for tumors to both assess their presence and to predict their type and subtype. Recently, thanks to the availability of miRNAs datasets, machine learning techniques have been successfully applied to tumor classification. The results, however, are difficult to assess and interpret by medical experts because the algorithms exploit information from thousands of miRNAs. In this work, we propose a novel technique that aims at reducing the necessary information to the smallest possible set of circulating miRNAs. The dimensionality reduction achieved reflects a very important first step in a potential, clinically actionable, circulating miRNA-based precision medicine pipeline. While it is currently under discussion whether this first step can be taken, we demonstrate here that it is possible to perform classification tasks by exploiting a recursive feature elimination procedure that integrates a heterogeneous ensemble of high-quality, state-of-the-art classifiers on circulating miRNAs. Heterogeneous ensembles can compensate inherent biases of classifiers by using different classification algorithms. Selecting features then further eliminates biases emerging from using data from different studies or batches, yielding more robust and reliable outcomes. The proposed approach is first tested on a tumor classification problem in order to separate 10 different types of cancer, with samples collected over 10 different clinical trials, and later is assessed on a cancer subtype classification task, with the aim to distinguish triple negative breast cancer from other subtypes of breast cancer. Overall, the presented methodology proves to be effective and compares favorably to other state-of-the-art feature selection methods.

Docking studies include many conformations to predict binding free energies (scoring functions) and to search (scoring sampling) for the most representative binding conformations. Therefore, several biological properties, from side chain residues to complete protein motions, have been included in docking studies to improve theoretical predictions.

Molecular Dynamics (MD) simulations is a computational method that employs Newton's laws to evaluate the motions of water, ions, small molecules, and macromolecules or more complex systems, for example, whole viruses, to reproduce the behavior of the biological environment, including water molecules and lipid membranes. Specifically, structural motions, such as those that are dependent of the temperature and solute/ solvent are very important to study the recognition pattern of ligandprotein or protein-protein complexes, in that sense, MD simulations are very useful because these motions can be modeled using this methodology. Furthermore, MD simulations for drug design provide insights into the structural cavities required to design novel structures with higher affinity to the target. Also, the employment of MD simulations to drug design can help to refine the three-dimensional (3D) structure of targets in order to obtain a better sampling of the binding poses and more reliable affinity values with better structural advantages, because they incorporate some biological conditions that include structural motions compared to traditional docking procedures. This work analyzes the concepts and applicability of MD simulations for drug design because molecular structural motions are considered, and these help to identify hot spots, decipher structural details in the reported protein sites, as well as to eliminate sites that could be structural artifacts which could be originated from the structural characterization conditions from MD. Moreover, better free energy values for protein ligand recognition can also be obtained, and these can be validated under experimental procedures due to the robustness of the MD simulation methods.

The worldwide health emergency of the SARS-CoV-2 pandemic and the absence of a specific treatment for this new coronavirus have led to the use of computational strategies (drug repositioning) to search for treatments. The aim of this work is to identify FDA (Food and Drug Administration)-approved drugs with the potential for binding to the spike structural glycoprotein at the hinge site, receptor binding motif (RBM), and fusion peptide (FP) using molecular docking simulations. Drugs that bind to amino acids are crucial for conformational changes, receptor recognition, and fusion of the viral membrane with the cell membrane. The results revealed some drugs that bind to hinge site amino acids (varenicline, or steroids such as betamethasone while other drugs bind to crucial amino acids in the RBM (naldemedine, atovaquone, cefotetan) or FP (azilsartan, maraviroc, and difluprednate); saquinavir binds both the RBM and the FP. Therefore, these drugs could inhibit spike glycoprotein and prevent viral entry as possible anti-COVID-19 drugs. Several drugs are in clinical studies; by focusing on other pharmacological agents (candesartan, atovaquone, losartan, maviroc and ritonavir) in this work we propose an additional target: the spike glycoprotein. These results can impact the proposed use of treatments that inhibit the first steps of the virus replication cycle.

Synaptic vesicle protein 2A (SV2A) is an integral membrane protein necessary for the proper function of the central nervous system and is associated to the physiopathology of epilepsy. SV2A is the molecular target of the anti-epileptic drug levetiracetam and its racetam analogs. The racetam binding site in SV2A and the non-covalent interactions between racetams and SV2A are currently unknown; therefore, an in silico study was performed to explore these issues. Since SV2A has not been structurally characterized with X-ray crystallography or nuclear magnetic resonance, a three-dimensional (3D) model was built. The model was refined by performing a molecular dynamics simulation (MDS) and the interactions of SV2A with the racetams were determined by docking studies. A reliable 3D model of SV2A was obtained; it reached structural equilibrium during the last 15 ns of the MDS (50 ns) with remaining structural motions in the N-terminus and long cytoplasmic loop. The docking studies revealed that hydrophobic interactions and hydrogen bonds participate importantly in ligand recognition within the binding site. Residues T456, S665, W666, D670 and L689 were important for racetam binding within the trans-membrane hydrophilic core of SV2A. Identifying the racetam binding site within SV2A should facilitate the synthesis of suitable radio-ligands to study treatment response and possibly epilepsy progression.

No abstract