Protein stability: a crystallographer's perspective

Inhibition of Chlamydial Infection by CRISPR/Cas9-SAM Mediated Enhancement of Human Peptidoglycan Recognition Proteins Gene Expression in HeLa Cells

The global problem of emerging resistance of microorganisms to antibiotics makes the search for new natural substances with antibacterial properties relevant. Such substances include peptidoglycan recognition proteins (PGLYRP), which are the components of the innate immunity of many organisms, including humans. These proteins have a unique mechanism of action that allows them to evade the resistance of bacteria to them, as well as to be active against both Gram-positive and Gram-negative bacteria. However, the use of antimicrobial recombinant proteins is not always advisable due to the complexity of local delivery of the proteins and their stability; in this regard it seems appropriate to activate the components of the innate immunity. The aim of this study was to increase the expression level of native peptidoglycan recognition protein genes in HeLa cells using genome-editing technology with synergistic activation mediators (CRISPR/Cas9-SAM) and evaluate antichlamydial effect of PGLYRP. We demonstrated activation of the chlamydial two-component gene system (ctcB-ctcC), which played a key role in the mechanism of action of the peptidoglycan recognition proteins. We generated the HeLa cell line transduced with lentiviruses encoding CRISPR/Cas9-SAM activation system with increased PGLYRP gene expression. It was shown that activation of the own peptidoglycan recognition proteins gene expression in the cell line caused inhibition of the chlamydial infection development. The proposed approach makes it possible to use the capabilities of innate immunity to combat infectious diseases caused by Gram-positive and Gram-negative bacteria.

In silico analysis of deleterious SNPs of human MTUS1 gene and their impacts on subsequent protein structure and function

The mitochondrial tumor suppressor 1 (MTUS1) gene acts as a crucial tumor suppressor by inhibiting growth and proliferation of eukaryotic cells including tumor cell lines. Down regulation of MTUS1 gene has been implicated in a wide range of cancers as well as various human diseases. Alteration through nsSNPs can potentially damage the structure and/or function of the protein. As characterization of functional SNPs in such disease linked genes is a major challenge, it is feasible to analyze putative functional SNPs prior to performing larger population studies. Hence, in this in silico study we differentiated the potentially harmful nsSNPs of the MTUS1 gene from the neutral ones by using various sequence and structure based bioinformatic tools. In a total of 215 nsSNPs, 9 were found to be most likely to exert deleterious effect using 7 prediction tools. From which, 5nsSNPs (S1259L, E960K, P503T, L1084V and L1143Q) were selected as potentially damaging due to their presence in the highly conserved region and ability to decrease protein stability. In fact, 2 nsSNPs (S1259L and E960K) among these 5 were found to be individually associated with two distinctive cancers named Stomach adenocarcinoma and Uterine corpus endometrial carcinoma. As this is the first comprehensive study analyzing the functional nsSNPs of MTUS1, the results of the current study would certainly be helpful in future prospects concerning large population-based studies as well as drug discovery, especially developing individualized medicine.

Heterologous Expression and Biochemical Characterization of the Human Zinc Transporter 1 (ZnT1) and Its Soluble C-Terminal Domain

Human zinc transporter 1 (hZnT1) belongs to the cation diffusion facilitator (CDF) family. It plays a major role in transporting zinc (Zn2+) from the cytoplasm across the plasma membrane and into the extracellular space thereby protecting cells from Zn2+ toxicity. Through homology with other CDF family members, ZnT1 is predicted to contain a transmembrane region and a soluble C-terminal domain though little is known about its biochemistry. Here, we demonstrate that human ZnT1 and a variant can be produced by heterologous expression in Saccharomyces cerevisiae cells and purified in the presence of detergent and cholesteryl hemisuccinate. We show that the purified hZnT1 variant has Zn2+/H+ antiporter activity. Furthermore, we expressed, purified and characterized the soluble C-terminal domain of hZnT1 (hZnT1-CTD) in a bacterial expression system. We found that the hZnT1-CTD melting temperature increases at acidic pH, thus, we used an acetate buffer at pH 4.5 for purifications and concentration of the protein up to 12 mg/mL. Small-angle X-ray scattering analysis of hZnT1-CTD is consistent with the formation of a dimer in solution with a V-shaped core.

Engineering of Ancestors as a Tool to Elucidate Structure, Mechanism, and Specificity of Extant Terpene Cyclase

Structural information is crucial for understanding catalytic mechanisms and to guide enzyme engineering efforts of biocatalysts, such as terpene cyclases. However, low sequence similarity can impede homology modeling, and inherent protein instability presents challenges for structural studies. We hypothesized that X-ray crystallography of engineered thermostable ancestral enzymes can enable access to reliable homology models of extant biocatalysts. We have applied this concept in concert with molecular modeling and enzymatic assays to understand the structure activity relationship of spiroviolene synthase, a class I terpene cyclase, aiming to engineer its specificity. Engineering a surface patch in the reconstructed ancestor afforded a template structure for generation of a high-confidence homology model of the extant enzyme. On the basis of structural considerations, we designed and crystallized ancestral variants with single residue exchanges that exhibited tailored substrate specificity and preserved thermostability. We show how the two single amino acid alterations identified in the ancestral scaffold can be transferred to the extant enzyme, conferring a specificity switch that impacts the extant enzyme's specificity for formation of the diterpene spiroviolene over formation of sesquiterpenes hedycaryol and farnesol by up to 25-fold. This study emphasizes the value of ancestral sequence reconstruction combined with enzyme engineering as a versatile tool in chemical biology.

Rational Design of Adenylate Kinase Thermostability through Coevolution and Sequence Divergence Analysis

Protein engineering is actively pursued in industrial and laboratory settings for high thermostability. Among the many protein engineering methods, rational design by bioinformatics provides theoretical guidance without time-consuming experimental screenings. However, most rational design methods either rely on protein tertiary structure information or have limited accuracies. We proposed a primary-sequence-based algorithm for increasing the heat resistance of a protein while maintaining its functions. Using adenylate kinase (ADK) family as a model system, this method identified a series of amino acid sites closely related to thermostability. Single- and double-point mutants constructed based on this method increase the thermal denaturation temperature of the mesophilic Escherichia coli (E. coli) ADK by 5.5 and 8.3 °C, respectively, while preserving most of the catalytic function at ambient temperatures. Additionally, the constructed mutants have improved enzymatic activity at higher temperature.

Membrane Protein Stabilization Strategies for Structural and Functional Studies

Accounting for nearly two-thirds of known druggable targets, membrane proteins are highly relevant for cell physiology and pharmacology. In this regard, the structural determination of pharmacologically relevant targets would facilitate the intelligent design of new drugs. The structural biology of membrane proteins is a field experiencing significant growth as a result of the development of new strategies for structure determination. However, membrane protein preparation for structural studies continues to be a limiting step in many cases due to the inherent instability of these molecules in non-native membrane environments. This review describes the approaches that have been developed to improve membrane protein stability. Membrane protein mutagenesis, detergent selection, lipid membrane mimics, antibodies, and ligands are described in this review as approaches to facilitate the production of purified and stable membrane proteins of interest for structural and functional studies.

Implication of post-translationally modified SOD1 in pathological aging

Why certain people relish healthy aging throughout their life span while others suffer pathological consequences? In this review, we focus on some of the dominant paradigms of pathological aging, such as amyotrophic lateral sclerosis (ALS), Alzheimer's disease (AD), and Parkinson's disease (PD), and predict that the antioxidant superoxide dismutase 1 (SOD1), when post-translationally modified by aging-associated oxidative stress, acts as a mechanism to accelerated aging in these age-related neurodegenerative diseases. Oxidative modifications of natively reduced SOD1 induce pathological confirmations such as misfolding, leading to a subsequent formation of monomeric, oligomeric, and multimeric aggregates. Misfolded SOD1 propagates like prions from cell to cell. These modified conformations are detected in brain tissues in ALS, AD, and PD, and are considered a contributing factor to their initial pathogenesis. We have also elaborated on oxidative stress-induced non-native modifications of SOD1 and offered a logistic argument on their global implication in accelerated or pathological aging in the context of ALS, AD, and PD.

Photo-induced radical thiol-ene chemistry: a versatile toolbox for peptide-based drug design

While the global market for peptide/protein-based therapeutics is witnessing significant growth, the development of peptide drugs remains challenging due to their low oral bioavailability, poor membrane permeability, and reduced metabolic stability. However, a toolbox of chemical approaches has been explored for peptide modification to overcome these obstacles. In recent years, there has been a revival of interest in photoinduced radical thiol-ene chemistry as a powerful tool for the construction of therapeutic peptides.

In-depth interrogation of protein thermal unfolding data with MoltenProt

Protein stability is a key factor in successful structural and biochemical research. However, the approaches for systematic comparison of protein stability are limited by sample consumption or compatibility with sample buffer components. Here we describe how miniaturized measurement of intrinsic tryptophan fluorescence (NanoDSF assay) in combination with a simplified description of protein unfolding can be used to interrogate the stability of a protein sample. We demonstrate that improved protein stability measures, such as apparent Gibbs free energy of unfolding, rather than melting temperature Tm , should be used to rank the results of thermostability screens. The assay is compatible with protein samples of any composition, including protein complexes and membrane proteins. Our data analysis software, MoltenProt, provides an easy and robust way to perform characterization of multiple samples. Potential applications of MoltenProt and NanoDSF include buffer and construct optimization for X-ray crystallography and cryo-electron microscopy, screening for small-molecule binding partners and comparison of effects of point mutations.

Exploring conformational preferences of proteins: ionic liquid effects on the energy landscape of avidin

In this work we experimentally investigate solvent and temperature induced conformational transitions of proteins and examine the role of ion–protein interactions in determining the conformational preferences of avidin, a homotetrameric glycoprotein, in choline-based ionic liquid (IL) solutions. Avidin was modified by surface cationisation and the addition of anionic surfactants, and the structural, thermal, and conformational stabilities of native and modified avidin were examined using dynamic light scattering, differential scanning calorimetry, and thermogravimetric analysis experiments. The protein-surfactant nanoconjugates showed higher thermostability behaviour compared to unmodified avidin, demonstrating distinct conformational ensembles. Small-angle X-ray scattering data showed that with increasing IL concentration, avidin became more compact, interpreted in the context of molecular confinement. To experimentally determine the detailed effects of IL on the energy landscape of avidin, differential scanning fluorimetry and variable temperature circular dichroism spectroscopy were performed. We show that different IL solutions can influence avidin conformation and thermal stability, and we provide insight into the effects of ILs on the folding pathways and thermodynamics of proteins. To further study the effects of ILs on avidin binding and correlate thermostability with conformational heterogeneity, we conducted a binding study. We found the ILs examined inhibited ligand binding in native avidin while enhancing binding in the modified protein, indicating ILs can influence the conformational stability of the distinct proteins differently. Significantly, this work presents a systematic strategy to explore protein conformational space and experimentally detect and characterise ‘invisible’ rare conformations using ILs.

Revealing solvent and temperature induced conformational transitions of proteins and the role of ion–protein interactions in determining the conformational preferences of avidin in ionic liquids.

Optimization and bioactivity verification of porcine recombinant visfatin with high expression and low endotoxin content using pig liver as template

In order to obtain the porcine recombinant visfatin protein with high expression and low endotoxin content, the current study aims to express and verify the biological activity of the purified porcine recombinant visfatin protein. Firstly, four different expression strains were successfully constructed. Then they were simultaneously induced at 37 °C for 4 h and 16 °C for 16 h. The results showed that Visfatin-pET28a-Transetta was the best strain with high protein expression and purity at 16 °C induction for 16 h. After that, endotoxin was reduced from the recombinant visfatin until the residual endotoxin was less than one endotoxin units per milliliter (EU/mL). Finally, the purified porcine recombinant visfatin protein was incubated with RAW264.7 cells. The results of cell counting kit-8 (CCK-8) showed the survival rate of the cells first increased and then decreased with the increase in visfatin concentration. When the concentration of visfatin was 700 ng/mL, the survival rate of the cells was the highest. Thereafter, control (PBS), Visfatin and Visfatin + PolymyxinB (Ploy.B) groups were incubated with the RAW264.7 cells for 6 h. Real-time quantitative polymerase chain reaction (RT-qPCR) and Enzyme Linked Immuno-Sorbent Assay (ELISA) results showed that, as compared to the control group, the expressions of interleukin (IL)-1β, tumor necrosis factor (TNF)-α and monocyte chemoattractant protein (MCP)-1 in Visfatin group were significantly increased (P < 0.05). However, there was no significant difference between the Visfatin and Visfatin + Poly.B groups, indicating that porcine recombinant visfatin protein promoted the inflammatory activity of RAW264.7 cells while the residual endotoxin did not play a role, suggesting biological activity of porcine recombinant visfatin protein.

Biochemical characterization of a glycosyltransferase Gtf3 from Mycobacterium smegmatis: a case study of improved protein solubilization

Glycosyltransferases (GTs) are widely present in several organisms. These enzymes specifically transfer sugar moieties to a range of substrates. The processes of bacterial glycosylation of the cell wall and their relations with host-pathogen interactions have been studied extensively, yet the majority of mycobacterial GTs involved in the cell wall synthesis remain poorly characterized. Glycopeptidolipids (GPLs) are major class of glycolipids present on the cell wall of various mycobacterial species. They play an important role in drug resistance and host-pathogen interaction virulence. Gtf3 enzyme performs a key step in the biosynthesis of triglycosylated GPLs. Here, we describe a general procedure to achieve expression, purification, and crystallization of recombinant protein Gtf3 from Mycobacterium smegmatis using an E. coli expression system. We reported also a combined bioinformatics and biochemical methods to predict aggregation propensity and improve protein solubilization of recombinant Gtf3. NVoy, a carbohydrate-based polymer reagent, was added to prevent protein aggregation by binding to hydrophobic protein surfaces of Gtf3. Using intrinsic tryptophan fluorescence quenching experiments, we also demonstrated that Gtf3-NVoy enzyme interacted with TDP and UDP nucleotide ligands. This case report proposes useful tools for the study of other glycosyltransferases which are rather difficult to characterize and crystallize.

BCrystal: an interpretable sequence-based protein crystallization predictor

MOTIVATION

X-ray crystallography has facilitated the majority of protein structures determined to date. Sequence-based predictors that can accurately estimate protein crystallization propensities would be highly beneficial to overcome the high expenditure, large attrition rate, and to reduce the trial-and-error settings required for crystallization.

RESULTS

In this study, we present a novel model, BCrystal, which uses an optimized gradient boosting machine (XGBoost) on sequence, structural and physio-chemical features extracted from the proteins of interest. BCrystal also provides explanations, highlighting the most important features for the predicted crystallization propensity of an individual protein using the SHAP algorithm. On three independent test sets, BCrystal outperforms state-of-the-art sequence-based methods by more than 12.5% in accuracy, 18% in recall and 0.253 in Matthew's correlation coefficient, with an average accuracy of 93.7%, recall of 96.63% and Matthew's correlation coefficient of 0.868. For relative solvent accessibility of exposed residues, we observed higher values to associate positively with protein crystallizability and the number of disordered regions, fraction of coils and tripeptide stretches that contain multiple histidines associate negatively with crystallizability. The higher accuracy of BCrystal enables it to accurately screen for sequence variants with enhanced crystallizability.

AVAILABILITY AND IMPLEMENTATION

Our BCrystal webserver is at https://machinelearning-protein.qcri.org/ and source code is available at https://github.com/raghvendra5688/BCrystal.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

In silico analysis predicting effects of deleterious SNPs of human RASSF5 gene on its structure and functions

Ras association domain-containing protein 5 (RASSF5), one of the prospective biomarkers for tumors, generally plays a crucial role as a tumor suppressor. As deleterious effects can result from functional differences through SNPs, we sought to analyze the most deleterious SNPs of RASSF5 as well as predict the structural changes associated with the mutants that hamper the normal protein-protein interactions. We adopted both sequence and structure based approaches to analyze the SNPs of RASSF5 protein. We also analyzed the putative post translational modification sites as well as the altered protein-protein interactions that encompass various cascades of signals. Out of all the SNPs obtained from the NCBI database, only 25 were considered as highly deleterious by six in silico SNP prediction tools. Among them, upon analyzing the effect of these nsSNPs on the stability of the protein, we found 17 SNPs that decrease the stability. Significant deviation in the energy minimization score was observed in P350R, F321L, and R277W. Besides this, docking analysis confirmed that P350R, A319V, F321L, and R277W reduce the binding affinity of the protein with H-Ras, where P350R shows the most remarkable deviation. Protein-protein interaction analysis revealed that RASSF5 acts as a hub connecting two clusters consisting of 18 proteins and alteration in the RASSF5 may lead to disassociation of several signal cascades. Thus, based on these analyses, our study suggests that the reported functional SNPs may serve as potential targets for different proteomic studies, diagnosis and therapeutic interventions.

Current pivotal strategies leading a difficult target protein to a sample suitable for crystallographic analysis

Crystallographic structural analysis is an essential method for the determination of protein structure. However, crystallization of a protein of interest is the most difficult process in the analysis. The process is often hampered during the sample preparation, including expression and purification. Even after a sample has been purified, not all candidate proteins crystallize. In this mini-review, the current methodologies used to overcome obstacles encountered during protein crystallization are sorted. Specifically, the strategy for an effective crystallization is compared with a pipeline where various expression hosts and constructs, purification and crystallization conditions, and crystallization chaperones as target-specific binder proteins are assessed by a precrystallization screening. These methodologies are also developed continuously to improve the process. The described methods are useful for sample preparation in crystallographic analysis and other structure determination techniques, such as cryo-electron microscopy.

Insight into the structural and functional analysis of the impact of missense mutation on cytochrome P450 oxidoreductase

Cytochrome P450 oxidoreductase (POR) is a steroidogenic and drug-metabolizing enzyme. It helps in the NADPH dependent transfer of electrons to cytochrome P450 (CYP) enzymes for their biological activity. In this study, we employed integrative computational approaches to decipher the impact of proline to leucine missense mutation at position 384 (P384L) in the connecting/hinge domain region which is essential for the catalytic activity of POR. Analysis of protein stability using DUET, MUpro, CUPSAT, I-Mutant2.0, iStable and SAAFEC servers predicted that mutation might alter the structural stability of POR. The significant conformational changes induced by the mutation to the POR structure were analyzed by long-range molecular dynamics simulation. The results revealed that missense mutation decreased the conformational stability of POR as compared to wild type (WT). The PCA based FEL analysis described the mutant-specific conformational alterations and dominant motions essential for the biological activity of POR. The LIGPLOT interaction analysis showed the different binding architecture of FMN, FAD, and NADPH as a result of mutation. The increased number of hydrogen bonds in the FEL conformation of WT proved the strong binding of cofactors in the binding pocket as compared to the mutant. The porcupine plot analysis associated with cross-correlation analysis depicted the high-intensity flexible motion exhibited by functionally important FAD and NADPH binding domain regions. The computational findings unravel the impact of mutation at the structural level, which could be helpful in understanding the molecular mechanism of drug metabolism.

Common mechanism of thermostability in small α- and β-proteins studied by molecular dynamics

Protein thermostability is important to evolution, diseases, and industrial applications. Proteins use diverse molecular strategies to achieve stability at high temperature, yet reducing the entropy of unfolding seems required. We investigated five small α-proteins and five β-proteins with known, distinct structures and thermostability (T_m ) using multi-seed molecular dynamics simulations at 300, 350, and 400 K. The proteins displayed diverse changes in hydrogen bonding, solvent exposure, and secondary structure with no simple relationship to T_m . Our dynamics were in good agreement with experimental B-factors at 300 K and insensitive to force-field choice. Despite the very distinct structures, the native-state (300 + 350 K) free-energy landscapes (FELs) were significantly broader for the two most thermostable proteins and smallest for the three least stable proteins in both the α- and β-group and with both force fields studied independently (tailed t-test, 95% confidence level). Our results suggest that entropic ensembles stabilize proteins at high temperature due to reduced entropy of unfolding, viz., ΔG = ΔH - TΔS. Supporting this mechanism, the most thermostable proteins were also the least kinetically stable, consistent with broader FELs, typified by villin headpiece and confirmed by specific comparison to a mesophilic ortholog of Thermus thermophilus apo-pyrophosphate phosphohydrolase. We propose that molecular strategies of protein thermostabilization, although diverse, tend to converge toward highest possible entropy in the native state consistent with the functional requirements. We speculate that this tendency may explain why many proteins are not optimally structured and why molten-globule states resemble native proteins so much.

MPTherm: database for membrane protein thermodynamics for understanding folding and stability

The functions of membrane proteins (MPs) are attributed to their structure and stability. Factors influencing the stability of MPs differ from globular proteins due to the presence of membrane spanning regions. Thermodynamic data of MPs aid to understand the relationship among their structure, stability and function. Although a wealth of experimental data on thermodynamics of MPs are reported in the literature, there is no database available explicitly for MPs. In this work, we have developed a database for MP thermodynamics, MPTherm, which contains more than 7000 thermodynamic data from about 320 MPs. Each entry contains protein sequence and structural information, membrane topology, experimental conditions, thermodynamic parameters such as melting temperature, free energy, enthalpy etc. and literature information. MPTherm assists users to retrieve the data by using different search and display options. We have also provided the sequence and structure visualization as well as cross-links to UniProt and PDB databases. MPTherm database is freely available at http://www.iitm.ac.in/bioinfo/mptherm/. It is implemented in HTML, PHP, MySQL and JavaScript, and supports the latest versions of major browsers, such as Firefox, Chrome and Opera. MPTherm would serve as an effective resource for understanding the stability of MPs, development of prediction tools and identifying drug targets for diseases associated with MPs.

In Situ Proteolysis Condition-Induced Crystallization of the XcpVWX Complex in Different Lattices

Although prevalent in the determination of protein structures; crystallography always has the bottleneck of obtaining high-quality protein crystals for characterizing a wide range of proteins; especially large protein complexes. Stable fragments or domains of proteins are more readily to crystallize; which prompts the use of in situ proteolysis to remove flexible or unstable structures for improving crystallization and crystal quality. In this work; we investigated the effects of in situ proteolysis by chymotrypsin on the crystallization of the XcpVWX complex from the Type II secretion system of Pseudomonas aeruginosa. Different proteolysis conditions were found to result in two distinct lattices in the same crystallization solution. With a shorter chymotrypsin digestion at a lower concentration; the crystals exhibited a P3 hexagonal lattice that accommodates three complex molecules in one asymmetric unit. By contrast; a longer digestion with chymotrypsin of a 10-fold higher concentration facilitated the formation of a compact P2₁2₁2₁ orthorhombic lattice with only one complex molecule in each asymmetric unit. The molecules in the hexagonal lattice have shown high atomic displacement parameter values compared with the ones in the orthorhombic lattice. Taken together; our results clearly demonstrate that different proteolysis conditions can result in the generation of distinct lattices in the same crystallization solution; which can be exploited in order to obtain different crystal forms of a better quality

Phytocystatins and their Potential Application in the Development of Drought Tolerance Plants in Soybeans (Glycine max L.)

Plant cystatins, also called phytocystatins constitute a family of specific cysteine protease inhibitors found in several monocots and dicots. In soybean, phytocystatins regulate several endogenous processes contributing immensely to this crop's tolerance to abiotic stress factors. Soybeans offer numerous nutritional, pharmaceutical and industrial benefits; however, their growth and yields is hampered by drought, which causes more than 10% yield losses recorded every harvest period worldwide. This review analyses the role of papain-like cysteine proteases and their inhibitors in soybean plant growth and development under drought stress. It also describes their localisation, regulation, target organs and tissues, and the overall impact of cystatins on generating drought tolerance soybean plants. These proteins have many functions that remain poorly characterized, particularly under abiotic stress. Although much information is available on the utilisation of proteases for industrial applications, very few reports have focused on the impact of proteases on plant stress responses. The exploitation of cystatins in plant engineering, as competitive proteases inhibitors is one of the means that will guarantee the continued utilisation of soybeans as an important oilseed crop.

Selection and screening strategies in directed evolution to improve protein stability

Protein stability is not only fundamental for experimental, industrial, and therapeutic applications, but is also the baseline for evolving novel protein functions. For decades, stability engineering armed with directed evolution has continued its rapid development and inevitably poses challenges. Generally, in directed evolution, establishing a reliable link between a genotype and any interpretable phenotype is more challenging than diversifying genetic libraries. Consequently, we set forth in a small picture to emphasize the screening or selection techniques in protein stability-directed evolution to secure the link. For a more systematic review, two main branches of these techniques, namely cellular or cell-free display and stability biosensors, are expounded with informative examples.

Multiple selective sweeps of ancient polymorphisms in and around LTα located in the MHC class III region on chromosome 6

BACKGROUND

Lymphotoxin-α (LTα), located in the Major Histocompatibility Complex (MHC) class III region on chromosome 6, encodes a cytotoxic protein that mediates a variety of antiviral responses among other biological functions. Furthermore, several genotypes at this gene have been implicated in the onset of a number of complex diseases, including myocardial infarction, autoimmunity, and various types of cancer. However, little is known about levels of nucleotide variation and linkage disequilibrium (LD) in and near LTα, which could also influence phenotypic variance. To address this gap in knowledge, we examined sequence variation across ~ 10 kilobases (kbs), encompassing LTα and the upstream region, in 2039 individuals from the 1000 Genomes Project originating from 21 global populations.

RESULTS

Here, we observed striking patterns of diversity, including an excess of intermediate-frequency alleles, the maintenance of multiple common haplotypes and a deep coalescence time for variation (dating > 1.0 million years ago), in global populations. While these results are generally consistent with a model of balancing selection, we also uncovered a signature of positive selection in the form of long-range LD on chromosomes with derived alleles primarily in Eurasian populations. To reconcile these findings, which appear to support different models of selection, we argue that selective sweeps (particularly, soft sweeps) of multiple derived alleles in and/or near LTα occurred in non-Africans after their ancestors left Africa. Furthermore, these targets of selection were predicted to alter transcription factor binding site affinity and protein stability, suggesting they play a role in gene function. Additionally, our data also showed that a subset of these functional adaptive variants are present in archaic hominin genomes.

CONCLUSIONS

Overall, this study identified candidate functional alleles in a biologically-relevant genomic region, and offers new insights into the evolutionary origins of these loci in modern human populations.

The Methyltransferase Region of Vesicular Stomatitis Virus L Polymerase Is a Target Site for Functional Intramolecular Insertion

The L-protein of vesicular stomatitis virus (VSV) is a single-chain multi-domain RNA-dependent RNA polymerase. Previously reported attempts of intramolecular insertions of fluorescent proteins into the L-protein resulted in temperature-sensitive and highly attenuated polymerase activity. Here, we describe a novel insertion site that was selected based on in silico prediction. Of five preselected locations, insertion of the fluorescent protein mCherry in the VSV polymerase between amino acids 1620 and 1621 preserved polymerase function even after extended passaging and showed only mild attenuation compared to wildtype VSV polymerase. High magnification fluorescence imaging revealed a corpuscular cytosolic pattern for the L-protein. To confirm that the insertion site tolerates inclusion of proteins others than mCherry, we cloned mWasabi into the same position in L, generating a VSV-LmWasabi, which was also functional. We also generated a functional dual-color-dual-insertion VSV construct with intramolecularly labeled P and L-proteins. Together, our data present an approach to tag VSV polymerase intramolecularly without perturbing enzymatic activity. This L fusion protein might enable future tracing studies to monitor intracellular location of the VSV transcription and replication machinery in real-time life-imaging studies.

Molecular modeling and MD-simulation studies: Fast and reliable tool to study the role of low-redox bacterial laccases in the decolorization of various commercial dyes

Synthetic dyes are toxic and carcinogenic in nature, which also causes environmental pollution. The present study was aimed to decolorize various commercial dyes using purified recombinant bacterial laccases. Laccase gene from Yersinia enterocolitica strain 8081 (yacK), Y. enterocolitica strain 7 (yacK) and Bacillus pumilus DSKK1 was cloned in vector pET28a and overproduced in host Escherichia coli BL21. The high yield of recombinant laccase protein resulted in the formation of inclusion bodies, which were further solubilized, refolded, and purified. The purified recombinant laccases were alkali-tolerant and thermostable, with pH optima at 7-8, temperature optima at 60-70 °C and low redox potential. For in silico studies, laccase protein models of B. pumilus DSKK1, Y. enterocolitica strain 7 and Y. enterocolitica strain 8081 were docked with commercial dyes. This is the first and foremost study where the stability of docked complexes of pathogenic and non-pathogenic microorganism has been explored via molecular dynamics (MD) simulations using Gromacs version 4.5.5 with the gromos96 43a force field. Finally, the in silico results were validated experimentally and it was found that purified laccases from B. pumilus DSKK1 and Y. enterocolitica strain 7 efficiently decolorized rose bengal (90.4%), malachite green (77.7%), and congo red (74.5%) dyes.

Accurate Calculation of Barnase and SNase Folding Energetics Using Short Molecular Dynamics Simulations and an Atomistic Model of the Unfolded Ensemble: Evaluation of Force Fields and Water Models

As proteins perform most cellular functions, quantitative understanding of protein energetics is required to gain control of biological phenomena. Accurate models of native proteins can be obtained experimentally, but the lack of equally fine models of unfolded ensembles impedes the calculation of protein folding energetics from first principles. Here, we show that an atomistic unfolded ensemble model, consisting of a few dozen conformations built from a protein sequence, can be used in conjunction with an X-ray structure of its native state to calculate accurately by difference the changes in enthalpy and heat capacity of the polypeptide upon folding. The calculation is done using molecular dynamics simulations, popular force fields, and water models, and for the two model proteins studied (barnase and SNase), the results agree within error or are very close to their experimentally determined properties. The enthalpy sampling of the unfolded ensemble is done through short 2 ns simulations that do not significantly modify the representative distribution of R_g of the starting conformations. The impressive accuracy obtained opens the possibility to investigate quantitatively systems or phenomena not amenable to experiment and paves the way for addressing the calculation of protein conformational stability (i.e., the change in Gibbs energy upon folding), a central goal of structural biology. So far, these calculated enthalpy and heat capacity changes, combined with the experimentally determined melting temperatures of the corresponding protein, allow us to reproduce the stability curves of both barnase and SNase.

Emerging Role of Migration and Invasion Enhancer 1 (MIEN1) in Cancer Progression and Metastasis

Tumor metastasis is a sequential event accounting for numerous cancer-related fatalities worldwide. The process of metastasis serially involves invasion, intravasation, extravasation, and tumor growth at the secondary site. Migration and invasion enhancer 1 (MIEN1) is a membrane associated protein overexpressed in various human cancers. Biological activity of MIEN1 is driven by geranylgeranyltransferase-I mediated prenylation at CAAX motif and methylation of the prenylated protein that anchors MIEN1 into the cellular membrane. Post-translationally modified MIEN1 interacts with Syk kinase and Annexin A2 protein; polymerizes G-actin and stabilizes F-actin filament; induces focal adhesion kinase phosphorylation and decrease cofilin phosphorylation implicated in both invasion and metastasis of different cancer types. In the present review, we discuss the structure, function, and involvement of MIEN1 in cancer progression. We also highlight the future prospects of MIEN1 as an emerging molecule and novel target in cancer cell invasion and metastasis.

Next Generation Techniques for Determination of Protein-Protein Interactions: Beyond the Crystal Structure

PURPOSE OF REVIEW

We discuss recent advancements in structural biology methods for investigating sites of protein-protein interactions. We will inform readers outside the field of structural biology about techniques beyond crystallography, and how these different technologies can be utilized for drug development.

RECENT FINDINGS

Advancements in cryo-electron microscopy (cryoEM) and micro-electron diffraction (microED) may change how we view atomic resolution structural biology, such that well-ordered macrocrystals of protein complexes are not required for interface identification. However, some drug discovery applications, such as lead peptide compound generation, may not require atomic resolution; mass spectrometry techniques can provide an expedited path to generation of lead compounds. New crosslinking compounds, more user-friendly data analysis, and novel protocols such as protein painting can advance drug discovery programs, even in the absence of atomic resolution structural data. Finally, artificial intelligence and machine learning methods, while never truly replacing experimental methods, may provide rational ways to stratify potential druggable regions identified with mass spectrometry into higher and lower priority candidates.

SUMMARY

Electron diffraction of nanocrystals combines the benefits of both x-ray diffraction and cryoEM, and may prove to be the next generation of atomic resolution protein-protein interface identification. However, in situations such as peptide drug discovery, mass spectrometry techniques supported by advancements in computational methods will likely prove sufficient to support drug discovery efforts. In addition, these methods can be significantly faster than any crystallographic or cryoEM methods for identification of interacting regions.

Structure Determination by Single-Particle Cryo-Electron Microscopy: Only the Sky (and Intrinsic Disorder) is the Limit

Traditionally, X-ray crystallography and NMR spectroscopy represent major workhorses of structural biologists, with the lion share of protein structures reported in protein data bank (PDB) being generated by these powerful techniques. Despite their wide utilization in protein structure determination, these two techniques have logical limitations, with X-ray crystallography being unsuitable for the analysis of highly dynamic structures and with NMR spectroscopy being restricted to the analysis of relatively small proteins. In recent years, we have witnessed an explosive development of the techniques based on Cryo-electron microscopy (Cryo-EM) for structural characterization of biological molecules. In fact, single-particle Cryo-EM is a special niche as it is a technique of choice for the structural analysis of large, structurally heterogeneous, and dynamic complexes. Here, sub-nanometer atomic resolution can be achieved (i.e., resolution below 10 Å) via single-particle imaging of non-crystalline specimens, with accurate 3D reconstruction being generated based on the computational averaging of multiple 2D projection images of the same particle that was frozen rapidly in solution. We provide here a brief overview of single-particle Cryo-EM and show how Cryo-EM has revolutionized structural investigations of membrane proteins. We also show that the presence of intrinsically disordered or flexible regions in a target protein represents one of the major limitations of this promising technique.

Genetically Engineered Proteins to Improve Biomass Conversion: New Advances and Challenges for Tailoring Biocatalysts

Protein engineering emerged as a powerful approach to generate more robust and efficient biocatalysts for bio-based economy applications, an alternative to ecologically toxic chemistries that rely on petroleum. On the quest for environmentally friendly technologies, sustainable and low-cost resources such as lignocellulosic plant-derived biomass are being used for the production of biofuels and fine chemicals. Since most of the enzymes used in the biorefinery industry act in suboptimal conditions, modification of their catalytic properties through protein rational design and in vitro evolution techniques allows the improvement of enzymatic parameters such as specificity, activity, efficiency, secretability, and stability, leading to better yields in the production lines. This review focuses on the current application of protein engineering techniques for improving the catalytic performance of enzymes used to break down lignocellulosic polymers. We discuss the use of both classical and modern methods reported in the literature in the last five years that allowed the boosting of biocatalysts for biomass degradation.

Hot CoFi Blot: A High-Throughput Colony-Based Screen for Identifying More Thermally Stable Protein Variants

Highly soluble and stable proteins are desirable for many different applications, from basic science to reaching a cancer patient in the form of a biological drug. For X-ray crystallography-where production of a protein crystal might take weeks and even months-a stable protein sample of high purity and concentration can greatly increase the chances of producing a well-diffracting crystal. For a patient receiving a specific protein drug, its safety, efficacy, and even cost are factors affected by its solubility and stability. Increased protein expression and protein stability can be achieved by randomly altering the coding sequence. As the number of mutants generated might be overwhelming, a powerful protein expression and stability screen is required. In this chapter, we describe a colony filtration technology, which allows us to screen random mutagenesis libraries for increased thermal stability-the Hot CoFi blot. We share how to create the random mutagenesis library, how to perform the Hot CoFi blot, and how to identify more thermally stable clones. We use the Tobacco Etch Virus protease as a target to exemplify the procedure.

In silico analysis of missense mutations in exons 1-5 of the F9 gene that cause hemophilia B

Background

Missense mutations in the first five exons of F9, which encodes factor FIX, represent 40% of all mutations that cause hemophilia B. To address the ongoing debate regarding in silico identification of disease-causing mutations at these exons, we analyzed 215 missense mutations from www.factorix.org using six in silico prediction tools, which are the most common used programs for analysis prediction of impact of mutations on the protein structure and function, with further advantage of using similar approaches. We developed different algorithms to integrate multiple predictions from such tools. In order to approach a structural analysis on FIX we performed a modeling of five selected pathogenic mutations.

Results

SIFT, PolyPhen-2 HumDiv, SNAP2, and MutationAssessor were the most successful in identifying true non-causative and causative mutations. A proposed function integrating these algorithms (wgP4) was the most sensitive (90.1%), specific (22.6%), and accurate (87%) than similar functions, and identified 187 variants as deleterious. Clinical phenotype was significantly associated with predicted causative mutations at all five exons. However, PolyPhen-2 HumDiv was more successful in linking clinical severity to specific exons, while functions that integrate 4–6 predictions were more successful in linking phenotype to genotypes at the light chain (exons 3–5). The most important value of integrating multiple predictions is the inclusion of scores derived from different approaches. Modeling of protein structure showed the effects of pathogenic nsSNPs on structure and function of FIX.

Conclusions

A simple function that integrates information from different in silico programs yields the best prediction of mutated phenotypes. However, the specificity, sensitivity, and accuracy of genotype-phenotype predictions depend on specific characteristics of the protein domain and the disease of interest as we validated by the structural analysis of selected pathogenic F9 mutations. The proposed function integrating algorithm (wgP4) might be useful for the analysis of nsSNPs impact on other genes.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-2919-x) contains supplementary material, which is available to authorized users.

Evaluation of potential molecular interaction between quorum sensing receptor, LuxP and grouper fatty acids: in-silico screening and simulation

Pathologically relevant behaviors of Vibrio, such as the expression of virulence factors, biofilm production, and swarming motility, have been shown to be controlled by quorum sensing. The autoinducer-2 quorum sensing receptor protein LuxP is one of the target proteins for drug development to suppress the virulence of Vibrio. Here, we reported the potential molecular interaction of fatty acids identified in vibriosis-resistant grouper with LuxP. Fatty acid, 4-oxodocosahexaenoic acid (4R8) showed significant binding affinity toward LuxP (-6.0 kcal/mol) based on molecular docking analysis. The dynamic behavior of the protein-ligand complex was illustrated by molecular dynamic simulations. The fluctuation of the protein backbone, the stability of ligand binding, and hydrogen bond interactions were assessed, suggesting 4R8 possesses potential interaction with LuxP, which was supported by the low binding free energy (-29.144 kJ/mol) calculated using the molecular mechanics Poisson-Boltzmann surface area.

Dual-Family Peptidylprolyl Isomerases (Immunophilins) of Select Monocellular Organisms

The dual-family peptidylprolyl cis-trans isomerases (immunophilins) represent a naturally occurring chimera of the classical FK506-binding protein (FKBP) and cyclophilin (CYN), connected by a flexible linker. They are found exclusively in monocellular organisms. The modular builds of these molecules represent two distinct types: CYN-(linker)-FKBP and FKBP-3TPR (tetratricopeptide repeat)-CYN. Abbreviated respectively as CFBP and FCBP, the two classes also exhibit distinct organism preference, the CFBP being found in prokaryotes, and the FCBP in eukaryotes. This review summarizes the mystery of these unique class of prolyl isomerases, focusing on their host organisms, potential physiological role, and likely routes of evolution.

Transcriptome and Proteome of Fish-Pathogenic Streptococcus agalactiae Are Modulated by Temperature

Streptococcus agalactiae is one of the most important pathogens associated with streptococcosis outbreaks in Nile tilapia farms worldwide. High water temperature (above 27°C) has been described as a predisposing factor for the disease in fish. At low temperatures (below 25°C), fish mortalities are not usually observed in farms. Temperature variation can modulate the expression of genes and proteins involved in metabolism, adaptation, and bacterial pathogenicity, thus increasing or decreasing the ability to infect the host. This study aimed to evaluate the transcriptome and proteome of a fish-pathogenic S. agalactiae strain SA53 subjected to in vitro growth at different temperatures using a microarray and label-free shotgun LC-HDMS^E approach. Biological triplicates of isolates were cultured in BHIT broth at 22 or 32°C for RNA and protein isolation and submitted for transcriptomic and proteomic analyses. In total, 1,730 transcripts were identified in SA53, with 107 genes being differentially expressed between the temperatures evaluated. A higher number of genes related to metabolism, mainly from the phosphotransferase system (PTS) and ATP-binding cassette (ABC) transport system, were upregulated at 32°C. In the proteome analysis, 1,046 proteins were identified in SA53, of which 81 were differentially regulated between 22 and 32°C. Proteins involved in defense mechanisms, lipid transport and metabolism, and nucleotide transport and metabolism were upregulated at 32°C. A higher number of interactions were observed in proteins involved in nucleotide transport and metabolism. We observed a low correlation between the transcriptome and proteome datasets. Our study indicates that the transcriptome and proteome of a fish-adapted S. agalactiae strain are modulated by temperature, particularly showing differential expression of genes/proteins involved in metabolism, virulence factors, and adaptation.

Carbon metabolism influenced for promoters and temperature used in the heterologous protein production using Pichia pastoris yeast

Nowadays, it is necessary to search for different high-scale production strategies to produce recombinant proteins of economic interest. Only a few microorganisms are industrially relevant for recombinant protein production: methylotrophic yeasts are known to use methanol efficiently as the sole carbon and energy source. Pichia pastoris is a methylotrophic yeast characterized as being an economical, fast and effective system for heterologous protein expression. Many factors can affect both the product and the production, including the promoter, carbon source, pH, production volume, temperature, and many others; but to control all of them most of the time is difficult and this depends on the initial selection of each variable. Therefore, this review focuses on the selection of the best promoter in the recombination process, considering different inductors, and the temperature as a culture medium variable in methylotrophic Pichia pastoris yeast. The goal is to understand the effects associated with different factors that influence its cell metabolism and to reach the construction of an expression system that fulfills the requirements of the yeast, presenting an optimal growth and development in batch, fed-batch or continuous cultures, and at the same time improve its yield in heterologous protein production.

The Goldilocks Approach: A Review of Employing Design of Experiments in Prokaryotic Recombinant Protein Production

The production of high yields of soluble recombinant protein is one of the main objectives of protein biotechnology. Several factors, such as expression system, vector, host, media composition and induction conditions can influence recombinant protein yield. Identifying the most important factors for optimum protein expression may involve significant investment of time and considerable cost. To address this problem, statistical models such as Design of Experiments (DoE) have been used to optimise recombinant protein production. This review examines the application of DoE in the production of recombinant proteins in prokaryotic expression systems with specific emphasis on media composition and culture conditions. The review examines the most commonly used DoE screening and optimisation designs. It provides examples of DoE applied to optimisation of media and culture conditions.

Role of simple descriptors and applicability domain in predicting change in protein thermostability

The melting temperature (Tm) of a protein is the temperature at which half of the protein population is in a folded state. Therefore, Tm is a measure of the thermostability of a protein. Increasing the Tm of a protein is a critical goal in biotechnology and biomedicine. However, predicting the change in melting temperature (dTm) due to mutations at a single residue is difficult because it depends on an intricate balance of forces. Existing methods for predicting dTm have had similar levels of success using generally complex models. We find that training a machine learning model with a simple set of easy to calculate physicochemical descriptors describing the local environment of the mutation performed as well as more complicated machine learning models and is 2-6 orders of magnitude faster. Importantly, unlike in most previous publications, we perform a blind prospective test on our simple model by designing 96 variants of a protein not in the training set. Results from retrospective and prospective predictions reveal the limited applicability domain of each model. This study highlights the current deficiencies in the available dTm dataset and is a call to the community to systematically design a larger and more diverse experimental dataset of mutants to prospectively predict dTm with greater certainty.

Structures of NHBA elucidate a broadly conserved epitope identified by a vaccine induced antibody

Neisserial heparin binding antigen (NHBA) is one of three main recombinant protein antigens in 4CMenB, a vaccine for the prevention of invasive meningococcal disease caused by Neisseria meningitidis serogroup B. NHBA is a surface-exposed lipoprotein composed of a predicted disordered N-terminal region, an arginine-rich region that binds heparin, and a C-terminal domain that folds as an anti-parallel β-barrel and that upon release after cleavage by human proteases alters endothelial permeability. NHBA induces bactericidal antibodies in humans, and NHBA-specific antibodies elicited by the 4CMenB vaccine contribute to serum bactericidal activity, the correlate of protection. To better understand the structural bases of the human antibody response to 4CMenB vaccination and to inform antigen design, we used X-ray crystallography to elucidate the structures of two C-terminal fragments of NHBA, either alone or in complex with the Fab derived from the vaccine-elicited human monoclonal antibody 5H2, and the structure of the unbound Fab 5H2. The structures reveal details on the interaction between an N-terminal β-hairpin fragment and the β-barrel, and explain how NHBA is capable of generating cross-reactive antibodies through an extensive conserved conformational epitope that covers the entire C-terminal face of the β-barrel. By providing new structural information on a vaccine antigen and on the human immune response to vaccination, these results deepen our molecular understanding of 4CMenB, and might also aid future vaccine design projects.

The Peptidisc, a simple method for stabilizing membrane proteins in detergent-free solution

Membrane proteins are difficult to work with due to their insolubility in aqueous solution and quite often their poor stability in detergent micelles. Here, we present the peptidisc for their facile capture into water-soluble particles. Unlike the nanodisc, which requires scaffold proteins of different lengths and precise amounts of matching lipids, reconstitution of detergent solubilized proteins in peptidisc only requires a short amphipathic bi-helical peptide (NSP_r) and no extra lipids. Multiple copies of the peptide wrap around to shield the membrane-exposed part of the target protein. We demonstrate the effectiveness of this ‘one size fits all’ method using five different membrane protein assemblies (MalFGK₂, FhuA, SecYEG, OmpF, BRC) during ‘on-column’, ‘in-gel’, and ‘on-bead’ reconstitution embedded within the membrane protein purification protocol. The peptidisc method is rapid and cost-effective, and it may emerge as a universal tool for high-throughput stabilization of membrane proteins to advance modern biological studies.

Surrounding every living cell is a biological membrane that is largely impermeable to water-soluble molecules. This hydrophobic (or “water-hating”) barrier preserves the contents of the cell and also regulates how the cell interacts with its environment. This latter function is critical and relies on a class of proteins that are embedded within the membrane and are also hydrophobic.

The hydrophobic nature of membrane proteins is however inconvenient for biochemical studies which usually take place in water-based solutions. Therefore, membrane proteins are under-represented in biological research compared to the water-soluble ones, even though roughly one quarter of a cell’s proteins are membrane proteins. Researchers have developed a few tricks to keep membrane proteins soluble after they have been extracted from the membrane. An old but popular technique makes use of detergents, which are chemicals with opposing hydrophobic and hydrophilic properties (hydrophilic literally means “water-loving”). However, even mild detergents can damage membrane proteins and will sometimes lead to experimental artifacts. More recent tricks to stabilize membrane proteins without detergents have been described but remain laborious, costly or difficult to perform.

To overcome these limitations, Carlson et al. developed a simple method to stabilize membrane proteins without detergent. Called the “peptidisc”, the method uses multiple copies of a unique peptide – a short sequence of the building blocks of protein – that had been redesigned to have optimal hydrophobic and hydrophilic properties. The idea was that the peptides would wrap around the hydrophobic parts of the membrane protein, and shield them from the watery solution. Indeed, when Carlson et al. mixed this peptide with five different membrane proteins from bacteria, all were perfectly soluble and functional without detergent. The ideal ratio of peptide needed to form a peptidisc around each membrane protein was reached automatically, without having to test many different conditions. This indicates that the peptidisc acts like a “one size fits all” scaffold.

The peptidisc is a new tool that will allow more researchers, including those who are not expert biochemists, to study membrane proteins. This will yield a better understanding of the structure of a cell’s membrane and how it interacts with the environment. Since the approach is both simple and easy to apply, more membrane proteins can now also be included in high-throughput searches for potential new drugs for various medical conditions.

Stability of allergens

For proteins to cause IgE-mediated allergic reactions, several common characteristics have to be defined, including small molecular size, solubility and stability to changing pH levels and enzymatic degradation. Nevertheless, these features are not unique for potent allergens, but are also observed in non-allergenic proteins. Due to the increasing awareness by regulatory authorities regarding the allergy pandemic, definition of characteristics unique to potent allergens would facilitate allergenicity assessment in the future. Despite major research efforts even to date the features unique for major allergens have not been elucidated so far. The route of allergen entry into the organism determines to a great extent these required characteristics. Especially orally ingested allergens are exposed to the harsh milieu of the gastrointestinal tract but might additionally be influenced by food processing. Depending on molecular properties such as disulphide bonds contributing to protein fold and formation of conformational IgE epitopes, posttranslational protein modification or protein food matrix interactions, enzymatic and thermal stability might differ between allergens. Moreover, also ligand binding influences structural stability. In the current review article, we aim at highlighting specific characteristics and molecular pattern contributing to a stabilized protein structure and overall allergenicity.

eModel-BDB: a database of comparative structure models of drug-target interactions from the Binding Database

Background

The structural information on proteins in their ligand-bound conformational state is invaluable for protein function studies and rational drug design. Compared to the number of available sequences, not only is the repertoire of the experimentally determined structures of holo-proteins limited, these structures do not always include pharmacologically relevant compounds at their binding sites. In addition, binding affinity databases provide vast quantities of information on interactions between drug-like molecules and their targets, however, often lacking structural data. On that account, there is a need for computational methods to complement existing repositories by constructing the atomic-level models of drug-protein assemblies that will not be determined experimentally in the near future.

Results

We created eModel-BDB, a database of  200,005 comparative models of drug-bound proteins based on   1,391,403 interaction data obtained from the Binding Database and the PDB library of 31 January 2017. Complex models in eModel-BDB were generated with a collection of the state-of-the-art techniques, including protein meta-threading, template-based structure modeling, refinement and binding site detection, and ligand similarity-based docking. In addition to a rigorous quality control maintained during dataset generation, a subset of weakly homologous models was selected for the retrospective validation against experimental structural data recently deposited to the Protein Data Bank. Validation results indicate that eModel-BDB contains models that are accurate not only at the global protein structure level but also with respect to the atomic details of bound ligands.

Conclusions

Freely available eModel-BDB can be used to support structure-based drug discovery and repositioning, drug target identification, and protein structure determination.

Egress and invasion machinery of malaria: an in-depth look into the structural and functional features of the flap dynamics of plasmepsin IX and X

Plasmepsins, a family of aspartic proteases expressed by Plasmodium falciparum parasite, have been identified as key mediators in the onset of lethal malaria. Precedence has been placed on this family of enzymes due their essential role in the virulence of the parasite, thus highlighting their importance as novel drug targets. A previously published study by our group proposed a set of parameters used to define the flap motion of aspartic proteases. These parameters were used in the study of Plm I-V and focused on the flap flexibility as well as structural dynamics. Recent studies have highlighted the essential role played by Plm IX and X in egress and invasion of the malarial parasite. This study aims to close the gap on the latter family, investigating the flap dynamics of Plms IX and X. Molecular dynamics simulations demonstrated an "open and close" mechanism at the region of the catalytic site. Further computation of the dihedral angles at the catalytic region revealed tractability at both the flap tip and flexible loop. This structural versatility enhances the interaction of variant ligand sizes, in comparison to other Plm family members. The results obtained from this study signify the essential role of structural flap dynamics and its resultant effect on the binding landscapes of Plm IX and X. We believe that this unique structural mechanism may be integrated in the design and development of effective anti-malarial drugs.

Deciphering RNA-Recognition Patterns of Intrinsically Disordered Proteins

Intrinsically disordered regions (IDRs) and protein (IDPs) are highly flexible owing to their lack of well-defined structures. A subset of such proteins interacts with various substrates; including RNA; frequently adopting regular structures in the final complex. In this work; we have analysed a dataset of protein–RNA complexes undergoing disorder-to-order transition (DOT) upon binding. We found that DOT regions are generally small in size (less than 3 residues) for RNA binding proteins. Like structured proteins; positively charged residues are found to interact with RNA molecules; indicating the dominance of electrostatic and cation-π interactions. However, a comparison of binding frequency shows that interface hydrophobic and aromatic residues have more interactions in only DOT regions than in a protein. Further; DOT regions have significantly higher exposure to water than their structured counterparts. Interactions of DOT regions with RNA increase the sheet formation with minor changes in helix forming residues. We have computed the interaction energy for amino acids–nucleotide pairs; which showed the preference of His–G; Asn–U and Ser–U at for the interface of DOT regions. This study provides insights to understand protein–RNA interactions and the results could also be used for developing a tool for identifying DOT regions in RNA binding proteins.

Expression of an 8R-Lipoxygenase From the Coral Plexaura homomalla

Methods are presented for the use of the coral 8R-lipoxygenase from the Caribbean sea whip coral Plexaura homomalla as a model enzyme for structural studies of animal lipoxygenases. The 8R-lipoxygenase is remarkably stable and can be stored at 4°C for 3 months with virtually no loss of activity. In addition, an engineered "pseudo wild-type" enzyme is soluble in the absence of detergents, which helps facilitate the preparation of enzyme:substrate complexes.

New Bacterial Phytase through Metagenomic Prospection

Alkaline phytases from uncultured microorganisms, which hydrolyze phytate to less phosphorylated myo-inositols and inorganic phosphate, have great potential as additives in agricultural industry. The development of metagenomics has stemmed from the ineluctable evidence that as-yet-uncultured microorganisms represent the vast majority of organisms in most environments on earth. In this study, a gene encoding a phytase was cloned from red rice crop residues and castor bean cake using a metagenomics strategy. The amino acid identity between this gene and its closest published counterparts is lower than 60%. The phytase was named PhyRC001 and was biochemically characterized. This recombinant protein showed activity on sodium phytate, indicating that PhyRC001 is a hydrolase enzyme. The enzymatic activity was optimal at a pH of 7.0 and at a temperature of 35 °C. β-propeller phytases possess great potential as feed additives because they are the only type of phytase with high activity at neutral pH. Therefore, to explore and exploit the underlying mechanism for β-propeller phytase functions could be of great benefit to biotechnology.

Impact of an N-terminal Polyhistidine Tag on Protein Thermal Stability

For years, the use of polyhistidine tags (His-tags) has been a staple in the isolation of recombinant proteins in immobilized metal affinity chromatography experiments. Their usage has been widely beneficial in increasing protein purity from crude cell lysates. For some recombinant proteins, a consequence of His-tag addition is that it can affect protein function and stability. Functional proteins are essential in the elucidation of their biological, kinetic, structural, and thermodynamic properties. In this study, we determine the effect of N-terminal His-tags on the thermal stability of select proteins using differential scanning fluorimetry and identify that the removal of the His-tag can have both beneficial and deleterious effects on their stability.