Faculty Publications
Permanent URI for this communityhttps://idr.nitk.ac.in/handle/123456789/18736
Publications by NITK Faculty
Browse
34 results
Search Results
Item Interstitial Voids and Resultant Density of Liquid Water: A First-Principles Molecular Dynamics Study(American Chemical Society service@acs.org, 2018) Biswas, S.; Chakraborty, D.; Mallik, B.S.Many anomalous properties of water can be explained on the basis of the coexistence of more than one density states: high-density water (HDW) and low-density water (LDW). We investigated these two phases of water molecules through first-principles molecular dynamics simulations using density functional theory (DFT) in conjunction with various van der Waals-corrected exchange and correlation functionals. Different density regions were found to exist due to the difference in short-range and long-range forces present in DFT potentials. These density regions were identified and analyzed on the basis of the distribution of molecules and voids present. We defined a local structure index to distinguish and find the probability of occurrence of these different states. HDW and LDW arise due to the presence of "interstitial water" molecules in between the first and second coordination shells. The population of interstitial water molecules is found to affect the overall dynamics of the system as they change the hydrogen bond pattern. © 2018 American Chemical Society.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 Hydrophilicity of the hydrophobic group: Effect of cosolvents and ions(Elsevier B.V., 2019) Dilip, H.N.; Chakraborty, D.Classical molecular dynamics simulations were performed to study the effect of cosolvents and ions on the solvation structure of zwitterionic glycine in liquid water. Simulations were carried out for 2 M and 1 M concentration of TMAO, Urea, KCl and LiCl solutions to observe the changes in liquid structure of water near the glycine molecule. Radial distribution functions and spatial distribution functions showed the presence of protective hydration layer near the C ? in presence of TMAO which gets reduced in case of urea, KCl and minimum in case of LiCl. LiCl is found to disrupt severely the solvation structure near the glycine molecule. For LiCl system, a small hydration layer is found near C ? unit at higher distances which is mainly due to the first hydration shell of lithium ion bonded to the carboxylate group. Presence of these hydration layers gives extra stabilization energy to the glycine water system. Stabilizing and destabilizing effect of water near the glycine molecule is calculated in terms of Potential Mean Force. The anomalous behaviour of lithium salts with respect to Group I cation salts in protein stabilization can be explained on the basis of this behaviour. We found maximum hydrogen bond lifetime for water molecules in presence of TMAO followed by LiCl, KCl and least in case of urea. The higher lifetimes in presence of ions are found mainly due to their electrostatic force. The stabilization of the hydrophobic part of the glycine molecule can be correlated with the stabilization of proteins in presence of these cosolvents. © 2019 Elsevier B.V.Item Carbohelicenes and thiahelicene from phthalaldehydes through Perkin approach(Elsevier B.V., 2019) Sarkar, P.; Das, B.K.; Chakraborty, D.; Muthamma, K.Synthesis and structural features of helical nanographene molecules comprising of seven benzene rings are examined. Thus dibutyl-dicarboxylate functional [7]helicene and its two regioisomers, dinaphtho[1,2–a:1?,2?–h]anthracene and naphtho[2,1–c]pentahelicene, have been synthesized in two steps through Perkin approach using napthalene-2-acetic acid and ortho- or meta-phthalaldehydes. The feasibility of this approach to construct sulfur doped twisted dithiaarenes is also investigated by using thiophene-3-acetic acid. While dithiaarenes from meta-phthalaldehyde remains challenging, synthesis and characterization of planar anthra[1,2–b:5,6–b']dithiophene and twisted 1,12-dithiapentahelicene is successful from ortho-phthalaldehyde. Conformational analysis with DFT calculation shows unique helicity preference in such doubly helical carbon nanostructures. Absorption and emission behavior of these ?-extended molecules shows enhanced conjugation. © 2019 Elsevier B.V.Item 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 Effect of cosolvents in the preferential binding affinity of water in aqueous solutions of amino acids and amides(Elsevier B.V., 2020) Dilip, H.N.; Chakraborty, D.Effects of two naturally occurring osmolytes, urea and trimethylamine-N-oxide (TMAO) on the solvation structure of hydrophobic moiety of alanine, glycine, N-methylacetamide and acetamide are investigated by classical molecular dynamics simulations. Our results are analysed in terms of site-site radial distribution functions (RDF), spatial distribution functions (SDF), number of hydrogen bonds, orientation profile, KB integrals, preferential binding coefficient and hydrogen bond dynamics. RDF and SDF showed presence of an extra hydration shell near the hydrophobic unit when TMAO is present in the solution. This hydration shell mainly consists of broken hydrogen bonds. In urea-water solution, intramolecular association is favoured compared to intermolecular association: which is in contrast to the TMAO-water solution. Alanine, glycine, NMA and acetamide showed preferred interactions with the water molecules in presence of TMAO compared to urea. Urea and TMAO both are found to be excluded from the alanine, glycine, NMA and acetamide surface but presence of urea was slightly favoured at higher distances in case of NMA and acetamide. The strong hydrogen bond between TMAO-water increases the hydrogen bond lifetime of other hydrogen bonds in the system. The preferential binding affinity of water with the protein molecules and strong hydrogen bonds are found to be the key reasons for stability in presence of TMAO. © 2019 Elsevier B.V.Item Structural and Thermophysical Anomalies of Liquid Water: A Tale of Molecules in the Instantaneous Low- And High-Density Regions(American Chemical Society service@acs.org, 2020) Priyadarsini, A.; Biswas, A.; Chakraborty, D.; Mallik, B.S.Water is believed to be a heterogeneous liquid comprising multiple density regions that arise because of the presence of interstitial molecules and can be differentiated by their structure as well as the existence of hydrogen-bonded pairs with varying strengths. First-principles molecular dynamics studies were performed at six different temperatures to investigate the effect of temperature on the thermophysical, structure, dynamics, and vibrational spectral properties of the water molecules using dispersion-corrected density functional theory. The variation of properties like density, cohesive energy, and compressibility with a change in temperature produces a trend that matches with the experiments and resembles the experimentally observed anomalous behavior. We explore the possibility of explaining the trends in calculated properties by analyzing the structure and dynamics of the water molecules in terms of instantaneous low- and instantaneous high-density regions that are found during the simulation time. The dynamics of these two types of water molecules were studied by calculating the lifetime from the proposed autocorrelation functions. The lifetime of formation of instantaneous low-density water is found to decrease with an increase in temperature, whereas the lifetime of instantaneous high-density water is found to be maximum at 298 K among all the considered temperatures. The presence of more interstitial water molecules is observed at this temperature. The signature of these water molecules is found in the radial distribution function, spatial distribution function, void distribution, configurational space, orientational dynamics, and spectral diffusion calculations. It is also found that around 298 K, these water molecules are present distinctively that mix up with the first and second solvation shells with the rise of the temperature. The outlook of the reported results can be extended to other thermodynamic conditions to explain some of the anomalous properties, which can be related to the presence of the interstitial molecules in water. © © 2020 American Chemical Society.Item In-silico epitope identification and design of Uricase mutein with reduced immunogenicity(Elsevier Ltd, 2020) Nelapati, A.K.; Das, B.K.; JagadeeshBabu, J.B.; Chakraborty, D.The clinical utilization of Uricase against gout is limited due to the immunogenicity. In the present article, we identified the antigenic determinants of Uricase and reduced their immunogenicity via in-silico mutagenesis. Multiple sequence alignment and motif analysis were carried out to identify the conserved residues in evolutionary process. Emini surface accessibility, Parker hydrophilicity, and Karplus & Schulz flexibility methods were employed to predict the linear B-cell epitopes of both Ag-Uricase and Bf-Uricase. Deimmunization approach identified T-cell epitopes and the hot spot residues. Reduced antigenic probability was obtained in case of T159W, D169C, N264W and Y203D mutations for Ag-Uricase, while S139 V, K215W, G216 F and I172 P mutations for Bf-Uricase. The binding affinity values of uric acid towards the catalytic pocket of Ag-Uricase and Bf-Uricase models were found to be -48.71 kcal/mol and -40.93 kcal/mol, respectively. This energy is further stabilized in the mutant model by -6.36 kcal/mol and -1.45 kcal/mol for Ag-Uricase and Bf-Uricase, respectively. About 100 ns molecular dynamics simulation was performed to evaluate the conformational stability of both native and mutated Uricase. Insights obtained from this study provide guidelines for experimental design of Uricase muteins with reduced antigenicity. © 2020 Elsevier LtdItem 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 Diverse interactions of aggregated insulin with selected coumarin dyes: Time dependent fluorogenicity, simulation studies and comparison with thioflavin T(Elsevier Ltd, 2021) Dalal, S.; Das, B.K.; Saini, M.; Chakraborty, D.; Sadhu, K.K.In this study, we have compared neutral coumarin based well-known commercially available probes C6, C7 and C545T for fluorogenic response from the aggregated insulin. The immediate fluorogenic responses were comparatively poor from all the three probes with respect to the previously reported response from thioflavin T (ThT) in the presence of aggregated insulin. Interestingly C6 among the three neutral coumarin derivative showed a significant steady increase of fluorescence intensity with time up to 6 h before reaching the saturation limit. Similar time dependent fluorogenic experiment with C7, C545T and ThT showed comparatively fast saturation within few minutes to 2 h. The molecular docking and simulation studies showed that these neutral probes could be stabilized in the aggregated form of the insulin predominantly by non-covalent weak interactions such as hydrogen bonding, ?-? and cation-? interactions. The probability distributions of the dihedral angles between two heterocyclic parts in C6 showed maximum probability of occurrence at 0° and 180°. These probability distributions of the dihedral angles between two heterocyclic parts within all the four fluorophores provided the justification of selective time dependent fluorescence enhancement from C6 in presence of insulin aggregate. The overall fluorogenic enhancement from C6 was comparable to the fluorogenic response from ThT and theoretical study confirmed distinctly different origin of this associated slow time dependent fluorogenic response. © 2020 Elsevier Ltd