Structure prediction and binding sites analysis of curcin protein of Jatropha curcas using computational approaches

Computational Methods for the Analysis and Prediction of EGFR-Mutated Lung Cancer Drug Resistance: Recent Advances in Drug Design, Challenges and Future Prospects

Lung cancer is a major cause of cancer deaths worldwide, and has a very low survival rate. Non-small cell lung cancer (NSCLC) is the largest subset of lung cancers, which accounts for about 85% of all cases. It has been well established that a mutation in the epidermal growth factor receptor (EGFR) can lead to lung cancer. EGFR Tyrosine Kinase Inhibitors (TKIs) are developed to target the kinase domain of EGFR. These TKIs produce promising results at the initial stage of therapy, but the efficacy becomes limited due to the development of drug resistance. In this paper, we provide a comprehensive overview of computational methods, for understanding drug resistance mechanisms. The important EGFR mutants and the different generations of EGFR-TKIs, with the survival and response rates are discussed. Next, we evaluate the role of important EGFR parameters in drug resistance mechanism, including structural dynamics, hydrogen bonds, stability, dimerization, binding free energies, and signaling pathways. Personalized drug resistance prediction models, drug response curve, drug synergy, and other data-driven methods are also discussed. Recent advancements in deep learning; such as AlphaFold2, deep generative models, big data analytics, and the applications of statistics and permutation are also highlighted. We explore limitations in the current methodologies, and discuss strategies to overcome them. We believe this review will serve as a reference for researchers; to apply computational techniques for precision medicine, analyzing structures of protein-drug complexes, drug discovery, and understanding the drug response and resistance mechanisms in lung cancer patients.

Prediction and validation of host-pathogen interactions by a versatile inference approach using Aspergillus fumigatus as a case study

Biological networks are characterized by diverse interactions and dynamics in time and space. Many regulatory modules operate in parallel and are interconnected with each other. Some pathways are functionally known and annotated accordingly, e.g., endocytosis, migration, or cytoskeletal rearrangement. However, many interactions are not so well characterized. For reconstructing the biological complexity in cellular networks, we combine here existing experimentally confirmed and analyzed interactions with a protein-interaction inference framework using as basis experimentally confirmed interactions from other organisms. Prediction scoring includes sequence similarity, evolutionary conservation of interactions, the coexistence of interactions in the same pathway, orthology as well as structure similarity to rank and compare inferred interactions. We exemplify our inference method by studying host-pathogen interactions during infection of Mus musculus (phagolysosomes in alveolar macrophages) with Aspergillus fumigatus (conidia, airborne, asexual spores). Three of nine predicted critical host-pathogen interactions could even be confirmed by direct experiments. Moreover, we suggest drugs that manipulate the host-pathogen interaction.

Structural Distinctive 26SK, a Ribosome-Inactivating Protein from Jatropha curcas and Its Biological Activities

Ribosome-inactivating proteins (RIPs) are a group of proteins exhibiting N-glycosidase activity leading to an inactivation of protein synthesis. Thirteen predicted Jatropha curcas RIP sequences could be grouped into RIP types 1 or 2. The expression of the RIP genes was detected in seed kernels, seed coats, and leaves. The full-length cDNA of two RIP genes (26SK and 34.7(A)SK) were cloned and studied. The 34.7(A)SK protein was successfully expressed in the host cells while it was difficult to produce even only a small amount of the 26SK protein. Therefore, the crude proteins were used from E. coli expressing 26SK and 34.7(A)SK constructs and they showed RIP activity. Only the cell lysate from 26SK could inhibit the growth of E. coli. In addition, the crude protein extracted from 26SK expressing cells displayed the effect on the growth of MDA-MB-231, a human breast cancer cell line. Based on in silico analysis, all 13 J. curcas RIPs contained RNA and ribosomal P2 stalk protein binding sites; however, the C-terminal region of the P2 stalk binding site was lacking in the 26SK structure. In addition, an amphipathic distribution between positive and negative potential was observed only in the 26SK protein, similar to that found in the anti-microbial peptide. These findings suggested that this 26SK protein structure might have contributed to its toxicity, suggesting potential uses against pathogenic bacteria in the future.

Molecular Characterization and Expression of Four Aquaporin Genes in Impatiens walleriana During Drought Stress and Recovery

Aquaporins comprise a large group of transmembrane proteins responsible for water transport, which is crucial for plant survival under stress conditions. Despite the vital role of aquaporins, nothing is known about this protein family in Impatiens walleriana, a commercially important horticultural plant, which is sensitive to drought stress. In the present study, attention is given to the molecular characterization of aquaporins in I. walleriana and their expression during drought stress and recovery. We identified four I. walleriana aquaporins: IwPIP1;4, IwPIP2;2, IwPIP2;7 and IwTIP4;1. All of them had conserved NPA motifs (Asparagine-Proline-Alanine), transmembrane helices (TMh), pore characteristics, stereochemical properties and tetrameric structure of holoprotein. Drought stress and recovery treatment affected the aquaporins expression in I. walleriana leaves, which was up- or downregulated depending on stress intensity. Expression of IwPIP2;7 was the most affected of all analyzed I. walleriana aquaporins. At 15% and 5% soil moisture and recovery from 15% and 5% soil moisture, IwPIP2;7 expression significantly decreased and increased, respectively. Aquaporins IwPIP1;4 and IwTIP4;1 had lower expression in comparison to IwPIP2;7, with moderate expression changes in response to drought and recovery, while IwPIP2;2 expression was of significance only in recovered plants. Insight into the molecular structure of I. walleriana aquaporins expanded knowledge about plant aquaporins, while its expression during drought and recovery contributed to I. walleriana drought tolerance mechanisms and re-acclimation.

Virtual screening for potential inhibitors of β(1,3)-D-glucan synthase as drug candidates against fungal cell wall

Background

To enhance the outcome in patients with invasive candidiasis, initiation of an efficient antifungal treatment in a suitable dosage is necessary. Echinocandins (e.g. caspofungin) inhibit the enzyme β(1,3)-D-glucan synthase of the fungal cell wall. Compared to azoles and other antifungal agents, echinocandins have lower adverse effects and toxicity in humans. Echinocandins are available in injectable (intravenous) form.

Methods

In this study, to identify the novel oral drug-like compounds that affect the fungal cell wall, downloaded oral drug-like compounds from the ZINC database were processed with a virtual screening procedure. The docking free energies were calculated and compared with the known inhibitor caspofungin. Four molecules were selected as the most potent ligands and subjected to hydrogen bonds analysis.

Results

Considering the hydrogen bond analysis, two compounds (ZINC71336662 and ZINC40910772) were predicted to better interact with the active site of β(1,3)-D-glucan synthase compared with caspofungin.

Conclusion

The introduced compound in this study may be valuable to analyze experimentally as a novel oral drug candidate targeting fungal cell walls.

Molecular Modeling of the Interaction Between Stem Cell Peptide and Immune Receptor in Plants

Molecular docking enables comprehensive exploration of interactions between chemical moieties and proteins. Modeling and docking approaches are useful to determine the three-dimensional (3D) structure of experimentally uncrystallized proteins and subsequently their interactions with various inhibitors and activators or peptides. Here, we describe a protocol for carrying out molecular modeling and docking of stem cell peptide CLV3p on plant innate immune receptor FLS2.

The functional impact of the C/N-terminal extensions of the mouse retinal IMPDH1 isoforms: a kinetic evaluation

Mutations in the retinal inosine monophosphate dehydrogenase1 (IMPDH1) gene is believed to be one cause of retinitis pigmentosa (RP). The main structural difference between the mutation-susceptible retinal isoforms with canonical one resides in the C- and N-terminal extensions. There are limited studies on the structure and function of terminal peptide extensions of the IMPDH1 retinal isoforms. Using recombinant murine IMPDH1 (mH1), we evaluated the kinetics of the retinal isoforms along with inhibition by some of the purine nucleotides. Molecular modeling tools were also applied to study the probable effect(s) of the terminal peptide tails on the function of the retinal isoforms. Molecular dynamic simulations indicated the possible impact of the end-terminal segments on the enzyme function through interactions with the enzyme's finger domain, affecting its critical pseudo barrel structure. The higher experimentally-determined K_m and K_i values of the retinal mIMPDH1 (546) and mIMPDH1 (603) relative to that of the canonical isoform, mIMPDH1 (514), might clearly be due to these interactions. Furthermore and despite of the canonical isoform, the retinal isoforms of mH1 exhibited no NAD⁺ substrate inhibition. The resent data would certainly provide the ground for future evaluation of the physiological significance of these variations.

Enhanced stability of manganese superoxide dismutase by amino acid replacement designed via molecular dynamics simulation

In order to improve manganese-SOD stability, three mutations were constructed via site-directed mutagenesis, and the root mean square fluctuation (RMSF) and root mean square deviation (RMSD) were used as stability assessment indexes. The amino acids of V140, E155 and E215 from wild-type mouse Mn-SOD was replaced to L140, W155 and W215, and a recombinant plasmid containing DNA segment coding wild-type and mutant Mn-SOD protein was transformed into Escherichia coli BL21 for expression. The highest enzyme activity of the mutations-MnSOD was 2050 U/mg. In addition, the recombinant protein, TM-MnSOD^{V140L, E155W, E215W} exhibited higher working temperature and improved stability compared with the wild-type Mn-SOD. Furthermore, CD spectrum analysis of the improved mutants and wild-type enzyme showed that there was no significant change in their secondary structures. This study not only expands the scope of the application of enzymes, but also helps us understand the relationship between protein structure and function.

Computational Analysis of High-Risk SNPs in Human DBY Gene Responsible for Male Infertility: A Functional and Structural Impact

BACKGROUND

DEAD-box helicase 3, Y-linked (DBY) is a candidate gene of the AZF region which is involved in spermatogenesis process. Mutations in the DBY gene may disrupt the spermatogenesis and lead to infertility in men. Identification of functionally neutral mutation from the disease-causing mutation is the biggest challenge in human genetic variation analysis. Owing to the importance of DBY in male infertility, functional analysis was carried out to reveal the association between genetic mutation and phenotypic variation through various in silico approaches.

METHODS

The present study analyzed the functional consequences of the nsSNPs in human DBY gene using SIFT, PolyPhen 2, PROVEAN, SNAP2, PMut, nsSNPAnalyzer, PhD-SNP and SNPs&GO along with stability analysis through I-Mutant2.0, MuPro and iPTREE-STAB. The conservational analysis of amino acid residues, biophysical properties and conserved domains of the DBY protein was analyzed using various computational tools. The 3D structure of the protein was generated using SPARKS-X and validated using RAMPAGE.

RESULTS

Out of 1130 SNPs reported in dbSNP, only one nsSNP (G300D) was found to have a functional effect on stability as well as the function of the DBY protein. The results showed the presence of G300 in the putative structure of DBY domain.

CONCLUSION

To the best of our knowledge, this is the first study to detect pathologically significant nsSNPs (G300D) through a computational approach in the DBY which can be useful for development in potent drug discovery studies.

Computational analysis, structural modeling and ligand binding site prediction of Plasmodium falciparum 1-deoxy-d-xylulose-5-phosphate synthase

Malaria remains one of the most serious infectious diseases in the world. There are five human species of the Plasmodium genus, of which Plasmodium falciparum is the most virulent and responsible for the vast majority of malaria related deaths. The unique biochemical processes that exist in Plasmodium falciparum provide a useful way to develop novel inhibitors. One such biochemical pathway is the methyl erythritol phosphate pathway (MEP), required to synthesize isoprenoid precursors. In the present study, a detailed computational analysis has been performed for 1-deoxy-d-xylulose-5-phosphate synthase, a key enzyme in MEP. The protein is found to be stable and residues from 825 to 971 are highly conserved across species. The homology model of the enzyme is developed using three web-based servers and Modeller software. It has twelve disordered regions indicating its druggability. Virtual screening of ZINC database identifies ten potential compounds in thiamine diphosphate binding region of the enzyme.

Structural modeling and in silico analysis of non-synonymous single nucleotide polymorphisms of human 3β-hydroxysteroid dehydrogenase type 2

Single-nucleotide polymorphisms (SNPs), a most common type of genetic mutations, result from single base pair alterations. Non-synonymous SNPs (nsSNP) occur in the coding regions of a gene and result in single amino acid substitution which might have the potential to affect the function as well as structure of the corresponding protein. In human the 3β-hydroxysteroid dehydrogenases/Δ^4,5-isomerase type 2 (HSD3B2) is an important membrane-bound enzyme involved in the dehydrogenation and Δ^4,5-isomerization of the Δ⁵-steroid precursors into their respective Δ⁴-ketosteroids in the biosynthesis of steroid hormones such as glucocorticoids, mineralocorticoids, progesterone, androgens, and estrogens in tissues such as adrenal gland, ovary, and testis. Most of the nsSNPs of HSD3B2 are still uncharacterized in terms of their disease causing potential. So, this study has been undertaken to explore and extend the knowledge related to the effect of nsSNPs on the stability and function of the HSD3B2. In this study sixteen nsSNP of HSD3B2 were subjected to in silico analysis using nine different algorithms: SIFT, PROVEAN, PolyPhen, MutPred, SNPeffect, nsSNP Analyzer, PhD SNP, stSNP, and I Mutant 2.0. The results obtained from the analysis revealed that the prioritization of diseases associated amino acid substitution as evident from possible alteration in structure–function relationship. Structural phylogenetic analysis using ConSurf revealed that the functional residues are highly conserved in human HSD3B2; and most of the disease associated nsSNPs are within these conserved residues. Structural theoritical models of HSD3B2 were created using HHPred, Phyre2 and RaptorX server. The predicted models were evaluated to get the best one for structural understanding of amino acid substitutions in three dimensional spaces.

Understanding the mechanism of atovaquone drug resistance in Plasmodium falciparum cytochrome b mutation Y268S using computational methods

The rapid appearance of resistant malarial parasites after introduction of atovaquone (ATQ) drug has prompted the search for new drugs as even single point mutations in the active site of Cytochrome b protein can rapidly render ATQ ineffective. The presence of Y268 mutations in the Cytochrome b (Cyt b) protein is previously suggested to be responsible for the ATQ resistance in Plasmodium falciparum (P. falciparum). In this study, we examined the resistance mechanism against ATQ in P. falciparum through computational methods. Here, we reported a reliable protein model of Cyt bc1 complex containing Cyt b and the Iron-Sulphur Protein (ISP) of P. falciparum using composite modeling method by combining threading, ab initio modeling and atomic-level structure refinement approaches. The molecular dynamics simulations suggest that Y268S mutation causes ATQ resistance by reducing hydrophobic interactions between Cyt bc1 protein complex and ATQ. Moreover, the important histidine contact of ATQ with the ISP chain is also lost due to Y268S mutation. We noticed the induced mutation alters the arrangement of active site residues in a fashion that enforces ATQ to find its new stable binding site far away from the wild-type binding pocket. The MM-PBSA calculations also shows that the binding affinity of ATQ with Cyt bc1 complex is enough to hold it at this new site that ultimately leads to the ATQ resistance.