Estimation of thermodynamic stability of human carbonic anhydrase IX from urea-induced denaturation and MD simulation studies

Integrated spectroscopic and MD simulation approach to decipher the effect of pH on the structure function of Staphylococcus aureus thymidine kinase

Staphylococcus aureus is a major human pathogen responsible for a variety of clinical infections, becoming increasingly resistant to antibiotics. To address this challenge, there is a need to identify new cellular targets and innovative approaches to expand treatment options. One such target is thymidine kinase (TK), a crucial enzyme in the pyrimidine salvage pathway, which plays a key role in the phosphorylation of thymidine, an essential component in DNA synthesis and repair. In this study, we have successfully cloned, expressed, and purified the TK protein. A comprehensive investigation into how different pH levels affect the structure and functional activity of TK, using a combination of spectroscopy, classical molecular dynamics simulations, and enzyme activity assays was conducted. Our study revealed that variation in pH disrupts secondary and tertiary structures of TK with noticeable aggregate formation at pH 5.0. Enzyme activity studies demonstrated that TK exhibited its maximum kinase activity within the physiological pH range. These findings strongly suggest a connection between structural changes and enzymatic activity, which was further supported by the agreement between the spectroscopic features we measured and the results of our MD simulations. Our study provides a deeper insight into the structural features of TK, which could potentially be harnessed for the development of therapeutic strategies aimed at combatting infectious diseases. Conformational dynamics plays an essential role in the design and development of effective inhibitors. Considering the effects of pH on the conformational dynamics of TK, our findings may be implicated in the development of potent and selective inhibitors.Communicated by Ramaswamy H. Sarma.

Structure and Dynamics of the Isozymes II and IX of Human Carbonic Anhydrase

Human carbonic anhydrases (HCAs) are responsible for the pH control and sensing in our body and constitute key components in the central pH paradigm connected to cancer therapeutics. However, little or no molecular level studies are available on the pH-dependent stability and functional dynamics of the known isozymes of HCA. The main objective of this Article is to report the first bench-marking study on the structure and dynamics of the two most efficient isozymes, HCA II and IX, at neutral pH using classical molecular dynamics (MD) and constant pH MD (CpHMD) simulations combined with umbrella sampling, transition path sampling, and Markov state models. Starting from the known crystal structures of HCA II and the monomeric catalytic domain of HCA IX (labeled as HCA IX-c), we have generated classical MD and CpHMD trajectories (of length 1 μs each). In all cases, the overall stability, RMSD, and secondary structure segments of the two isozymes are found to be quite similar. Functionally important dynamics of these two enzymes have been probed in terms of active site hydration, coordination of the Zn(II) ion to a transient excess water, and the formation of putative proton transfer paths. The most important difference between the two isozymes is observed for the side-chain fluctuations of His-64 that is expected to shuttle an excess proton out of the active site as a part of the rate-determining intramolecular proton transfer reaction. The relative stability of the stable inward and outward conformations of the His-64 side-chain and the underlying free energy surfaces are found to depend strongly on the isozyme. In each case, a lower free energy barrier is detected between predominantly inward conformations from predominantly outward ones when simulated under constant pH conditions. The kinetic rate constants of interconversion between different free energy basins are found to span 10⁷-10⁸ s^-1 with faster conformational transitions predicted at constant pH condition. The estimated rate constants and free energies are expected to validate if the fluctuation of the His-64 side-chain in HCA IX may have a significance similar to that known in the multistep catalytic cycle of HCA II.

Discovery of Natural Compounds as Potential Inhibitors of Human Carbonic Anhydrase II: An Integrated Virtual Screening, Docking, and Molecular Dynamics Simulation Study

Carbonic anhydrase II (CAII) is one of the zinc metalloenzymes that catalyze the reversible hydration of carbon dioxide, leading to the formation of bicarbonate and proton. CAII plays a significant role in health and disease. For example, CAII helps to maintain eye pressure while regulating the pH of the tumor microenvironment, and by extension, contributing to cancer progression. Owing to its remarkable role in cancer, visual health, and other human diseases, CAII can serve as an attractive therapeutic target. We report an original study based on high-throughput virtual screening of natural compounds from the ZINC database in search of potential inhibitors of CAII. We selected the hits based on the physicochemical, absorption, distribution, metabolism, excretion, and toxicity (ADMET) properties, pan-assay interference compound (PAINS) patterns, and interaction analysis. Importantly, two natural compounds were identified, ZINC08918123 and ZINC00952700, bearing considerable affinity and specific interactions to the residues of the CAII-binding pocket with well-organized conformational fitting compatibility. We investigated the conformational dynamics of CAII in complex with the identified compounds through molecular dynamics simulation, which revealed the formation of a stable complex preserved throughout the 100 ns trajectories. The stability of the protein/ligand complexes is maintained by significant numbers of noncovalent interactions throughout the simulations. In conclusion, natural compounds identified in the present study specifically and computer-assisted drug design broadly offer a reliable resource and strategy to discover potential promising therapeutic inhibitors of CAII to cure various cancers and glaucoma after further experimental validation and clinical studies.

Mechanistic insights into the urea-induced denaturation of human sphingosine kinase 1

Sphingosine kinase 1 (SphK1) plays a significant role in various cellular processes, including cell proliferation, apoptosis, and angiogenesis. SphK1 is considered as an attractive target for drug development owing to its connection with several diseases, including cancer. In the current work, the urea-induced unfolding of SphK1 was performed at pH 8.0 and 25 °C using CD and fluorescence spectroscopy. SphK1 follows a biphasic unfolding transition (N ⇌ I ⇌ D) with an intermediate (I) state populated around 4.0 M urea concentration. The circular dichroism ([θ]₂₂₂) and fluorescence emission spectra (λ_max) of SphK1 with increasing concentrations of urea were analyzed to calculate Gibbs free energy (ΔG⁰) for both the transitions (N ⇌ I and I ⇌ D). A significant overlap of both the transitions obtained by two spectroscopic properties ([θ]₂₂₂ and λ_max) was observed, indicating that both N ⇌ I and I ⇌ D transition follow two-step equilibrium unfolding pattern. Also, we performed 100 ns molecular dynamics (MD) simulations to get atomistic insights into the structural changes in SphK1 with increasing urea concentrations. Our results showed a consistent pattern of the SphK1 unfolding with increasing urea concentrations. Together, spectroscopic and MD simulation findings provide deep insights into the unfolding mechanism and conformational features of SphK1.

Unravelling the unfolding pathway of human Fas-activated serine/threonine kinase induced by urea

Fas-activated serine/threonine kinase (FASTK) is a mitochondria-associated nuclear protein that inhibits Fas- and UV-induced apoptosis. This protein is generally activated during Fas-mediated apoptosis by phosphorylating a nuclear RNA-binding protein T-cell intracellular antigen-1 and thus considered as a modulator of apoptosis. In the present study, we have examined the equilibrium unfolding and conformational stability of the kinase domain of FASTK (FASTK_353-444). The kinase domain of FASTK_353-444 was cloned, expressed, and purified. The folding ↔ unfolding transitions of urea-induced denaturation was monitored with the help of circular dichroism, intrinsic fluorescence, and UV absorption spectroscopies. Analysis of transition curves obtained from different probes revealed a coincidence of denaturation curves, suggesting that folding/unfolding of FASTK follows a two-state process with the midpoint (C_m) value at 3.50 ± 0.1 M. Urea-induced denaturation curves were further analyzed to estimate change in the Gibbs free energy in the absence of urea (ΔG_D⁰) associated with the equilibrium of denaturation. To get atomistic insights into the urea-induced denaturation of FASTK, we performed an all-atom molecular dynamics simulation for 100 ns. A close agreement was noticed between experimental and computational studies. This study will help to understand the unfolding mechanism and structural stability of the kinase domain of FASTK.Communicated by Ramaswamy H. Sarma.

Understanding Thermostability Factors of Barley Limit Dextrinase by Molecular Dynamics Simulations

Limit dextrinase (LD) is the only endogenous starch-debranching enzyme in barley (Hordeum vulgare, Hv), which is the key factor affecting the production of a high degree of fermentation. Free LD will lose its activity in the mashing process at high temperature in beer production. However, there remains a lack of understanding on the factor affecting the themostability of HvLD at the atomic level. In this work, the molecular dynamics simulations were carried out for HvLD to explore the key factors affecting the thermal stability of LD. The higher value of root mean square deviation (RMSD), radius of gyration (R_g), and surface accessibility (SASA) suggests the instability of HvLD at high temperatures. Intra-protein hydrogen bonds and hydrogen bonds between protein and water decrease at high temperature. Long-lived hydrogen bonds, salt bridges, and hydrophobic contacts are lost at high temperature. The salt bridge interaction analysis suggests that these salt bridges are important for the thermostability of HvLD, including E568–R875, D317–R378, D803–R884, D457–R214, D468–R395, D456–R452, D399–R471, and D541–R542. Root mean square fluctuation (RMSF) analysis identified the thermal-sensitive regions of HvLD, which will facilitate enzyme engineering of HvLD for enhanced themostability.

Effect of pH on the structure and function of pyruvate dehydrogenase kinase 3: Combined spectroscopic and MD simulation studies

Pyruvate dehydrogenase kinase-3 (PDK3) plays important role in the glucose metabolism and is associated with cancer progression, and thus being considered as an attractive target for cancer therapy. In this study, we employed spectroscopic techniques to study the structural and conformational changes in the PDK3 at varying pH conditions ranging from pH 2.0 to 12.0. UV/Vis, fluorescence and circular dichroism spectroscopic measurements revealed that PDK3 maintains its native-like structure (both secondary and tertiary) in the alkaline conditions (pH 7.0-12.0). However, a significant loss in the structure was observed under acidic conditions (pH 2.0-6.0). The propensity of aggregate formation at pH 4.0 was estimated by thioflavin T fluorescence measurements. To further complement structural data, kinase activity assay was performed, and maximum activity of PDK3 was observed at pH 7.0-8.0 range; whereas, its activity was lost under acidic pH. To further see conformational changes at atomistic level we have performed all-atom molecular dynamics at different pH conditions for 150 ns. A well defined correlation was observed between experimental and computational studies. This work highlights the significance of structural dependence of pH for wide implications in protein-protein interaction, biological function and drug design procedures.

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α.

Unravelling the unfolding mechanism of human integrin linked kinase by GdmCl-induced denaturation

Integrin-linked kinase (ILK) is a ubiquitously expressed Ser/Thr kinase which plays significant role in the cell-matrix interactions and growth factor signalling. In this study, guanidinium chloride (GdmCl)-induced unfolding of kinase domain of ILK (ILK_193-446) was carried out at pH 7.5 and 25 °C. Eventually, denaturation curves of mean residue ellipticity at 222 nm ([θ]₂₂₂) and fluorescence emission spectrum were analysed to estimate stability parameters. The optical properties maximum emission (λ_max) and difference absorption coefficient at 292 nm (Δε₂₉₂) were analysed. The denaturation curve was measured only in the GdmCl molar concentration ranging 3.0-4.2 M because protein was aggregating below 3.0 M of GdmCl concentrations. The denaturation process of ILK_193-446 was found as reversible at [GdmCl] ≥ 3.0 M. Moreover, a coincidence of normalized denaturation curves of optical properties ([θ]₂₂₂, Δε₂₉₂ and λ_max) suggesting that GdmCl-induced denaturation of ILK_193-446 is a two-state process. In addition, 100 ns molecular dynamics simulations were performed to see the effects of GdmCl on the structure and stability of ILK_193-446. Both the spectroscopic and molecular dynamics approaches provided clear insights into the stability and conformational properties of ILK_193-446.

Molecular dynamics perspective on the thermal stability of mandelate racemase

Mandelate racemase from Pseudomonas putida is a promising candidate for the dynamic kinetic resolution of α-hydroxy carboxylic acids. In the present study, the thermal stability of mandelate racemase was investigated through molecular dynamics simulations in the temperature range of 303-363 K, which can guide the design of mandelate racemase with higher stability. The basic features such as radius of gyration, surface accessibility, and secondary structure content suggested the instability of mandelate racemase at high temperatures. With increase in temperature, α-helix content reduced significantly, especially the α-helices exposed to the environment. At the simulation time scale considered, intra-protein hydrogen bonds, hydrogen bonds between protein and water decreased at 363 K, while the number of salt-bridges increased. The long-distance networks remarkably changed at 363 K. A considerable number of long-lived (percentage existence time higher than 90%) hydrogen bonds and Cα contacts were lost. Root mean square fluctuation analysis revealed regions with high fluctuation, which should be helpful in the reengineering of mandelate racemase for enhanced thermal stability.