Insights into structural dynamics of allosteric binding sites in HCV RNA-dependent RNA polymerase

Molecular analysis of dUTPase of Helicobacter pylori for identification of novel inhibitors using in silico studies

The human gastric pathogen Helicobacter pylori chronically affects the gastric mucosal layer of approximately half of world's population. The emergence of resistant strains urges the need for identification of novel and selective drug against new molecular targets. A ubiquitous enzyme, Deoxyuridine 5'-triphosphate nucleotidohydrolase (dUTPase), is considered as first line of defense against uracil mis-incorporation into DNA, and essential for genome integrity. Lack of dUTPase triggers an elevated recombination frequency, DNA breaks and ultimately cell death. Hence, dUTPase can be considered as a promising target for development of novel lead inhibitor compounds in H. pylori treatment. Herein, we report the generation of three-dimensional model of the target protein using comparative modelling and its validation. To identify dUTPase inhibitors, a high throughput virtual screening approach utilizing Knowledge-based inhibitors and DrugBank database was implemented. Top ranked compounds were scrutinized based on investigations of the protein-ligand interaction fingerprints, molecular interaction maps and binding affinities and the drug potentiality. The best ligands were studied further for complex stability and intermolecular interaction profiling with respect to time under 100 ns classical molecular dynamic stimulation, establishing significant stability in dynamic states as observed from RMSD and RMSF parameters and interactions with the catalytic site residues. The binding free energy calculation computed using MM-GBSA method from the MD simulation trajectories demonstrated that our molecules possess strong binding affinity towards the Helicobacter pylori dUTPase protein. We conclude that our proposed molecules may be potential lead molecules for effective inhibition against the H. pylori dUTPase protein subject to experimental validation.Communicated by Ramaswamy H. Sarma.

Discovery of non-nucleoside oxindole derivatives as potent inhibitors against dengue RNA-dependent RNA polymerase

A series of thiazole linked Oxindole-5-Sulfonamide (OSA) derivatives were designed as inhibitors of RNA-dependent RNA polymerase (RdRp) activity of Dengue virus. These were synthesized and then evaluated for their efficacy in ex-vivo virus replication assay using human cell lines. Among 20 primary compounds in the series, OSA-15 was identified as a hit. A series of analogues were synthesized by replacing the difluoro benzyl group of OSA-15 with different substituted benzyl groups. The efficacy of OSA-15derivatives was less than that of the parent compound, except OSA-15-17, which has shown improved efficacy than OSA-15. The further optimization was carried out by adding dimethyl (DM) groups to both the sulfonamide and oxindole NH's to produce OSA-15-DM and OSA-15-17-DM. These two compounds were showing no detectable cytotoxicity and the latter was more efficacious. Further, both these compounds were tested for inhibition in all the serotypes of the Dengue virus using an ex-vivo assay. The EC₅₀ of OSA-15-17-DM was observed in a low micromolar range between 2.5 and 5.0 µg/ml. Computation docking and molecular dynamics simulation studies confirmed the binding of identified hits to DENV RdRp. OSA15-17-DM blocks the RNA entrance and elongation site for their biological activity with high binding affinity. Overall, the identified oxindole derivatives are novel compounds that can inhibit Dengue replication, working as non-nucleoside inhibitors (NNI) to explore as anti-viral RdRp activity.

Conformational Characterization of the Co-Activator Binding Site Revealed the Mechanism to Achieve the Bioactive State of FXR

FXR bioactive states are responsible for the regulation of metabolic pathways, which are modulated by agonists and co-activators. The synergy between agonist binding and ‘co-activator’ recruitment is highly conformationally driven. The characterization of conformational dynamics is essential for mechanistic and therapeutic understanding. To shed light on the conformational ensembles, dynamics, and structural determinants that govern the activation process of FXR, molecular dynamic (MD) simulation is employed. Atomic insights into the ligand binding domain (LBD) of FXR revealed significant differences in inter/intra molecular bonding patterns, leading to structural anomalies in different systems of FXR. The sole presence of an agonist or ‘co-activator’ fails to achieve the essential bioactive conformation of FXR. However, the presence of both establishes the bioactive conformation of FXR as they modulate the internal wiring of key residues that coordinate allosteric structural transitions and their activity. We provide a precise description of critical residue positioning during conformational changes that elucidate the synergy between its binding partners to achieve an FXR activation state. Our study offers insights into the associated modulation occurring in FXR at bound and unbound forms. Thereafter, we also identified hot-spots that are critical to arrest the activation mechanism of FXR that would be helpful for the rational design of its agonists.

Molecular Dynamics Simulations Reveal the Interaction Fingerprint of Remdesivir Triphosphate Pivotal in Allosteric Regulation of SARS-CoV-2 RdRp

The COVID-19 pandemic has now strengthened its hold on human health and coronavirus' lethal existence does not seem to be going away soon. In this regard, the optimization of reported information for understanding the mechanistic insights that facilitate the discovery towards new therapeutics is an unmet need. Remdesivir (RDV) is established to inhibit RNA-dependent RNA polymerase (RdRp) in distinct viral families including Ebola and SARS-CoV-2. Therefore, its derivatives have the potential to become a broad-spectrum antiviral agent effective against many other RNA viruses. In this study, we performed comparative analysis of RDV, RMP (RDV monophosphate), and RTP (RDV triphosphate) to undermine the inhibition mechanism caused by RTP as it is a metabolically active form of RDV. The MD results indicated that RTP rearranges itself from its initial RMP-pose at the catalytic site towards NTP entry site, however, RMP stays at the catalytic site. The thermodynamic profiling and free-energy analysis revealed that a stable pose of RTP at NTP entrance site seems critical to modulate the inhibition as its binding strength improved more than its initial RMP-pose obtained from docking at the catalytic site. We found that RTP not only occupies the residues K545, R553, and R555, essential to escorting NTP towards the catalytic site, but also interacts with other residues D618, P620, K621, R624, K798, and R836 that contribute significantly to its stability. From the interaction fingerprinting it is revealed that the RTP interact with basic and conserved residues that are detrimental for the RdRp activity, therefore it possibly perturbed the catalytic site and blocked the NTP entrance site considerably. Overall, we are highlighting the RTP binding pose and key residues that render the SARS-CoV-2 RdRp inactive, paving crucial insights towards the discovery of potent inhibitors.

Binding mode characterization of 13b in the monomeric and dimeric states of SARS-CoV-2 main protease using molecular dynamics simulations

The main protease, M^pro/3CL^pro, plays an essential role in processing polyproteins translated from viral RNA to produce functional viral proteins and therefore serve as an attractive target for discovering COVID-19 therapeutics. The availability of both monomer and dimer crystal bound with a common ligand, '13b' (α-ketoamide inhibitor), opened up opportunities to understand the M^pro mechanism of action. A comparative analysis of both forms of M^pro was carried out to elucidate the binding site architectural differences in the presence and absence of '13b'. Molecular dynamics simulations suggest that the presence of '13b' enhances the stability of M^pro than the unbound APO form. The N- and C- terminals of both the protomers stabilize each other, and making it's interface essential for the active form of M^pro. In comparison to monomer, the relatively high affinity of '13b' is gained in dimer pocket due to the high stability of the pocket by the interaction of S1 residue of chain B with residues F140, E166 and H172 of chain A, which is absent in monomer. The comprehensive essential dynamics, protein structure network analysis and thermodynamic profiling highlight the hot-spots, pivotal in molecular recognition process at protein-ligand and protein-protein interaction levels, cross-validated through computational alanine scanning study. A comparative description of '13b' binding mechanism in both forms illustrates valuable insights into the inhibition mechanism and the selection of critical residues suitable for the structure-based approaches for the identification of more potent M^pro inhibitors.Communicated by Ramaswamy H. Sarma.

Interplay among Structural Stability, Plasticity, and Energetics Determined by Conformational Attuning of Flexible Loops in PD-1

The dynamics and plasticity of the PD-1/PD-L1 axis are the bottlenecks for the discovery of small-molecule antagonists to perturb this interaction interface significantly. Understanding the process of this protein-protein interaction (PPI) is of fundamental biological interest in structure-based drug designing. Food and Drug Administration (FDA)-approved anti-PD-1 monoclonal antibodies (mAbs) are the first-in-class with distinct binding modes to access this axis clinically; however, their mechanistic aspects remain elusive. Here, we have unveiled the interactive interfaces with PD-L1 and mAbs to investigate the native plasticity of PD-1 at global (structural and dynamical) and local (residue side-chain orientations) levels. We found that the structural stability and coordinated C_α movements are increased in the presence of PD-1's binding partners. The rigorous analysis of these PPIs using computational biophysical approaches revealed PD-1's intrinsic plasticity, its concerted loops' movement (BC, FG, and CC'), distal side-chain motions, and the thermodynamic landscape, which are perturbed remarkably from its unbound to bound states. Based on intra-/inter-residues' contact networks and energetics, the hot-spots have been identified that were found to be essential to arrest the dynamical motions of PD-1 significantly for the rational design of therapeutic agents by mimicking the mAbs mechanism.

Tracing the footsteps of autophagy in computational biology

Autophagy plays a crucial role in maintaining cellular homeostasis through the degradation of unwanted materials like damaged mitochondria and misfolded proteins. However, the contribution of autophagy toward a healthy cell environment is not only limited to the cleaning process. It also assists in protein synthesis when the system lacks the amino acids’ inflow from the extracellular environment due to diet consumptions. Reduction in the autophagy process is associated with diseases like cancer, diabetes, non-alcoholic steatohepatitis, etc., while uncontrolled autophagy may facilitate cell death. We need a better understanding of the autophagy processes and their regulatory mechanisms at various levels (molecules, cells, tissues). This demands a thorough understanding of the system with the help of mathematical and computational tools. The present review illuminates how systems biology approaches are being used for the study of the autophagy process. A comprehensive insight is provided on the application of computational methods involving mathematical modeling and network analysis in the autophagy process. Various mathematical models based on the system of differential equations for studying autophagy are covered here. We have also highlighted the significance of network analysis and machine learning in capturing the core regulatory machinery governing the autophagy process. We explored the available autophagic databases and related resources along with their attributes that are useful in investigating autophagy through computational methods. We conclude the article addressing the potential future perspective in this area, which might provide a more in-depth insight into the dynamics of autophagy.

Identification of potential anti-TMPRSS2 natural products through homology modelling, virtual screening and molecular dynamics simulation studies

Recent outbreak of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has led to a pandemic of COVID-19. The absence of a therapeutic drug and vaccine is causing severe loss of life and economy worldwide. SARS-CoV and SARS-CoV-2 employ the host cellular serine protease TMPRSS2 for spike (S) protein priming for viral entry into host cells. A potential way to reduce the initial site of SARS-CoV-2 infection may be to inhibit the activity of TMPRSS2. In the current study, the three-dimensional structure of TMPRSS2 was generated by homology modelling and subsequently validated with a number of parameters. The structure-based virtual screening of Selleckchem database was performed through 'Virtual Work Flow' (VSW) to find out potential lead-like TMPRSS2 inhibitors. Camostat and bromhexine are known TMPRSS2 inhibitor drugs, hence these were used as control molecules throughout the study. Based on better dock score, binding-free energy and binding interactions compared to the control molecules, six molecules (Neohesperidin, Myricitrin, Quercitrin, Naringin, Icariin, and Ambroxol) were found to be promising against the TMPRSS2. Binding interactions analysis revealed a number of significant binding interactions with binding site amino residues of TMPRSS2. The all-atoms molecular dynamics (MD) simulation study indicated that all proposed molecules retain inside the receptor in dynamic states. The binding energy calculated from the MD simulation trajectories also favour the strong affinity of the molecules towards the TMPRSS2. Proposed molecules belong to the bioflavonoid class of phytochemicals and are reported to possess antiviral activity, our study indicates their possible potential for application in COVID-19.Communicated by Ramaswamy H. Sarma.

Bile acids mediated potential functional interaction between FXR and FATP5 in the regulation of Lipid Metabolism

Perturbation in lipid homeostasis is one of the major bottlenecks in metabolic diseases, especially Non-alcoholic Fatty Liver Disease (NAFLD), which has emerged as a leading global cause of chronic liver disease. The bile acids (BAs) and their derivatives exert a variety of metabolic effects through complex and intertwined pathways, thus becoming the attractive target for metabolic syndrome treatment. To modulate the lipid homeostasis, the role of BAs, turn out to be paramount as it is essential for the absorption, transport of dietary lipids, regulation of metabolic enzymes and transporters that are essential for lipid modulation, flux, and excretion. The synthesis and transport of BAs (conjugated and unconjugated) is chiefly controlled by nuclear receptors and the uptake of long-chain fatty acids (LCFA) and BA conjugation via transporters. Among them, from in-vivo studies, farnesoid X receptor (FXR) and liver-specific fatty acid transport protein 5 (FATP5) have shown convincing evidence for their key roles in lipid homeostasis and reversal of fatty liver disease substantially. BAs have a wider range of biological effects as they are identified as modulators for FXR and FATP5 both and therefore hold a significant promise for altering the lipid content in the treatment of a metabolic disorder. BAs also have received noteworthy interest in drug delivery research due to its peculiar physicochemical properties and biocompatibility. Here, we are highlighting the connecting possibility of BAs as an agonist for FXR and antagonist for FATP5, paving an avenue to target them for designing synthetic small molecules for lipid homeostasis.

Fluoresceinated Aminohexanol Tethered Inositol Hexakisphosphate: Studies on Arabidopsis thaliana and Drosophila melanogaster and Docking with 2P1M Receptor

Inositol hexakisphosphate (InsP6; phytic acid) is considered as the second messenger and plays a very important role in plants, animals, and human beings. It is the principal storage form of phosphorus in many plant tissues, especially in dry fruits, bran, and seeds. The resulting anion is a colorless species that plays a critical role in nutrition and is believed to cure many diseases. A fluoresceinated aminohexanol tethered inositol hexakisphosphate (III) had been synthesized earlier involving many complicated steps. We describe here a simple two-step synthesis of (III) and its characterization using different techniques such as matrix-assisted laser desorption ionization mass spectrometry, tandem mass spectrometry, and Fourier transform infrared, ultraviolet-visible, ultraviolet-fluorescence, 1H nuclear magnetic resonance (NMR), and two-dimensional NMR spectroscopies. The effect of (III) has been investigated in the model systems, Arabidopsis thaliana and Drosophila melanogaster. Using Schrodinger software, computational studies on the binding of (III) with the protein 2P1M (Auxin-receptor TIR1-adaptor ASK1 complex) has revealed strong binding propensity with this compound. These studies on the fluoresceinated tethered phytic acid could have far reaching implications on its efficacy for human health and treatment of diseases (cancer/tumor and glioblastoma) and for understanding phosphorous recycling in the environment, especially for plant systems.

Unique and exclusive peptide signatures directly identify intrinsically disordered proteins from sequences without structural information

Intrinsically disordered proteins are now widely accepted to play crucial roles in biological functions. Identification of signatures of intrinsic disorder is one of the key steps towards building a proper repertoire for their occurrence in proteomes. In this work, systematic computational synthesis of a library of all possible (3368400) dipeptides, tripeptides, tetrapeptides and pentapeptides using the natural 20 amino acids allowed us to identify 36 unique tetrapeptides present exclusively in intrinsically disordered proteins and absent in the complete primary sequence space of naturally occurring structured proteins. Further, out of more than 530000 known naturally occurring primary sequences without any structural information, 1349 sequences contain the above identified unique signatures of intrinsic disorder. These sequences, having cellular functions varying from housekeeping to metabolic to transport, more than double the number of the currently known intrinsically disordered proteins. On similar lines, we report that 26577 pentapeptide signatures exclusive to intrinsically disordered proteins, and absent in naturally occurring structured proteins, identify ∼50% of more than half-a-million curated protein sequences without structural information to be intrinsically disordered. The results reported are a major leap forward in exploring functional manifestations of intrinsically disordered proteins.Communicated by Ramaswamy H. Sarma.