In silico screening, molecular docking, and molecular dynamics studies of SNP-derived human P5CR mutants

Mechanistic and structural insight into R2R3-MYB transcription factor in plants: molecular dynamics based binding free energy analysis

The MYB transcription factor (TF) family is essential for various plant growth and development processes, including responses to biotic and abiotic stresses. This study investigated the R2R3-MYB protein structure from five plants, including cereal crops. The R2R3-MYB protein structure was docked with the DNA structure, and the best complexes were selected for two runs of molecular dynamics (MD) simulations to investigate the key interacting residues and conformational changes in the R2R3-MYB proteins caused by DNA binding. The MM/PBSA method calculated the binding free energy for each R2R3-MYB protein-DNA complex, showing strong interaction. Hydrophobic and hydrogen bonds significantly stabilized the R2R3-MYB protein-DNA complexes. The principal component analysis showed high restrictions on the movement of protein atoms in the phase space. A similar MD simulation analysis was performed using the crystal structure of the R2R3-MYB protein-DNA complex from Arabidopsis thaliana, and the generated complexes resembled the X-ray crystal structure. This is the first detailed study on the R2R3-MYB protein-DNA complex in cereal crops, providing a cost-effective solution to identify the key interacting residues and analyze the conformational changes in the MYB domain before and after DNA binding.Communicated by Ramaswamy H. Sarma.

Exploring the prognostic significance of PKCε variants in cervical cancer

BACKGROUND

Protein Kinase C-epsilon (PKCε) is a member of the novel subfamily of PKCs (nPKCs) that plays a role in cancer development. Studies have revealed that its elevated expression levels are associated with cervical cancer. Previously, we identified pathogenic variations in its different domains through various bioinformatics tools and molecular dynamic simulation. In the present study, the aim was to find the association of its variants rs1553369874 and rs1345511001 with cervical cancer and to determine the influence of these variants on the protein-protein interactions of PKCε, which can lead towards cancer development and poor survival rates.

METHODS

The association of the variants with cervical cancer and its clinicopathological features was determined through genotyping analysis. Odds ratio and relative risk along with Fisher exact test were calculated to evaluate variants significance and disease risk. Protein-protein docking was performed and docked complexes were subjected to molecular dynamics simulation to gauge the variants impact on PKCε's molecular interactions.

RESULTS

This study revealed that genetic variants rs1553369874 and rs1345511001 were associated with cervical cancer. Smad3 interacts with PKCε and this interaction promotes cervical cancer angiogenesis; therefore, Smad3 was selected for protein-protein docking. The analysis revealed PKCε variants promoted aberrant interactions with Smad3 that might lead to the activation of oncogenic pathways. The data obtained from this study suggested the prognostic significance of PRKCE gene variants rs1553369874 and rs1345511001.

CONCLUSION

Through further in vitro and in vivo validation, these variants can be used at the clinical level as novel prognostic markers and therapeutic targets against cervical cancer.

Computational insight into in silico analysis and molecular dynamics simulation of the dimer interface residues of ALS-linked hSOD1 forms in apo/holo states: a combined experimental and bioinformatic perspective

The aggregation of misfolded SOD1 proteins in neurodegenerative illnesses is a key pathological hallmark in amyotrophic lateral sclerosis (ALS). SOD1 is stabilized and enzymatically activated after binding to Cu/Zn and forming intramolecular disulfide. SOD1 aggregation/oligomerization is triggered by the dissociation of Cu and/or Zn ions. Therefore, we compared the possible effects of ALS-associated point mutations of the holo/apo forms of WT/I149T/V148G SOD1 variants located at the dimer interface to determine structural characterization using spectroscopic methods, computational approaches as well as molecular dynamics (MD) simulations. Predictive results of computational analysis of single-nucleotide polymorphisms (SNPs) suggested that mutant SOD1 has a deleterious effect on activity and structure destabilization. MD data analysis indicated that changes in flexibility, stability, hydrophobicity of the protein as well as increased intramolecular interactions of apo-SOD1 were more than holo-SOD1. Furthermore, a decrease in enzymatic activity in apo-SOD1 was observed compared to holo-SOD1. Comparative intrinsic and ANS fluorescence results of holo/apo-WT-hSOD1 and mutants indicated structural alterations in the local environment of tryptophan residue and hydrophobic patches, respectively. Experimental and MD data supported that substitution effect and metal deficiency of mutants (apo forms) in the dimer interface may promote the tendency to protein mis-folding and aggregation, consequently disrupting the dimer-monomer equilibrium and increased propensity to dissociation dimer into SOD-monomer ultimately leading to loss of stability and function. Overall, data analysis of apo/holo SOD1 forms on protein structure and function using computational and experimental studies will contribute to a better understanding of ALS pathogenicity.

Comprehensive in silico analysis and molecular dynamics of the superoxide dismutase 1 (SOD1) variants related to amyotrophic lateral sclerosis

Amyotrophic Lateral Sclerosis (ALS) is the most frequent motor neuron disorder, with a significant social and economic burden. ALS remains incurable, and the only drugs approved for its treatments confers a survival benefit of a few months for the patients. Missense mutations in superoxide dismutase 1 (SOD1), a major cytoplasmic antioxidant enzyme, has been associated with ALS development, accounting for 23% of its familial cases and 7% of all sporadic cases. This work aims to characterize in silico the structural and functional effects of SOD1 protein variants. Missense mutations in SOD1 were compiled from the literature and databases. Twelve algorithms were used to predict the functional and stability effects of these mutations. ConSurf was used to estimate the evolutionary conservation of SOD1 amino-acids. GROMACS was used to perform molecular dynamics (MD) simulations of SOD1 wild-type and variants A4V, D90A, H46R, and I113T, which account for approximately half of all ALS-SOD1 cases in the United States, Europe, Japan, and United Kingdom, respectively. 233 missense mutations in SOD1 protein were compiled from the databases and literature consulted. The predictive analyses pointed to an elevated rate of deleterious and destabilizing predictions for the analyzed variants, indicating their harmful effects. The ConSurf analysis suggested that mutations in SOD1 mainly affect conserved and possibly functionally essential amino acids. The MD analyses pointed to flexibility and essential dynamics alterations at the electrostatic and metal-binding loops of variants A4V, D90A, H46R, and I113T that could lead to aberrant interactions triggering toxic protein aggregation. These alterations may have harmful implications for SOD1 and explain their association with ALS. Understanding the effects of SOD1 mutations on protein structure and function facilitates the design of further experiments and provides relevant information on the molecular mechanism of pathology, which may contribute to improvements in existing treatments for ALS.

Comprehensive in silico mutational-sensitivity analysis of PTEN establishes signature regions implicated in pathogenesis of Autism Spectrum Disorders

An extensively studied cancer and Autism Spectrum Disorders (ASD) gene like PTEN provided an exclusive opportunity to map its mutational-landscape, compare and establish plausible genotypic predictors of ASD-associated phenotypic outcomes. Our exhaustive in silico analysis on 4252 SNPs using >30 tools identified increased mutational-density in exon7. Phosphatase domain, although evolutionarily conserved, had the most nsSNPs localised within signature regions. The evolutionarily variable C-terminal side contained the highest truncating-SNPs outside signature regions of C2 domain and most PTMs within C-tail site which displayed maximum intolerance to polymorphisms, and permitted benign but destabilising nsSNPs that enhanced its intrinsically-disordered nature. ASD-associated SNPs localised within ATP-binding motifs and Nuclear-Localising-Sequences were the most potent triggers of ASD manifestation. These, along with variations within P, WPD and TI loops, M1 within phosphatase domain, M2 and MoRFs of C2 domain, caused severe long-range conformational fluctuations altering PTEN's dynamic stability- not observed in variations outside signature regions. 3'UTR-SNPs affected 44 strong miRNA brain-specific targets; several 5' UTR-SNPs targeted transcription-factor POLR2A and 10 pathogenic Splice-Affecting-Variants were identified.

In Silico Tools and Approaches for the Prediction of Functional and Structural Effects of Single-Nucleotide Polymorphisms on Proteins: An Expert Review

Single-nucleotide polymorphisms (SNPs) are single-base variants that contribute to human biological variation and pathogenesis of many human diseases. Among all SNP types, nonsynonymous single-nucleotide polymorphisms (nsSNPs) can alter many structural, biochemical, and functional features of a protein such as folding characteristics, charge distribution, stability, dynamics, and interactions with other proteins/nucleotides. These modifications in the protein structure can lead nsSNPs to be closely associated with many multifactorial diseases such as cancer, diabetes, and neurodegenerative diseases. Predicting structural and functional effects of nsSNPs with experimental approaches can be time-consuming and costly; hence, computational prediction tools and algorithms are being widely and increasingly utilized in biology and medical research. This expert review examines the in silico tools and algorithms for the prediction of functional or structural effects of SNP variants, in addition to the description of the phenotypic effects of nsSNPs on protein structure, association between pathogenicity of variants, and functional or structural features of disease-associated variants. Finally, case studies investigating the functional and structural effects of nsSNPs on selected protein structures are highlighted. We conclude that creating a consistent workflow with a combination of in silico approaches or tools should be considered to increase the performance, accuracy, and precision of the biological and clinical predictions made in silico.

Anti-HIV drug repurposing against SARS-CoV-2

A novel severe acute respiratory syndrome human coronavirus (SARS HCoV) was identified from respiratory illness patients (named SARS-CoV-2 by ICTV) in December 2019 and has recently emerged as a serious threat to world public health. However, no approved drugs have been found to effectively inhibit the virus. Since it has been reported that HIV protease inhibitors can be used as anti-SARS drugs by targeting SARS-CoV-1 3CLpro, we chose six approved anti-HIV drugs and investigated their binding interactions with 3CLpro to evaluate their potential to become clinical drugs for the new coronavirus pneumonia (COVID-19) caused by SARS-CoV-2 infection. The molecular docking results indicate that the 3CLpro of SARS-CoV-2 has a higher binding affinity for all the studied inhibitors than does SARS-CoV-1. Two docking complexes (indinavir and darunavir) with high docking scores were further subjected to MM-PBSA binding free energy calculations to detail the molecular interactions between these two protease inhibitors and SARS HCoV 3CLpro. Our results show that, among the inhibitors tested, darunavir has the highest binding affinity with SARS-CoV-2 and SARS-CoV-1 3CLpro, indicating that it may have the potential to be used as an anti-COVID-19 clinical drug. The mechanism behind the increased binding affinity of HIV protease inhibitors toward SARS-CoV-2 3CLpro (as compared to SARS-CoV-1) was investigated by MD simulations. Our study provides insight into the possible role of structural flexibility during interactions between SARS HCoV 3CLpro and inhibitors and sheds light on structure-based design of anti-COVID-19 drugs targeting SARS-CoV-2 3CLpro.

Deciphering the structural basis for glucocorticoid resistance caused by missense mutations in the ligand binding domain of glucocorticoid receptor

The glucocorticoid resistance hereditary condition may emerge from the occurrence of point mutations in the glucocorticoid receptor (GR), which could impair its functionality. Because the main feature of such pathology is the resistance of the hypothalamic-pituitary-adrenal axis to the hormone cortisol, we used the GR ligand binding domain three-dimensional structure to perform computational analysis for eight variants known to cause this clinical condition (I559 N, V571A, D641V, G679S, F737L, I747 M, L753F and L773P), aiming to understand, on the atom scale, how they cause glucocorticoid resistance. We observed that the mutations generated a reduced affinity to cortisol and they alter some loop conformations, which could be a consequence from changes in protein motion, which in turn could result from the reduced stability of mutant GR structures. Therefore, the analyzed mutations compromise the GR ligand binding domain structure and cortisol binding, which could characterize the glucocorticoid resistance phenotype.

A profound computational study to prioritize the disease-causing mutations in PRPS1 gene

Charcot-Marie-Tooth disease (CMT) is one of the most commonly inherited congenital neurological disorders, affecting approximately 1 in 2500 in the US. About 80 genes were found to be in association with CMT. The phosphoribosyl pyrophosphate synthetase 1 (PRPS1) is an essential enzyme in the primary stage of de novo and salvage nucleotide synthesis. The mutations in the PRPS1 gene leads to X-linked Charcot-Marie-Tooth neuropathy type 5 (CMTX5), PRS super activity, Arts syndrome, X-linked deafness-1, breast cancer, and colorectal cancer. In the present study, we obtained 20 missense mutations from UniProt and dbSNP databases and applied series of comprehensive in silico prediction methods to assess the degree of pathogenicity and stability. In silico tools predicted four missense mutations (D52H, M115 T, L152P, and D203H) to be potential disease causing mutations. We further subjected the four mutations along with native protein to 50 ns molecular dynamics simulation (MDS) using Gromacs package. The resulting trajectory files were analyzed to understand the stability differences caused by the mutations. We used the Root Mean Square Deviation (RMSD), Radius of Gyration (Rg), solvent accessibility surface area (SASA), Covariance matrix, Principal Component Analysis (PCA), Free Energy Landscape (FEL), and secondary structure analysis to assess the structural changes in the protein upon mutation. Our study suggests that the four mutations might affect the PRPS1 protein function and stability of the structure. The proposed study may serve as a platform for drug repositioning and personalized medicine for diseases that are caused by the PRPS1 deficiency.

Resolving the cofactor-binding site in the proline biosynthetic enzyme human pyrroline-5-carboxylate reductase 1

Pyrroline-5-carboxylate reductase (PYCR) is the final enzyme in proline biosynthesis, catalyzing the NAD(P)H-dependent reduction of Δ1-pyrroline-5-carboxylate (P5C) to proline. Mutations in the PYCR1 gene alter mitochondrial function and cause the connective tissue disorder cutis laxa. Furthermore, PYCR1 is overexpressed in multiple cancers, and the PYCR1 knock-out suppresses tumorigenic growth, suggesting that PYCR1 is a potential cancer target. However, inhibitor development has been stymied by limited mechanistic details for the enzyme, particularly in light of a previous crystallographic study that placed the cofactor-binding site in the C-terminal domain rather than the anticipated Rossmann fold of the N-terminal domain. To fill this gap, we report crystallographic, sedimentation-velocity, and kinetics data for human PYCR1. Structures of binary complexes of PYCR1 with NADPH or proline determined at 1.9 Å resolution provide insight into cofactor and substrate recognition. We see NADPH bound to the Rossmann fold, over 25 Å from the previously proposed site. The 1.85 Å resolution structure of a ternary complex containing NADPH and a P5C/proline analog provides a model of the Michaelis complex formed during hydride transfer. Sedimentation velocity shows that PYCR1 forms a concentration-dependent decamer in solution, consistent with the pentamer-of-dimers assembly seen crystallographically. Kinetic and mutational analysis confirmed several features seen in the crystal structure, including the importance of a hydrogen bond between Thr-238 and the substrate as well as limited cofactor discrimination.

Effect of R119G Mutation on Human P5CR1 Dynamic Property and Enzymatic Activity

Pyrroline-5-carboxylate reductase (P5CR1) is a universal housekeeping enzyme that catalyzes the reduction of Δ1-pyrroline-5-carboxylate (P5C) to proline with concomitant oxidation of NAD(P)H to NAD(P)+. The enzymatic cycle between P5C and proline is important for function in amino acid metabolism, apoptosis, and intracellular redox potential balance in mitochondria. Autosomal recessive cutis laxa (ARCL) results from a mutation in P5CR1 encoded by PYCR1. Specifically, the R119G mutation is reported to be linked to ARCL although it has not yet been characterized. We synthesized R119G P5CR1 and compared it to WT P5CR1. Foldx prediction of WT and R119G mutant P5CR1 protein stability suggests that the R119G mutation could significantly reduce protein stability. We also performed enzymatic activity assays to determine how the mutation impacts P5CR1 enzymatic function. The results of these experiments show that mutagenesis of R119 to G decreases P5CR1 catalytic efficiency for 3,4-dehydro-L-proline relative to WT. Mutagenesis and kinetic studies reveal that the activity of the mutant decreases as temperature increases from 5°C to 37°C, with almost no activity at 37°C, indicating that this mutation impairs P5CR1 function in vivo. Conversely, WT P5CR1 retains its activity after incubation at 37°C and has essentially no remaining activity at 75°C. Taken together, our experimental results indicate the R119G mutation could be an involving pathomechanism for ARCL.