Our group has previously reported the isolation, partial characterisation, and application of a Galbeta1-3GalNAc- and IgA1-reactive lectin from the seeds of champedak (Artocarpus integer). In the present study, we have subjected the purified lectin to reverse-phase high performance liquid chromatography and sequenced its subunits. Determination of the N-terminal sequence of the first 47 residues of the large subunit demonstrated at least 95% homology to the N-terminal sequence of the alpha chains of a few other galactose-binding Artocarpus lectins. The two smaller subunits of the lectin, each comprised of 21 amino acid residues, demonstrated minor sequence variability. Their sequences were generally comparable to the beta chains of the other galactose-binding Artocarpus lectins. When used to probe human serum glycopeptides that were separated by two-dimensional gel electrophoresis, the lectin demonstrated strong apparent interactions with glycopeptides of IgA1, hemopexin, alpha2-HS glycoprotein, alpha1-antichymotrypsin, and a few unknown glycoproteins. Immobilisation of the lectin to Sepharose generated an affinity column that may be used to isolate the O-glycosylated serum glycoproteins.
The arginine repressor (ArgR) of Escherichia coli binds to six L-arginine molecules that act as its co-repressor in order to bind to DNA. The binding of L-arginine molecules as well as its structural analogues is compared by means of computational docking. A grid-based energy evaluation method combined with a Monte Carlo simulated annealing process was used in the automated docking. For all ligands, the docking procedure proposed more than one binding site in the C-terminal domain of ArgR (ArgRc). Interaction patterns of ArgRc with L-arginine were also observed for L-canavanine and L-citrulline. L-lysine and L-homoarginine, on the other hand, were shown to bind poorly at the binding site. Figure A general overview of the sites found from docking the various ligands into ArgRc ( grey ribbons). Red coloured sticks: residues in binding site H that was selected for docking
Protein-ligand binding free energy values of wild-type and mutant C-terminal domain of Escherichia coli arginine repressor (ArgRc) protein systems bound to L-arginine or L-citrulline molecules were calculated using the linear interaction energy (LIE) method by molecular dynamics (MD) simulation. The binding behaviour predicted by the dissociation constant (K(d)) calculations from the binding free energy values showed preferences for binding of L-arginine to the wild-type ArgRc but not to the mutant ArgRc(D128N). On the other hand, L-citrulline do not favour binding to wild-type ArgRc but prefer binding to mutant ArgRc(D128N). The dissociation constant for the wild-type ArgRc-L-arginine complex obtained in this study is in agreement with reported experimental results. Our results also support the experimental data for the binding of L-citrulline to the mutant ArgRc(D128N). These showed that LIE method for protein-ligand binding free energy calculation could be applied to the wild-type and the mutant E. coli ArgRc-L-arginine and ArgRc-L-citrulline protein-ligand complexes and possibly to other transcriptional repressor-co-repressor systems as well.
Filamentous phage random peptide libraries were used to identify the epitopes of Burkholderia pseudomallei protease by panning against IgG polyclonal sera that exhibited protease neutralizing properties. The isolated fusion peptides presented a consensus peptide sequence, TKSMALSG, which closely resembles part of the active site sequence, 435GTSMATPHVAG445, of B. pseudomallei serine metalloprotease. By comparing the consensus sequence, TKSMALSG, with the predicted three-dimensional molecular model of B. pseudomallei serine metalloprotease, it appears that the potential antibody binding epitope was buried within the molecule. This active site was conformational whereby one continuous sub-region (SMA) was located between two discontinuous sub-regions, supplied by the flanking residues in the same polypeptide. All phages selected from the biopanning with IgG polyclonal sera showed good binding towards the polyclonal antibodies when compared to the negative control. In addition, these peptide-bearing phages showed competitive inhibition of B. pseudomallei serine metalloprotease binding to the polyclonal IgG.
A non-redundant database of 4536 structural domains, comprising more than 790,000 residues, has been used for the calculation of their solvent accessibility in the native protein environment and then in the isolated domain environment. Nearly 140,000 (18%) residues showed a change in accessible surface area in the above two conditions. General features of this change under these two circumstances have been pointed out. Propensities of these interfacing amino acid residues have been calculated and their variation for different secondary structure types has been analyzed. Actual amount of surface area lost by different secondary structures is higher in the case of helix and strands compared to coil and other conformations. Overall change in surface area in hydrophobic and uncharged residues is higher than that in charged residues. An attempt has been made to know the predictability of interface residues from sequence environments. This analysis and prediction results have significant implications towards determining interacting residues in proteins and for the prediction of protein-protein, protein-ligand, protein-DNA and similar interactions.
Bacterial proteases are an important group of enzymes that have very diverse biochemical and cellular functions. Proteases from prokaryotic sources also have a wide range of uses, either in medicine as pathogenic factors or in industry and therapeutics. ProLysED (Prokaryotic Lysis Enzymes Database), our meta-server integrated database of bacterial proteases, is a useful, albeit very niche, resource. The features include protease classification browsing and searching, organism-specific protease browsing, molecular information and visualisation of protease structures from the Protein Data Bank (PDB) as well as predicted protease structures.
In this study, the hypothesis was tested that behaviour of rats under the open field test condition and effects of subsequent acute stress relate to conformational properties of the main plasma carrier protein, albumin.To evaluate albumin properties, fluorescence intensity of a molecular probe CAPIDAN (N-carboxyphenylimide of dimethylaminonaphthalic acid) at N (at pH 7.4) and F (at pH 4.2) albumin conformations was measured and the N-F signal ratio was calculated. The data obtained showed that CAPIDAN fluoresces selectively from albumin in rat serum and its fluorescence is sensitive to binding of fatty acids and some other ligands to albumin. Behaviour of 78 Wistar male rats was characterized from the fraction of time taken for exploratory and ambulatory activity during the open field test. In rats not subjected to stress (n = 40), a negative correlation was revealed between open field activity and CAPIDAN N-to-F ratio for albumin (r = - 0.55, p < 0.0005). In the group of rats subjected to acute stress (immobilization plus stochastic electrocutaneous stimulation) the correlation between behavioural activity and the albumin conformational properties was significantly positive (r = 0.59, p < 0.0001): the CAPIDAN albumin fluorescence ratio increased in the highly active rats and decreased in the low-activity rats. The mechanisms of the observed effects may involve differences in nonesterified fatty acid production during stress.
Previous molecular dynamic simulations have reported elongation of the existing beta-sheet in prion proteins. Detailed examination has shown that these elongations do not extend beyond the proline residues flanking these beta-sheets. In addition, proline has also been suggested to possess a possible structural role in preserving protein interaction sites by preventing invasion of neighboring secondary structures. In this work, we have studied the possible structural role of the flanking proline residues by simulating mutant structures with alternate substitution of the proline residues with valine. Simulations showed a directional inhibition of elongation, with the elongation progressing in the direction of valine including evident inhibition of elongation by existing proline residues. This suggests that the flanking proline residues in prion proteins may have a containment role and would confine the beta-sheet within a specific length.
Curculin isolated from Curculigo latifolia, a plant grown in Malaysia, has an intriguing property of modifying sour taste into sweet taste. In addition to this taste-modifying activity, curculin itself elicits a sweet taste. Although these activities have been attributed to the heterodimeric isoform and not homodimers of curculin, the underlying mechanisms for the dual action of this protein have been largely unknown. To identify critical sites for these activities, we performed a mutational and structural study of recombinant curculin. Based on the comparison of crystal structures of curculin homo- and heterodimers, a series of mutants was designed and subjected to tasting assays. Mapping of amino acid residues on the three-dimensional structure according to their mutational effects revealed that the curculin heterodimer exhibits sweet-tasting and taste-modifying activities through its partially overlapping but distinct molecular surfaces. These findings suggest that the two activities of the curculin heterodimer are expressed through its two different modes of interactions with the T1R2-T1R3 heterodimeric sweet taste receptor.
The continuing rise in tuberculosis incidence and the problem of drug resistance strains have prompted the research on new drug candidates and the mechanism of drug resistance. Molecular docking and molecular dynamics simulation (MD) were performed to study the binding of isoniazid onto the active site of Mycobacterium tuberculosis enoyl-acyl carrier protein reductase (InhA) in an attempt to address the mycobacterial resistance against isoniazid. Results show that isonicotinic acyl-NADH (INADH) has an extremely high binding affinity toward the wild type InhA by forming stronger interactions compared to the parent drug (isoniazid) (INH). Due to the increase of hydrophobicity and reduction in the side chain's volume of A94 of mutant type InhA, both INADH and the mutated protein become more mobile. Due to this reason, the molecular interactions of INADH with mutant type are weaker than that observed with the wild type. However, the reduced interaction caused by the fluctuation of INADH and the mutant protein only inflected minor resistance in the mutant strain as inferred from free energy calculation. MD results also showed there exists a water-mediated hydrogen bond between INADH and InhA. However, the bridged water molecule is only present in the INADH-wild type complex, reflecting the putative role of the water molecule in the binding of INADH to the wild type protein. The results support the assumption that the conversion of prodrug isoniazid into its active form INADH is mediated by KatG as a necessary step prior to target binding on InhA. Our findings also contribute to a better understanding of INH resistance in mutant type.
Molecular Dynamics (MD) simulations have been used to understand how protein structure, dynamics, and flexibility are affected by adaptation to high temperature for several years. We report here the results of the high temperature MD simulations of Bacillus stearothermophilus L1 (L1 lipase). We found that the N-terminal moiety of the enzyme showed a high flexibility and dynamics during high temperature simulations which preceded and followed by clear structural changes in two specific regions; the small domain and the main catalytic domain or core domain of the enzyme. These two domains interact with each other through a Zn(2+)-binding coordination with Asp-61 and Asp-238 from the core domain and His-81 and His-87 from the small domain. Interestingly, the His-81 and His-87 were among the highly fluctuated and mobile residues at high temperatures. The results appear to suggest that tight interactions of Zn(2+)-binding coordination with specified residues became weak at high temperature which suggests the contribution of this region to the thermostability of the enzyme.
Using 100-fold molar excess of succinic anhydride, about 99% of lysine residues of hen egg white lysozyme (HEWL) were modified. Succinylated (S(99)) HEWL showed both charge and size homogeneity as judged by PAGE and gel filtration, respectively. Hydrodynamic parameters such as Stokes radius and frictional ratio (f/f(o)) showed more expanded conformation of S(99) HEWL compared to native HEWL as evident from the increase in Stokes radius (from 1.36 to 1.86 nm) and f/f(o) (from 0.86 to 1.15) values. Guanidine hydrochloride (GdnHCl) denaturation studies using fluorescence spectroscopy connoted a marked decrease in conformational stability of HEWL upon succinylation. Complete denaturation of S(99) HEWL was achieved at lower GdnHCl concentration ( approximately 3.8 M) compared to native HEWL ( approximately 5 M). Furthermore, free energy of stabilization (DeltaG(D)(H(2)O)) value also showed a notable decrease from 8,559 and 7,956 cal/mol (for native HEWL) to 4,404 and 4,669 cal/mol (for succinylated HEWL) using excitation at 280 and 295 nm, respectively. Both expanded conformation and decreased DeltaG(D)(H(2)O) can be attributed to the increase in the net negative charge on the protein upon succinylation. All these results manifested the importance of positively charged lysine residues in maintaining the conformational stability of HEWL through electrostatic interactions.
Matched MeSH terms: Protein Conformation/drug effects
Twenty percent of genes that encode for hypothetical proteins from Klebsiella pneumoniae MGH78578 were identified, leading to KPN00728 and KPN00729 after bioinformatics analysis. Both open reading frames showed high sequence homology to Succinate dehydrogenase Chain C (SdhC) and D (SdhD) from Escherichia coli. Recently, KPN00729 was assigned as SdhD. KPN00728 thus remains of particular interest as no annotated genes from the complete genome sequence encode for SdhC. We discovered KPN00728 has a missing region with conserved residues important for ubiquinone (UQ) and heme group binding. Structure and function prediction of KPN00728 coupled with secondary structure analysis and transmembrane topology showed KPN00728 adopts SDH-(subunit C)-like structure. To further probe its functionality, UQ was docked on the built model (consisting KPN00728 and KPN00729) and formation of hydrogen bonds between UQ and Ser27, Arg31 (KPN00728) and Tyr84 (KPN00729) further reinforces and supports that KPN00728 is indeed SDH. This is the first report on the structural and function prediction of KPN00728 of K. pneumoniae MGH78578 as SdhC.
T1 lipase is isolated from the palm Geobacillus zalihae strain T1 in Malaysia, functioning as a secreted protein responsible for the catalyzing hydrolysis of long-chain triglycerides into fatty acids and glycerol at high temperatures. In the current study, using 30 ns molecular dynamics simulations at different temperatures, an aqueous activation was detected for T1 lipase. This aqueous activation in T1 lipase was mainly caused by a double-flap movement mechanism. The double flaps were constituted by the hydrophobic helices 6 and 9. Helix 6 employed two major components with the hydrophilic part at the surface and the hydrophobic part inside. In the aqueous solution, the hydrophobic part could provide enough power for helix 6 to move away, driving the protein into an open configuration and exposing the catalytic triad. Our findings could provide structural evidence to support the double-flap movement, revealing the catalytic mechanism for T1 lipase.
An all-atom level MD simulation in explicit solvent at high temperature is a powerful technique to increase our knowledge about the structurally important regions modulating thermal stability in thermenzymes. In this respect, two large-sized thermoalkalophilic enzymes from Bacillus stearothermophilus L1 (L1 lipase) and Geobacillus zalihae strain T1 (T1 lipase) are well-established representatives. In this paper, comparative results from temperature-induced MD simulations of both model systems at 300 K, 400 K and 500 K are presented and discussed with respect to identification of highly flexible regions critical to thermostability. From our MD simulation results, specific regions along the L1 lipase and T1 lipase polypeptide chain including the small domain and the main catalytic domain or core domain of both enzymes show a marked increase in fluctuations and dynamics followed by clear structural changes. Overall, the N-terminal moiety of both enzymes and their small domains exhibit hyper-sensitivity to thermal stress. The results appear to propose that these regions are critical in determining of the overall thermal stability of both organisms.
Newcastle disease virus (NDV) is an avian paramyxovirus that has gained a lot of interest in cancer viro-therapeutic applications because of its ability to selectively induce apoptosis in human cancer cells. However, the underlying mechanisms by which NDV induces apoptosis in human cancer cells are still not entirely understood.
Japanese encephalitis virus (JEV) is the most important cause of epidemic encephalitis worldwide but its origin is unknown. Epidemics of encephalitis suggestive of Japanese encephalitis (JE) were described in Japan from the 1870s onwards. Four genotypes of JEV have been characterised and representatives of each genotype have been fully sequenced. Based on limited information, a single isolate from Malaysia is thought to represent a putative fifth genotype. We have determined the complete nucleotide and amino acid sequence of Muar strain and compared it with other fully sequenced JEV genomes. Muar was the least similar, with nucleotide divergence ranging from 20.2 to 21.2% and amino acid divergence ranging from 8.5 to 9.9%. Phylogenetic analysis of Muar strain revealed that it does represent a distinct fifth genotype of JEV. We elucidated Muar signature amino acids in the envelope (E) protein, including E327 Glu on the exposed lateral surface of the putative receptor binding domain which distinguishes Muar strain from the other four genotypes. Evolutionary analysis of full-length JEV genomes revealed that the mean evolutionary rate is 4.35 × 10(-4) (3.4906 × 10(-4) to 5.303 × 10(-4)) nucleotides substitutions per site per year and suggests JEV originated from its ancestral virus in the mid 1500s in the Indonesia-Malaysia region and evolved there into different genotypes, which then spread across Asia. No strong evidence for positive selection was found between JEV strains of the five genotypes and the E gene has generally been subjected to strong purifying selection.