Faculty Publications
Permanent URI for this communityhttps://idr.nitk.ac.in/handle/123456789/18736
Publications by NITK Faculty
Browse
9 results
Search Results
Item Computational insights into factor affecting the potency of diaryl sulfone analogs as Escherichia coli dihydropteroate synthase inhibitors(Elsevier Ltd, 2019) Das, B.K.; PV, P.; Chakraborty, D.Dihydropteroate synthase (DHPS) is an alluring target for designing novel drug candidates to prevent infections caused by pathogenic Escherichia coli strains. Diaryl Sulfone (SO) compounds are found to inhibit DHPS competitively with respect to the substrate pABA (p-aminobenzoate). The extra aromatic ring of diaryl sulfone compounds found to stabilize them in highly flexible pABA binding loops. In this present study, a statistically significant 3D-QSAR model was developed using a data set of diaryl sulfone compounds. The favourable and unfavourable contributions of substitutions in sulfone compounds were illustrated by contour plot obtained from the developed 3D-QSAR model. Molecular docking calculations were performed to investigate the putative binding mode of diaryl sulfone compounds at the catalytic pocket. DFT calculations were carried out using SCF approach, B3LYP- 6-31 G (d) basis set to compute the HOMO, LUMO energies and their respective location at pABA binding pocket. Further, the developed model was validated by FEP (Free Energy Perturbation) calculations. The calculated relative free energy of binding between the highly potent and less potent sulfone compound was found to be ?3.78 kcal/ mol which is comparable to the experimental value of ?5.85 kcal/mol. A 10 ns molecular dynamics simulation of inhibitor and DHPS confirmed its stability at pABA catalytic site. Outcomes of the present work provide deeper insight in designing novel drug candidates for pathogenic Escherichia coli strains. © 2018 Elsevier LtdItem Effect of hydrophobic and hydrogen bonding interactions on the potency of ß-alanine analogs of G-protein coupled glucagon receptor inhibitors(John Wiley and Sons Inc. P.O.Box 18667 Newark NJ 07191-8667, 2020) Venugopal, P.P.; Das, B.K.; Soorya, E.; Chakraborty, D.G-protein coupled glucagon receptors (GCGRs) play an important role in glucose homeostasis and pathophysiology of Type-II Diabetes Mellitus (T2DM). The allosteric pocket located at the trans-membrane domain of GCGR consists of hydrophobic (TM5) and hydrophilic (TM7) units. Hydrophobic interactions with the amino acid residues present at TM5, found to facilitate the favorable orientation of antagonist at GCGR allosteric pocket. A statistically robust and highly predictive 3D-QSAR model was developed using 58 ?-alanine based GCGR antagonists with significant variation in structure and potency profile. The correlation coefficient (R2) and cross-validation coefficient (Q2) of the developed model were found to be 0.9981 and 0.8253, respectively at the PLS factor of 8. The analysis of the favorable and unfavorable contribution of different structural features on the glucagon receptor antagonists was done by 3D-QSAR contour plots. Hydrophobic and hydrogen bonding interactions are found to be main dominating non-bonding interactions in docking studies. Presence of highest occupied molecular orbital (HOMO) in the polar part and lowest unoccupied molecular orbital (LUMO) in the hydrophobic part of antagonists leads to favorable protein-ligand interactions. Molecular mechanics/generalized born surface area (MM/GBSA) calculations showed that van der Waals and nonpolar solvation energy terms are crucial components for thermodynamically stable binding of the inhibitors. The binding free energy of highly potent compound was found to be ?63.475 kcal/mol; whereas the least active compound exhibited binding energy of ?41.097 kcal/mol. Further, five 100 ns molecular dynamics simulation (MD) simulations were done to confirm the stability of the inhibitor-receptor complex. Outcomes of the present study can serve as the basis for designing improved GCGR antagonists. © 2019 Wiley Periodicals, Inc.Item Epitope-Based Potential Vaccine Candidate for Humoral and Cell-Mediated Immunity to Combat Severe Acute Respiratory Syndrome Coronavirus 2 Pandemic(American Chemical Society, 2020) Das, B.K.; Chakraborty, D.The emergence of severe acute respiratory syndrome from novel Coronavirus (SARS-CoV-2) has put an immense pressure worldwide where vaccination is believed to be an efficient way for developing hard immunity. Herein, we employ immunoinformatic tools to identify B-cell, T-cell epitopes associated with the spike protein of SARS-CoV-2, which is important for genome release. The results showed that the highly immunogenic epitopes located at the stalk part are mostly conserved compared to the receptor binding domain (RDB). Further, two vaccine candidates were computationally modeled from the linear B-cell, T-cell epitopes. Molecular docking reveals the crucial interactions of the vaccines with immune-receptors, and their stability is assessed by MD simulation studies. The chimeric vaccines showed remarkable binding affinity toward the immune cell receptors computed by the MM/PBSA method. van der Waals and electrostatic interactions are found to be the dominant factors for the stability of the complexes. The molecular-level interaction obtained from this study may provide deeper insight into the process of vaccine development against the pandemic of COVID-19. © 2020 American Chemical Society.Item Molecular mechanism of inhibition of COVID-19 main protease by ?-adrenoceptor agonists and adenosine deaminase inhibitors using in silico methods(Taylor and Francis Ltd., 2022) Venugopal, P.P.; Chakraborty, D.Novel coronavirus (COVID-19) responsible for viral pneumonia which emerged in late 2019 has badly affected the world. No clinically proven drugs are available yet as the targeted therapeutic agents for the treatment of this disease. The viral main protease which helps in replication and transcription inside the host can be an effective drug target. In the present study, we aimed to discover the potential of ?-adrenoceptor agonists and adenosine deaminase inhibitors which are used in asthma and cancer/inflammatory disorders, respectively, as repurposing drugs against protease inhibitor by ligand-based and structure-based virtual screening using COVID-19 protease-N3 complex. The AARRR pharmacophore model was used to screen a set of 22,621 molecules to obtain hits, which were subjected to high-throughput virtual screening. Extra precision docking identified four top-scored molecules such as +/?-fenoterol, FR236913 and FR230513 with lower binding energy from both categories. Docking identified three major hydrogen bonds with Gly143, Glu166 and Gln189 residues. 100 ns MD simulation was performed for four top-scored molecules to analyze the stability, molecular mechanism and energy requirements. MM/PBSA energy calculation suggested that van der Waals and electrostatic energy components are the main reasons for the stability of complexes. Water-mediated hydrogen bonds between protein-ligand and flexibility of the ligand are found to be responsible for providing extra stability to the complexes. The insights gained from this combinatorial approach can be used to design more potent and bio-available protease inhibitors against novel coronavirus. Communicated by Ramaswamy H. Sarma. © 2020 Informa UK Limited, trading as Taylor & Francis Group.Item Deciphering the competitive inhibition of dihydropteroate synthase by 8 marcaptoguanine analogs: enhanced potency in phenylsulfonyl fragments(Taylor and Francis Ltd., 2022) Das, B.K.; Chakraborty, D.The emergence of sulfa-drug resistance and reduced efficacy of pterin-based analogs towards Dihydropteroate synthase (DHPS) inhibition dictate a pressing need of developing novel antimicrobial agents for immune-compromised patients. Recently, a series of 8-Marcaptoguanin (8-MG) derivatives synthesized for 6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (experimental KD ∼ 100–.0.36) showed remarkable homology with the pteroic-acid and serve as a template for product antagonism in DHPS. The present work integrates ligand-based drug discovery techniques with structure-based docking, enhanced MD simulation, and MM/PBSA techniques to demonstrate the essential features of 8-MG analogs which make it a potent inhibitor for DHPS. The delicate balance in hydrophilic, hydrophobic substitutions on the 8-MG core is the crucial signature for DHPS inhibition. It is found that the dynamic interactions of active compounds are mainly dominated by consistent hydrogen bonding network with Asp 96, Asn 115, Asp 185, Ser 222, Arg 255 and π-π stacking, π-cation interactions with Phe 190, Lys 221. Further, two new 8-MG compounds containing N-phenylacetamide (compound S1, ΔGbind-eff = –62.03 kJ/mol) and phenylsulfonyl (compound S3, ΔGbind-eff = −71.29 kJ/mol) fragments were found to be the most potent inhibitor of DHPS, which stabilize the flexible pABA binding loop, thereby increasing their binding affinity. MM/PBSA calculation shows electrostatic energy contribution to be the principal component in stabilizing the inhibitors in the binding pocket. This fact is further confirmed by the higher energy barrier obtained in umbrella sampling for this class of inhibitors. © 2021 Informa UK Limited, trading as Taylor & Francis Group.Item Exploring the Barriers in the Aggregation of a Hexadecameric Human Prion Peptide through the Markov State Model(American Chemical Society, 2023) Das, B.K.; Singh, O.; Chakraborty, D.The prefibrillar aggregation kinetics of prion peptides are still an enigma. In this perspective, we employ atomistic molecular dynamics (MD) simulations of the shortest human prion peptide (HPP) (127GYMLGS132) at various temperatures and peptide concentrations and apply the Markov state model to determine the various intermediates and lag phases. Our results reveal that the natural mechanism of prion peptide self-assembly in the aqueous phase is impeded by two significant kinetic barriers with oligomer sizes of 6-9 and 12-13 peptides, respectively. The first one is the aggregation of unstructured lower-order oligomers, and the second is fibril nucleation, which impedes the further growth of prion aggregates. Among these two activation barriers, the second one is found to be dominant irrespective of the increase in temperature and peptide concentration. These lag phases are captured in all three different force-field parameters, namely, GROMOS-54a7, AMBER-99SB-ILDN, and CHARMMS 36m, at different concentrations. The GROMOS-54a7 and AMBER-99SB-ILDN force fields showed a comparatively higher percentage of β-sheet formation in the metastable aggregate that evolved during the aggregation process. In contrast, the CHARMM-36m force field showed mostly coil or turn conformations. The addition of a novel catecholamine derivative (naphthoquinone dopamine (NQDA)) arrests the aggregation process between the lag phases by increasing the activation barrier for the Lag1 and Lag2 phases in all of the force fields, which further validates the existence of these lag phases. The preferential binding of NQDA with the peptides increases the hydration of peptides and eventually disrupts the organized morphology of prefibrillar aggregates. It reduces the dimer dissociation energy by −24.34 kJ/mol. © 2023 American Chemical Society.Item Effect of peptide hydrophilicity on membrane curvature and permeation(American Institute of Physics, 2024) Mathath, A.V.; Chakraborty, D.Using a well-developed reaction coordinate in umbrella sampling, we studied the single peptide permeation through a model cancerous cell membrane, varying the hydrophilicity and the charge of the peptides. Two peptides, melittin and pHD108, were studied. The permeation mechanism differs from a barrel-stave-like mechanism to toroidal pore and vesicle formation based on the number and the placement of the hydrophilic amino acids in the peptide. Membrane curvature changes dynamically as the permeation process occurs. In the case of vesicles, the peptide traverses along a smooth, homogenous pathway, whereas a rugged, steep pathway was found when lipid molecules did not line up along the wall of the membrane (barrel-stave-like mechanism). A mechanism similar to a toroidal pore consists of multiple minima. Higher free energy was found for the permeating terminal containing charged amino acid residues. Vesicle formation was found for pHD108 peptide N-terminal with a maximum membrane thinning effect of 54.4% with free energy cost of 8.20 ± 0.10 kcal mol?1 and pore radius of 2.33 ± 0.07 nm. Insights gained from this study can help to build synthetic peptides for drug delivery. © 2024 Author(s).Item Evaluating the impact of the membrane thickness on the function of the intramembrane protease GlpG(Elsevier B.V., 2024) Engberg, O.; Mathath, A.V.; Döbel, V.; Frie, C.; Lemberg, M.K.; Chakraborty, D.; Huster, D.Cellular membranes exhibit a huge diversity of lipids and membrane proteins that differ in their properties and chemical structure. Cells organize these molecules into distinct membrane compartments characterized by specific lipid profiles and hydrophobic thicknesses of the respective domains. If a hydrophobic mismatch occurs between a membrane protein and the surrounding lipids, there can be functional consequences such as reduced protein activity. This phenomenon has been extensively studied for single-pass transmembrane proteins, rhodopsin, and small polypeptides such as gramicidin. Here, we investigate the E. coli rhomboid intramembrane protease GlpG as a model to systematically explore the impact of membrane thickness on GlpG activity. We used fully saturated 1,2-dilauroyl-sn-glycero-3-phosphocholine (DLPC) and 1,2-dimyristoyl-sn-glycero-3-phosphocholine(DMPC) model lipids and altered membrane thickness by varying the cholesterol content. Physical membrane parameters were determined by 2H and 31P NMR spectroscopy and correlated with GlpG activity measurements in the respective host membranes. Differences in bulk and annular lipids as well as alterations in protein structure in the respective host membranes were determined using molecular dynamics simulations. Our findings indicate that GlpG can influence the membrane thickness in DLPC/cholesterol membranes but not in DMPC/cholesterol membranes. Moreover, we observe that GlpG protease activity is reduced in DLPC membranes at low cholesterol content, which was not observed for DMPC. While a change in GlpG activity can already be due to smallest differences in the lipid environment, potentially enabling allosteric regulation of intramembrane proteolysis, there is no overall correlation to cholesterol-mediated lipid bilayer organization and phase behavior. Additional factors such as the influence of cholesterol on membrane bending rigidity and curvature energy need to be considered. In conclusion, the functionality of ?-helical membrane proteins such as GlpG relies not only on hydrophobic matching but also on other membrane properties, specific lipid interaction, and the composition of the annular layer. © 2024 The Author(s)Item A new reaction coordinate to study the translocation pathway of cell-penetrating peptides across lipid bilayers: The cases of transportan-10 and penetratin(American Institute of Physics, 2025) Mathath, A.V.; Chakraborty, D.Translocation pathway of cell-penetrating peptides remains elusive, as it is hard to observe by experimental and theoretical studies, which limits their effective use. Furthermore, lipid dynamics influence the translocation pathway, which is often overlooked due to its slow timescale. Current studies lack the effect of multiple peptides on the translocation process. Therefore, in this work, we employ the umbrella sampling technique with a preferential lipid–peptide interaction term in the reaction coordinate to explore the translocation activity of penetratin and transportan-10 (TP10) peptides in a heterogeneous membrane. In experiments, they follow different pathways according to their concentration, but the cause of this difference is unknown. We considered single and multiple (two and four) peptide translocation processes to understand the differences. Self-aggregation process is taken into account for multiple peptides. The interaction between peptides and peptides–lipids is found to be important for a proper overview of the translocation process. Peptide translocation was found to be related to the dynamics of the lipids, which change during the translocation process, making the system complex to study. In the case of multiple penetratin translocation, the anionic lipids were found to aggregate on the positive curvature of the upper leaflet, helping fold the membrane. Lipid composition of the TP10 multiple peptide case was found random. The increased mass and size of the solute in this case helped attain a radius more than the threshold value, leading to pore formation. Free energy barriers of single TP10 and penetratin are found to be 45.4 ± 2 and 33.7 ± 0.8 kJ mol?1, respectively. © 2025 Author(s).