Synthesis, Characterization, Antiglycation Evaluation, Molecular Docking, and ADMET Studies of 4-Thiazolidinone Derivatives.

The design and development of new small-molecule glycation inhibitors are essential for preventing various chronic diseases, including diabetes mellitus, immunoinflammation, cardiovascular, and neurodegenerative diseases. 4-Thiazolidinone or thiazolidine-4-one is a well-known heterocyclic compound with the potential to inhibit the formation of advanced glycation end products. In the present work, we report the synthesis and characterization of four new 5-arylidene 3-cyclopropyl-2-(phenylimino)thiazolidin-4-one (1-4) compounds and their human serum albumin glycation inhibitory activity. One of the compounds 5-(2H-1,3-benzodioxol-5-ylmethylidene)-3-cyclopropyl-2-(phenylimino)-1,3-thiazolidin-4-one (3) showed potent inhibition in the synthesis of initial, intermediary, and final products of glycation reactions. Besides, conformational changes in the α-helix and ß-sheet (due to hyperglycemia) were also found to be reversed upon the addition of (3). Experimental findings were complemented by computational [molecular docking, ADME/Tox, and density functional theory (DFT)] studies. The docking scores of the compounds were in order 1 > 3 > 2 > 4, indicating the importance of the polar group at the 5-arylidene moiety. The results of ADME/Tox and DFT calculations revealed the safe nature of the compounds with high drug-likeness and stability. Overall, we speculate that the results of this study could provide valuable insights into the biological activity of 4-thiazolidinones.

In silico evaluation of 4-thiazolidinone-based inhibitors against the receptor for advanced glycation end products (RAGE).

Non-enzymatic glycation of biomolecules by reducing sugars led to several products, including the advanced glycation end products (AGEs), the accumulation of which has been linked to various life-threatening diseases. The binding of AGEs to their respective protein receptors for advanced glycation end products (RAGE) can initiate a cascade of reactions, which may alter physiological conditions. The present work investigates the potential of 4-thiazolidinones as RAGE inhibitors. We performed an extensive computational study to identify the structural requirements needed to act as RAGE inhibitors. To achieve this goal, 4-thiazolidinone-based compounds available in PubChem, ZINC15, ChEMBL, and ChEBI databases were screened against RAGE (PDB: 4LP5), leading to the identification of top five drug-like candidates with a high binding affinity to RAGE V-domain catalytic region. Drug likeness, absorption, distribution, metabolism, excretion, and toxicity (ADMET) of the top-scoring compounds have been studied and discussed. Global molecular descriptors, chemical reactivity, hardness, softness, etc., have been estimated. Finally, molecular dynamics (MD) simulations at 100 ns were carried out to check the stability and other properties. Overall, we believe that the identified compounds can potentially attenuate RAGE-AGE interactions.Communicated by Ramaswamy H. Sarma.

In silico studies of alkaloids and their derivatives against N-acetyltransferase EIS protein from Mycobacterium tuberculosis.

Antibiotic resistance against Mycobacterium tuberculosis (M.tb.) has been a significant cause of death worldwide. The Enhanced intracellular survival (EIS) protein of the bacteria is an acetyltransferase that multiacetylates aminoglycoside antibiotics, preventing them from binding to the bacterial ribosome. To overcome the EIS-mediated antibiotics resistance of M.tb., we compiled 888 alkaloids and derivatives from five different databases and virtually screened them against the EIS receptor. The compound library was filtered down to 87 compounds, which underwent additional analysis and filtration. Moreover, the top 15 most prominent phytocompounds were obtained after the drug-likeness prediction and ADMET screening. Out of 15, nine compounds confirmed the maximum number of hydrogen bond interactions and reliable binding energies during molecular docking. Additionally, the Molecular dynamics (MD) simulation of nine compounds showed the three most stable complexes, further verified by re-docking with mutated protein. The density functional theory (DFT) calculation was performed to identify the HOMO-LUMO energy gaps of the selected three potential compounds. Finally, our selected top lead compounds i.e., Alkaloid AQC2 (PubChem85634496), Nobilisitine A (ChEbi68116), and N-methylcheilanthifoline (ChEbi140673) demonstrated more favourable outcomes when compared with reference compounds (i.e., 39b and 2i) in all parameters used in this study. Therefore, we anticipate that our findings will help to explore and develop natural compound therapy against multi and extensively drug-resistant strains of M.tb.Communicated by Ramaswamy H. Sarma.

Identification of Novel Tau-Tubulin Kinase 2 Inhibitors Using Computational Approaches.

Tau tubulin kinase 2 (TTBK2) associated with multiple diseases is one of the kinases which phosphorylates tau and tubulin. Numerous efforts have been made to understand the role of TTBK2 in protein folding mechanisms and misfolding behavior. The misfolded protein intermediates form polymers with unwanted aggregation properties that initiate several diseases, including Alzheimer's. The availability of TTBK2 inhibitors can enhance the understanding of the molecular mechanism of action of the kinase and assist in developing novel therapeutics. In the quest for TTBK2 inhibitors, this study focuses on screening two chemical libraries (ChEMBL and ZINC-FDA). The molecular docking, RO5/absorption, distribution, metabolism, and excretion/toxicity, density functional theory, molecular dynamics (MD) simulations, and molecular mechanics with generalized Born and surface area solvation techniques enabled shortlisting of the four most active compounds, namely, ChEMBL1236395, ChEMBL2104398, ChEMBL3427435, and ZINC000000509440. Moreover, 500 ns MD simulation was performed for each complex, which provided valuable insights into the structural changes in the complexes. The relative fluctuation, solvent accessible surface area, atomic gyration, compactness covariance, and free energy landscapes revealed that the compounds could stabilize the TTBK2 protein. Overall, this study would be valuable for the researchers targeting the development of novel TTBK2 inhibitors.

Lead optimization, pharmacophore development and scaffold design of protein kinase CK2 inhibitors as potential COVID-19 therapeutics.

Therapeutic agents being designed against COVID-19 have targeted either the virus directly or the host cellular machinery. A particularly attractive host target is the ubiquitous and constitutively active serine-threonine kinase, Protein kinase CK2 (CK2). CK2 enhances viral protein synthesis by inhibiting the sequestration of host translational machinery as stress granules and assists in viral egression via association with the N-protein at filopodial protrusions of the infected cell. CK2 inhibitors such as Silmitasertib have been proposed as possible therapeutic candidates in COVID-19 infections. The present study aims to optimize Silmitasertib, develop pharmacophore models and design unique scaffolds to modulate CK2. The lead optimization phase involved the generation of compounds structurally similar to Silmitasertib via bioisostere replacement followed by a multi-stage docking approach to identify drug-like candidates. Molecular dynamics (MD) simulations were performed for two promising candidates (ZINC-43206125 and PC-57664175) to estimate their binding stability and interaction. Top scoring candidates from the lead optimization phase were utilized to build ligand-based pharmacophore models. These models were then merged with structure-based pharmacophores (e-pharmacophores) to build a hybrid hypothesis. This hybrid hypothesis was validated against a decoy set and used to screen a diverse kinase inhibitors library to identify favored chemical features in the retrieved actives. These chemical features include; an anion, an aromatic ring and an H-bond acceptor. Based on the knowledge of these features; de-novo scaffold design was carried out which identified phenindiones, carboxylated steroids, macrocycles and peptides as novel scaffolds with the potential to modulate CK2.Communicated by Ramaswamy H. Sarma.

Strand 1A variant in neuroserpin shows increased aggregation and no loss of inhibition: implication in ameliorating polymerization to retain activity.

Neuroserpin (NS) is predominantly expressed in the brain and is the primary inhibitor of tissue plasminogen activator (tPA). NS variants are associated with the neurogenerative disease termed familial encephalopathy with neuroserpin inclusion bodies (FENIB). The disease is characterized by variable age of onset and severity. The reactive center loop (RCL) insertion-based inhibitory mechanism of NS requires a coordinated conformational change leading to a shift in the strands of the ß-sheet A and movement of helix F. Strand 1A is connected to the helix F at its C terminal end and with the strand 2A at its N terminal, both these domain move for accommodating the inserting loop; therefore, a variant that influences their movement may alter the inhibition rates. A molecular dynamic simulation analysis of a H138C NS variant from strand 1A showed a large decrease in conformational fluctuations as compared with wild-type NS. H138 was mutated, expressed, purified and a native-PAGE and transmission electron microscopy (TEM) analysis showed that this variant forms large molecular weight aggregates on a slight increase in temperature. However, a circular dichroism analysis showed its secondary structure to be largely conserved. Surprisingly, its tPA inhibition activity and complex formation remain unhindered even after the site-specific labeling of H138C with Alexa fluor C5 maleimide. Further, a helix F-strand 1A (W154C-H138C) double variant still shows appreciable inhibitory activity. Increasingly, it appears that aggregation and not loss of inhibition is the more likely cause of shutter region-based variants phenotypes, indicating that hindering polymer formation using small molecules may retain inhibitory activity in pathological variants of NS.

Identification of PARP12 Inhibitors By Virtual Screening and Molecular Dynamics Simulations.

Poly [adenosine diphosphate (ADP)-ribose] polymerases (PARPs) are members of a family of 17 enzymes that performs several fundamental cellular processes. Aberrant activity (mutation) in PARP12 has been linked to various diseases including inflammation, cardiovascular disease, and cancer. Herein, a large library of compounds (ZINC-FDA database) has been screened virtually to identify potential PARP12 inhibitor(s). The best compounds were selected on the basis of binding affinity scores and poses. Molecular dynamics (MD) simulation and binding free energy calculation (MMGBSA) were carried out to delineate the stability and dynamics of the resulting complexes. To this end, root means deviations, relative fluctuation, atomic gyration, compactness, covariance, residue-residue contact map, and free energy landscapes were studied. These studies have revealed that compounds ZINC03830332, ZINC03830554, and ZINC03831186 are promising agents against mutated PARP12.

Anti-Fungal Drug Anidulafungin Inhibits SARS-CoV-2 Spike-Induced Syncytia Formation by Targeting ACE2-Spike Protein Interaction.

Drug repositioning continues to be the most effective, practicable possibility to treat COVID-19 patients. The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus enters target cells by binding to the ACE2 receptor via its spike (S) glycoprotein. We used molecular docking-based virtual screening approaches to categorize potential antagonists, halting ACE2-spike interactions by utilizing 450 FDA-approved chemical compounds. Three drug candidates (i.e., anidulafungin, lopinavir, and indinavir) were selected, which show high binding affinity toward the ACE2 receptor. The conformational stability of selected docked complexes was analyzed through molecular dynamics (MD) simulations. The MD simulation trajectories were assessed and monitored for ACE2 deviation, residue fluctuation, the radius of gyration, solvent accessible surface area, and free energy landscapes. The inhibitory activities of the selected compounds were eventually tested in-vitro using Vero and HEK-ACE2 cells. Interestingly, besides inhibiting SARS-CoV-2 S glycoprotein induced syncytia formation, anidulafungin and lopinavir also blocked S-pseudotyped particle entry into target cells. Altogether, anidulafungin and lopinavir are ranked the most effective among all the tested drugs against ACE2 receptor-S glycoprotein interaction. Based on these findings, we propose that anidulafungin is a novel potential drug targeting ACE2, which warrants further investigation for COVID-19 treatment.

Insights into the structure and dynamics of SARS-CoV-2 spike glycoprotein double mutant L452R-E484Q.

The Receptor Binding Domain (RBD) of SARS-CoV-2, located on the S1 subunit, plays a vital role in the virus binding and its entry into the host cell through angiotensin-converting enzyme 2 (ACE2) receptor. Therefore, understanding the dynamic effects of mutants on the SARS-CoV-2 RBD is essential for discovering drugs to inhibit the virus binding and disrupt its entry into the host cells. A recent study reported a double mutant of SARS-CoV-2, L452R-E484Q, located in the RBD region. Thus, this study employed various computational algorithms and methods to understand the structural impact of both individual variants L452R, E484Q, and the double mutant L452R-E484Q on the native RBD of spike glycoprotein. The effects of the mutations on native RBD structure were predicted by in silico algorithms, which predicted changes in the protein structure and function upon the mutations. Subsequently, molecular dynamics (MD) simulations were employed to understand the conformational stability and functional changes on the RBD upon the mutations. The comparative results of MD simulation parameters displayed that the double mutant induces significant conformational changes in the spike glycoprotein RBD, which may alter its biological functions. Supplementary Information: The online version contains supplementary material available at 10.1007/s13205-022-03151-0.

Insights into the structural and dynamical changes of spike glycoprotein mutations associated with SARS-CoV-2 host receptor binding.

Novel Coronavirus or SARS-CoV-2 has received worldwide attention due to the COVID-19 pandemic, which originated in Wuhan, China leading to thousands of deaths to date. The SARS-CoV-2 Spike glycoprotein protein is one of the main focus of COVID-19 related research as it is a structural protein that facilitates its attachment, entry, and infection to the host cells. We have focused our work on mutations in two of the several functional domains in the virus spike glycoprotein, namely, receptor-binding domain (RBD) and heptad repeat 1 (HR1) domain. These domains are majorly responsible for the stability of spike glycoprotein and play a key role in the host cell attachment and infection. In our study, several mutations like R408I, L455Y, F486L, Q493N, Q498Y, N501T of RBD (319-591), and A930V, D936Y of HR1 (912-984) have been studied to examine its role on the spike glycoprotein native structure. Comparisons of MD simulations in the WT and mutants revealed a significant de-stabilization effect of the mutations on RBD and HR1 domains. We have investigated the impact of mapped mutations on the stability of the spike glycoprotein, before binding to the receptor, which may be consequential to its binding properties to the receptor and other ligands.Communicated by Ramaswamy H. Sarma.

Targeting SARS-CoV-2 nucleocapsid oligomerization: Insights from molecular docking and molecular dynamics simulations.

The outbreak of COVID-19 caused by SARS-CoV-2 virus continually led to infect a large population worldwide. Currently, there is no specific viral protein-targeted therapeutics. The Nucleocapsid (N) protein of the SARS-CoV-2 virus is necessary for viral RNA replication and transcription. The C-terminal domain of N protein (CTD) involves in the self-assembly of N protein into a filament that is packaged into new virions. In this study, the CTD (PDB ID: 6WJI) was targeted for the identification of possible inhibitors of oligomerization of N protein. Herein, multiple computational approaches were employed to explore the potential mechanisms of binding and inhibitor activity of five antiviral drugs toward CTD. The five anti-N drugs studied in this work are 4E1RCat, Silmitasertib, TMCB, Sapanisertib, and Rapamycin. Among the five drugs, 4E1RCat displayed highest binding affinity (-10.95 kcal/mol), followed by rapamycin (-8.91 kcal/mol), silmitasertib (-7.89 kcal/mol), TMCB (-7.05 kcal/mol), and sapanisertib (-6.14 kcal/mol). Subsequently, stability and dynamics of the protein-drug complex were examined with molecular dynamics (MD) simulations. Overall, drug binding increases the stability of the complex with maximum stability observed in the case of 4E1RCat. The CTD-drug complex systems behave differently in terms of the free energy landscape and showed differences in population distribution. Overall, the MD simulation parameters like RMSD, RMSF, Rg, hydrogen bonds analysis, PCA, FEL, and DCCM analysis indicated that 4E1RCat and TMCB complexes were more stable as compared to silmitasertib and sapanisertib and thus could act as effective drug compounds against CTD.Communicated by Ramaswamy H. Sarma.

Structural stability predictions and molecular dynamics simulations of RBD and HR1 mutations associated with SARS-CoV-2 spike glycoprotein.

The COVID-19 pandemic is caused by human transmission and infection of Severe Acute Respiratory Syndrome Corona Virus-2 (SARS-CoV-2). There is no trusted drug against the virus; hence, efforts are on discovering novel inhibitors for the virus. The entry of a SARS-CoV-2 virus particle into a host cell is initiated by its spike glycoprotein and host Angiotensin-Converting Enzyme 2 (ACE2) receptor interaction. Spike glycoprotein domains, namely, the Receptor Binding Domain (RBD) and Heptad Repeat (HR) domains, are essential for this activity. We have studied the impact of mutations such as A348T, N354D, D364Y, G476S, V483A, S494D in the RBD (319-591), and S939F, S940T, T941A, S943P (912-984) in the HR1 domains of spike glycoprotein. Summarily, we utilized the computational screening algorithms to rank the deleterious, damaging and disease-associated spike glycoprotein mutations. Subsequently, to understand the changes in conformation, flexibility and function of the spike glycoprotein mutants, Molecular Dynamics (MD) simulations were performed. The computational predictions and analysis of the MD trajectories suggest that the RBD and HR1 mutations induce significant phenotypic effects on the pre-binding spike glycoprotein structure, which are presumably consequential to its binding to the receptor and provides lead to design inhibitors against the binding.Communicated by Ramaswamy H. Sarma.

Targeting cathepsins: A potential link between COVID-19 and associated neurological manifestations.

Many studies have shown that the lysosomal cathepsins, especially cathepsins B/L (CTSB/L) are required for SARS-CoV-2 entry into host cells. Lysosomal proteases, cathepsins are indispensable for normal health and are involved in several brain disorders occurring at different development age periods. On the other hand, it has been well known that COVID-19 infection is largely associated with several neurological disorders. Taken together these findings and given the high levels of expression of CTSB/L in the brain, we here proposed a reasonable hypothesis about the involvement of CTSB/L in the neurological manifestations linked to COVID-19. Pharmacological inhibitions of the CTSB/L could be a potential therapeutic target to block the virus entry as well as to mitigate the brain disorders. To this end, we utilized the network-based drug repurposing analyses to identify the possible drugs that can target CTSB/L. This study identifies the molecules like cyclosporine, phenytoin, and paclitaxel as potential drugs with binding ability to the CTSB/L. Further, we have performed molecular docking and all-atom molecular dynamics (MD) simulations to investigate the stability of CTSL-drug complexes. The results showed strong and stable binding of drugs with CTSL.

The structural, functional, and dynamic effect of Tau tubulin kinase1 upon a mutation: A neuro-degenerative hotspot.

Alzheimer's disease (AD) is a progressive disorder that causes brain cells to degenerate and die. AD is one of the common causes of dementia that leads to a decline in thinking, behavioral and social skills that disrupts a person's ability to function independently. Tau-tubulin kinase1 (TTBK1) is a crucial disease regulating AD protein, which is majorly responsible for the phosphorylation and accumulation of tau protein at specific Serine/Threonine residues found in paired helical filaments, suggesting its role in tauopathy. TTBK1 involvement in many diseases and the restricted expression of TTBK1 to the central nervous system (CNS) makes TTBK1 an attractive therapeutic target for tauopathies. The genetic variations in TTBK1 are primarily involved in the TTBK1 pathogenesis. This study highlighted the destabilizing, damaging and deleterious effect of the mutation R142Q on TTBK1 structure through computational predictions and molecular dynamics simulations. The protein deviation, fluctuations, conformational dynamics, solvent accessibility, hydrogen bonding, and the residue-residue mapping confirmed the mutant effect to cause structural aberrations, suggesting overall destabilization due to the protein mutation. The presence of well-defined free energy minima was observed in TTBK1-wild type, as opposed to that in the R142Q mutant, reflecting structural deterioration. The overall findings from the study reveal that the presence of R142Q mutation on TTBK1 is responsible for the structural instability, leading to disruption of its biological functions. The mutation could be used as future diagnostic markers in treating AD.

Atomic Resolution Homology Models and Molecular Dynamics Simulations of Plasmodium falciparum Tubulins.

Microtubules are tubulin polymers present in the eukaryotic cytoskeleton essential for structural stability and cell division that are also roadways for intracellular transport of vesicles and organelles. In the human malaria parasite Plasmodium falciparum, apart from providing structural stability and cell division, microtubules also facilitate important biological activities crucial for parasite survival in hosts, such as egression and motility. Hence, parasite structures and processes involving microtubules are among the most important drug targets for discovering much-needed novel Plasmodium inhibitors. The current study aims to construct reliable and high-quality 3D models of α-, ß-, and γ-tubulins using various modeling techniques. We identified a common binding pocket specific to Plasmodium α-, ß-, and γ-tubulins. Molecular dynamics simulations confirmed the stability of the Plasmodium tubulin 3D structures. The models generated in the present study may be used for protein-protein and protein-drug interaction investigations targeted toward designing malaria parasite tubulin-specific inhibitors.

Detection of truncated isoforms of human neuroserpin lacking the reactive center loop: Implications in noninhibitory role.

Neuroserpin is a serine protease inhibitor expressed mainly in the brain and at low levels in other tissues like the kidney, testis, heart, and spinal cord. It is involved in the inhibition of tissue plasminogen activator (tPA), plasmin, and to a lesser extent, urokinase-type plasminogen (uPA). Neuroserpin has also been shown to plays noninhibitory roles in the regulation of N-cadherin-mediated cell adhesion. It is involved in neuroprotection from seizure and stroke through tPA-mediated inhibition and also through its other protease targets. Mutations in critical domains of neuroserpin lead to its polymerization and neuronal death. In this study, a novel truncated isoform of human neuroserpin was identified in the brain and liver, which was confirmed by reverse transcriptase-PCR and DNA sequencing using exon-specific primers. Structural characterization of novel isoform using MD simulations studies indicated that it lacks the reactive center loop (RCL) but largely maintains its secondary structure fold. The novel truncated variant was cloned, expressed, and purified. A comparative intrinsic fluorescence and 4,4'-bis-1-anilino naphthalene 8-sulfonate studies revealed a decrease in fluorescence emission intensity and a more exposed hydrophobic surface as compared to the reported isoform. However, the novel isoform has lost its ability for tPA inhibition and complex formation. The absence of RCL indicates a noninhibitory role for the truncated isoform, prompting a detailed search and identification of two smaller isoforms in the human brain. With indications of the noninhibitory role of neuroserpin, identifying novel isoforms that appear to be without the tPA recognition domain is significant.

Simultaneous Inhibition of SARS-CoV-2 Entry Pathways by Cyclosporine.

The COVID-19 pandemic caused by SARS-CoV-2 represents a global public health emergency. The entry of SARS-CoV-2 into host cells requires the activation of its spike protein by host cell proteases. The serine protease, TMPRSS2, and cysteine proteases, Cathepsins B/L, activate spike protein and enable SARS-CoV-2 entry to the host cell through two completely different and independent pathways. Therefore, inhibiting either TMPRSS2 or cathepsin B/L may not sufficiently block the virus entry. We here hypothesized that simultaneous targeting of both the entry pathways would be more efficient to block the virus entry rather than targeting the entry pathways individually. To this end, we utilized the network-based drug repurposing analyses to identify the possible common drugs that can target both the entry pathways. This study, for the first time, reports the molecules like cyclosporine, calcitriol, and estradiol as candidate drugs with the binding ability to the host proteases, TMPRSS2, and cathepsin B/L. Next, we analyzed drug-gene and gene-gene interaction networks using 332 human targets of SARS-CoV-2 proteins. The network results indicate that, out of 332 human proteins, cyclosporine interacts with 216 (65%) proteins. Furthermore, we performed molecular docking and all-atom molecular dynamics (MD) simulations to explore the binding of drug with TMPRSS2 and cathepsin L. The molecular docking and MD simulation results showed strong and stable binding of cyclosporine A (CsA) with TMPRSS2 and CTSL genes. The above results indicate cyclosporine as a potential drug molecule, as apart from interacting with SARS-CoV-2 entry receptors, it also interacts with most of SARS-CoV-2 target host genes; thus it could potentially interfere with functions of SARS-CoV-2 proteins in human cells. We here also suggest that these antiviral drugs alone or in combination can simultaneously target both the entry pathways and thus can be considered as a potential treatment option for COVID-19.

In-silico profiling and structural insights into the impact of nSNPs in the P. falciparum acetyl-CoA transporter gene to understand the mechanism of drug resistance in malaria.

The continuous emergence of resistance to the available drugs poses major constraints in the development of effective therapeutics against malaria. Malaria drug resistance has been attributed to be the manifestation of numerous factors. For example, mutations in the parasite transporter protein acetyl-CoA transporter (Pfact) can remarkably affect its uptake affinity for a drug molecule against malaria, and hence enhance its susceptibility to resistance. To identify major contributors to its loss of functionality, we have thoroughly scrutinized eight such recently reported resistant mutants, via in-silico tools in terms of alterations in different properties. We performed molecular dynamics simulations of the selected Pfact mutants to gain deeper insights into the structural perturbation and dynamicity. Comparison of residue interaction network map of mutants with that of Wild type (WT) protein suggests structural changes as a result of the mutation(s) that translate into the weakening of intra-protein interactions, especially around the drug binding pocket. This, in turn, diminishes the affinity of drug molecules towards the binding site, which was validated by docking analysis. Finally, collating all the observations, we have delineated R108K mutant to deviate the most from WT protein, which, intriguingly suggests that replacing an amino acid with another of similar nature can even translate into greater functional effects as those with dissimilar substitutions.Communicated by Ramaswamy H. Sarma.

Investigation of conformational dynamics of Tyr89Cys mutation in protection of telomeres 1 gene associated with familial melanoma.

Protection of telomeres 1 (POT1) is a component of the shelterin complex which is crucial for the regulation of telomere length and maintenance. Many naturally occurring mutations in the POT1 gene have been found to be associated with cardiac angiosarcoma, glioma, familial melanoma, and chronic lymphocytic leukemia. In particular, Y89C is a naturally occurring mutation of POT1, responsible for familial melanoma, and the molecular basis of this mutation is unexplored. In this study, we have extensively analyzed the structure of WT and Y89C mutant of POT1 to see the change in the conformational dynamics, free energy landscape, molecular motions and configurational frustration using molecular dynamics (MD) and other bioinformatics approaches. Y89C mutation shows a significant change in the backbone orientation, compactness, residual fluctuation, solvent accessibility, and hydrogen bonding, suggesting an overall destabilization of the protein structure. Besides, essential dynamics, conformation, magnitude, direction of motion and frustration analysis further suggesting the structural loss in POT1 due to Y89C mutation. Free energy landscape analysis also indicates the presence of a single well-defined free-energy minima in case of WT compared to multiple wells defined free energy minima observed in Y89C, clearly suggesting that this mutation leads to reduce the stability of POT1. This study possibly provides a valuable path to understand the molecular basis of Y89C-mediated development of familial melanoma.Communicated by Ramaswamy H. Sarma.

A molecular journey to check the conformational dynamics of tau tubulin kinase 2 mutations associated with Alzheimer's disease.

Proteins are one of the most vital components of biological functions. Proteins have evolutionarily conserved structures as the shape and folding pattern predominantly determine their function. Considerable research efforts have been made to study the protein folding mechanism. The misfolding of protein intermediates of large groups form polymers with unwanted aggregates that may initiate various diseases. Amongst the diseases caused by misfolding of proteins, Alzheimer's disease (AD) is one of the most prevalent neuro-disorders which has a worldwide impact on human health. The disease is associated with several vital proteins and single amino acid mutations. Tau tubulin kinase 2 (TTBK2) is one of the kinases which is known to phosphorylate tau and tubulin. The literature strongly supports that the mutations-K50E, D163A, R181E, A184E and K143E are associated with multiple important cellular processes of TTBK2. In this study, to understand the molecular basis of the functional effects of the mutations, we have performed structural modeling for TTBK2 and its mutations, using computational prediction algorithms and Molecular Dynamics (MD) simulations. The MD simulations highlighted the impact of the mutations on the Wild Type (WT) by the conformational dynamics, Free Energy Landscape (FEL) and internal molecular motions, indicating the structural de-stabilization which may lead to the disruption of its biological functions. The destabilizing effect of TTBK2 upon mutations provided valuable information about individuals carrying this mutant which could be used as a diagnostic marker in AD.