Template-based protein structure prediction in CASP11 and retrospect of I-TASSER in the last decade

Design of a novel multi-epitopes vaccine against Escherichia fergusonii: a pan-proteome based in- silico approach

Escherichia fergusonii a gram-negative rod-shaped bacterium in the Enterobacteriaceae family, infect humans, causing serious illnesses such as urinary tract infection, cystitis, biliary tract infection, pneumonia, meningitis, hemolytic uremic syndrome, and death. Initially treatable with penicillin, antibiotic misuse led to evolving resistance, including resistance to colistin, a last-resort drug. With no licensed vaccine, the study aimed to design a multi-epitope vaccine against E. fergusonii. The study started with the retrieval of the complete proteome of all known strains and proceeded to filter the surface exposed virulent proteins. Seventeen virulent proteins (4 extracellular, 4 outer membranes, 9 periplasmic) with desirable physicochemical properties were identified from the complete proteome of known strains. Further, these proteins were processed for B-cell and T-cell epitope mapping. Obtained epitopes were evaluated for antigenicity, allergenicity, solubility, MHC-binding, and toxicity and the filtered epitopes were fused by specific linkers and an adjuvant into a vaccine construct. Structure of the vaccine candidate was predicted and refined resulting in 78.1% amino acids in allowed regions and VERIFY3D score of 81%. Vaccine construct was docked with TLR-4, MHC-I, and MHC-II, showing binding energies of -1040.8 kcal/mol, -871.4 kcal/mol, and -1154.6 kcal/mol and maximum interactions. Further, molecular dynamic simulation of the docked complexes was carried out resulting in a significant stable nature of the docked complexes (high B-factor and deformability values, lower Eigen and high variance values) in terms of intermolecular binding conformation and interactions. The vaccine was also reported to stimulate a variety of immunological pathways after administration. In short, the designed vaccine revealed promising predictions about its immune protective potential against E. fergusonii infections however experimental validation is needed to validate the results.

Rapid Generation and Molecular Docking Analysis of Single-Chain Fragment Variable (scFv) Antibody Selected by Ribosome Display Targeting Cholecystokinin B Receptor (CCK-BR) for Reduction of Chronic Neuropathic Pain

A robust cell-free platform technology, ribosome display in combination with cloning, expression, and purification was utilized to develop single chain Fragment variable (scFv) antibody variants as pain therapy directed at the mouse cholecystokinin B (CCK-B) receptor. Three effective CCK-B peptide-specific scFvs were generated through ribosomal display technology. Soluble expression and ELISA analysis showed that one antibody, scFv77-2 had the highest binding and could be purified from bacterial cells in large quantities. Octet measurements further revealed that the CCK-B scFv77-2 antibody had binding kinetics of K_D = 1.794 × 10^-8 M. Molecular modeling and docking analyses suggested that the scFv77-2 antibody shaped a proper cavity to embed the whole CCK-B peptide molecule and that a steady-state complex was formed relying on intermolecular forces, including hydrogen bonding, electrostatic force, and hydrophobic interactions. Thus, the scFv antibody can be applied for mechanistic intermolecular interactions and functional in vivo studies of CCK-BR. The high affinity scFv77-2 antibody showed good efficacy with binding to CCK-BR tested in a chronic pain model. In vivo studies validated the efficacy of the CCK-B receptor (CCK-BR) scFv77-2 antibody as a potential therapy for chronic trigeminal nerve injury-induced pain. Mice were given a single dose of the CCK-B receptor (CCK-BR) scFv antibody 3 weeks after induction of a chronic trigeminal neuropathic pain model, during the transition from acute to chronic pain. The long-term effectiveness for the reduction of mechanical hypersensitivity was evident, persisting for months. The anxiety- and depression-related behaviors typically accompanying persisting hypersensitivity subsequently never developed in the mice given CCK-BR scFv. The effectiveness of the antibody is the basis for further development of the lead CCK-BR scFv as a promising non-opioid therapeutic for chronic pain and the long-term reduction of chronic pain- and anxiety-related behaviors.

Application of Molecular Simulation Methods in Food Science: Status and Prospects

Molecular simulation methods, such as molecular docking, molecular dynamic (MD) simulation, and quantum chemical (QC) calculation, have become popular as characterization and/or virtual screening tools because they can visually display interaction details that in vitro experiments can not capture and quickly screen bioactive compounds from large databases with millions of molecules. Currently, interdisciplinary research has expanded molecular simulation technology from computer aided drug design (CADD) to food science. More food scientists are supporting their hypotheses/results with this technology. To understand better the use of molecular simulation methods, it is necessary to systematically summarize the latest applications and usage trends of molecular simulation methods in the research field of food science. However, this type of review article is rare. To bridge this gap, we have comprehensively summarized the principle, combination usage, and application of molecular simulation methods in food science. We also analyzed the limitations and future trends and offered valuable strategies with the latest technologies to help food scientists use molecular simulation methods.

Combination of Docking-Based and Pharmacophore-Based Virtual Screening Identifies Novel Agonists That Target the Urotensin Receptor

The urotensin receptor (UT receptor), a G-protein-coupled receptor mediating urotensin-II and urotensin-II-related peptide signaling in the urotensinergic system, has multiple pharmacological activities. However, there is no drug targeting the UT receptor currently in clinical use, and the discovery of new leads is still important. The complete crystal structure of the UT receptor has not yet been resolved and a screening strategy combining multiple methods can improve the accuracy and efficiency of drug screening. This study aimed to identify novel UT receptor agonists using a combination of docking-based, pharmacophore-based, and cell-based drug screening. First, the three-dimensional structures of the UT receptor were constructed through single-template, multi-template homologous modeling and threading strategies. After structure evaluation and ligand enrichment analysis, a model from the threading modeling was selected for docking-based virtual screening based on stepwise filtering, and 1368 positive compounds were obtained from our compound library. Second, the pharmacophore models were constructed using known ligands targeting the UT receptor for pharmacophore-based virtual screening. A model was selected after model validation, and 300 positive compounds were retrieved. Then, after intersecting the results of two different virtual screening methods with 570 compound entities from our primary screening, 14 compounds were obtained. Finally, three hits were obtained after in vitro confirmation. Furthermore, preliminary evaluation of the hits showed that they influenced glucose consumption. In summary, by integrating docking-based, pharmacophore-based, and in vitro drug screening, three new agonists targeting the UT receptor were identified which may serve as promising therapeutic agents for urotensinergic system disorders.

A Deep-Learning Proteomic-Scale Approach for Drug Design

Computational approaches have accelerated novel therapeutic discovery in recent decades. The Computational Analysis of Novel Drug Opportunities (CANDO) platform for shotgun multitarget therapeutic discovery, repurposing, and design aims to improve their efficacy and safety by employing a holistic approach that computes interaction signatures between every drug/compound and a large library of non-redundant protein structures corresponding to the human proteome fold space. These signatures are compared and analyzed to determine if a given drug/compound is efficacious and safe for a given indication/disease. In this study, we used a deep learning-based autoencoder to first reduce the dimensionality of CANDO-computed drug-proteome interaction signatures. We then employed a reduced conditional variational autoencoder to generate novel drug-like compounds when given a target encoded "objective" signature. Using this approach, we designed compounds to recreate the interaction signatures for twenty approved and experimental drugs and showed that 16/20 designed compounds were predicted to be significantly (p-value ≤ 0.05) more behaviorally similar relative to all corresponding controls, and 20/20 were predicted to be more behaviorally similar relative to a random control. We further observed that redesigns of objectives developed via rational drug design performed significantly better than those derived from natural sources (p-value ≤ 0.05), suggesting that the model learned an abstraction of rational drug design. We also show that the designed compounds are structurally diverse and synthetically feasible when compared to their respective objective drugs despite consistently high predicted behavioral similarity. Finally, we generated new designs that enhanced thirteen drugs/compounds associated with non-small cell lung cancer and anti-aging properties using their predicted proteomic interaction signatures. his study represents a significant step forward in automating holistic therapeutic design with machine learning, enabling the rapid generation of novel, effective, and safe drug leads for any indication.

Structural Insights on the SARS-CoV-2 Variants of Concern Spike Glycoprotein: A Computational Study With Possible Clinical Implications

Coronavirus disease 2019 (COVID-19) pandemic has been attributed to SARS-CoV-2 (SARS2) and, consequently, SARS2 has evolved into multiple SARS2 variants driving subsequent waves of infections. In particular, variants of concern (VOC) were identified to have both increased transmissibility and virulence ascribable to mutational changes occurring within the spike protein resulting to modifications in the protein structural orientation which in-turn may affect viral pathogenesis. However, this was never fully elucidated. Here, we generated spike models of endemic HCoVs (HCoV 229E, HCoV OC43, HCoV NL63, HCoV HKU1, SARS CoV, MERS CoV), original SARS2, and VOC (alpha, beta, gamma, delta). Model quality check, structural superimposition, and structural comparison based on RMSD values, TM scores, and contact mapping were all performed. We found that: 1) structural comparison between the original SARS2 and VOC whole spike protein model have minor structural differences (TM > 0.98); 2) the whole VOC spike models putatively have higher structural similarity (TM > 0.70) to spike models from endemic HCoVs coming from the same phylogenetic cluster; 3) original SARS2 S1-CTD and S1-NTD models are structurally comparable to VOC S1-CTD (TM = 1.0) and S1-NTD (TM > 0.96); and 4) endemic HCoV S1-CTD and S1-NTD models are structurally comparable to VOC S1-CTD (TM > 0.70) and S1-NTD (TM > 0.70) models belonging to the same phylogenetic cluster. Overall, we propose that structural similarities (possibly ascribable to similar conformational epitopes) may help determine immune cross-reactivity, whereas, structural differences (possibly associated with varying conformational epitopes) may lead to viral infection (either reinfection or breakthrough infection).

AlphaFold 2: Why It Works and Its Implications for Understanding the Relationships of Protein Sequence, Structure, and Function

AlphaFold 2 (AF2) was the star of CASP14, the last biannual structure prediction experiment. Using novel deep learning, AF2 predicted the structures of many difficult protein targets at or near experimental resolution. Here, we present our perspective of why AF2 works and show that it is a very sophisticated fold recognition algorithm that exploits the completeness of the library of single domain PDB structures. It has also learned local side chain packing rearrangements that enable it to refine proteins to high resolution. The benefits and limitations of its ability to predict the structures of many more proteins at or close to atomic detail are discussed.

A practical guide to large-scale docking

Structure-based docking screens of large compound libraries have become common in early drug and probe discovery. As computer efficiency has improved and compound libraries have grown, the ability to screen hundreds of millions, and even billions, of compounds has become feasible for modest-sized computer clusters. This allows the rapid and cost-effective exploration and categorization of vast chemical space into a subset enriched with potential hits for a given target. To accomplish this goal at speed, approximations are used that result in undersampling of possible configurations and inaccurate predictions of absolute binding energies. Accordingly, it is important to establish controls, as are common in other fields, to enhance the likelihood of success in spite of these challenges. Here we outline best practices and control docking calculations that help evaluate docking parameters for a given target prior to undertaking a large-scale prospective screen, with exemplification in one particular target, the melatonin receptor, where following this procedure led to direct docking hits with activities in the subnanomolar range. Additional controls are suggested to ensure specific activity for experimentally validated hit compounds. These guidelines should be useful regardless of the docking software used. Docking software described in the outlined protocol (DOCK3.7) is made freely available for academic research to explore new hits for a range of targets.

Possible therapeutic targets and promising drugs based on unsymmetrical hetaryl-substituted porphyrins to combat SARS-CoV-2

Coronavirus disease 2019 is a serious disease that causes acute respiratory syndrome and negatively affects the central nervous system. SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) crosses the blood-brain barrier due to the S protein on the surface of the viral particles. Thus, it is important to develop compounds that not only have an inhibitory effect but are also capable of completely deactivating the S-protein function. This study describes the purposeful modification of porphyrins and proposes compounds, asymmetrically hetaryl-substituted porphyrins with benzothiazole, benzoxazole, and N-methylbenzimidazole residues, to deactivate the S-protein functions. Molecular docking of SARS-CoV-2 proteins with hetaryl-substituted porphyrins showed that the viral proteins S-, N-, and nsp13 exhibited the highest binding affinity. Hetaryl-substituted porphyrins form strong complexes (13-14 kcal/mol) with the receptor-binding domain of the S protein, while the distance from the porphyrins to the receptor-binding motif (RBM) does not exceed 20 Å; therefore, RBM can be oxidized by ¹O₂, which is generated by porphyrin. Hetaryl-substituted porphyrins interact with the nucleocapsid protein in the serine/arginine-rich region, and a number of vulnerable amino acid residues are located in the photooxidation zone. This damage complicates the packaging of viral RNA into new virions. High-energy binding of hetaryl-substituted porphyrins with the N- and C-terminal domains of nsp13 was observed. This binding blocks the action of nsp13 as an enzyme of exoribonuclease and methyltransferase, thereby preventing RNA replication and processing. A procedure for the synthesis of hetaryl-substituted porphyrins was developed, new compounds were obtained, their structures were identified, and their photocatalytic properties were studied.

Deep Learning-Based Advances in Protein Structure Prediction

Obtaining an accurate description of protein structure is a fundamental step toward understanding the underpinning of biology. Although recent advances in experimental approaches have greatly enhanced our capabilities to experimentally determine protein structures, the gap between the number of protein sequences and known protein structures is ever increasing. Computational protein structure prediction is one of the ways to fill this gap. Recently, the protein structure prediction field has witnessed a lot of advances due to Deep Learning (DL)-based approaches as evidenced by the success of AlphaFold2 in the most recent Critical Assessment of protein Structure Prediction (CASP14). In this article, we highlight important milestones and progresses in the field of protein structure prediction due to DL-based methods as observed in CASP experiments. We describe advances in various steps of protein structure prediction pipeline viz. protein contact map prediction, protein distogram prediction, protein real-valued distance prediction, and Quality Assessment/refinement. We also highlight some end-to-end DL-based approaches for protein structure prediction approaches. Additionally, as there have been some recent DL-based advances in protein structure determination using Cryo-Electron (Cryo-EM) microscopy based, we also highlight some of the important progress in the field. Finally, we provide an outlook and possible future research directions for DL-based approaches in the protein structure prediction arena.

Prediction of LncRNA-encoded small peptides in glioma and oligomer channel functional analysis using in silico approaches

Glioma is a lethal malignant brain cancer, and many reports have shown that abnormalities in the behavior of water and ion channels play an important role in regulating tumor proliferation, migration, apoptosis, and differentiation. Recently, new studies have suggested that some long noncoding RNAs containing small open reading frames can encode small peptides and form oligomers for water or ion regulation. However, because the peptides are difficult to identify, their functional mechanisms are far from being clearly understood. In this study, we used bioinformatics methods to identify and evaluate lncRNAs, which may encode small transmembrane peptides in gliomas. Combining ab initio homology modeling, molecular dynamics simulations, and free energy calculations, we constructed a predictive model and predicted the oligomer channel activity of peptides by identifying the lncRNA ORFs. We found that one key hub lncRNA, namely, DLEU1, which contains two smORFs (ORF1 and ORF8), encodes small peptides that form pentameric channels. The mechanics of water and ion (Na⁺ and Cl^-) transport through this pentameric channel were simulated. The potential mean force of the H₂O molecules along the two ORF-encoded peptide channels indicated that the energy barrier was different between ORF1 and ORF8. The ORF1-encoded peptide pentamer acted as a self-assembled water channel but not as an ion channel, and the ORF8 permeated neither ions nor water. This work provides new methods and theoretical support for further elucidation of the function of lncRNA-encoded small peptides and their role in cancer. Additionally, this study provides a theoretical basis for drug development.

Structural Comparison of the SARS CoV 2 Spike Protein Relative to Other Human-Infecting Coronaviruses

Coronaviruses (CoV) are enveloped positive-stranded RNA viruses and, historically, there are seven known human-infecting CoVs with varying degrees of virulence. CoV attachment to the host is the first step of viral pathogenesis and mainly relies on the spike glycoprotein located on the viral surface. Among the human-infecting CoVs, only the infection of SARS CoV 2 (SARS2) among humans resulted to a pandemic which would suggest that the protein structural conformation of SARS2 spike protein is distinct as compared to other human-infecting CoVs. Surprisingly, the possible differences and similarities in the protein structural conformation between the various human-infecting CoV spike proteins have not been fully elucidated. In this study, we utilized a computational approach to generate models and analyze the seven human-infecting CoV spike proteins, namely: HCoV 229E, HCoV OC43, HCoV NL63, HCoV HKU1, SARS CoV, MERS CoV, and SARS2. Model quality assessment of all CoV models generated, structural superimposition of the whole protein model and selected S1 domains (S1-CTD and S1-NTD), and structural comparison based on RMSD values, Tm scores, and contact mapping were all performed. We found that the structural orientation of S1-CTD is a potential structural feature associated to both the CoV phylogenetic cluster and lineage. Moreover, we observed that spike models in the same phylogenetic cluster or lineage could potentially have similar protein structure. Additionally, we established that there are potentially three distinct S1-CTD orientation (Pattern I, Pattern II, Pattern III) among the human-infecting CoVs. Furthermore, we postulate that human-infecting CoVs in the same phylogenetic cluster may have similar S1-CTD and S1-NTD structural orientation. Taken together, we propose that the SARS2 spike S1-CTD follows a Pattern III orientation which has a higher degree of similarity with SARS1 and some degree of similarity with both OC43 and HKU1 which coincidentally are in the same phylogenetic cluster and lineage, whereas, the SARS2 spike S1-NTD has some degree of similarity among human-infecting CoVs that are either in the same phylogenetic cluster or lineage.

Improving homology modeling from low-sequence identity templates in Rosetta: A case study in GPCRs

As sequencing methodologies continue to advance, the availability of protein sequences far outpaces the ability of structure determination. Homology modeling is used to bridge this gap but relies on high-identity templates for accurate model building. G-protein coupled receptors (GPCRs) represent a significant target class for pharmaceutical therapies in which homology modeling could fill the knowledge gap for structure-based drug design. To date, only about 17% of druggable GPCRs have had their structures characterized at atomic resolution. However, modeling of the remaining 83% is hindered by the low sequence identity between receptors. Here we test key inputs in the model building process using GPCRs as a focus to improve the pipeline in two critical ways: Firstly, we use a blended sequence- and structure-based alignment that accounts for structure conservation in loop regions. Secondly, by merging multiple template structures into one comparative model, the best possible template for every region of a target can be used expanding the conformational space sampled in a meaningful way. This optimization allows for accurate modeling of receptors using templates as low as 20% sequence identity, which accounts for nearly the entire druggable space of GPCRs. A model database of all non-odorant GPCRs is made available at www.rosettagpcr.org. Additionally, all protocols are made available with insights into modifications that may improve accuracy at new targets.

Computational Analysis of IDH1, IDH2, and TP53 Mutations in Low-Grade Gliomas Including Oligodendrogliomas and Astrocytomas

Introduction:

The emergence of new omics approaches, such as genomic algorithms to identify tumor mutations and molecular modeling tools to predict the three-dimensional structure of proteins, has facilitated the understanding of the dynamic mechanisms involved in the pathogenesis of low-grade gliomas including oligodendrogliomas and astrocytomas.

Methods:

In this study, we targeted known mutations involved in low-grade gliomas, starting with the sequencing of genomic regions encompassing exon 4 of isocitrate dehydrogenase 1 (IDH1) and isocitrate dehydrogenase 2 (IDH2) and the four exons (5-6 and 7-8) of TP53 from 32 samples, followed by computational analysis to study the impact of these mutations on the structure and function of 3 proteins IDH1, IDH2, and p53.

Results:

We obtain a mutation that has an effect on the catalytic site of the protein IDH1 as R132H and on the catalytic site of the protein IDH2 as R172M. Other mutations at p53 have been identified as K305N, which is a pathogenic mutation; R175 H, which is a benign mutation; and R158G, which disrupts the structural conformation of the tumor suppressor protein.

Conclusion:

In low-grade gliomas, mutations in IDH1, IDH2, and TP53 may be the key to tumor progression because they have an effect on the function of the protein such as mutations R132H in IDH1 and R172M in IDH2, which change the function of the enzyme alpha-ketoglutarate, or R158G in TP53, which affects the structure of the generated protein, thus their importance in understanding gliomagenesis and for more accurate diagnosis complementary to the anatomical pathology tests.

Structural and functional divergence of GDAP1 from the glutathione S-transferase superfamily

Mutations in ganglioside-induced differentiation-associated protein 1 (GDAP1) alter mitochondrial morphology and result in several subtypes of the inherited peripheral neuropathy Charcot-Marie-Tooth disease; however, the mechanism by which GDAP1 functions has remained elusive. GDAP1 contains primary sequence homology to the GST superfamily; however, the question of whether GDAP1 is an active GST has not been clearly resolved. Here, we present biochemical evidence, suggesting that GDAP1 has lost the ability to bind glutathione without a loss of substrate binding activity. We have revealed that the α-loop, located within the H-site motif is the primary determinant for substrate binding. Using structural data of GDAP1, we have found that critical residues and configurations in the G-site which canonically interact with glutathione are altered in GDAP1, rendering it incapable of binding glutathione. Last, we have found that the overexpression of GDAP1 in HeLa cells results in a mitochondrial phenotype which is distinct from oxidative stress-induced mitochondrial fragmentation. This phenotype is dependent on the presence of the transmembrane domain, as well as a unique hydrophobic domain that is not found in canonical GSTs. Together, we data point toward a non-enzymatic role for GDAP1, such as a sensor or receptor.

Ranking near-native candidate protein structures via random forest classification

Background

In ab initio protein-structure predictions, a large set of structural decoys are often generated, with the requirement to select best five or three candidates from the decoys. The clustered central structures with the most number of neighbors are frequently regarded as the near-native protein structures with the lowest free energy; however, limitations in clustering methods and three-dimensional structural-distance assessments make identifying exact order of the best five or three near-native candidate structures difficult.

Results

To address this issue, we propose a method that re-ranks the candidate structures via random forest classification using intra- and inter-cluster features from the results of the clustering. Comparative analysis indicated that our method was better able to identify the order of the candidate structures as comparing with current methods SPICKR, Calibur, and Durandal. The results confirmed that the identification of the first model were closer to the native structure in 12 of 43 cases versus four for SPICKER, and the same as the native structure in up to 27 of 43 cases versus 14 for Calibur and up to eight of 43 cases versus two for Durandal.

Conclusions

In this study, we presented an improved method based on random forest classification to transform the problem of re-ranking the candidate structures by an binary classification. Our results indicate that this method is a powerful method for the problem and the effect of this method is better than other methods.

An overview of human leptospirosis vaccine design and future perspectives

Introduction: It's been 20 years since the first report of a recombinant vaccine that protected against leptospirosis. Since then, numerous recombinant vaccines have been evaluated; however, no recombinant vaccine candidate has advanced to clinical trials. With the ever-increasing burden of leptospirosis, there is an urgent need for a universal vaccine against leptospirosis.Areas covered: This review covers the most promising vaccine candidates that induced significant, reproducible, protection and how advances in the field of bioinformatics has led to the discovery of hundreds of novel protein targets. The authors also discuss the most recent findings regarding the innate immune response and host-pathogen interactions and their impact on the discovery of novel vaccine candidates. In addition, the authors have identified what they believe are the most challenging problems for the discovery and development of a universal vaccine and their potential solutions.Expert opinion: A universal vaccine for leptospirosis will likely only be achieved using a recombinant vaccine as the bacterins are of limited use due to the lack of a cross-protective immune response. Although there are hundreds of novel targets, due to the lack of immune correlates and the need for more research into the basic microbiology of Leptospira spp., a universal vaccine is 10-15 years away.

Comparative evaluation of structure and characteristic of peptidyl-prolyl cis-trans isomerase proteins and their function in Salmonella Typhimurium stress responses and virulence

Peptidyl-prolyl cis-trans isomerases (PPIase) exhibit chaperone activity and assist in protein folding by increasing the rate of cis-trans transition on proline-peptide bonds. The current study aimed to identify and characterize three genes, ppiA, ppiB, and ppiC, which encode proteins of the PPIase family in the bacterium Salmonella enterica serovar Typhimurium. Salmonella Typhimurium is a facultative intracellular zoonotic pathogen that causes food- and water-borne gastroenteritis in humans (leading to bacteremia in immune-compromised subjects). Recombinant clones for the three genes were constructed and sequenced and the sequences submitted to NCBI GenBank. Three-dimensional structures for the corresponding proteins were predicted by comparative modeling. A maximum-likelihood phylogenetic gene tree constructed for the three genes showed a low evolutionary mean diversity, indicating strong evolutionary conservation. Further, single-gene deletion mutant strains, generated for the respective genes, were observed to be more susceptible to the stationary phase of growth and heat stress conditions and showed reduced survival within macrophage cells line. The present study thus indicates that ppiA, ppiB, and ppiC genes are conserved among Salmonella genome, are critical for the growth of Salmonella Typhimurium in the examined stress conditions, and may play a role in its responses and virulence.

Simulations of CYP51A from Aspergillus fumigatus in a model bilayer provide insights into triazole drug resistance

Azole antifungal drugs target CYP51A in Aspergillus fumigatus by binding with the active site of the protein, blocking ergosterol biosynthesis. Resistance to azole antifungal drugs is now common, with a leucine to histidine amino acid substitution at position 98 the most frequent, predominantly conferring resistance to itraconazole, although cross-resistance has been reported in conjunction with other mutations. In this study, we create a homology model of CYP51A using a recently published crystal structure of the paralog protein CYP51B. The derived structures, wild type, and L98H mutant are positioned within a lipid membrane bilayer and subjected to molecular dynamics simulations in order improve the accuracy of both models. The structural analysis from our simulations suggests a decrease in active site surface from the formation of hydrogen bonds between the histidine substitution and neighboring polar side chains, potentially preventing the binding of azole drugs. This study yields a biologically relevant structure and set of dynamics of the A. fumigatus Lanosterol 14 alpha-demethylase enzyme and provides further insight into azole antifungal drug resistance.

Repository of Enriched Structures of Proteins Involved in the Red Blood Cell Environment (RESPIRE)

The Red Blood Cell (RBC) is a metabolically-driven cell vital for processes such a gas transport and homeostasis. RBC possesses at its surface exposing antigens proteins that are critical in blood transfusion. Due to their importance, numerous studies address the cell function as a whole but more and more details of RBC structure and protein content are now studied using massive state-of-the art characterisation techniques. Yet, the resulting information is frequently scattered in many scientific articles, in many databases and specialized web servers. To provide a more compendious view of erythrocytes and of their protein content, we developed a dedicated database called RESPIRE that aims at gathering a comprehensive and coherent ensemble of information and data about proteins in RBC. This cell-driven database lists proteins found in erythrocytes. For a given protein entry, initial data are processed from external portals and enriched by using state-of-the-art bioinformatics methods. As structural information is extremely useful to understand protein function and predict the impact of mutations, a strong effort has been put on the prediction of protein structures with a special treatment for membrane proteins. Browsing the database is available through text search for reference gene names or protein identifiers, through pre-defined queries or via hyperlinks. The RESPIRE database provides valuable information and unique annotations that should be useful to a wide audience of biologists, clinicians and structural biologists.

Database URL:http://www.dsimb.inserm.fr/respire

In silico prediction of prolactin molecules as a tool for equine genomics reproduction

The prolactin hormone is involved in several biological functions, although its main role resides on reproduction. As it interferes on fertility changes, studies focused on human health have established a linkage of this hormone to fertility losses. Regarding animal research, there is still a lack of information about the structure of prolactin. In case of horse breeding, prolactin has a particular influence; once there is an individualization of these animals and equines are known for presenting several reproductive disorders. As there is no molecular structure available for the prolactin hormone and receptor, we performed several bioinformatics analyses through prediction and refinement softwares, as well as manual modifications. Aiming to elucidate the first computational structure of both molecules and analyse structural and functional aspects related to these proteins, here we provide the first known equine model for prolactin and prolactin receptor, which obtained high global quality scores in diverse software's for quality assessment. QMEAN overall score obtained for ePrl was (- 4.09) and QMEANbrane for ePrlr was (- 8.45), which proves the structures' reliability. This study will implement another tool in equine genomics in order to give light to interactions of these molecules, structural and functional alterations and therefore help diagnosing fertility problems, contributing in the selection of a high genetic herd.

Modeling of Disordered Protein Structures Using Monte Carlo Simulations and Knowledge-Based Statistical Force Fields

The description of protein disordered states is important for understanding protein folding mechanisms and their functions. In this short review, we briefly describe a simulation approach to modeling protein interactions, which involve disordered peptide partners or intrinsically disordered protein regions, and unfolded states of globular proteins. It is based on the CABS coarse-grained protein model that uses a Monte Carlo (MC) sampling scheme and a knowledge-based statistical force field. We review several case studies showing that description of protein disordered states resulting from CABS simulations is consistent with experimental data. The case studies comprise investigations of protein⁻peptide binding and protein folding processes. The CABS model has been recently made available as the simulation engine of multiscale modeling tools enabling studies of protein⁻peptide docking and protein flexibility. Those tools offer customization of the modeling process, driving the conformational search using distance restraints, reconstruction of selected models to all-atom resolution, and simulation of large protein systems in a reasonable computational time. Therefore, CABS can be combined in integrative modeling pipelines incorporating experimental data and other modeling tools of various resolution.

Advances in Protein Super-Secondary Structure Prediction and Application to Protein Structure Prediction

Due to the advancement in various sequencing technologies, the gap between the number of protein sequences and the number of experimental protein structures is ever increasing. Community-wide initiatives like CASP have resulted in considerable efforts in the development of computational methods to accurately model protein structures from sequences. Sequence-based prediction of super-secondary structure has direct application in protein structure prediction, and there have been significant efforts in the prediction of super-secondary structure in the last decade. In this chapter, we first introduce the protein structure prediction problem and highlight some of the important progress in the field of protein structure prediction. Next, we discuss recent methods for the prediction of super-secondary structures. Finally, we discuss applications of super-secondary structure prediction in structure prediction/analysis of proteins. We also discuss prediction of protein structures that are composed of simple super-secondary structure repeats and protein structures that are composed of complex super-secondary structure repeats. Finally, we also discuss the recent trends in the field.

Experimental accuracy in protein structure refinement via molecular dynamics simulations

Refinement is the last step in protein structure prediction pipelines to convert approximate homology models to experimental accuracy. Protocols based on molecular dynamics (MD) simulations have shown promise, but current methods are limited to moderate levels of consistent refinement. To explore the energy landscape between homology models and native structures and analyze the challenges of MD-based refinement, eight test cases were studied via extensive simulations followed by Markov state modeling. In all cases, native states were found very close to the experimental structures and at the lowest free energies, but refinement was hindered by a rough energy landscape. Transitions from the homology model to the native states require the crossing of significant kinetic barriers on at least microsecond time scales. A significant energetic driving force toward the native state was lacking until its immediate vicinity, and there was significant sampling of off-pathway states competing for productive refinement. The role of recent force field improvements is discussed and transition paths are analyzed in detail to inform which key transitions have to be overcome to achieve successful refinement.

Recent advances in sequence-based protein structure prediction

The most accurate characterizations of the structure of proteins are provided by structural biology experiments. However, because of the high cost and labor-intensive nature of the structural experiments, the gap between the number of protein sequences and solved structures is widening rapidly. Development of computational methods to accurately model protein structures from sequences is becoming increasingly important to the biological community. In this article, we highlight some important progress in the field of protein structure prediction, especially those related to free modeling (FM) methods that generate structure models without using homologous templates. We also provide a short synopsis of some of the recent advances in FM approaches as demonstrated in the recent Computational Assessment of Structure Prediction competition as well as recent trends and outlook for FM approaches in protein structure prediction.

PREFMD: a web server for protein structure refinement via molecular dynamics simulations

Summary

Refinement of protein structure models is a long-standing problem in structural bioinformatics. Molecular dynamics-based methods have emerged as an avenue to achieve consistent refinement. The PREFMD web server implements an optimized protocol based on the method successfully tested in CASP11. Validation with recent CASP refinement targets shows consistent and more significant improvement in global structure accuracy over other state-of-the-art servers.

Availability and implementation

PREFMD is freely available as a web server at http://feiglab.org/prefmd. Scripts for running PREFMD as a stand-alone package are available at https://github.com/feiglab/prefmd.git.

Contact

feig@msu.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

Be positive: optimizing pramlintide from microcanonical analysis of amylin isoforms

Amylin, or human islet amyloid polypeptide (hIAPP), is a 37-residue hormone synergistic to insulin and co-secreted with it by β-cells in the pancreas. The deposition of its cytotoxic amyloid fibrils is strongly related to the progression of Type II diabetes (T2D) and islet graft failures. Notably, isoforms from some mammalian species, such as rats (rIAPP) and porcine (pIAPP), present a few key mutations preventing aggregation. This has lead to biotechnological development of drugs for adjunct therapies of T2D, such as pramlintide, a variant of hIAPP inspired by rIAPP whose proline substitutions have β-strand fibril-breaking properties. Ideally, such a drug should be formulated with insulin and co-administered, but this has been prevented by a poor solubility profile at the appropriate pH. Hopefully, this could be improved with appropriate point mutations, increasing the molecular net charge. Despite experimental progress, preliminary screening during rational drug design can greatly benefit from thermodynamic insight derived from molecular simulations. So we introduce microcanonical thermostatistics analysis of multicanonical (MUCA) simulations of wild-type amylin isoforms as a systematic assessment of protein thermostability. As a consequence of this comprehensive investigation, the most suitable single-point mutations able to optimize pramlintide are located among the wild-type amylin isoforms. In particular, we find that aggregation inhibition and increased solubility are inherited by pramlintide through further S20R substitution typical of pIAPP. Thus, we provide a consistent thermostatistical methodology to aid the design of improved adjunct therapies for T2D according to current clinical knowledge.

Identification of coding sequence and its use for functional and structural characterization of catalase from Phyllanthus emblica

Catalase is an essential antioxidant enzyme that is well characterized from microbial and animal sources. The structure of plant catalase is unknown. Therefore, it is of interest to understand the functional and structural characteristics of catalase from an Indian gooseberry, Phyllanthus emblica (or Emblica officinalis). Hence, catalase from P. emblica was cloned in pUC18 plasmid, sequenced and submitted to GenBank with the accession numbers "MF979112" and "ATO98311.1". InterProScan showed that the coding sequence is monofunctional and haem-dependent catalase-like superfamily. Multiple sequence alignment (MSA) followed by phylogenetic analysis showed that the P. emblica catalase groups with soybean catalase. We further report the characteristics of structural model of the enzyme for functional characterization.

High-resolution structure prediction of β-barrel membrane proteins

[Formula: see text]-Barrel membrane proteins ([Formula: see text]MPs) play important roles, but knowledge of their structures is limited. We have developed a method to predict their 3D structures. We predict strand registers and construct transmembrane (TM) domains of [Formula: see text]MPs accurately, including proteins for which no prediction has been attempted before. Our method also accurately predicts structures from protein families with a limited number of sequences and proteins with novel folds. An average main-chain rmsd of 3.48 Å is achieved between predicted and experimentally resolved structures of TM domains, which is a significant improvement ([Formula: see text]3 Å) over a recent study. For [Formula: see text]MPs with NMR structures, the deviation between predictions and experimentally solved structures is similar to the difference among the NMR structures, indicating excellent prediction accuracy. Moreover, we can now accurately model the extended [Formula: see text]-barrels and loops in non-TM domains, increasing the overall coverage of structure prediction by [Formula: see text]%. Our method is general and can be applied to genome-wide structural prediction of [Formula: see text]MPs.

Structural Analysis of Recombinant Human Preproinsulins by Structure Prediction, Molecular Dynamics, and Protein-Protein Docking

More effective production of human insulin is important, because insulin is the main medication that is used to treat multiple types of diabetes and because many people are suffering from diabetes. The current system of insulin production is based on recombinant DNA technology, and the expression vector is composed of a preproinsulin sequence that is a fused form of an artificial leader peptide and the native proinsulin. It has been reported that the sequence of the leader peptide affects the production of insulin. To analyze how the leader peptide affects the maturation of insulin structurally, we adapted several in silico simulations using 13 artificial proinsulin sequences. Three-dimensional structures of models were predicted and compared. Although their sequences had few differences, the predicted structures were somewhat different. The structures were refined by molecular dynamics simulation, and the energy of each model was estimated. Then, protein-protein docking between the models and trypsin was carried out to compare how efficiently the protease could access the cleavage sites of the proinsulin models. The results showed some concordance with experimental results that have been reported; so, we expect our analysis will be used to predict the optimized sequence of artificial proinsulin for more effective production.

Template-based and free modeling of I-TASSER and QUARK pipelines using predicted contact maps in CASP12

We develop two complementary pipelines, "Zhang-Server" and "QUARK", based on I-TASSER and QUARK pipelines for template-based modeling (TBM) and free modeling (FM), and test them in the CASP12 experiment. The combination of I-TASSER and QUARK successfully folds three medium-size FM targets that have more than 150 residues, even though the interplay between the two pipelines still awaits further optimization. Newly developed sequence-based contact prediction by NeBcon plays a critical role to enhance the quality of models, particularly for FM targets, by the new pipelines. The inclusion of NeBcon predicted contacts as restraints in the QUARK simulations results in an average TM-score of 0.41 for the best in top five predicted models, which is 37% higher than that by the QUARK simulations without contacts. In particular, there are seven targets that are converted from non-foldable to foldable (TM-score >0.5) due to the use of contact restraints in the simulations. Another additional feature in the current pipelines is the local structure quality prediction by ResQ, which provides a robust residue-level modeling error estimation. Despite the success, significant challenges still remain in ab initio modeling of multi-domain proteins and folding of β-proteins with complicated topologies bound by long-range strand-strand interactions. Improvements on domain boundary and long-range contact prediction, as well as optimal use of the predicted contacts and multiple threading alignments, are critical to address these issues seen in the CASP12 experiment.

SURPASS Low-Resolution Coarse-Grained Protein Modeling

Coarse-grained modeling of biomolecules has a very important role in molecular biology. In this work we present a novel SURPASS (Single United Residue per Pre-Averaged Secondary Structure fragment) model of proteins that can be an interesting alternative for existing coarse-grained models. The design of the model is unique and strongly supported by the statistical analysis of structural regularities characteristic for protein systems. Coarse-graining of protein chain structures assumes a single center of interactions per residue and accounts for preaveraged effects of four adjacent residue fragments. Knowledge-based statistical potentials encode complex interaction patterns of these fragments. Using the Replica Exchange Monte Carlo sampling scheme and a generic version of the SURPASS force field we performed test simulations of a representative set of single-domain globular proteins. The method samples a significant part of conformational space and reproduces protein structures, including native-like, with surprisingly good accuracy. Future extension of the SURPASS model on large biomacromolecular systems is briefly discussed.

An overview of comparative modelling and resources dedicated to large-scale modelling of genome sequences

Computational modelling of proteins has been a major catalyst in structural biology. Bioinformatics groups have exploited the repositories of known structures to predict high-quality structural models with high efficiency at low cost. This article provides an overview of comparative modelling, reviews recent developments and describes resources dedicated to large-scale comparative modelling of genome sequences. The value of subclustering protein domain superfamilies to guide the template-selection process is investigated. Some recent cases in which structural modelling has aided experimental work to determine very large macromolecular complexes are also cited.

Computational protein structure refinement: Almost there, yet still so far to go

Protein structures are essential in modern biology yet experimental methods are far from being able to catch up with the rapid increase in available genomic data. Computational protein structure prediction methods aim to fill the gap while the role of protein structure refinement is to take approximate initial template-based models and bring them closer to the true native structure. Current methods for computational structure refinement rely on molecular dynamics simulations, related sampling methods, or iterative structure optimization protocols. The best methods are able to achieve moderate degrees of refinement but consistent refinement that can reach near-experimental accuracy remains elusive. Key issues revolve around the accuracy of the energy function, the inability to reliably rank multiple models, and the use of restraints that keep sampling close to the native state but also limit the degree of possible refinement. A different aspect is the question of what exactly the target of high-resolution refinement should be as experimental structures are affected by experimental conditions and different biological questions require varying levels of accuracy. While improvement of the global protein structure is a difficult problem, high-resolution refinement methods that improves local structural quality such as favorable stereochemistry and the avoidance of atomic clashes are much more successful.

Identify High-Quality Protein Structural Models by Enhanced K-Means

Background. One critical issue in protein three-dimensional structure prediction using either ab initio or comparative modeling involves identification of high-quality protein structural models from generated decoys. Currently, clustering algorithms are widely used to identify near-native models; however, their performance is dependent upon different conformational decoys, and, for some algorithms, the accuracy declines when the decoy population increases. Results. Here, we proposed two enhanced K-means clustering algorithms capable of robustly identifying high-quality protein structural models. The first one employs the clustering algorithm SPICKER to determine the initial centroids for basic K-means clustering (SK-means), whereas the other employs squared distance to optimize the initial centroids (K-means++). Our results showed that SK-means and K-means++ were more robust as compared with SPICKER alone, detecting 33 (59%) and 42 (75%) of 56 targets, respectively, with template modeling scores better than or equal to those of SPICKER. Conclusions. We observed that the classic K-means algorithm showed a similar performance to that of SPICKER, which is a widely used algorithm for protein-structure identification. Both SK-means and K-means++ demonstrated substantial improvements relative to results from SPICKER and classical K-means.

Why Is There a Glass Ceiling for Threading Based Protein Structure Prediction Methods?

Despite their different implementations, comparison of the best threading approaches to the prediction of evolutionary distant protein structures reveals that they tend to succeed or fail on the same protein targets. This is true despite the fact that the structural template library has good templates for all cases. Thus, a key question is why are certain protein structures threadable while others are not. Comparison with threading results on a set of artificial sequences selected for stability further argues that the failure of threading is due to the nature of the protein structures themselves. Using a new contact map based alignment algorithm, we demonstrate that certain folds are highly degenerate in that they can have very similar coarse grained fractions of native contacts aligned and yet differ significantly from the native structure. For threadable proteins, this is not the case. Thus, contemporary threading approaches appear to have reached a plateau, and new approaches to structure prediction are required.

Mapping the Geometric Evolution of Protein Folding Motor

Polypeptide chain has an invariant main-chain and a variant side-chain sequence. How the side-chain sequence determines fold in terms of its chemical constitution has been scrutinized extensively and verified periodically. However, a focussed investigation on the directive effect of side-chain geometry may provide important insights supplementing existing algorithms in mapping the geometrical evolution of protein chains and its structural preferences. Geometrically, folding of protein structure may be envisaged as the evolution of its geometric variables: ϕ, and ψ dihedral angles of polypeptide main-chain directed by χ1, and χ2 of side chain. In this work, protein molecule is metaphorically modelled as a machine with 4 rotors ϕ, ψ, χ1 and χ2, with its evolution to the functional fold is directed by combinations of its rotor directions. We observe that differential rotor motions lead to different secondary structure formations and the combinatorial pattern is unique and consistent for particular secondary structure type. Further, we found that combination of rotor geometries of each amino acid is unique which partly explains how different amino acid sequence combinations have unique structural evolution and functional adaptation. Quantification of these amino acid rotor preferences, resulted in the generation of 3 substitution matrices, which later on plugged in the BLAST tool, for evaluating their efficiency in aligning sequences. We have employed BLOSUM62 and PAM30 as standard for primary evaluation. Generation of substitution matrices is a logical extension of the conceptual framework we attempted to build during the development of this work. Optimization of matrices following the conventional routines and possible application with biologically relevant data sets are beyond the scope of this manuscript, though it is a part of the larger project design.

Performance of protein-structure predictions with the physics-based UNRES force field in CASP11

Participating as the Cornell-Gdansk group, we have used our physics-based coarse-grained UNited RESidue (UNRES) force field to predict protein structure in the 11th Community Wide Experiment on the Critical Assessment of Techniques for Protein Structure Prediction (CASP11). Our methodology involved extensive multiplexed replica exchange simulations of the target proteins with a recently improved UNRES force field to provide better reproductions of the local structures of polypeptide chains. All simulations were started from fully extended polypeptide chains, and no external information was included in the simulation process except for weak restraints on secondary structure to enable us to finish each prediction within the allowed 3-week time window. Because of simplified UNRES representation of polypeptide chains, use of enhanced sampling methods, code optimization and parallelization and sufficient computational resources, we were able to treat, for the first time, all 55 human prediction targets with sizes from 44 to 595 amino acid residues, the average size being 251 residues. Complete structures of six single-domain proteins were predicted accurately, with the highest accuracy being attained for the T0769, for which the CαRMSD was 3.8 Å for 97 residues of the experimental structure. Correct structures were also predicted for 13 domains of multi-domain proteins with accuracy comparable to that of the best template-based modeling methods. With further improvements of the UNRES force field that are now underway, our physics-based coarse-grained approach to protein-structure prediction will eventually reach global prediction capacity and, consequently, reliability in simulating protein structure and dynamics that are important in biochemical processes.

AVAILABILITY AND IMPLEMENTATION

Freely available on the web at http://www.unres.pl/ CONTACT: has5@cornell.edu.

ORION: a web server for protein fold recognition and structure prediction using evolutionary hybrid profiles

Protein structure prediction based on comparative modeling is the most efficient way to produce structural models when it can be performed. ORION is a dedicated webserver based on a new strategy that performs this task. The identification by ORION of suitable templates is performed using an original profile-profile approach that combines sequence and structure evolution information. Structure evolution information is encoded into profiles using structural features, such as solvent accessibility and local conformation -with Protein Blocks-, which give an accurate description of the local protein structure. ORION has recently been improved, increasing by 5% the quality of its results. The ORION web server accepts a single protein sequence as input and searches homologous protein structures within minutes. Various databases such as PDB, SCOP and HOMSTRAD can be mined to find an appropriate structural template. For the modeling step, a protein 3D structure can be directly obtained from the selected template by MODELLER and displayed with global and local quality model estimation measures. The sequence and the predicted structure of 4 examples from the CAMEO server and a recent CASP11 target from the 'Hard' category (T0818-D1) are shown as pertinent examples. Our web server is accessible at http://www.dsimb.inserm.fr/ORION/.

Identification of the High-affinity Substrate-binding Site of the Multidrug and Toxic Compound Extrusion (MATE) Family Transporter from Pseudomonas stutzeri

Multidrug and toxic compound extrusion (MATE) transporters exist in all three domains of life. They confer multidrug resistance by utilizing H(+) or Na(+) electrochemical gradients to extrude various drugs across the cell membranes. The substrate binding and the transport mechanism of MATE transporters is a fundamental process but so far not fully understood. Here we report a detailed substrate binding study of NorM_PS, a representative MATE transporter from Pseudomonas stutzeri Our results indicate that NorM_PS is a proton-dependent multidrug efflux transporter. Detailed binding studies between NorM_PS and 4',6-diamidino-2-phenylindole (DAPI) were performed by isothermal titration calorimetry (ITC), differential scanning calorimetry (DSC), and spectrofluorometry. Two exothermic binding events were observed from ITC data, and the high-affinity event was directly correlated with the extrusion of DAPI. The affinities are about 1 μm and 0.1 mm for the high and low affinity binding, respectively. Based on our homology model of NorM_PS, variants with mutations of amino acids that are potentially involved in substrate binding, were constructed. By carrying out the functional characterization of these variants, the critical amino acid residues (Glu-257 and Asp-373) for high-affinity DAPI binding were determined. Taken together, our results suggest a new substrate-binding site for MATE transporters.

CASP11--An Evaluation of a Modular BCL::Fold-Based Protein Structure Prediction Pipeline

In silico prediction of a protein's tertiary structure remains an unsolved problem. The community-wide Critical Assessment of Protein Structure Prediction (CASP) experiment provides a double-blind study to evaluate improvements in protein structure prediction algorithms. We developed a protein structure prediction pipeline employing a three-stage approach, consisting of low-resolution topology search, high-resolution refinement, and molecular dynamics simulation to predict the tertiary structure of proteins from the primary structure alone or including distance restraints either from predicted residue-residue contacts, nuclear magnetic resonance (NMR) nuclear overhauser effect (NOE) experiments, or mass spectroscopy (MS) cross-linking (XL) data. The protein structure prediction pipeline was evaluated in the CASP11 experiment on twenty regular protein targets as well as thirty-three 'assisted' protein targets, which also had distance restraints available. Although the low-resolution topology search module was able to sample models with a global distance test total score (GDT_TS) value greater than 30% for twelve out of twenty proteins, frequently it was not possible to select the most accurate models for refinement, resulting in a general decay of model quality over the course of the prediction pipeline. In this study, we provide a detailed overall analysis, study one target protein in more detail as it travels through the protein structure prediction pipeline, and evaluate the impact of limited experimental data.