Identification of deleterious missense variants of human Piwi like RNA-mediated gene silencing 1 gene and their impact on PAZ domain structure, stability, flexibility and dimension: in silico analysis

Computational study of the potential impact of AURKC missense SNPs on AURKC-INCENP interaction and their correlation to macrozoospermia

Aurora Kinase C (AURKC) is considered an important element in Chromosome Passenger Complex (CPC), its interaction with Inner Centromere Protein (INCENP) plays a critical role in the establishment and the recruitment of a stable CPC during spermatogenesis. Genetic variations of AURKC gene are susceptible to impact AURKC-INCENP interaction, which may affect CPC stability and predispose male subjects to macrozoospermia. In this study, we systematically applied computational approaches using different bioinformatic tools to predict the effect of missense SNPs reported on AURKC gene, we selected the deleterious ones and we introduced their corresponding amino acid substitutions on AURKC protein structure. Then we did a protein-protein docking between AURKC variants and INCENP followed by a structural assessment of each resulting complex using PRODIGY server, Yassara view, Ligplot + and we choose the complexes of the most impactful variants for molecular dynamics (MD) simulation study. Seventeen missense SNPs of AURKC were identified as deleterious between all reported ones. All of them were located on relatively conserved positions on AURKC protein according to Consurf server. Only the four missense SNPs; E91K, D166V, D221Y and G235V were ranked as the most impactful ones and were chosen for MD simulation. D221Y and G235V were responsible for the most remarkable changes on AURKC-INCENP structural stability, therefore, they were selected as the most deleterious ones. Experimental studies are recommended to test the actual effect of these two variants and their actual impact on the morphology of sperm cells.Communicated by Ramaswamy H. Sarma.

Exploring the dysregulated mRNAs-miRNAs-lncRNAs interactions associated to idiopathic non-obstructive azoospermia

Non-obstructive azoospermia (NOA) is the most clinical problem in case of infertility. About 70% of NOA patients are idiopathic with uncharacterized molecular mechanisms. This study aimed to analyze the possible pathogenic miRNA-target gene interaction and lncRNA-miRNA association involved in NOA. In the current study, differentially expressed (DE) nRNAs, miRNAs and lncRNAs were determined using the microarray dataset and statistical software R. miRNAs-mRNA and miRNA-lncRNA interactions were identified and the base-pair binding between the seed region of miRNAs and complementary nucleotides in 3' UTR of mRNAs were analyzed. The influence of the validated single nucleotide polymorphisms (SNPs) was described by calculating the minimum free energy (MFE) of the interaction. A total of 74 mRNAs, 14 miRNAs, and 10 lncRNAs were identified to have significant differential expression in testicular tissue between patients and the fertile group. Four of the DE-mRNAs and all of the reported DE-miRNAs were upregulated. In addition, all of the represented DE-lncRNAs were showed to be downregulated. miR-509-5p and miR-27b-3p were found to interact with target gene polo-like kinase 1 (PLK1) and Cysteine-rich secretory protein2 (CRISP2), respectively. Rs550967205 (A > G) positioned at 3' UTR CRISP2 and rs544604911 (T > C) located at 3' UTR PLK1, with lowest MFE in miRNA-mRNA interaction, were assumed to have possible pathogenic roles linked to spermatogenesis arrest. The results of the study provide new clues to understand the regulatory roles of miRNAs and lncRNAs in the pathogenesis and diagnosis of idiopathic azoospermia. Communicated by Ramaswamy H. Sarma.

Dominant ACO2 mutations are a frequent cause of isolated optic atrophy

Biallelic mutations in ACO2, encoding the mitochondrial aconitase 2, have been identified in individuals with neurodegenerative syndromes, including infantile cerebellar retinal degeneration and recessive optic neuropathies (locus OPA9). By screening European cohorts of individuals with genetically unsolved inherited optic neuropathies, we identified 61 cases harbouring variants in ACO2, among whom 50 carried dominant mutations, emphasizing for the first time the important contribution of ACO2 monoallelic pathogenic variants to dominant optic atrophy. Analysis of the ophthalmological and clinical data revealed that recessive cases are affected more severely than dominant cases, while not significantly earlier. In addition, 27% of the recessive cases and 11% of the dominant cases manifested with extraocular features in addition to optic atrophy. In silico analyses of ACO2 variants predicted their deleterious impacts on ACO2 biophysical properties. Skin derived fibroblasts from patients harbouring dominant and recessive ACO2 mutations revealed a reduction of ACO2 abundance and enzymatic activity, and the impairment of the mitochondrial respiration using citrate and pyruvate as substrates, while the addition of other Krebs cycle intermediates restored a normal respiration, suggesting a possible short-cut adaptation of the tricarboxylic citric acid cycle. Analysis of the mitochondrial genome abundance disclosed a significant reduction of the mitochondrial DNA amount in all ACO2 fibroblasts. Overall, our data position ACO2 as the third most frequently mutated gene in autosomal inherited optic neuropathies, after OPA1 and WFS1, and emphasize the crucial involvement of the first steps of the Krebs cycle in the maintenance and survival of retinal ganglion cells.

MicroRNAs Responding to Space Radiation

High-energy and high-atom-number (HZE) space radiation poses an inevitable potential threat to astronauts on deep space exploration missions. Compared with low-LET radiation, high-energy and high-LET radiation in space is more efficient in inducing clustered DNA damage with more serious biological consequences, such as carcinogenesis, central nervous system injury and degenerative disease. Space radiation also causes epigenetic changes in addition to inducing damage at the DNA level. Considering the important roles of microRNAs in the regulation of biological responses of radiation, we systematically reviewed both expression profiling and functional studies relating to microRNAs responding to space radiation as well as to space compound environment. Finally, the directions for improvement of the research related to microRNAs responding to space radiation are proposed. A better understanding of the functions and underlying mechanisms of the microRNAs responding to space radiation is of significance to both space radiation risk assessment and therapy development for lesions caused by space radiation.

Potential inhibitors of SARS-cov-2 RNA dependent RNA polymerase protein: molecular docking, molecular dynamics simulations and MM-PBSA analyses

The SARS-cov-2 RNA dependent RNA polymerase (nsp12) is a crucial viral enzyme that catalyzes the replication of RNA from RNA templates. The fixation of some ligands in the active site may alter the viral life cycle. The aim of the present study is to identify the conservation level of nsp12 motifs (A-G), using consurf server, and discover their interactions with rifabutin, rifampicin, rifapentin, sorangicin A, streptolydigin, myxopyronin B, VXR and VRX using AutoDockTools-1.5.6, Gromacs 2018.2 and g-mmpbsa. Thus, the most of amino acids residues located in nsp12 protein Motifs (A-G) were predicted as highly conserved. The binding energies of streptolydigin, VXR, rifabutin, rifapentine, VRX, sorangicin A, myxopyronin B and rifampicin with nsp12 protein are -8.11, -8.23, -7.14, -6.94, -6.55, -5.46, -5.33 and -5.26 kcal/mol, respectively. In the other hand, the binding energies of ligand in the same order with nsp7-nsp8-nsp12 complex are -7.23, -7.08, -7.21, -7, -6.59, -8.73, -5.52, -5.87 kcal/mol, respectively. All ligands interact with at least two nsp12 motifs. The molecular dynamics simulation of nsp12-streptolydigin and nsp12-VXR complexes shows that these two complexes are stable and the number of hydrogen bonds as a function of time, after 30 ns of simulation, varies between 0 and 6 for nsp12-streptolydigin complex and between 0 and 4 for nsp12-VXR complex. The average of free binding energies obtained using g_mmpbsa, after 30 ns of simulation, is -191.982 Kj/mol for nsp12-streptolydigin complex and -153.583 Kj/mol for nsp12-VXR complex. Our results suggest that these ligands may be used as inhibitors of SARS-cov-2 nsp12 protein.Communicated by Ramaswamy H. Sarma.