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.
Tuberculosis (TB), an epidemic disease, affects the world with death rate of two million people every year. The bacterium Mycobacterium tuberculosis was found to be a more potent and disease-prolonged bacterium among the world due to multi-drug resistance. Emergence of new drug targets is needed to overcome the bacterial resistance that leads to control epidemic tuberculosis. The pathway thiamine biosynthesis was targeting M. tuberculosis due to its role in intracellular growth of the bacterium. The screening of enzymes involved in thiamin biosynthesis showed novel target thiazole synthase (ThiG) involved in catalysis of rearrangement of 1-deoxy-D-xylulose 5-phosphate (DXP) to produce the thiazole phosphate moiety of thiamine. We carried out homology modeling for ThiG to understand the structure-function relationship, and the model was refined with MD simulations. The results showed that the model predicted with (α + β)8-fold of synthase family proteins. Molecular docking of ThiG model with substrate DXP showed binding mode and key residues ARG46, ASN69, THR41, and LYS96 involved in the catalysis. First-line anti-tuberculosis drugs were docked with ThiG to identify the inhibition. The report showed the anti-tuberculosis drugs interact well with ThiG which may lead to block thiamin biosynthesis pathway.
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.
A new subfamily of glycosyl hydrolase family GH13 was recently proposed for α-amylases from Anoxybacillus species (ASKA and ADTA), Geobacillus thermoleovorans (GTA, Pizzo, and GtamyII), Bacillus aquimaris (BaqA), and 95 other putative protein homologues. To understand this new GH13 subfamily, we report crystal structures of truncated ASKA (TASKA). ASKA is a thermostable enzyme capable of producing high levels of maltose. Unlike GTA, biochemical analysis showed that Ca(2+) ion supplementation enhances the catalytic activities of ASKA and TASKA. The crystal structures reveal the presence of four Ca(2+) ion binding sites, with three of these binding sites are highly conserved among Anoxybacillus α-amylases. This work provides structural insights into this new GH13 subfamily both in the apo form and in complex with maltose. Furthermore, structural comparison of TASKA and GTA provides an overview of the conformational changes accompanying maltose binding at each subsite.
The low presentation efficiency of Fab (fragment antigen binding) fragments during phage display is largely due to the complexity of disulphide bond formation. This can result in the presentation of Fab fragments devoid of a light chain during phage display. Here we propose the use of a coplasmid system encoding several molecular chaperones (DsbA, DsbC, FkpA, and SurA) to improve Fab packaging. A comparison was done using the Fab fragment from IgG and IgD. We found that the use of the coplasmid during phage packaging was able to improve the presentation efficiency of the Fab fragment on phage surfaces. A modified version of panning using the coplasmid system was evaluated and was successful at enriching Fab binders. Therefore, the coplasmid system would be an attractive alternative for improved Fab presentation for phage display.
Polyhydroxyalkanoates (PHAs) are biopolymers synthesized by a wide range of bacteria, which serve as a promising candidate in replacing some conventional petrochemical-based plastics. PHA synthase (PhaC) is the key enzyme in the polymerization of PHA, and the crystal structures were successfully determined using the catalytic domain of PhaC from Cupriavidus necator (PhaCCn-CAT) and Chromobacterium sp. USM2 (PhaCCs-CAT). Here, we review the beneficial mutations discovered in PhaCs from a structural perspective. The structural comparison of the residues involved in beneficial mutation reveals that the residues are near to the catalytic triad, but not inside the catalytic pocket. For instance, Ala510 of PhaCCn is near catalytic His508 and may be involved in the open-close regulation, which presumably play an important role in substrate specificity and activity. In the class II PhaC1 from Pseudomonas sp. 61-3 (PhaC1Ps), Ser325 stabilizes the catalytic cysteine through hydrogen bonding. Another residue, Gln508 of PhaC1Ps is located in a conserved hydrophobic pocket which is next to the catalytic Asp and His. A class I, II-conserved Phe420 of PhaCCn is one of the residues involved in dimerization and its mutation to serine greatly reduced the lag phase. The current structural analysis shows that the Phe362 and Phe518 of PhaC from Aeromonas caviae (PhaCAc) are assisting the dimer formation and maintaining the integrity of the core beta-sheet, respectively. The structure-function relationship of PhaCs discussed in this review will serve as valuable reference for future protein engineering works to enhance the performance of PhaCs and to produce novel biopolymers.
Matched MeSH terms: Protein Conformation, beta-Strand
Clostridial neurotoxins, including tetanus and botulinum neurotoxins, generally target vertebrates. We show here that this family of toxins has a much broader host spectrum, by identifying PMP1, a clostridial-like neurotoxin that selectively targets anopheline mosquitoes. Isolation of PMP1 from Paraclostridium bifermentans strains collected in anopheline endemic areas on two continents indicates it is widely distributed. The toxin likely evolved from an ancestral form that targets the nervous system of similar organisms, using a common mechanism that disrupts SNARE-mediated exocytosis. It cleaves the mosquito syntaxin and employs a unique receptor recognition strategy. Our research has an important impact on the study of the evolution of clostridial neurotoxins and provides the basis for the use of P. bifermentans strains and PMP1 as innovative, environmentally friendly approaches to reduce malaria through anopheline control.
Mycobacterium tuberculosis (Mtb) 16.3 kDa heat shock protein 16.3 (HSP16.3) is a latency-associated antigen that can be targeted for latent tuberculosis (TB) diagnostic and therapeutic development. We have previously developed human VH domain antibodies (dAbs; clone E3 and F1) specific against HSP16.3. In this work, we applied computational methods to optimise and design the antibodies in order to improve the binding affinity with HSP16.3. The VH domain antibodies were first docked to the dimer form of HSP16.3 and further sampled using molecular dynamics simulation. The calculated binding free energy of the HSP16.3-dAb complexes showed non-polar interactions were responsible for the antigen-antibody association. Per-residue free energy decomposition and computational alanine scanning have identified one hotspot residue for E3 (Y391) and 4 hotspot residues for F1 (M394, Y396, R397 and M398). These hotspot residues were then mutated and evaluated by binding free energy calculations. Phage ELISA assay was carried out on the potential mutants (E3Y391W, F1M394E, F1R397N and F1M398Y). The experimental assay showed improved binding affinities of E3Y391W and F1M394E against HSP16.3 compared with the wild type E3 and F1. This case study has thus showed in silico methods are able to assist in optimisation or improvement of antibody-antigen binding.
A number of living creatures in the Antarctic region have developed characteristic adaptation of cold weather by producing antifreeze proteins (AFP). Antifreeze peptide (Afp1m) fragment have been designed in the sequence of strings from native proteins. The objectives of this study were to assess the properties of Afp1m to cryopreserve skin graft at the temperature of -10 °C and -20 °C and to assess sub-zero injuries in Afp1m cryopreserved skin graft using light microscopic techniques. In the present study, a process was developed to cryopreserve Sprague-Dawley (SD) rat skin grafts with antifreeze peptide, Afp1m, α-helix peptide fragment derived from Glaciozyma antractica yeast. Its viability assessed by different microscopic techniques. This study also described the damages caused by subzero temperatures (-10 and -20 °C) on tissue cryopreserved in different concentrations of Afp1m (0.5, 1, 2, 5 and 10 mg/mL) for 72 h. Histological scores of epidermis, dermis and hypodermis of cryopreserved skin grafts showed highly significant difference (p
Matched MeSH terms: Protein Conformation, alpha-Helical
The potentiality of bioactive phenolic compounds may result in plant extracts having multiple biological activities. The aim of this study was to investigate into the biological activities of the methanolic, ethyl acetate, and water extracts of Tchihatchewia isatidea Boiss, an endemic medicinal plant of Turkey. The phenolic compositions of the extracts were confirmed using RP-HPLC. Extracts were screened for their potential antioxidant through a panoply of assays; their anti-diabetic potential, and plausible inhibitory activity against tyrosinase and acetylcholinesterase. Molecular modelling methods were also used to assess the docking properties of phenolic compounds on tyrosinase. The major and most abundant compounds were rosmarinic acid (570 ± 14 μg/g extract in the methanolic extract), ferrulic acid (336 ± 6 μg/g extract in the methanolic extract), (+)-catechin (340 ± 4 μg/g extract in the water extract), apigenin (182 ± 4 μg/g extract in the methanolic extract), and epicatechin (188 ± 12 μg/g extract in the water extract). Radical scavenging, reducing capacity, and metal chelating activities were detected in the extracts, with preponderance activity observed in the methanolic extract. In conclusion, the potential clinical applications observed during this study may provide new insights into the molecular aspect particularly for neuroprotective and anti-diabetic mechanisms involving oxidative stress.
Pulmonary tuberculosis, caused by Mycobacterium tuberculosis, is one of the most persistent diseases leading to death in humans. As one of the key targets during the latent/dormant stage of M. tuberculosis, isocitrate lyase (ICL) has been a subject of interest for new tuberculosis therapeutics. In this work, the cleavage of the isocitrate by M. tuberculosis ICL was studied using quantum mechanics/molecular mechanics method at M06-2X/6-31+G(d,p): AMBER level of theory. The electronic embedding approach was applied to provide a better depiction of electrostatic interactions between MM and QM regions. Two possible pathways (pathway I that involves Asp108 and pathway II that involves Glu182) that could lead to the metabolism of isocitrate was studied in this study. The results suggested that the core residues involved in isocitrate catalytic cleavage mechanism are Asp108, Cys191 and Arg228. A water molecule bonded to Mg2+ acts as the catalytic base for the deprotonation of isocitrate C(2)-OH group, while Cys191 acts as the catalytic acid. Our observation suggests that the shuttle proton from isocitrate hydroxyl group C(2) atom is favourably transferred to Asp108 instead of Glu182 with a lower activation energy of 6.2 kcal/mol. Natural bond analysis also demonstrated that pathway I involving the transfer of proton to Asp108 has a higher intermolecular interaction and charge transfer that were associated with higher stabilization energy. The QM/MM transition state stepwise catalytic mechanism of ICL agrees with the in vitro enzymatic assay whereby Asp108Ala and Cys191Ser ICL mutants lost their isocitrate cleavage activities.
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.
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.
DNA polymerase β is a 39 kDa enzyme that is a major component of Base Excision Repair in human cells. The enzyme comprises two major domains, a 31 kDa domain responsible for the polymerase activity and an 8 kDa domain, which bind ssDNA and has a deoxyribose phosphate (dRP) lyase activity. DNA polymerase β was shown to be phosphorylated in vitro with protein kinase C (PKC) at serines 44 and 55 (S44 and S55), resulting in loss of its polymerase enzymic activity, but not its ability to bind ssDNA. In this study, we investigate the potential phosphorylation-induced structural changes for DNA polymerase β using molecular dynamics simulations. The simulations show drastic conformational changes of the polymerase structure as a result of S44 phosphorylation. Phosphorylation-induced conformational changes transform the closed (active) enzyme structure into an open one. Further analysis of the results points to a key hydrogen bond and newly formed salt bridges as potential drivers of these structural fluctuations. The changes observed with S55/44 and S55 phosphorylation were less dramatic and the integrity of the H-bond was not compromised. Thus the phosphorylation of S44 is the major contributor to structural fluctuations that lead to loss of enzymatic activity.
Effective novel peptide inhibitors which targeted the domain III of the dengue envelope (E) protein by blocking dengue virus (DENV) entry into target cells, were identified. The binding affinities of these peptides towards E-protein were evaluated by using a combination of docking and explicit solvent molecular dynamics (MD) simulation methods. The interactions of these complexes were further investigated by using the Molecular Mechanics-Poisson Boltzmann Surface Area (MMPBSA) and Molecular Mechanics Generalized Born Surface Area (MMGBSA) methods. Free energy calculations of the peptides interacting with the E-protein demonstrated that van der Waals (vdW) and electrostatic interactions were the main driving forces stabilizing the complexes. Interestingly, calculated binding free energies showed good agreement with the experimental dissociation constant (Kd) values. Our results also demonstrated that specific residues might play a crucial role in the effective binding interactions. Thus, this study has demonstrated that a combination of docking and molecular dynamics simulations can accelerate the identification process of peptides as potential inhibitors of dengue virus entry into host cells.
Matched MeSH terms: Protein Conformation, beta-Strand
Schistosomes are parasitic blood flukes that infect approximately 250 million people worldwide. The disease known as schistosomiasis, is the second most significant tropical parasitic disease after malaria. Praziquantel is the only effective drug currently licensed for schistosomiasis and there are concerns about resistance to the drug. There has been much effort to develop vaccines against schistosomiasis to produce long-term protection in endemic regions. Surface-associated proteins, and in particular, those expressed in the body wall, or tegument, have been proposed as potential vaccine targets. Of these, annexins are thought to be of integral importance for the stability of this apical membrane system. Here, we present the structural and immunobiochemical characterization of four homologous annexins namely annexin B30, annexin B5a, annexin B7a and annexin B5b from S. mansoni. Bioinformatics analysis showed that there was no signal peptide predicted for any annexin in this study. Further analysis showed that each of all four annexin protein possesses a primary structure consisting of a short but variable N-terminal region and a long C-terminal core containing four homologous annexin repeats (I-IV), which contain five alpha-helices. The life cycle expression profile of each annexin was assessed using quantitative PCR. The results showed that the overall transcript levels of the each of four homologous annexins were relatively low in the egg stage, but increased gradually after the transition of cercariae (the invasive schistosome larvae) to schistosomula (the post-invasive larvae). Circular dichroism (CD) demonstrated that rAnnexin B30, rAnnexin B5a and rAnnexin 7a were folded, showing a secondary structure content rich in alpha-helices. The membrane binding affinity was enhanced when rAnnexin B30, rAnnexin B5a and rAnnexin 7a was incubated in the presence of Ca2+. All annexin members evaluated in this study were immunolocalized to the tegument, with immunoreactivity also occurring in cells and in muscle of adult parasites. All four recombinant annexins were immunoreactive and they were recognized by the sera of mice infected with S. mansoni. In conclusion, the overall results present the molecular characterization of annexin B30, annexin B5a, annexin B7a and annexin B5b from S. mansoni in host-parasite interactions and strongly suggest that the molecules could be useful candidates for vaccine or diagnostic development.
Matched MeSH terms: Protein Conformation, alpha-Helical
Dehalogenases are of high interest due to their potential applications in bioremediation and in synthesis of various industrial products. DehL is an L-2-haloacid dehalogenase (EC 3.8.1.2) that catalyses the cleavage of halide ion from L-2-halocarboxylic acid to produce D-2-hydroxycarboxylic acid. Although DehL utilises the same substrates as the other L-2-haloacid dehalogenases, its deduced amino acid sequence is substantially different (<25%) from those of the rest L-2-haloacid dehalogenases. To date, the 3D structure of DehL is not available. This limits the detailed understanding of the enzyme's reaction mechanism. The present work predicted the first homology-based model of DehL and defined its active site. The monomeric unit of the DehL constitutes α/β structure that is organised into two distinct structural domains: main and subdomains. Despite the sequence disparity between the DehL and other L-2-haloacid dehalogenases, its structural model share similar fold as the experimentally solved L-DEX and DehlB structures. The findings of the present work will play a crucial role in elucidating the molecular details of the DehL functional mechanism.
Matched MeSH terms: Protein Conformation, alpha-Helical
Acid denaturation of champedak galactose-binding (CGB) lectin was studied in the pH range, 7.0-1.0 using intrinsic fluorescence and ANS fluorescence measurements. The lectin remained stable up to pH 5.0 and showed local disordering in the vicinity of the protein fluorophores within the pH range, 5.0-3.5. Decrease in the pH from pH 3.5 to pH 2.5 led to structural transition, marked by the decrease in the intrinsic fluorescence and increase in the ANS fluorescence signals. This can be ascribed to the dissociation of the tetrameric lectin into monomeric forms. Further decrease in the pH up to pH 1.5 produced another transition, which specified the unfolding of monomers as reflected from the decrease in both intrinsic fluorescence and ANS fluorescence signals. Characterization of the conformational states obtained at pH 7.0, pH 2.5 and pH 1.5 based on intrinsic and ANS fluorescence spectra, gel chromatographic behavior and thermal denaturation confirmed the existence of folded monomeric forms at pH 2.5 and unfolded states at pH 1.5. However, the aciddenatured state of CGB lectin at pH 1.5 retained significant residual structure, as evident from the greater loss of both secondary and tertiary structures in the presence of 6 M guanidine hydrochloride at low pH values. Anion-induced refolding below pH 1.5 was also seen using ANS fluorescence measurements.