Mohd. Shahbaaz

Large scale analysis of the mutational landscape in β-glucuronidase: A major player of mucopolysaccharidosis type VII

The lysosomal storage disorders are a group of 50 unique inherited diseases characterized by unseemly lipid storage in lysosomes. These malfunctions arise due to genetic mutations that result in deficiency or reduced activities of the lysosomal enzymes, which are responsible for catabolism of biological macromolecules. Sly syndrome or mucopolysaccharidosis type VII is a lysosomal storage disorder associated with the deficiency of β-glucuronidase (EC 3.2.1.31) that catalyzes the hydrolysis of β-D-glucuronic acid residues from the non-reducing terminal of glycosaminoglycan. The effects of the disease causing mutations on the framework of the sequences and structure of β-glucuronidase (GUSBp) were analyzed utilizing a variety of bioinformatic tools. These analyses showed that 211 mutations may result in alteration of the biological activity of GUSBp, including previously experimentally validated mutations. Finally, we refined 90 disease causing mutations, which presumably cause a significant impact on the structure, function, and stability of GUSBp. Stability analyses showed that mutations p.Phe208Pro, p.Phe539Gly, p.Leu622Gly, p.Ile499Gly and p.Ile586Gly caused the highest impact on GUSBp stability and function because of destabilization of the protein structure. Furthermore, structures of wild type and mutant GUSBp were subjected to molecular dynamics simulation to examine the relative structural behaviors in the explicit conditions of water. In a broader view, the use of in silico approaches provided a useful understanding of the effect of single point mutations on the structure-function relationship of GUSBp.

Functional Annotation of Conserved Hypothetical Proteins from Haemophilus influenzae Rd KW20

Haemophilus influenzae is a Gram negative bacterium that belongs to the family Pasteurellaceae, causes bacteremia, pneumonia and acute bacterial meningitis in infants. The emergence of multi-drug resistance H. influenzae strain in clinical isolates demands the development of better/new drugs against this pathogen. Our study combines a number of bioinformatics tools for function predictions of previously not assigned proteins in the genome of H. influenzae. This genome was extensively analyzed and found 1,657 functional proteins in which function of 429 proteins are unknown, termed as hypothetical proteins (HPs). Amino acid sequences of all 429 HPs were extensively annotated and we successfully assigned the function to 296 HPs with high confidence. We also characterized the function of 124 HPs precisely, but with less confidence. We believed that sequence of a protein can be used as a framework to explain known functional properties. Here we have combined the latest versions of protein family databases, protein motifs, intrinsic features from the amino acid sequence, pathway and genome context methods to assign a precise function to hypothetical proteins for which no experimental information is available. We found these HPs belong to various classes of proteins such as enzymes, transporters, carriers, receptors, signal transducers, binding proteins, virulence and other proteins. The outcome of this work will be helpful for a better understanding of the mechanism of pathogenesis and in finding novel therapeutic targets for H. influenzae.

Identification of Functional Candidates amongst Hypothetical Proteins of Treponema pallidum ssp. pallidum

Syphilis is a globally occurring venereal disease, and its infection is propagated through sexual contact. The causative agent of syphilis, Treponema pallidum ssp. pallidum, a Gram-negative sphirochaete, is an obligate human parasite. Genome of T. pallidum ssp. pallidum SS14 strain (RefSeq NC_010741.1) encodes 1,027 proteins, of which 444 proteins are known as hypothetical proteins (HPs), i.e., proteins of unknown functions. Here, we performed functional annotation of HPs of T. pallidum ssp. pallidum using various database, domain architecture predictors, protein function annotators and clustering tools. We have analyzed the sequences of 444 HPs of T. pallidum ssp. pallidum and subsequently predicted the function of 207 HPs with a high level of confidence. However, functions of 237 HPs are predicted with less accuracy. We found various enzymes, transporters, binding proteins in the annotated group of HPs that may be possible molecular targets, facilitating for the survival of pathogen. Our comprehensive analysis helps to understand the mechanism of pathogenesis to provide many novel potential therapeutic interventions.

PKR-inhibitor binds efficiently with human microtubule affinity-regulating kinase 4

MAP/microtubule affinity-regulating kinase 4 (MARK4) plays a central role in the cellular physiology, and it is inseparably linked with many human diseases including cancer, diet induced obesity, type2 diabetes and neurodegenerative disorders. Here, we studied the interaction of PKR-inhibitor with two variants of human MARK4. One variant is named as MARK4-F1 which has 59 N-terminal residues along with kinase domain while another variant is MARK4-F2 which has kinase domain only. Molecular-docking, molecular dynamics (MD) simulation and fluorescence-binding studies were undertaken to understand the role of N-terminal 59-residues in the binding of substrate/inhibitors. Molecular docking studies revealed that the PKR-inhibitor binds in the large hydrophobic cavity of the kinase domain of MARK4 through several hydrophobic and hydrogen-bonded interactions. Furthermore, MD simulation showed a stable parameters for the complexes of both MARK4-F1 and MARK4-F2 to PKR-inhibitor with marginal difference in their binding affinities. A significant decrease in the fluorescence intensity of MARK4 was observed on successive addition of the PKR-inhibitor. Using fluorescence data we have calculated the binding-affinity and the number of binding site of PKR-inhibitor to the MARK4. A significantly high binding affinity was observed for the PKR-inhibitor to the MARK4 variants. However, there is no any significant difference in the binding affinity of PKR-inhibitor to the MARK4 variants was observed, indicating that 59 N-terminal residues of MARK4-F1 are not playing a crucial role in the ligand binding. The present study will provide an insights into designing of new PKR-inhibitor derivative as potent and selective therapeutic agent against many life threatening diseases which are associated with MARK4.

Urea-induced denaturation of human calcium/calmodulin-dependent protein kinase IV: a combined spectroscopic and MD simulation studies

Calcium/calmodulin-dependent protein kinase IV (CaMKIV) is a multifunctional enzyme which belongs to the Ser/Thr kinase family. CaMKIV plays important role in varieties of biological processes such as gene expression regulation, memory consolidation, bone growth, T-cell maturation, sperm motility, regulation of microtubule dynamics, cell-cycle progression, and apoptosis. To measure stability parameters, urea-induced denaturation of CaMKIV was carried out at pH 7.4 and 25°C, using three different probes, namely far-UV CD, near-UV absorption, and tryptophan fluorescence. A coincidence of normalized denaturation curves of these optical properties suggests that urea-induced denaturation is a two-state process. Analysis of these denaturation curves gave values of 4.20 ± 0.12 kcal mol−1, 2.95 ± 0.15 M, and 1.42 ± 0.06 kcal mol−1 M−1 for (Gibbs free energy change (ΔGD) in the absence of urea), Cm (molar urea concentration ([urea]) at the midpoint of the denaturation curve), and m (=∂ΔGD/∂[urea]), respectively. All these experimental observations have been fully supported by 30 ns molecular dynamics simulation studies.

Effect of pH on the structure, function, and stability of human calcium/calmodulin-dependent protein kinase IV: combined spectroscopic and MD simulation studies

Human calcium/calmodulin-dependent protein kinase IV (CAMKIV) is a member of Ser/Thr protein kinase family. It is regulated by the calcium-calmodulin dependent signal through a secondary messenger, Ca(2+), which leads to the activation of its autoinhibited form. The over-expression and mutation in CAMKIV as well as change in Ca(2+) concentration is often associated with numerous neurodegenerative diseases and cancers. We have successfully cloned, expressed, and purified a functionally active kinase domain of human CAMKIV. To observe the effect of different pH conditions on the structural and functional properties of CAMKIV, we have used spectroscopic techniques such as circular diachroism (CD) absorbance and fluorescence. We have observed that within the pH range 5.0-11.5, CAMKIV maintained both its secondary and tertiary structures, along with its function, whereas significant aggregation was observed at acidic pH (2.0-4.5). We have also performed ATPase activity assays under different pH conditions and found a significant correlation between the structure and enzymatic activities of CAMKIV. In-silico validations were further carried out by modeling the 3-dimensional structure of CAMKIV and then subjecting it to molecular dynamics (MD) simulations to understand its conformational behavior in explicit water conditions. A strong correlation between spectroscopic observations and the output of molecular dynamics simulation was observed for CAMKIV.

Effect of pH on structure, function, and stability of mitochondrial carbonic anhydrase VA

Mitochondrial carbonic anhydrase VA (CAVA) catalyzes the hydration of carbon dioxide to produce proton and bicarbonate which is primarily expressed in the mitochondrial matrix of liver, and involved in numerous physiological processes including lipogenesis, insulin secretion from pancreatic cells, ureagenesis, gluconeogenesis, and neuronal transmission. To understand the effect of pH on the structure, function, and stability of CAVA, we employed spectroscopic techniques such as circular dichroism, fluorescence, and absorbance measurements in wide range of pH (from pH 2.0 to pH 11.5). CAVA showed an aggregation at acidic pH range from pH 2.0 to pH 5.0. However, it remains stable and maintains its secondary structure in the pH range, pH 7.0–pH 11.5. Furthermore, this enzyme has an appreciable activity at more than pH 7.0 (7.0 < pH ≤ 11.5) with maximum activity at pH 9.0. The maximal values of kcat and kcat/Km at pH 9.0 are 3.7 × 106 s−1 and 5.5 × 107 M−1 s−1, respectively. However, this enzyme loses its activity in the acidic pH range. We further performed 20-ns molecular dynamics simulation of CAVA to see the dynamics at different pH values. An excellent agreement was observed between in silico and in vitro studies. This study provides an insight into the activity of CAVA in the pH range of subcellular environment.

Structure-based functional annotation of putative conserved proteins having lyase activity from Haemophilus influenzae

Haemophilus influenzae is a small pleomorphic Gram-negative bacteria which causes several chronic diseases, including bacteremia, meningitis, cellulitis, epiglottitis, septic arthritis, pneumonia, and empyema. Here we extensively analyzed the sequenced genome of H. influenzae strain Rd KW20 using protein family databases, protein structure prediction, pathways and genome context methods to assign a precise function to proteins whose functions are unknown. These proteins are termed as hypothetical proteins (HPs), for which no experimental information is available. Function prediction of these proteins would surely be supportive to precisely understand the biochemical pathways and mechanism of pathogenesis of Haemophilus influenzae. During the extensive analysis of H. influenzae genome, we found the presence of eight HPs showing lyase activity. Subsequently, we modeled and analyzed three-dimensional structure of all these HPs to determine their functions more precisely. We found these HPs possess cystathionine-β-synthase, cyclase, carboxymuconolactone decarboxylase, pseudouridine synthase A and C, D-tagatose-1,6-bisphosphate aldolase and aminodeoxychorismate lyase-like features, indicating their corresponding functions in the H. influenzae. Lyases are actively involved in the regulation of biosynthesis of various hormones, metabolic pathways, signal transduction, and DNA repair. Lyases are also considered as a key player for various biological processes. These enzymes are critically essential for the survival and pathogenesis of H. influenzae and, therefore, these enzymes may be considered as a potential target for structure-based rational drug design. Our structure–function relationship analysis will be useful to search and design potential lead molecules based on the structure of these lyases, for drug design and discovery.

In silico approaches for the identification of virulence candidates amongst hypothetical proteins of Mycoplasma pneumoniae 309

Mycoplasma pneumoniae type 2a strain 309 is a simplest known bacterium and is the primary cause of community acquired pneumonia in the children. It mainly causes severe atypical pneumonia as well as several other non-pulmonary manifestations such as neurological, hepatic, hemolytic anemia, cardiac diseases and polyarthritis. The size of M. pneumoniae genome (Accession number: NC_016807.1) is relatively smaller as compared to other bacteria and contains 707 functional proteins, in which 204 are classified as hypothetical proteins (HPs) because of the unavailability of experimentally validated functions. The functions of the HPs were predicted by integrating a variety of protein classification systems, motif discovery tools as well as methods that are based on characteristic features obtained from the protein sequence and metabolic pathways. The probable functions of 83HPs were predicted successfully. The accuracy of the diverse tools used in the adopted pipeline was evaluated on the basis of statistical techniques of Receiver Operating Characteristic (ROC), which indicated the reliability of the functional predictions. Furthermore, the virulent HPs present in the set of 83 functionally annotated proteins were predicted by using the Bioinformatics tools and the conformational behaviours of the proteins with highest virulence scores were studied by using the molecular dynamics (MD) simulations. This study will facilitate in the better understanding of various drug resistance and pathogenesis mechanisms present in the M. pneumoniae and can be utilized in designing of better therapeutic agents.

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

Carbonic anhydrase IX (CAIX) is a transmembrane glycoprotein, overexpressed in cancer cells under hypoxia condition. In cancerous cells, CAIX plays an important role to combat the deleterious effects of a high rate of glycolytic metabolism. In order to favor tumor survival, CAIX maintains intracellular pH neutral or slightly alkaline and extracellular acidic pH. The equilibrium unfolding and conformational stability of CAIX were measured in the presence of increasing urea concentrations to understand its structural features under stressed conditions. Two different spectroscopic techniques were used to follow urea-induced denaturation and observed that urea induces a reversible denaturation of CAIX. Coincidence of the normalized transition curves of both optical properties suggesting that denaturation of CAIX is a two-state process, i.e., native state ↔ denatured state. Each denaturation curve was analyzed to estimate thermodynamic parameters, ΔGD0,value of Gibbs free energy change (ΔGD) associated with the urea-induced denaturation, Cm (midpoint of denaturation) and m (=δΔGD/δ[urea]). We further performed molecular dynamics simulation of CAIX for 50ns to see the dynamics of protein structure in the presence of different urea concentrations. An excellent agreement was observed between in silico and in vitro studies.