Evaluation of predictions in the CASP10 model refinement category

Design, synthesis, and bioevaluation of novel unsaturated cyanoacetamide derivatives: In vitro and in silico exploration

In this study, we synthesized novel α,β-unsaturated 2-cyanoacetamide derivatives (1-5) using microwave-assisted Knoevenagel condensation. Characterization of these compounds was carried out using FTIR and 1H NMR spectroscopy. We then evaluated their in vitro antibacterial activity against both gram-positive and gram-negative pathogenic bacteria. Additionally, we employed in silico methods, including ADMET prediction and density functional theory (DFT) calculations of molecular orbital properties, to investigate these cyanoacetamide derivatives (1-5). Molecular docking was used to assess the binding interactions of these derivatives (1-5) with seven target proteins (5MM8, 4NZZ, 7FEQ, 5NIJ, ITM2, 6SE1, and 5GVZ) and compared them to the reference standard tyrphostin AG99. Notably, derivative 5 exhibited the most favorable binding affinity, with a binding energy of -7.7 kcal mol-1 when interacting with the staphylococcus aureus (PDB:5MM8), while also meeting all drug-likeness criteria. Additionally, molecular dynamics simulations were carried out to evaluate the stability of the interaction between the protein and ligand, utilizing parameters such as Root-Mean-Square Deviation (RMSD), Root-Mean-Square Fluctuation (RMSF), Radius of Gyration (Rg), and Principal Component Analysis (PCA). A 50 nanosecond molecular dynamics (MD) simulation was performed to investigate stability further, incorporating RMSD and RMSF analyses on compound 5 within the active binding site of the modeled protein across different temperatures (300, 305, 310, and 320 K). Among these temperatures, compound 5 exhibited an RMSD value ranging from approximately 0.2 to 0.3 nm at 310 K (body temperature) with the 5MM8 target, which differed from the other temperature conditions. The in silico results suggest that compound 5 maintained significant conformational stability throughout the 50 ns simulation period. It is consistent with its low docking energy and in vitro findings concerning α,β-unsaturated cyanoacetamides. Key insights from this study include:•The creation of innovative α,β-unsaturated 2-cyanoacetamide derivatives (1-5) employing cost-effective, licensed, versatile, and efficient software for both in silico and in vitro assessment of antibacterial activity.•Utilization of FTIR and NMR techniques for characterizing compounds 1-5.

Immunoinformatic Identification of Multiple Epitopes of gp120 Protein of HIV-1 to Enhance the Immune Response against HIV-1 Infection

Acquired Immunodeficiency Syndrome is caused by the Human Immunodeficiency Virus (HIV), and a significant number of fatalities occur annually. There is a dire need to develop an effective vaccine against HIV-1. Understanding the structural proteins of viruses helps in designing a vaccine based on immunogenic peptides. In the current experiment, we identified gp120 epitopes using bioinformatic epitope prediction tools, molecular docking, and MD simulations. The Gb-1 peptide was considered an adjuvant. Consecutive sequences of GTG, GSG, GGTGG, and GGGGS linkers were used to bind the B cell, Cytotoxic T Lymphocytes (CTL), and Helper T Lymphocytes (HTL) epitopes. The final vaccine construct consisted of 315 amino acids and is expected to be a recombinant protein of approximately 35.49 kDa. Based on docking experiments, molecular dynamics simulations, and tertiary structure validation, the analysis of the modeled protein indicates that it possesses a stable structure and can interact with Toll-like receptors. The analysis demonstrates that the proposed vaccine can provoke an immunological response by activating T and B cells, as well as stimulating the release of IgA and IgG antibodies. This vaccine shows potential for HIV-1 prophylaxis. The in-silico design suggests that multiple-epitope constructs can be used as potentially effective immunogens for HIV-1 vaccine development.

Leveraging immunoinformatics for developing a multi-epitope subunit vaccine against Helicobacter pylori and Fusobacterium nucleatum

Gastric ulcers caused by Helicobacter pylori and Fusobacterium nucleatum remain a significant global health concern without an established vaccine. In this study, we utilized immunoinformatics methods to design a multi-epitope vaccine targeting these pathogens. Outer membrane proteins from H. pylori and F. nucleatum were scrutinized to identify high antigenic T-cell and B-cell epitopes. The resulting vaccine comprised carefully analyzed and evaluated epitopes, including cytotoxic T-lymphocytes, helper T-lymphocytes, and linear B-lymphocytes epitopes. This vaccine exhibited notable antigenicity, suitable immunogenicity, and demonstrated non-allergenicity and non-toxicity. It displayed favorable physiochemical characteristics and high solubility. In interaction studies, the vaccine exhibited robust binding to toll-like receptor 4 (TLR4). Molecular dynamic simulations revealed cohesive structural integrity and stable attachment. Codon adaptation utilizing Escherichia coli K12 host yielded a vaccine with elevated Codon Adaptation Index (CAI) and optimal GC content. In silico cloning into the pET28+(a) vector demonstrated efficient expression. Immune simulations indicated the vaccine's ability to initiate immune responses in humans, mirroring real-life scenarios. Based on these comprehensive findings, we propose that our developed vaccine has the potential to confer robust immunity against H. pylori and F. nucleatum infections.Communicated by Ramaswamy H. Sarma.

A computational approach to design a polyvalent vaccine against human respiratory syncytial virus

Human Respiratory Syncytial Virus (RSV) is one of the leading causes of lower respiratory tract infections (LRTI), responsible for infecting people from all age groups-a majority of which comprises infants and children. Primarily, severe RSV infections are accountable for multitudes of deaths worldwide, predominantly of children, every year. Despite several efforts to develop a vaccine against RSV as a potential countermeasure, there has been no approved or licensed vaccine available yet, to control the RSV infection effectively. Therefore, through the utilization of immunoinformatics tools, a computational approach was taken in this study, to design a multi-epitope polyvalent vaccine against two major antigenic subtypes of RSV, RSV-A and RSV-B. Potential predictions of the T-cell and B-cell epitopes were followed by extensive tests of antigenicity, allergenicity, toxicity, conservancy, homology to human proteome, transmembrane topology, and cytokine-inducing ability. The peptide vaccine was modeled, refined, and validated. Molecular docking analysis with specific Toll-like receptors (TLRs) revealed excellent interactions with suitable global binding energies. Additionally, molecular dynamics (MD) simulation ensured the stability of the docking interactions between the vaccine and TLRs. Mechanistic approaches to imitate and predict the potential immune response generated by the administration of vaccines were determined through immune simulations. Subsequent mass production of the vaccine peptide was evaluated; however, there remains a necessity for further in vitro and in vivo experiments to validate its efficacy against RSV infections.

Computational designing of a novel subunit vaccine for human cytomegalovirus by employing the immunoinformatics framework

Human cytomegalovirus (HCMV) is a widespread virus that can cause serious and irreversible neurological damage in newborns and even death in children who do not have the access to much-needed medications. While some vaccines and drugs are found to be effective against HCMV, their extended use has given rise to dose-limiting toxicities and the development of drug-resistant mutants among patients. Despite half a century's worth of research, the lack of a licensed HCMV vaccine heightens the need to develop newer antiviral therapies and vaccine candidates with improved effectiveness and reduced side effects. In this study, the immunoinformatics approach was utilized to design a potential polyvalent epitope-based vaccine effective against the four virulent strains of HCMV. The vaccine was constructed using seven CD8+ cytotoxic T lymphocytes epitopes, nine CD4+ helper T lymphocyte epitopes, and twelve linear B-cell lymphocyte epitopes that were predicted to be antigenic, non-allergenic, non-toxic, fully conserved, and non-human homologous. Subsequently, molecular docking study, protein-protein interaction analysis, molecular dynamics simulation (including the root mean square fluctuation (RMSF) and root mean square deviation (RMSD)), and immune simulation study rendered promising results assuring the vaccine to be stable, safe, and effective. Finally, in silico cloning was conducted to develop an efficient mass production strategy of the vaccine. However, further in vitro and in vivo research studies on the proposed vaccine are required to confirm its safety and efficacy.Communicated by Ramaswamy H. Sarma.

Using Local Protein Model Quality Estimates to Guide a Molecular Dynamics-Based Refinement Strategy

The refinement of predicted 3D models aims to bring them closer to the native structure by fixing errors including unusual bonds and torsion angles and irregular hydrogen bonding patterns. Refinement approaches based on molecular dynamics (MD) simulations using different types of restraints have performed well since CASP10. ReFOLD, developed by the McGuffin group, was one of the many MD-based refinement approaches, which were tested in CASP 12. When the performance of the ReFOLD method in CASP12 was evaluated, it was observed that ReFOLD suffered from the absence of a reliable guidance mechanism to reach consistent improvement for the quality of predicted 3D models, particularly in the case of template-based modelling (TBM) targets. Therefore, here we propose to utilize the local quality assessment score produced by ModFOLD6 to guide the MD-based refinement approach to further increase the accuracy of the predicted 3D models. The relative performance of the new local quality assessment guided MD-based refinement protocol and the original MD-based protocol ReFOLD are compared utilizing many different official scoring methods. By using the per-residue accuracy (or local quality) score to guide the refinement process, we are able to prevent the refined models from undesired structural deviations, thereby leading to more consistent improvements. This chapter will include a detailed analysis of the performance of the local quality assessment guided MD-based protocol versus that deployed in the original ReFOLD method.

Immunoinformatics-guided designing and in silico analysis of epitope-based polyvalent vaccines against multiple strains of human coronavirus (HCoV)

OBJECTIVES

The group of human coronaviruses (HCoVs) consists of some highly pathogenic viruses that have caused several outbreaks in the past. The newly emerged strain of HCoV, the SARS-CoV-2 is responsible for the recent global pandemic that has already caused the death of hundreds of thousands of people due to the lack of effective therapeutic options.

METHODS

In this study, immunoinformatics methods were used to design epitope-based polyvalent vaccines which are expected to be effective against four different pathogenic strains of HCoV i.e., HCoV-OC43, HCoV-SARS, HCoV-MERS, and SARS-CoV-2.

RESULTS

The constructed vaccines consist of highly antigenic, non-allergenic, nontoxic, conserved, and non-homologous T-cell and B-cell epitopes from all the four viral strains. Therefore, they should be able to provide strong protection against all these strains. Protein-protein docking was performed to predict the best vaccine construct. Later, the MD simulation and immune simulation of the best vaccine construct also predicted satisfactory results. Finally, in silico cloning was performed to develop a mass production strategy of the vaccine.

CONCLUSION

If satisfactory results are achieved in further in vivo and in vitro studies, then the vaccines designed in this study might be effective as preventative measures against the selected HCoV strains.

Designing a vaccine-based therapy against Epstein-Barr virus-associated tumors using immunoinformatics approach

Epstein-Barr virus (EBV) is widely known due to its role in the etiology of infectious mononucleosis. However, it is the first oncovirus that was identified and has been implicated in the etiology of several types of cancers. Globally, EBV infection is associated with more than 200, 000 new cancer cases and 150, 000 deaths yearly. A prophylactic or therapeutic vaccine targeting tumors associated with EBV infection is currently lacking. Therefore, this study aimed to develop a multiepitope-based polyvalent vaccine against EBV-associated tumors using immunoinformatics approach. The latency-associated proteins (LAP) of three strains of the virus were used in this study. Potential epitopes predicted from the proteins were analyzed and selected based on several predicted properties. Thirty viable B-cell and T-cell epitopes were selected and conjugated using various linkers alongside beta-defensin 3 as an adjuvant and pan HLA DR-binding epitope (PADRE) sequence to improve the immunogenicity of the vaccine construct. Molecular docking studies of the vaccine construct against toll-like receptors (TLRs) showed it is capable of inducing immune response via recognition by TLRs while immune simulation studies showed it could induce both cellular and humoral immune responses. Furthermore, molecular dynamics study of the complex formed by the vaccine candidate and TLR-4 showed that the complex was stable. Ultimately, the designed vaccine showed desirable properties based on in silico evaluation; however, experimental studies are needed to validate the efficacy of the vaccine against EBV-associated tumors.

Immunoinformatics Approach to Design Novel Subunit Vaccine against the Epstein-Barr Virus

Epstein-Barr virus (EBV) is a lymphotropic virus responsible for numerous epithelial and lymphoid cell malignancies, including gastric carcinoma, Hodgkin's lymphoma, nasopharyngeal carcinoma, and Burkitt lymphoma. Hundreds of thousands of people worldwide get infected with this virus, and in most cases, this viral infection leads to cancer. Although researchers are trying to develop potential vaccines and drug therapeutics, there is still no effective vaccine to combat this virus. In this study, the immunoinformatics approach was utilized to develop a potential multiepitope subunit vaccine against the two most common subtypes of EBV, targeting three of their virulent envelope glycoproteins. Eleven cytotoxic T lymphocyte (CTL) epitopes, 11 helper T lymphocyte (HTL) epitopes, and 10 B-cell lymphocyte (BCL) epitopes were predicted to be antigenic, nonallergenic, nontoxic, and fully conserved among the two subtypes, and nonhuman homologs were used for constructing the vaccine after much analysis. Later, further validation experiments, including molecular docking with different immune receptors (e.g., Toll-like receptors [TLRs]), molecular dynamics simulation analyses (including root means square deviation [RMSD], root mean square fluctuation [RMSF], radius of gyration [Rg], principal-component analysis [PCA], dynamic cross-correlation [DCC], definition of the secondary structure of proteins [DSSP], and Molecular Mechanics Poisson-Boltzmann Surface Area [MM-PBSA]), and immune simulation analyses generated promising results, ensuring the safe and stable response of the vaccine with specific immune receptors after potential administration within the human body. The vaccine's high binding affinity with TLRs was revealed in the docking study, and a very stable interaction throughout the simulation proved the potential high efficacy of the proposed vaccine. Further, in silico cloning was also conducted to design an efficient mass production strategy for future bulk industrial vaccine production. IMPORTANCE Epstein-Barr virus (EBV) vaccines have been developing for over 30 years, but polyphyletic and therapeutic vaccines have failed to get licensed. Our vaccine surpasses the limitations of many such vaccines and remains very promising, which is crucial because the infection rate is higher than most viral infections, affecting a whopping 90% of the adult population. One of the major identifications covers a holistic analysis of populations worldwide, giving us crucial information about its effectiveness for everyone's unique immunological system. We targeted three glycoproteins that enhance the virulence of the virus to design an epitope-based polyvalent vaccine against two different strains of EBV, type 1 and 2. Our methodology in this study is nonconventional yet swift to show effective results while designing vaccines.

Core Proteomics and Immunoinformatic Approaches to Design a Multiepitope Reverse Vaccine Candidate against Chagas Disease

Chagas disease is a tropical ailment indigenous to South America and caused by the protozoan parasite Trypanosoma cruzi, which has serious health consequences globally. Insect vectors transmit the parasite and, due to the lack of vaccine availability and limited treatment options, we implemented an integrated core proteomics analysis to design a reverse vaccine candidate based on immune epitopes for disease control. Firstly, T. cruzi core proteomics was used to identify immunodominant epitopes. Therefore, we designed the vaccine sequence to be non-allergic, antigenic, immunogenic, and to have better solubility. After predicting the tertiary structure, docking and molecular dynamics simulation (MDS) were performed with TLR4, MHC-I, and MHC-II receptors to discover the binding affinities. The final vaccine design demonstrated significant hydrogen bond interactions upon docking with TLR4, MHC-I, and MHC-II receptors. This indicated the efficacy of the vaccine candidate. A server-based immune simulation approach was generated to predict the efficacy. Significant structural compactness and binding stability were found based on MDS. Finally, by optimizing codons on Escherichia coli K12, a high GC content and CAI value were obtained, which were then incorporated into the cloning vector pET2+ (a). Thus, the developed vaccine sequence may be a viable therapy option for Chagas disease.

Application of reverse vaccinology to design a multi-epitope subunit vaccine against a new strain of Aeromonas veronii

Background

Aeromonas veronii is one of the most common pathogens of freshwater fishes that cause sepsis and ulcers. There are increasing numbers of cases showing that it is a significant zoonotic and aquatic agent. Epidemiological studies have shown that A. veronii virulence and drug tolerance have both increased over the last few years as a result of epidemiological investigations. Cadaverine reverse transporter (CadB) and maltoporin (LamB protein) contribute to the virulence of A. veronii TH0426. TH0426 strain is currently showing severe cases on fish species, and its resistance against therapeutic has been increasing. Despite these devastating complications, there is still no effective cure or vaccine for this strain of A.veronii.

Results

In this regard, an immunoinformatic method was used to generate an epitope-based vaccine against this pathogen. The immunodominant epitopes were identified using the CadB and LamB protein of A. veronii. The final constructed vaccine sequence was developed to be immunogenic, non-allergenic as well as have better solubility. Molecular dynamic simulation revealed significant binding stability and structural compactness. Finally, using Escherichia coli K12 as a model, codon optimization yielded ideal GC content and a higher CAI value, which was then included in the cloning vector pET2+ (a).

Conclusion

Altogether, our outcomes imply that the proposed peptide vaccine might be a good option for A. veronii TH0426 prophylaxis.

Immunoinformatic Design of a Multivalent Peptide Vaccine Against Mucormycosis: Targeting FTR1 Protein of Major Causative Fungi

Mucormycosis is a potentially fatal illness that arises in immunocompromised people due to diabetic ketoacidosis, neutropenia, organ transplantation, and elevated serum levels of accessible iron. The sudden spread of mucormycosis in COVID-19 patients engendered massive concern worldwide. Comorbidities including diabetes, cancer, steroid-based medications, long-term ventilation, and increased ferritin serum concentration in COVID-19 patients trigger favorable fungi growth that in turn effectuate mucormycosis. The necessity of FTR1 gene-encoded ferrous permease for host iron acquisition by fungi has been found in different studies recently. Thus, targeting the transit component could be a potential solution. Unfortunately, no appropriate antifungal vaccine has been constructed as of yet. To date, mucormycosis has been treated with antiviral therapy and surgical treatment only. Thus, in this study, the FTR1 protein has been targeted to design a convenient and novel epitope-based vaccine with the help of immunoinformatics against four different virulent fungal species. Furthermore, the vaccine was constructed using 8 CTL, 2 HTL, and 1 LBL epitopes that were found to be highly antigenic, non-allergenic, non-toxic, and fully conserved among the fungi under consideration. The vaccine has very reassuring stability due to its high pI value of 9.97, conclusive of a basic range. The vaccine was then subjected to molecular docking, molecular dynamics, and immune simulation studies to confirm the biological environment’s safety, efficacy, and stability. The vaccine constructs were found to be safe in addition to being effective. Finally, we used in-silico cloning to develop an effective strategy for vaccine mass production. The designed vaccine will be a potential therapeutic not only to control mucormycosis in COVID-19 patients but also be effective in general mucormycosis events. However, further in vitro, and in vivo testing is needed to confirm the vaccine’s safety and efficacy in controlling fungal infections. If successful, this vaccine could provide a low-cost and effective method of preventing the spread of mucormycosis worldwide.

Use of Integrated Core Proteomics, Immuno-Informatics, and In Silico Approaches to Design a Multiepitope Vaccine against Zoonotic Pathogen Edwardsiella tarda

Multidrug-resistant Edwardsiella tarda has been reported as the main causative agent for massive fish mortality. The pathogen is well-known for causing hemorrhagic septicemia in fish and has been linked to gastrointestinal infections in humans. Formalin-inactivated Edwardsiella vaccination has previously been found to be ineffective in aquaculture species. Therefore, based on E. tarda’s integrated core complete sequenced genomes, the study aimed to design a subunit vaccine based on T and B cell epitopes employing immunoinformatics approach. Initially, the top immunodominant and antigenic epitopes were predicted from the core complete sequenced genomes of the E. tarda genome and designed the vaccine by using linkers and adjuvant. In addition, vaccine 3D structure was predicted followed by refinement, and molecular docking was performed for the analysis of interacting residues between vaccines with TLR5, MHC-I, and MHC-II, respectively. The final vaccine constructs demonstrated strong hydrogen bond interactions. Molecular dynamic simulation of vaccine-TLR5 receptor complex showed a stable structural binding and compactness. Furthermore, E. coli used as a model organism for codon optimization proved optimal GC content and CAI value, which were subsequently cloned in vector pET2+ (a). Overall, the findings of the study imply that the designed epitope vaccine might be a good option for prophylaxis for E. tarda.

Designing a novel mRNA vaccine against Vibrio harveyi infection in fish: an immunoinformatics approach

Vibrio harveyi belongs to the family Vibrionaceae of class Gammaproteobacteria. Around 12 Vibrio species can cause gastroenteritis (gastrointestinal illness) in humans. A large number of bacterial particles can be found in the infected cells, which may cause death. Despite these devastating complications, there is still no cure or vaccine for the bacteria. As a result, we used an immunoinformatics approach to develop a multi-epitope vaccine against the most pathogenic hemolysin gene of V. harveyi. The immunodominant T- and B-cell epitopes were identified using the hemolysin protein. We developed a vaccine employing three possible epitopes: cytotoxic T-lymphocytes, helper T-lymphocytes, and linear B-lymphocyte epitopes, after thorough testing. The vaccine was developed to be antigenic, immunogenic, and non-allergenic, as well as have a better solubility. Molecular dynamics simulation revealed significant structural stiffness and binding stability. In addition, the immunological simulation generated by computers revealed that the vaccination might elicit immune reactions Escherichia coli K12 as a model, codon optimization yielded ideal GC content and a higher codon adaptation index value, which was then included in the cloning vector pET2+ (a). Altogether, our experiment implies that the proposed peptide vaccine might be a good option for vibriosis prophylaxis.

Structure‐conditioned amino‐acid couplings: how contact geometry affects pairwise sequence preferences

Relating a protein's sequence to its conformation is a central challenge for both structure prediction and sequence design. Statistical contact potentials, as well as their more descriptive versions that account for side‐chain orientation and other geometric descriptors, have served as simplistic but useful means of representing second‐order contributions in sequence–structure relationships. Here we ask what happens when a pairwise potential is conditioned on the fully defined geometry of interacting backbones fragments. We show that the resulting structure‐conditioned coupling energies more accurately reflect pair preferences as a function of structural contexts. These structure‐conditioned energies more reliably encode native sequence information and more highly correlate with experimentally determined coupling energies. Clustering a database of interaction motifs by structure results in ensembles of similar energies and clustering them by energy results in ensembles of similar structures. By comparing many pairs of interaction motifs and showing that structural similarity and energetic similarity go hand‐in‐hand, we provide a tangible link between modular sequence and structure elements. This link is applicable to structural modeling, and we show that scoring CASP models with structured‐conditioned energies results in substantially higher correlation with structural quality than scoring the same models with a contact potential. We conclude that structure‐conditioned coupling energies are a good way to model the impact of interaction geometry on second‐order sequence preferences.

Physics-based protein structure refinement in the era of artificial intelligence

Protein structure refinement is the last step in protein structure prediction pipelines. Physics-based refinement via molecular dynamics (MD) simulations has made significant progress during recent years. During CASP14, we tested a new refinement protocol based on an improved sampling strategy via MD simulations. MD simulations were carried out at an elevated temperature (360 K). An optimized use of biasing restraints and the use of multiple starting models led to enhanced sampling. The new protocol generally improved the model quality. In comparison with our previous protocols, the CASP14 protocol showed clear improvements. Our approach was successful with most initial models, many based on deep learning methods. However, we found that our approach was not able to refine machine-learning models from the AlphaFold2 group, often decreasing already high initial qualities. To better understand the role of refinement given new types of models based on machine-learning, a detailed analysis via MD simulations and Markov state modeling is presented here. We continue to find that MD-based refinement has the potential to improve AI predictions. We also identified several practical issues that make it difficult to realize that potential. Increasingly important is the consideration of inter-domain and oligomeric contacts in simulations; the presence of large kinetic barriers in refinement pathways also continues to present challenges. Finally, we provide a perspective on how physics-based refinement could continue to play a role in the future for improving initial predictions based on machine learning-based methods.

Assessing the utility of CASP14 models for molecular replacement

The assessment of CASP models for utility in molecular replacement is a measure of their use in a valuable real‐world application. In CASP7, the metric for molecular replacement assessment involved full likelihood‐based molecular replacement searches; however, this restricted the assessable targets to crystal structures with only one copy of the target in the asymmetric unit, and to those where the search found the correct pose. In CASP10, full molecular replacement searches were replaced by likelihood‐based rigid‐body refinement of models superimposed on the target using the LGA algorithm, with the metric being the refined log‐likelihood‐gain (LLG) score. This enabled multi‐copy targets and very poor models to be evaluated, but a significant further issue remained: the requirement of diffraction data for assessment. We introduce here the relative‐expected‐LLG (reLLG), which is independent of diffraction data. This reLLG is also independent of any crystal form, and can be calculated regardless of the source of the target, be it X‐ray, NMR or cryo‐EM. We calibrate the reLLG against the LLG for targets in CASP14, showing that it is a robust measure of both model and group ranking. Like the LLG, the reLLG shows that accurate coordinate error estimates add substantial value to predicted models. We find that refinement by CASP groups can often convert an inadequate initial model into a successful MR search model. Consistent with findings from others, we show that the AlphaFold2 models are sufficiently good, and reliably so, to surpass other current model generation strategies for attempting molecular replacement phasing.

Capabilities of bioinformatics tools for optimizing physicochemical features of proteins used in Nano biosensors: A short overview of the tools related to bioinformatics

Protein-protein ligand is one of the most detection methods used in Nano biosensors. Based on the advantage of specific docking between two special 3D structures, they have become a potent candidate in bioanalysis and Nanodiagnostic tools. These tools lease users to do a simple, fast, cost-effective, sensitive, and specific detection of molecular biomarkers in real samples. Recent advantages of using protein-protein ligand Nano-biosensors application is remarkable due to its special docking that refers to each protein unique 3D conformation. However, it challenges different problems such as low rate of docking and hard process for fixation on the basic layer. These challenges make developers to optimize the structure and functions of proteins. The process has different Nano scale calculation that could be done with algorithms and solutions are available as bioinformatics tools. This article aimed to have a short overview of the abilities of bioinformatics tools for modeling and optimization of physiochemical features of proteins in Nano scale.

•

Nano biosensors use different strategies which based on docking between two molecules to detect and identify different proteins.

•

Molecular docking between transducer in Nano biosensors and proteins rely on physicochemical features of transducer, protein and docking strategy.

•

Nano bioinformatics use bioinformatics tools and algorithms as a collective solution for developing functional structure in Nano scale.

•

Nano bioinformatics use different bioinformatics tools to optimize physicochemical features of proteins as a new approach in Nano biosensors and drug discovery.

DeepRefiner: high-accuracy protein structure refinement by deep network calibration

The DeepRefiner webserver, freely available at http://watson.cse.eng.auburn.edu/DeepRefiner/, is an interactive and fully configurable online system for high-accuracy protein structure refinement. Fuelled by deep learning, DeepRefiner offers the ability to leverage cutting-edge deep neural network architectures which can be calibrated for on-demand selection of adventurous or conservative refinement modes targeted at degree or consistency of refinement. The method has been extensively tested in the Critical Assessment of Techniques for Protein Structure Prediction (CASP) experiments under the group name 'Bhattacharya-Server' and was officially ranked as the No. 2 refinement server in CASP13 (second only to 'Seok-server' and outperforming all other refinement servers) and No. 2 refinement server in CASP14 (second only to 'FEIG-S' and outperforming all other refinement servers including 'Seok-server'). The DeepRefiner web interface offers a number of convenient features, including (i) fully customizable refinement job submission and validation; (ii) automated job status update, tracking, and notifications; (ii) interactive and interpretable web-based results retrieval with quantitative and visual analysis and (iv) extensive help information on job submission and results interpretation via web-based tutorial and help tooltips.

Immunoinformatics-guided design of a multi-epitope vaccine based on the structural proteins of severe acute respiratory syndrome coronavirus 2

Coronavirus disease 2019 (COVID-19) is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), resulting in a contagious respiratory tract infection that has become a global burden since the end of 2019. Notably, fewer patients infected with SARS-CoV-2 progress from acute disease onset to death compared with the progression rate associated with two other coronaviruses, SARS-CoV and Middle East respiratory syndrome coronavirus (MERS-CoV). Several research organizations and pharmaceutical industries have attempted to develop successful vaccine candidates for the prevention of COVID-19. However, increasing evidence indicates that the SARS-CoV-2 genome undergoes frequent mutation; thus, an adequate analysis of the viral strain remains necessary to construct effective vaccines. The current study attempted to design a multi-epitope vaccine by utilizing an approach based on the SARS-CoV-2 structural proteins. We predicted the antigenic T- and B-lymphocyte responses to four structural proteins after screening all structural proteins according to specific characteristics. The predicted epitopes were combined using suitable adjuvants and linkers, and a secondary structure profile indicated that the vaccine shared similar properties with the native protein. Importantly, the molecular docking analysis and molecular dynamics simulations revealed that the constructed vaccine possessed a high affinity for toll-like receptor 4 (TLR4). In addition, multiple descriptors were obtained from the simulation trajectories, including the root-mean-square deviation (RMSD), root-mean-square fluctuation (RMSF), solvent-accessible surface area (SASA), and radius of gyration (R g), demonstrating the rigid nature and inflexibility of the vaccine and receptor molecules. In addition, codon optimization, based on Escherichia coli K12, was used to determine the GC content and the codon adaptation index (CAI) value, which further followed for the incorporation into the cloning vector pET28+(a). Collectively, these findings suggested that the constructed vaccine could be used to modulate the immune reaction against SARS-CoV-2.

Molecular free energy optimization on a computational graph

Free energy is arguably the most important property of molecular systems. Despite great progress in both its efficient estimation by scoring functions/potentials and more rigorous computation based on extensive sampling, we remain far from accurately predicting and manipulating biomolecular structures and their interactions. There are fundamental limitations, including accuracy of interaction description and difficulty of sampling in high dimensional space, to be tackled. Computational graph underlies major artificial intelligence platforms and is proven to facilitate training, optimization and learning. Combining autodifferentiation, coordinates transformation and generalized solvation free energy theory, we construct a computational graph infrastructure to realize seamless integration of fully trainable local free energy landscape with end to end differentiable iterative free energy optimization. This new framework drastically improves efficiency by replacing local sampling with differentiation. Its specific implementation in protein structure refinement achieves superb efficiency and competitive accuracy when compared with state of the art all-atom mainstream methods.

Improved Sampling Strategies for Protein Model Refinement Based on Molecular Dynamics Simulation

Protein structures provide valuable information for understanding biological processes. Protein structures can be determined by experimental methods such as X-ray crystallography, nuclear magnetic resonance spectroscopy, or cryogenic electron microscopy. As an alternative, in silico methods can be used to predict protein structures. These methods utilize protein structure databases for structure prediction via template-based modeling or for training machine-learning models to generate predictions. Structure prediction for proteins distant from proteins with known structures often results in lower accuracy with respect to the true physiological structures. Physics-based protein model refinement methods can be applied to improve model accuracy in the predicted models. Refinement methods rely on conformational sampling around the predicted structures, and if structures closer to the native states are sampled, improvements in the model quality become possible. Molecular dynamics simulations have been especially successful for improving model qualities but although consistent refinement can be achieved, the improvements in model qualities are still moderate. To extend the refinement performance of a simulation-based protocol, we explored new schemes that focus on optimized use of biasing functions and the application of increased simulation temperatures. In addition, we tested the use of alternative initial models so that the simulations can explore the conformational space more broadly. Based on the insights of this analysis, we are proposing a new refinement protocol that significantly outperformed previous state-of-the-art molecular dynamics simulation-based protocols in the benchmark tests described here.

Engineering a novel subunit vaccine against SARS-CoV-2 by exploring immunoinformatics approach

As the number of infections and deaths caused by the recent COVID-19 pandemic is increasing dramatically day-by-day, scientists are rushing towards developing possible countermeasures to fight the deadly virus, SARS-CoV-2. Although many efforts have already been put forward for developing potential vaccines; however, most of them are proved to possess negative consequences. Therefore, in this study, immunoinformatics methods were exploited to design a novel epitope-based subunit vaccine against the SARS-CoV-2, targeting four essential proteins of the virus i.e., spike glycoprotein, nucleocapsid phosphoprotein, membrane glycoprotein, and envelope protein. The highly antigenic, non-allergenic, non-toxic, non-human homolog, and 100% conserved (across other isolates from different regions of the world) epitopes were used for constructing the vaccine. In total, fourteen CTL epitopes and eighteen HTL epitopes were used to construct the vaccine. Thereafter, several in silico validations i.e., the molecular docking, molecular dynamics simulation (including the RMSF and RMSD studies), and immune simulation studies were also performed which predicted that the designed vaccine should be quite safe, effective, and stable within the biological environment. Finally, in silico cloning and codon adaptation studies were also conducted to design an effective mass production strategy of the vaccine. However, more in vitro and in vivo studies are required on the predicted vaccine to finally validate its safety and efficacy.

Toxoplasma gondii ROP38 protein: Bioinformatics analysis for vaccine design improvement against toxoplasmosis

Rhoptry proteins (ROPs) play a significant role in various stages of Toxoplasma gondii (T. gondii) life cycle, being critical for both invasion and intracellular survival. ROP38 is a key manipulator of host gene expression and has a function in tachyzoite to bradyzoite conversion. In this study, we've employed various bioinformatics online tools for immunogenicity prediction of ROP38 protein, comprising physico-chemical, antigenic and allergenic profiles, transmembrane domain, subcellular localization, post-translational modification (PTM) sites, secondary and 3D structure, B-cell, MHC-binding and cytotoxic T-lymphocyte (CTL) epitopes. The findings showed 54 PTM sites without a transmembrane domain. Also, ROP38 was proved a non-allergenic and antigenic protein. The protein had Sec signal peptide (Sec/SPI) with 0.8762 likelihood. The secondary structure included 52.68% random coil, 29.57% alpha helix and 17.74% extended strand. Based on Ramachandran plot output for refined model, 95.3%, 3.4%, and 1.4% of amino acid residues were incorporated in the favored, allowed, and outlier regions, respectively. B-cell epitopes TFPGDDIQTSS (67-72) and KAKNKWGRTRYTLQG (207-221) as well as T-cell epitope LSPVGFFTAL (6-15) possessed the highest antigenic index in the protein sequence. This paper is a premise for further researches, and provides insights for the development of a suitable vaccine against toxoplasmosis. More empirical studies are required using the ROP38 alone or in combination with other antigens/epitopes in the future.

In-depth computational analysis of calcium-dependent protein kinase 3 of Toxoplasma gondii provides promising targets for vaccination

Purpose

The Toxoplasma gondii calcium-dependent protein kinase-3 (CDPK3) is a key enzyme for parasite egress, control of calcium-dependent permeabilization in parasitophorous vacuole membrane and tissue cyst formation. In this study, we comprehensively explored the bioinformatics features of this protein to improve vaccine design against T. gondii.

Materials and Methods

Various web servers were employed for the analysis of physico-chemical properties, post-translational modifications, localization in the subcellular milieu, secondary and tertiary structures, as well as B-cell, major histocompatibility complex (MHC)-binding and cytotoxic T-lymphocyte (CTL) epitopes.

Results

This protein was a 537 amino acid antigenic and non-allergenic molecule with a molecular weight of 60.42 kDa, a grand average of hydropathicity score of −0.508, and aliphatic index of 79.50. There exists 46.74% alpha helix, 12.48% extended strand, and 40.78% random coil in the secondary structure. Ramachandran plot of the refined model demonstrated 99.3%, 0.7%, and 0.0% of residues in the favored, allowed and outlier areas, respectively. Besides, various potential B-cell (continuous and conformational), MHC-binding and CTL epitopes were predicted for Toxoplasma CDPK3 protein.

Conclusion

This article provides a foundation for further investigations, and laid a theoretical basis for the development of an appropriate vaccine against T. gondii infection.

Designing a multi-epitope vaccine against SARS-CoV-2: an immunoinformatics approach

Ongoing COVID-19 outbreak has raised a drastic challenge to global public health security. Most of the patients with COVID-19 suffer from mild flu-like illnesses such as cold and fever; however, few percentages of the patients progress from severe illness to death, mostly in an immunocompromised individual. The causative agent of COVID-19 is an RNA virus known as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Despite these debilitating conditions, no medication to stop the disease progression or vaccination is available till now. Therefore, we aimed to formulate a multi-epitope vaccine against SARS-CoV-2 by utilizing an immunoinformatics approach. For this purpose, we used the SARS-CoV-2 spike glycoprotein to determine the immunodominant T- and B-cell epitopes. After rigorous assessment, we designed a vaccine construct using four potential epitopes from each of the three epitope classes such as cytotoxic T-lymphocytes, helper T-lymphocyte, and linear B-lymphocyte epitopes. The designed vaccine was antigenic, immunogenic, and non-allergenic with suitable physicochemical properties and has higher solubility. More importantly, the predicted vaccine structure was similar to the native protein. Further investigations indicated a strong and stable binding interaction between the vaccine and the toll-like receptor (TLR4). Strong binding stability and structural compactness were also evident in molecular dynamics simulation. Furthermore, the computer-generated immune simulation showed that the vaccine could trigger real-life-like immune responses upon administration into humans. Finally, codon optimization based on Escherichia coli K12 resulted in optimal GC content and higher CAI value followed by incorporating it into the cloning vector pET28+(a). Overall, these results suggest that the designed peptide vaccine can serve as an excellent prophylactic candidate against SARS-CoV-2.

Communicated by Ramaswamy H. Sarma

refineD: improved protein structure refinement using machine learning based restrained relaxation

MOTIVATION

Protein structure refinement aims to bring moderately accurate template-based protein models closer to the native state through conformational sampling. However, guiding the sampling towards the native state by effectively using restraints remains a major issue in structure refinement.

RESULTS

Here, we develop a machine learning based restrained relaxation protocol that uses deep discriminative learning based binary classifiers to predict multi-resolution probabilistic restraints from the starting structure and subsequently converts these restraints to be integrated into Rosetta all-atom energy function as additional scoring terms during structure refinement. We use four restraint resolutions as adopted in GDT-HA (0.5, 1, 2 and 4 Å), centered on the Cα atom of each residue that are predicted by ensemble of four deep discriminative classifiers trained using combinations of sequence and structure-derived features as well as several energy terms from Rosetta centroid scoring function. The proposed method, refineD, has been found to produce consistent and substantial structural refinement through the use of cumulative and non-cumulative restraints on 150 benchmarking targets. refineD outperforms unrestrained relaxation strategy or relaxation that is restrained to starting structures using the FastRelax application of Rosetta or atomic-level energy minimization based ModRefiner method as well as molecular dynamics (MD) simulation based FG-MD protocol. Furthermore, by adjusting restraint resolutions, the method addresses the tradeoff that exists between degree and consistency of refinement. These results demonstrate a promising new avenue for improving accuracy of template-based protein models by effectively guiding conformational sampling during structure refinement through the use of machine learning based restraints.

AVAILABILITY AND IMPLEMENTATION

http://watson.cse.eng.auburn.edu/refineD/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Proteome-wide screening for designing a multi-epitope vaccine against emerging pathogen Elizabethkingia anophelis using immunoinformatic approaches

Elizabethkingia anophelis is an emerging human pathogen causing neonatal meningitis, catheter-associated infections and nosocomial outbreaks with high mortality rates. Besides, they are resistant to most antibiotics used in empirical therapy. In this study, therefore, we used immunoinformatic approaches to design a prophylactic peptide vaccine against E. anophelis as an alternative preventive measure. Initially, cytotoxic T-lymphocyte (CTL), helper T-lymphocyte (HTL), and linear B-lymphocyte (LBL) epitopes were predicted from the highest antigenic protein. The CTL and HTL epitopes together had a population coverage of 99.97% around the world. Eventually, six CTL, seven HTL, and two LBL epitopes were selected and used to construct a multi-epitope vaccine. The vaccine protein was found to be highly immunogenic, non-allergenic, and non-toxic. Codon adaptation and in silico cloning were performed to ensure better expression within E. coli K12 host system. The stability of the vaccine structure was also improved by disulphide bridging. In addition, molecular docking and dynamics simulation revealed strong and stable binding affinity between the vaccine and toll-like receptor 4 (TLR4) molecule. The immune simulation showed higher levels of T-cell and B-cell activities which was in coherence with actual immune response. Repeated exposure simulation resulted in higher clonal selection and faster antigen clearance. Nevertheless, experimental validation is required to ensure the immunogenic potency and safety of this vaccine to control E. anophelis infection in the future.Communicated by Ramaswamy H. Sarma.

Critical assessment of methods of protein structure prediction (CASP)-Round XIII

CASP (critical assessment of structure prediction) assesses the state of the art in modeling protein structure from amino acid sequence. The most recent experiment (CASP13 held in 2018) saw dramatic progress in structure modeling without use of structural templates (historically "ab initio" modeling). Progress was driven by the successful application of deep learning techniques to predict inter-residue distances. In turn, these results drove dramatic improvements in three-dimensional structure accuracy: With the proviso that there are an adequate number of sequences known for the protein family, the new methods essentially solve the long-standing problem of predicting the fold topology of monomeric proteins. Further, the number of sequences required in the alignment has fallen substantially. There is also substantial improvement in the accuracy of template-based models. Other areas-model refinement, accuracy estimation, and the structure of protein assemblies-have again yielded interesting results. CASP13 placed increased emphasis on the use of sparse data together with modeling and chemical crosslinking, SAXS, and NMR all yielded more mature results. This paper summarizes the key outcomes of CASP13. The special issue of PROTEINS contains papers describing the CASP13 assessments in each modeling category and contributions from the participants.

Evaluation of model refinement in CASP13

Performance in the model refinement category of the 13th round of Critical Assessment of Structure Prediction (CASP13) is assessed, showing that some groups consistently improve most starting models whereas the majority of participants continue to degrade the starting model on average. Using the ranking formula developed for CASP12, it is shown that only 7 of 32 groups perform better than a “naïve predictor” who just submits the starting model. Common features in their approaches include a dependence on physics‐based force fields to judge alternative conformations and the use of molecular dynamics to relax models to local minima, usually with some restraints to prevent excessively large movements. In addition to the traditional CASP metrics that focus largely on the quality of the overall fold, alternative metrics are evaluated, including comparisons of the main‐chain and side‐chain torsion angles, and the utility of the models for solving crystal structures by the molecular replacement method. It is proposed that the introduction of these metrics, as well as consideration of the accuracy of coordinate error estimates, would improve the discrimination between good and very good models.

Applications of Molecular Dynamics Simulation in Structure Prediction of Peptides and Proteins

Compared with rapid accumulation of protein sequences from high-throughput DNA sequencing, obtaining experimental 3D structures of proteins is still much more difficult, making protein structure prediction (PSP) potentially very useful. Currently, a vast majority of PSP efforts are based on data mining of known sequences, structures and their relationships (informatics-based). However, if closely related template is not available, these methods are usually much less reliable than experiments. They may also be problematic in predicting the structures of naturally occurring or designed peptides. On the other hand, physics-based methods including molecular dynamics (MD) can utilize our understanding of detailed atomic interactions determining biomolecular structures. In this mini-review, we show that all-atom MD can predict structures of cyclic peptides and other peptide foldamers with accuracy similar to experiments. Then, some notable successes in reproducing experimental 3D structures of small proteins through MD simulations (some with replica-exchange) of the folding were summarized. We also describe advancements of MD-based refinement of structure models, and the integration of limited experimental or bioinformatics data into MD-based structure modeling.

Driven to near-experimental accuracy by refinement via molecular dynamics simulations

Protein model refinement has been an essential part of successful protein structure prediction. Molecular dynamics simulation-based refinement methods have shown consistent improvement of protein models. There had been progress in the extent of refinement for a few years since the idea of ensemble averaging of sampled conformations emerged. There was little progress in CASP12 because conformational sampling was not sufficiently diverse due to harmonic restraints. During CASP13, a new refinement method was tested that achieved significant improvements over CASP12. The new method intended to address previous bottlenecks in the refinement problem by introducing new features. Flat-bottom harmonic restraints replaced harmonic restraints, sampling was performed iteratively, and a new scoring function and selection criteria were used. The new protocol expanded conformational sampling at reduced computational costs. In addition to overall improvements, some models were refined significantly to near-experimental accuracy.

ReFOLD: a server for the refinement of 3D protein models guided by accurate quality estimates

ReFOLD is a novel hybrid refinement server with integrated high performance global and local Accuracy Self Estimates (ASEs). The server attempts to identify and to fix likely errors in user supplied 3D models of proteins via successive rounds of refinement. The server is unique in providing output for multiple alternative refined models in a way that allows users to quickly visualize the key residue locations, which are likely to have been improved. This is important, as global refinement of a full chain model may not always be possible, whereas local regions, or individual domains, can often be much improved. Thus, users may easily compare the specific regions of the alternative refined models in which they are most interested e.g. key interaction sites or domains. ReFOLD was used to generate hundreds of alternative refined models for the CASP12 experiment, boosting our group's performance in the main tertiary structure prediction category. Our successful refinement of initial server models combined with our built-in ASEs were instrumental to our second place ranking on Template Based Modeling (TBM) and Free Modeling (FM)/TBM targets. The ReFOLD server is freely available at: http://www.reading.ac.uk/bioinf/ReFOLD/.

Methods for the Refinement of Protein Structure 3D Models

The refinement of predicted 3D protein models is crucial in bringing them closer towards experimental accuracy for further computational studies. Refinement approaches can be divided into two main stages: The sampling and scoring stages. Sampling strategies, such as the popular Molecular Dynamics (MD)-based protocols, aim to generate improved 3D models. However, generating 3D models that are closer to the native structure than the initial model remains challenging, as structural deviations from the native basin can be encountered due to force-field inaccuracies. Therefore, different restraint strategies have been applied in order to avoid deviations away from the native structure. For example, the accurate prediction of local errors and/or contacts in the initial models can be used to guide restraints. MD-based protocols, using physics-based force fields and smart restraints, have made significant progress towards a more consistent refinement of 3D models. The scoring stage, including energy functions and Model Quality Assessment Programs (MQAPs) are also used to discriminate near-native conformations from non-native conformations. Nevertheless, there are often very small differences among generated 3D models in refinement pipelines, which makes model discrimination and selection problematic. For this reason, the identification of the most native-like conformations remains a major challenge.

Higher Accuracy Achieved in the Simulations of Protein Structure Refinement, Protein Folding, and Intrinsically Disordered Proteins Using Polarizable Force Fields

The accuracy of molecular mechanics force fields is of vital importance in biomolecular simulations. However, the admittedly more accurate polarizable force fields were recently reported to be less able to reproduce the experimental properties in comparison to additive force fields in some cases. Here, we perform long-time-scale molecular dynamics simulations to systematically evaluate the effect of explicit electronic polarization in polarizable force fields. The results show that the inclusion of electrostatic polarization effect in polarizable force fields can improve their accuracies in protein structure refinement and generate conformational ensembles more approximate to experiments for intrinsically disordered proteins. In contrast, it is difficult for polarizable force fields to approach the native structure, let alone to predict the native state when it is unknown a priori in the real protein structure predictions. We speculate that these effects might be attributed to the preference of protein-water interactions in polarizable force fields.

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.

Refining protein structures using enhanced sampling techniques with restraints derived from an ensemble-based model

This paper reports a method for high-accuracy protein structural refinement, which is a direct extension of the method in our recent publication (Zang, J Chem Phys 2018; 149:072319). It combines a parallel continuous simulated tempering (PCST) method with a temperature-dependent restraint and a blind model selection scheme. In this work, a single-reference-based restraint in previous work was changed to an ensemble-based model (EBM), in which the non-bonded Lennard-Jones term for each contacting atomic pair in previous restraining potential was replaced by a multi-Gaussian function whose parameters are derived from an ensemble of structures such as the ones from various CASP participating groups. The purpose of EBM is to take advantage of partial "correctness" distributed among members of the structural ensemble. Totally 18 targets were refined from the refinement category of CASP10, CASP11 and CASP12. In Top-1 group, 11 out of 18 targets had better models (greater GDT_TS scores) than the CASPR participants. In Top-5 group, nine out of 18 were better. Our results show that PCST-EBM method can considerably improve the low-accuracy structures.

Predicting improved protein conformations with a temporal deep recurrent neural network

Accurate protein structure prediction from amino acid sequence is still an unsolved problem. The most reliable methods centre on template based modelling. However, the accuracy of these models entirely depends on the availability of experimentally resolved homologous template structures. In order to generate more accurate models, extensive physics based molecular dynamics (MD) refinement simulations are performed to sample many different conformations to find improved conformational states. In this study, we propose a deep recurrent network model, called DeepTrajectory, that is able to identify these improved conformational states, with high precision, from a variety of different MD based sampling protocols. The proposed model learns the temporal patterns of features computed from MD trajectory data in order to classify whether each recorded simulation snapshot is an improved quality conformational state, decreased quality conformational state or whether there is no perceivable change in state with respect to the starting conformation. The model was trained and tested on 904 trajectories from 42 different protein systems with a cumulative number of more than 1.7 million snapshots. We show that our model outperforms other state of the art machine-learning algorithms that do not consider temporal dependencies. To our knowledge, DeepTrajectory is the first implementation of a time-dependent deep-learning protocol that is re-trainable and able to adapt to any new MD based sampling procedure, thereby demonstrating how a neural network can be used to learn the latter part of the protein folding funnel.

Improving low-accuracy protein structures using enhanced sampling techniques

In this paper, we report results of using enhanced sampling and blind selection techniques for high-accuracy protein structural refinement. By combining a parallel continuous simulated tempering (PCST) method, previously developed by Zang et al. [J. Chem. Phys. 141, 044113 (2014)], and the structure based model (SBM) as restraints, we refined 23 targets (18 from the refinement category of the CASP10 and 5 from that of CASP12). We also designed a novel model selection method to blindly select high-quality models from very long simulation trajectories. The combined use of PCST-SBM with the blind selection method yielded final models that are better than initial models. For Top-1 group, 7 out of 23 targets had better models (greater global distance test total scores) than the critical assessment of structure prediction participants. For Top-5 group, 10 out of 23 were better. Our results justify the crucial position of enhanced sampling in protein structure prediction and refinement and demonstrate that a considerable improvement of low-accuracy structures is achievable with current force fields.

Protein structure prediction

Predicting 3D structure of protein from its amino acid sequence is one of the most important unsolved problems in biophysics and computational biology. This paper attempts to give a comprehensive introduction of the most recent effort and progress on protein structure prediction. Following the general flowchart of structure prediction, related concepts and methods are presented and discussed. Moreover, brief introductions are made to several widely-used prediction methods and the community-wide critical assessment of protein structure prediction (CASP) experiments.

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.

What makes it difficult to refine protein models further via molecular dynamics simulations?

Protein structure refinement remains a challenging yet important problem as it has the potential to bring already accurate template-based models to near-native resolution. Refinement based on molecular dynamics simulations has been a highly promising approach and the performance of MD-based refinement in the Feig group during CASP12 is described here. During CASP12, sampling was extended well into the microsecond scale, an improved force field was applied, and new protocol variations were tested. Progress over previous rounds of CASP was found to be limited which is analyzed in terms of the quality of the initial models and dependency on the amount of sampling and refinement protocol variations. As current MD-based refinement protocols appear to be reaching a plateau, detailed analysis is presented to provide new insight into the major challenges towards more extensive structure refinement, focusing in particular on sampling with and without restraints.

Assessment of the model refinement category in CASP12

We here report on the assessment of the model refinement predictions submitted to the 12th Experiment on the Critical Assessment of Protein Structure Prediction (CASP12). This is the fifth refinement experiment since CASP8 (2008) and, as with the previous experiments, the predictors were invited to refine selected server models received in the regular (nonrefinement) stage of the CASP experiment. We assessed the submitted models using a combination of standard CASP measures. The coefficients for the linear combination of Z-scores (the CASP12 score) have been obtained by a machine learning algorithm trained on the results of visual inspection. We identified eight groups that improve both the backbone conformation and the side chain positioning for the majority of targets. Albeit the top methods adopted distinctively different approaches, their overall performance was almost indistinguishable, with each of them excelling in different scores or target subsets. What is more, there were a few novel approaches that, while doing worse than average in most cases, provided the best refinements for a few targets, showing significant latitude for further innovation in the field.

Simultaneous refinement of inaccurate local regions and overall structure in the CASP12 protein model refinement experiment

Advances in protein model refinement techniques are required as diverse sources of protein structure information are available from low-resolution experiments or informatics-based computations such as cryo-EM, NMR, homology models, or predicted residue contacts. Given semi-reliable or incomplete structural information, structure quality of a protein model has to be improved by ab initio methods such as energy-based simulation. In this study, we describe a new automatic refinement server method designed to improve locally inaccurate regions and overall structure simultaneously. Locally inaccurate regions may occur in protein structures due to non-convergent or missing information in template structures used in homology modeling or due to intrinsic structural flexibilities not resolved by experimental techniques. However, such variable or dynamic regions often play important functional roles by participating in interactions with other biomolecules or in transitions between different functional states. The new refinement method introduced here utilizes diverse types of geometric operators which drive both local and global changes, and the effect of structure changes and relaxations are accumulated. This resulted in consistent refinement of both local and global structural features. Performance of this method in CASP12 is discussed.