Dengue virus (DENV) presents a significant threat to global public health with more than 500,000 hospitalizations and 25,000 deaths annually. Currently, there is no clinically approved antiviral drug to treat DENV infection. The envelope (E) glycoprotein of DENV is a promising target for drug discovery as the E protein is important for viral attachment and fusion. Understanding the structure and function of DENV E protein has led to the exploration of structure-based drug discovery of antiviral compounds and peptides against DENV infections. This review summarizes the structural information of the DENV E protein with regards to DENV attachment and fusion. The information enables the development of antiviral agents through structure-based approaches. In addition, this review compares the potency of antivirals targeting the E protein with the antivirals targeting DENV multifunctional enzymes, repurposed drugs and clinically approved antiviral drugs. None of the current DENV antiviral candidates possess potency similar to the approved antiviral drugs which indicates that more efforts and resources must be invested before an effective DENV drug materializes.
Certain proteins have demonstrated proficient human immunodeficiency virus (HIV-1) life cycle disturbance. Recently, the ankyrin repeat protein targeting the HIV-1 capsid, AnkGAG1D4, showed a negative effect on the viral assembly of the HIV-1NL4-3 laboratory strain. To extend its potential for future clinical application, the activity of AnkGAG1D4 in the inhibition of other HIV-1 circulating strains was evaluated. Chimeric NL4-3 viruses carrying patient-derived Gag/PR-coding regions were generated from 131 antiretroviral drug-naïve HIV-1 infected individuals in northern Thailand during 2001⁻2012. SupT1, a stable T-cell line expressing AnkGAG1D4 and ankyrin non-binding control (AnkA32D3), were challenged with these chimeric viruses. The p24CA sequences were analysed and classified using the K-means clustering method. Among all the classes of virus classified using the p24CA sequences, SupT1/AnkGAG1D4 demonstrated significantly lower levels of p24CA than SupT1/AnkA32D3, which was found to correlate with the syncytia formation. This result suggests that AnkGAG1D4 can significantly interfere with the chimeric viruses derived from patients with different sequences of the p24CA domain. It supports the possibility of ankyrin-based therapy as a broad alternative therapeutic molecule for HIV-1 gene therapy in the future.
beta-D-N-Acetylhexosaminidase, a family 20 glycosyl hydrolase, catalyzes the removal of β-1,4-linked N-acetylhexosamine residues from oligosaccharides and their conjugates. We constructed phylogenetic tree of β-hexosaminidases to analyze the evolutionary history and predicted functions of plant hexosaminidases. Phylogenetic analysis reveals the complex history of evolution of plant β-hexosaminidase that can be described by gene duplication events. The 3D structure of tomato β-hexosaminidase (β-Hex-Sl) was predicted by homology modeling using 1now as a template. Structural conformity studies of the best fit model showed that more than 98% of the residues lie inside the favoured and allowed regions where only 0.9% lie in the unfavourable region. Predicted 3D structure contains 531 amino acids residues with glycosyl hydrolase20b domain-I and glycosyl hydrolase20 superfamily domain-II including the (β/α)8 barrel in the central part. The α and β contents of the modeled structure were found to be 33.3% and 12.2%, respectively. Eleven amino acids were found to be involved in ligand-binding site; Asp(330) and Glu(331) could play important roles in enzyme-catalyzed reactions. The predicted model provides a structural framework that can act as a guide to develop a hypothesis for β-Hex-Sl mutagenesis experiments for exploring the functions of this class of enzymes in plant kingdom.
Effect of 1,1,1,3,3,3-hexafluoroisopropanol (HFIP) on acid-denatured Bacillus licheniformis α -amylase (BLA) at pH 2.0 was investigated by far-UV CD, intrinsic fluorescence, and ANS fluorescence measurements. Addition of increasing HFIP concentrations led to an increase in the mean residue ellipticity at 222 nm (MRE 222 nm) up to 1.5 M HFIP concentration beyond which it sloped off. A small increase in the intrinsic fluorescence and a marked increase in the ANS fluorescence were also observed up to 0.4 M HFIP concentration, both of which decreased thereafter. Far- and near-UV CD spectra of the HFIP-induced state observed at 0.4 M HFIP showed significant retention of the secondary structures closer to native BLA but a disordered tertiary structure. Increase in the ANS fluorescence intensity was also observed with the HFIP-induced state, suggesting exposure of the hydrophobic clusters to the solvent. Furthermore, thermal denaturation of HFIP-induced state showed a non-cooperative transition. Taken together, all these results suggested that HFIP-induced state of BLA represented a molten globule-like state at pH 2.0.
Polygonum minus is an aromatic plant, which contains high abundance of terpenoids, especially the sesquiterpenes C15H24. Sesquiterpenes were believed to contribute to the many useful biological properties in plants. This study aimed to functionally characterize a full length sesquiterpene synthase gene from P. minus. P. minus sesquiterpene synthase (PmSTS) has a complete open reading frame (ORF) of 1689 base pairs encoding a 562 amino acid protein. Similar to other sesquiterpene synthases, PmSTS has two large domains: the N-terminal domain and the C-terminal metal-binding domain. It also consists of three conserved motifs: the DDXXD, NSE/DTE, and RXR. A three-dimensional protein model for PmSTS built clearly distinguished the two main domains, where conserved motifs were highlighted. We also constructed a phylogenetic tree, which showed that PmSTS belongs to the angiosperm sesquiterpene synthase subfamily Tps-a. To examine the function of PmSTS, we expressed this gene in Arabidopsis thaliana. Two transgenic lines, designated as OE3 and OE7, were further characterized, both molecularly and functionally. The transgenic plants demonstrated smaller basal rosette leaves, shorter and fewer flowering stems, and fewer seeds compared to wild type plants. Gas chromatography-mass spectrometry analysis of the transgenic plants showed that PmSTS was responsible for the production of β -sesquiphellandrene.
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.
The limited sequence similarity of protein sequences with known structures has led to an indispensable need for computational technology to predict their structures. Structural bioinformatics (SB) has become integral in elucidating the sequence-structure-function relationship of a protein. This report focuses on the applications of SB within the context of protein engineering including its limitation and future challenges.
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.
The stability of enzymes is critical for their application in industrial processes, which generally require different conditions from the natural enzyme environment. Both rational and random protein engineering approaches have been used to increase stability, with the latter requiring extensive experimental effort for the screening of variants. Moreover, some general rules addressing the molecular origin of protein thermostability have been established. Herein, we demonstrate the use of molecular dynamics simulations to gain molecular level understanding of protein thermostability and to engineer stabilizing mutations. Carbonic anhydrase (CA) is an enzyme with a high potential for biotechnological carbon capture applications, provided it can be engineered to withstand the high temperature process environments, inevitable in most gas treatment units. In this study, we used molecular dynamics simulations at 343, 353, and 363 K to study the relationship between structure flexibility and thermostability in bacterial α-CAs and applied this knowledge to the design of mutants with increased stability. The most thermostable α-CA known, TaCA from Thermovibrio ammonificans, had the most rigid structure during molecular dynamics simulations, but also showed regions with high flexibility. The most flexible amino acids in these regions were identified from root mean square fluctuation (RMSF) studies, and stabilizing point mutations were predicted based on their capacity to improve the calculated free energy of unfolding. Disulfide bonds were also designed at sites with suitable geometries and selected based on their location at flexible sites, assessed by B-factor calculation. Molecular dynamics simulations allowed the identification of five mutants with lower RMSF of the overall structure at 400 K, compared to wild-type TaCA. Comparison of free-energy landscapes between wild-type TaCA and the most promising mutants, Pro165Cys-Gln170Cys and Asn140Gly, showed an increased conformational stability of the mutants at 400 K.
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.
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.
ST50, an outer-membrane component of the multi-drug efflux system from Salmonella enterica serovar Typhi, is an obligatory diagnostic antigen for typhoid fever. ST50 is an excellent and unique diagnostic antigen with 95% specificity and 90% sensitivity and is used in the commercial diagnosis test kit (TYPHIDOT(TM)). The crystal structure of ST50 at a resolution of 2.98 Å reveals a trimer that forms an α-helical tunnel and a β-barrel transmembrane channel traversing the periplasmic space and outer membrane. Structural investigations suggest significant conformational variations in the extracellular loop regions, especially extracellular loop 2. This is the location of the most plausible antibody-binding domain that could be used to target the design of new antigenic epitopes for the development of better diagnostics or drugs for the treatment of typhoid fever. A molecule of the detergent n-octyl-β-D-glucoside is observed in the D-cage, which comprises three sets of Asp361 and Asp371 residues at the periplasmic entrance. These structural insights suggest a possible substrate transport mechanism in which the substrate first binds at the periplasmic entrance of ST50 and subsequently, via iris-like structural movements to open the periplasmic end, penetrates the periplasmic domain for efflux pumping of molecules, including poisonous metabolites or xenobiotics, for excretion outside the pathogen.
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.
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.
Antibodies have different chemical properties capable of targeting a diverse nature of antigens. Traditionally, immune antibody libraries are perceived to be disease-specific with a skewed repertoire. The complexity during the generation of a combinatorial antibody library allows for a skewed but diverse repertoire to be generated. Strongyloides stercoralis is a parasite that causes strongyloidiasis, a potentially life-threatening disease with a complex diagnosis that impedes effective control and treatment of the disease. This study describes the isolation of monoclonal antibodies against S. stercoralis NIE recombinant protein using an immune antibody phage display library derived from lymphatic filaria-infected individuals. The isolated antibody clones showed both lambda and kappa light chains gene usage, with diverse amino acid distributions. Structural analysis showed that electropositivity and the interface area could determine the binding affinity of the clones with NIE. The successful identification of S. stercoralis antibodies from the filarial immune library highlights the breadth of antibody gene diversification in an immune antibody library that can be applied for closely related infections.
Polyhydroxyalkanoate (PHA) is a promising candidate for use as an alternative bioplastic to replace petroleum-based plastics. Our understanding of PHA synthase PhaC is poor due to the paucity of available three-dimensional structural information. Here we present a high-resolution crystal structure of the catalytic domain of PhaC from Chromobacterium sp. USM2, PhaC Cs -CAT. The structure shows that PhaC Cs -CAT forms an α/β hydrolase fold comprising α/β core and CAP subdomains. The active site containing Cys291, Asp447 and His477 is located at the bottom of the cavity, which is filled with water molecules and is covered by the partly disordered CAP subdomain. We designated our structure as the closed form, which is distinct from the recently reported catalytic domain from Cupriavidus necator (PhaC Cn -CAT). Structural comparison showed PhaC Cn -CAT adopting a partially open form maintaining a narrow substrate access channel to the active site, but no product egress. PhaC Cs -CAT forms a face-to-face dimer mediated by the CAP subdomains. This arrangement of the dimer is also distinct from that of the PhaC Cn -CAT dimer. These findings suggest that the CAP subdomain should undergo a conformational change during catalytic activity that involves rearrangement of the dimer to facilitate substrate entry and product formation and egress from the active site.
MYB proteins are highly conserved DNA-binding domains (DBD) and mutations in MYB oncoproteins have been reported to cause aberrant and augmented cancer progression. Identification of MYB molecular biomarkers predictive of cancer progression can be used for improving cancer management. To address this, a biomarker discovery pipeline was employed in investigating deleterious non-synonymous single nucleotide polymorphisms (nsSNPs) in predicting damaging and potential alterations on the properties of proteins. The nsSNP of the MYB family; MYB, MYBL1, and MYBL2 was extracted from the NCBI database. Five in silico tools (PROVEAN, SIFT, PolyPhen-2, SNPs&GO and PhD-SNP) were utilized to investigate the outcomes of nsSNPs. A total of 45 nsSNPs were predicted as high-risk and damaging, and were subjected to PMut and I-Mutant 2.0 for protein stability analysis. This resulted in 32 nsSNPs with decreased stability with a DDG score lower than - 0.5, indicating damaging effect. G111S, N183S, G122S, and S178C located within the helix-turn-helix (HTH) domain were predicted to be conserved, further posttranslational modifications and 3-D protein analysis indicated these nsSNPs to shift DNA-binding specificity of the protein thus altering the protein function. Findings from this study would help in the field of pharmacogenomic and cancer therapy towards better intervention and management of cancer.
AlyQ from Persicobacter sp. CCB-QB2 is an alginate lyase with three domains - a carbohydrate-binding domain modestly resembling family 16 carbohydrate-binding module (CBM16), a family 32 CBM (CBM32) domain, and an alginate lyase domain belonging to polysaccharide lyase family 7 (PL7). Although AlyQ can also act on polyguluronate (poly-G) and polymannuronate (poly-M), it is most active on alginate. Studies with truncated AlyQ showed that the CBM32 domain did not contribute to enhancing AlyQ's activity under the assayed conditions. Nevertheless, it could bind to cleaved but not intact alginate, indicating that the CBM32 domain recognises alginate termini. The crystal structure containing both CBM32 and catalytic domains show that they do not interact with one another. The CBM32 domain contains a conserved Arg that may bind to the carboxyl group of alginate. The catalytic domain, meanwhile, shares a conserved substrate-binding groove, and the presence of two negatively charged Asp residues may dictate substrate specificity especially at subsite +1. As Persicobacter sp. CCB-QB2 was unable to utilise alginate, AlyQ may function to help the bacterium degrade cell walls more efficiently.
Owing to their ability to maintain a thermodynamically stable fold at extremely high temperatures, thermophilic proteins (TTPs) play a critical role in basic research and a variety of applications in the food industry. As a result, the development of computation models for rapidly and accurately identifying novel TTPs from a large number of uncharacterized protein sequences is desirable. In spite of existing computational models that have already been developed for characterizing thermophilic proteins, their performance and interpretability remain unsatisfactory. We present a novel sequence-based thermophilic protein predictor, termed SCMTPP, for improving model predictability and interpretability. First, an up-to-date and high-quality dataset consisting of 1853 TPPs and 3233 non-TPPs was compiled from published literature. Second, the SCMTPP predictor was created by combining the scoring card method (SCM) with estimated propensity scores of g-gap dipeptides. Benchmarking experiments revealed that SCMTPP had a cross-validation accuracy of 0.883, which was comparable to that of a support vector machine-based predictor (0.906-0.910) and 2-17% higher than that of commonly used machine learning models. Furthermore, SCMTPP outperformed the state-of-the-art approach (ThermoPred) on the independent test dataset, with accuracy and MCC of 0.865 and 0.731, respectively. Finally, the SCMTPP-derived propensity scores were used to elucidate the critical physicochemical properties for protein thermostability enhancement. In terms of interpretability and generalizability, comparative results showed that SCMTPP was effective for identifying and characterizing TPPs. We had implemented the proposed predictor as a user-friendly online web server at http://pmlabstack.pythonanywhere.com/SCMTPP in order to allow easy access to the model. SCMTPP is expected to be a powerful tool for facilitating community-wide efforts to identify TPPs on a large scale and guiding experimental characterization of TPPs.
Crystal structures of holo vitamin D receptor (VDR) revealed a canonical conformation in which the ligand is entrapped in a hydrophobic cavity buried in the ligand-binding domain (LBD). The mousetrap model postulates that helix 12 is positioned away from the domain to expose the interior cavity. However, the extended form of helix 12 is likely due to artifacts during crystallization. In this study, we set out to investigate conformational dynamics of apo VDR using molecular dynamics simulation on microsecond timescale. Here we show the neighboring backbones of helix 2-helix 3n and beta strand 2-helix 6 of LBD, instead of the helix 12, undergo large-scale motion, possibly gating the entrance of ligand to the ligand binding domain. Docking analysis to the simulated open structure of VDR with the estimated free energy of -37.0 kJ/mol, would emphasise the role of H2-H3n and S2-H6 in facilitating the entrance of calcitriol to the LBD of VDR.
Matched MeSH terms: Protein Conformation; Protein Conformation, alpha-Helical; Protein Conformation, beta-Strand