Effects of Pro1266Leu mutation on structure and function of glycoprotein Ib binding domain of von Willebrand factor.

Glycoprotein Ibα (GpIbα) binding ability of A1 domain of von Willebrand factor (vWF) facilitates platelet adhesion that plays a crucial role in maintaining hemostasis and thrombosis at the site of vascular damage. There are both "loss as well as gain of function" mutations observed in this domain. Naturally occurring "gain of function" mutations leave self-activating impacts on the A1 domain which turns the normal binding to characteristic constitutive binding with GPIbα. These "gain of function" mutations are associated with the von Willebrand disease type 2B. In recent years, studies focused on understanding the mechanism and conformational patterns attached to these phenomena have been conducted, but the conformational pathways leading to such binding patterns are poorly understood as of now. To obtain a microscopic picture of such events for the better understanding of pathways, we used molecular dynamics (MD) simulations along with principal component analysis and normal mode analysis to study the effects of Pro1266Leu (Pro503Leu in structural context) mutation on the structure and function of A1 domain of vWF. MD simulations have provided atomic-level details of intermolecular motions as a function of time to understand the dynamic behavior of A1 domain of vWF. Comparative analysis of the trajectories obtained from MD simulations of both the wild type and Pro503Leu mutant suggesting appreciable conformational changes in the structure of mutant which might provide a basis for assuming the "gain of function" effects of these mutations on the A1 domain of vWF, resulting in the constitutive binding with GpIbα.

Investigating the role of transcription factors of pancreas development in pancreatic cancer.

Pancreatic cancer (PC) is the seventh most common cause of cancer-related deaths worldwide that kills more than 300,000 people every year. Prognosis of PC is very poor with a five-year survival rate about 5%. The most common and highly observed type of PC is pancreatic ductal adenocarcinoma (PDAC). It is preceded by the progression of precursor lesions such as Pancreatic Intraepithelial Neoplasia (PanIN), Intraductal Papillary Neoplasm (IPMN) and Mucinous Cystic Neoplasm (MCN). PanIN is the most common among these premalignant lesions. Genes orchestrating the origin and differentiation of cells during organogenesis have the tendency to produce tumor cells in response to activating or inactivating mutations. Based on the following premise, we discuss the role of transcription factors (TFs) of pancreas development and cell fate differentiation in PC. Pancreas/duodenum homeobox protein 1 (PDX1), Pancreas transcription factor 1 subunit alpha (PTF1A), Nuclear receptor subfamily 5 group A member 2 (NR5A2), Hepatocyte nuclear factor 1-alpha (HNF1A) and Hepatocyte nuclear factor 1-beta (HNF1B) play vital role in the development and differentiation of pancreatic precursor cells. Mutated KRAS induces abnormalities in the regular function of these TFs which in turn cause abnormal cell growth and proliferation that leads to cancer. Thus, these TFs are highly susceptible for the origin of PC. Therefore, we propose that these TFs can be treated as therapeutic targets for the development of anticancer drugs.

Identification of functional candidates amongst hypothetical proteins of Mycobacterium leprae Br4923, a causative agent of leprosy.

Mycobacterium leprae is an intracellular obligate parasite that causes leprosy in humans, and it leads to the destruction of peripheral nerves and skin deformation. Here, we report an extensive analysis of the hypothetical proteins (HPs) from M. leprae strain Br4923, assigning their functions to better understand the mechanism of pathogenesis and to search for potential therapeutic interventions. The genome of M. leprae encodes 1604 proteins, of which the functions of 632 are not known (HPs). In this paper, we predicted the probable functions of 312 HPs. First, we classified all HPs into families and subfamilies on the basis of sequence similarity, followed by domain assignment, which provides many clues for their possible function. However, the functions of 320 proteins were not predicted because of low sequence similarity with proteins of known function. Annotated HPs were categorized into enzymes, binding proteins, transporters, and proteins involved in cellular processes. We found several novel proteins whose functions were unknown for M. leprae. These proteins have a requisite association with bacterial virulence and pathogenicity. Finally, our sequence-based analysis will be helpful for further validation and the search for potential drug targets while developing effective drugs to cure leprosy.

Pan-cancer analysis of Chromobox (CBX) genes for prognostic significance and cancer classification.

Polycomb group of proteins play a significant role in chromatin remodelling essential for epigenetic regulation of transcription. Chromobox (CBX) gene family is an important part of canonical polycomb repressive complex 1 (PRC1), belonging to the polycomb group involved in chromatin remodelling. Aberrations in CBX expression are linked to various cancers. To assess their biomarker significance, we performed a pan-cancer analysis of CBX mRNA levels in 18 cancer types. We also performed cancer classification using CBX genes as distinctive features for machine learning model development. Logistic regression (L.R.), support vector machine (SVM), random forest (R.F.), decision tree (D.T.), and XGBoost (XGB) algorithms for model training and classification. The expression of CBX genes was significantly changed in four cancer types, i.e., cholangiocarcinoma (CHOL), colon adenocarcinoma (COAD), lung adenocarcinoma (LUAD), and lung squamous cell carcinoma (LUSC). The fold change (FC) values suggest that CBX2 was significantly upregulated in CHOL (FC = 1.639), COAD (FC = 1.734), and LUSC (FC = 1.506). On the other hand, CBX7 was found downregulated in COAD (FC = -1.209), LUAD (FC = -1.190), and LUSC (FC = -1.214). The performance of machine learning models for classification was excellent. L.R., R.F., SVM, and XGB obtained a prediction accuracy of 100 % for most cancers. However, D.T. performed comparatively poorly in prediction accuracy. The results suggest that CBX expression is significantly altered in all the cancers studied; therefore, they might be treated as potential biomarkers for therapeutic intervention of these cancers.

Insight into the Conformational Transitions of Serine Acetyl Transferase Isoforms in E. histolytica: Implications for Structural and Functional Balance.

Serine acetyl transferase (SAT) is one of the crucial enzymes in the cysteine biosynthetic pathway and an essential enzyme for the survival of Entamoeba histolytica, the causative agent of amoebiasis. E. histolytica expresses three isoforms of SAT, where SAT1 and SAT2 are inhibited by the final product cysteine, while SAT3 is not inhibited. SAT3 has a slightly elongated C-terminus compared to SAT1. To understand the stability and conformational transition between two secondary structures of proteins, we measured the effect of urea, a chemical denaturant, on two isoforms of SAT (SAT1 and SAT3) of E. histolytica. The effect of urea on the structure and stability of SAT1 and SAT3 was determined by measuring changes in their far-UV circular dichroism (CD), Trp fluorescence, and near-UV absorption spectra. The urea-induced normal transition curves suggested that the structural transition is reversible and follows a two-state process. Analysis of the urea-induced transition of all optical properties for the stability parameters ΔG D° (Gibbs free energy change (ΔG D) in the absence of urea), m (dependence of ΔG D on urea concentration), and C m (midpoint of urea transition) suggested that SAT1 is more stable than SAT3. Characterization of the end product of the urea-induced transition of both proteins by the far-UV CD and Trp-fluorescence and near-UV absorbance suggested that urea causes α-helix to ß-sheet transition and burial of Trp residues, respectively. To support the in vitro findings, 100 ns molecular dynamics simulations (in silico study) were performed. Both the spectroscopic and molecular dynamics approaches clearly indicated that SAT1 is more stable than SAT3. SAT3 has evolved to escape the feedback inhibition to keep producing cysteine, but in the process, it compromises its structural stability relative to SAT1.

Measuring Structural Changes in Cytochrome c under Crowded Conditions Using In Vitro and In Silico Approaches.

It is known from in vitro studies that macromolecular crowding in the cell effects protein structure, stability and function; but predictive studies are relatively unexplored. There are few reports where the effect of various crowder mixtures has been exploited to discern their combined effect on the structural stability of proteins. These studies are more significant because their effect can mimicked with in vivo conditions, where the environment is heterogeneous. Effects of two crowders, polyethylene glycol (PEG 400 Da), and its monomer ethylene glycol (EG) alone and in mixture on the structural stability of cytochrome c (cyt c) were determined using various spectroscopic and bioinformatics tools. The main conclusions of our study are (i) the monomer EG has a kosmotropic effect on the protein (stabilizes the protein), and has no significant effect on the tertiary structure; (ii) PEG 400 destabilizes the structure as well as the stability of the protein; and (iii) EG counteracts the destabilizing effect of PEG 400. From this investigation, it seems evident that proteins may fold or unfold in the crowded environment of the cell where various interactions assist them to maintain their structure for their functions. Bioinformatics approaches were also used to support all of the in vitro observations. Cyt c is functional protein; if the structure of the protein is modulated due to change in the environment its nature of function will also change. Our research addresses the question by modulating the environment around the protein, and the macromolecule (protein) conformation dynamics and interaction study via in vitro and in silico approaches which indirectly compares with that of the environment in-cellular milieu, which is highly crowded.

Impact of glioblastoma multiforme associated mutations on the structure and function of MAP/microtubule affinity regulating kinase 4.

MAP/Microtubule affinity regulating kinase 4 (MARK4) plays an important role in the regulation of microtubule dynamics by phosphorylation of tau protein. A higher expression of MARK4 is observed in the glioblastoma multiforme (GBM) cell lines. We identified eight synonymous and non-synonymous mutations in the MARK4 gene related to GBM in The Cancer Genome Atlas (TCGA) consortium. Out of these, three non-synonymous mutations were found in the catalytic domain of the protein (Lys231Asn, Tyr247His and Arg265Gln), were selected to see the possible deleterious effects on the structure and function using the cutting-edge in-silico tools. In addition, molecular dynamics simulation, principal component analysis, dynamic cross correlation matrix analysis and correlation network analysis were performed to gain insights into the conformation of the MARK4 and its mutants. We found that, Tyr247His shows a maximum deleterious impact, reflected from structural deviation in comparison to Lys231Asn and Arg265Gln. In conclusion, Tyr247His mutant of MARK4 has relatively higher chances of affecting the structure and function of the protein thus leading to abnormal MARK4 activity which is associated to GBM.Communicated by Ramaswamy H. Sarma.

Investigating host-virus interaction mechanism and phylogenetic analysis of viral proteins involved in the pathogenesis.

Since the emergence of yellow fever in the Americas and the devastating 1918 influenza pandemic, biologists and clinicians have been drawn to human infecting viruses to understand their mechanisms of infection better and develop effective therapeutics against them. However, the complex molecular and cellular processes that these viruses use to infect and multiply in human cells have been a source of great concern for the scientific community since the discovery of the first human infecting virus. Viral disease outbreaks, such as the recent COVID-19 pandemic caused by a novel coronavirus, have claimed millions of lives and caused significant economic damage worldwide. In this study, we investigated the mechanisms of host-virus interaction and the molecular machinery involved in the pathogenesis of some common human viruses. We also performed a phylogenetic analysis of viral proteins involved in host-virus interaction to understand the changes in the sequence organization of these proteins during evolution for various strains of viruses to gain insights into the viral origin's evolutionary perspectives.

Evaluation of pyrazolopyrimidine derivatives as microtubule affinity regulating kinase 4 inhibitors: Towards therapeutic management of Alzheimer's disease.

Microtubule affinity regulating kinase 4 (MARK4) plays essential role in the tau-assisted regulation of microtubule dynamics. Over expression of MARK4 causes early phosphorylation of Ser262 of tau protein which is essential for microtubule binding. Hyperphosphorylation of tau protein causes the formation of paired helical fragments and neurofibrillary tangles, the hallmarks of Alzheimer's disease. Targeting the modulation of MARK4 activity is an effective strategy for therapeutic intervention of Alzheimer's and other MARK4 associated neurodegenerative diseases. Having role of pyrazolopyrimidine derivatives in the therapeutic management of neurodegenerative diseases, we have tried to estimate their binding affinity with the MARK4. We performed in silico screening of 59 pyrazolopyrimidine derivatives against MARK4 and obtained a few best possible inhibitors. Molecular docking-based interaction analysis suggested five potential leads that were further analyzed using molecular dynamics simulations, MM/PBSA, principal component analysis and graph theory based dynamic network analysis to observe structural changes caused due to ligand binding. All these computational analyses suggested that compounds with PubChem IDs: 91895678, 91895679, 91895692, 91145515 and 90794095 may be further exploited to address Alzheimer's and other neurodegenerative diseases.Communicated by Ramaswamy H. Sarma.

Insights into SARS-CoV-2 genome, structure, evolution, pathogenesis and therapies: Structural genomics approach.

The sudden emergence of severe respiratory disease, caused by a novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has recently become a public health emergency. Genome sequence analysis of SARS-CoV-2 revealed its close resemblance to the earlier reported SARS-CoV and Middle East respiratory syndrome coronavirus (MERS-CoV). However, initial testing of the drugs used against SARS-CoV and MERS-CoV has been ineffective in controlling SARS-CoV-2. The present study highlights the genomic, proteomic, pathogenesis, and therapeutic strategies in SARS-CoV-2 infection. We have carried out sequence analysis of potential drug target proteins in SARS-CoV-2 and, compared them with SARS-CoV and MERS viruses. Analysis of mutations in the coding and non-coding regions, genetic diversity, and pathogenicity of SARS-CoV-2 has also been done. A detailed structural analysis of drug target proteins has been performed to gain insights into the mechanism of pathogenesis, structure-function relationships, and the development of structure-guided therapeutic approaches. The cytokine profiling and inflammatory signalling are different in the case of SARS-CoV-2 infection. We also highlighted possible therapies and their mechanism of action followed by clinical manifestation. Our analysis suggests a minimal variation in the genome sequence of SARS-CoV-2, may be responsible for a drastic change in the structures of target proteins, which makes available drugs ineffective.

Investigating the structural features of chromodomain proteins in the human genome and predictive impacts of their mutations in cancers.

Epigenetic readers are specific proteins which recognize histone marks and represents the underlying mechanism for chromatin regulation. Histone H3 lysine methylation is a potential epigenetic code for the chromatin organization and transcriptional control. Recognition of histone methylation is achieved by evolutionary conserved reader modules known as chromodomain, identified in several proteins, and is involved in transcriptional silencing and chromatin remodelling. Genetic perturbations within the structurally conserved chromodomain could potentially mistarget the reader protein and impair their regulatory pathways, ultimately leading to cellular chaos by setting the stage for tumor development and progression. Here, we report the structural conservations associated with diverse functions, prognostic significance and functional consequences of mutations within chromodomain of human proteins in distinct cancers. We have extensively analysed chromodomain containing human proteins in terms of their structural-functional ability to act as a molecular switch in the recognition of methyl-lysine recognition. We further investigated the combinatorial potential, target promiscuity and binding specificity associated with their underlying mechanisms. Indeed, the molecular mechanism of epigenetic silencing significantly underlies a newer cancer therapy approach. We hope that a critical understanding of chromodomains will pave the way for novel paths of research providing newer insights into the designing of effective anti-cancer therapies.

Advancements in Docking and Molecular Dynamics Simulations Towards Ligand-receptor Interactions and Structure-function Relationships.

Protein-ligand interaction is an imperative subject in structure-based drug design and protein function prediction process. Molecular docking is a computational method which predicts the binding of a ligand molecule to the particular receptor. It predicts the binding pose, strength and binding affinity of the molecules using various scoring functions. Molecular docking and molecular dynamics simulations are widely used in combination to predict the binding modes, binding affinities and stability of different protein-ligand systems. With advancements in algorithms and computational power, molecular dynamics simulation is now a fundamental tool to investigative bio-molecular assemblies at atomic level. These methods in association with experimental support have been of great value in modern drug discovery and development. Nowadays, it has become an increasingly significant method in drug discovery process. In this review, we focus on protein-ligand interactions using molecular docking, virtual screening and molecular dynamics simulations. Here, we cover an overview of the available methods for molecular docking and molecular dynamics simulations, and their advancement and applications in the area of modern drug discovery. The available docking software and their advancement including application examples of different approaches for drug discovery are also discussed. We have also introduced the physicochemical foundations of molecular docking and simulations, mainly from the perception of bio-molecular interactions.

Correction: Identification of Functional Candidates amongst Hypothetical Proteins of Treponema pallidum ssp. pallidum.

[This corrects the article DOI: 10.1371/journal.pone.0124177.].

Sequence Analysis of Hypothetical Proteins from Helicobacter pylori 26695 to Identify Potential Virulence Factors.

Helicobacter pylori is a Gram-negative bacteria that is responsible for gastritis in human. Its spiral flagellated body helps in locomotion and colonization in the host environment. It is capable of living in the highly acidic environment of the stomach with the help of acid adaptive genes. The genome of H. pylori 26695 strain contains 1,555 coding genes that encode 1,445 proteins. Out of these, 340 proteins are characterized as hypothetical proteins (HP). This study involves extensive analysis of the HPs using an established pipeline which comprises various bioinformatics tools and databases to find out probable functions of the HPs and identification of virulence factors. After extensive analysis of all the 340 HPs, we found that 104 HPs are showing characteristic similarities with the proteins with known functions. Thus, on the basis of such similarities, we assigned probable functions to 104 HPs with high confidence and precision. All the predicted HPs contain representative members of diverse functional classes of proteins such as enzymes, transporters, binding proteins, regulatory proteins, proteins involved in cellular processes and other proteins with miscellaneous functions. Therefore, we classified 104 HPs into aforementioned functional groups. During the virulence factors analysis of the HPs, we found 11 HPs are showing significant virulence. The identification of virulence proteins with the help their predicted functions may pave the way for drug target estimation and development of effective drug to counter the activity of that protein.

Identification of functional candidates amongst hypothetical proteins of Treponema pallidum ssp. pallidum.

Syphilis is a globally occurring venereal disease, and its infection is propagated through sexual contact. The causative agent of syphilis, Treponema pallidum ssp. pallidum, a Gram-negative sphirochaete, is an obligate human parasite. Genome of T. pallidum ssp. pallidum SS14 strain (RefSeq NC_010741.1) encodes 1,027 proteins, of which 444 proteins are known as hypothetical proteins (HPs), i.e., proteins of unknown functions. Here, we performed functional annotation of HPs of T. pallidum ssp. pallidum using various database, domain architecture predictors, protein function annotators and clustering tools. We have analyzed the sequences of 444 HPs of T. pallidum ssp. pallidum and subsequently predicted the function of 207 HPs with a high level of confidence. However, functions of 237 HPs are predicted with less accuracy. We found various enzymes, transporters, binding proteins in the annotated group of HPs that may be possible molecular targets, facilitating for the survival of pathogen. Our comprehensive analysis helps to understand the mechanism of pathogenesis to provide many novel potential therapeutic interventions.