Basic protein structure prediction for the biologist: A review

Structure and dynamics of whole-sequence homology model of ORF3a protein of SARS-CoV-2: An insight from microsecond molecular dynamics simulations

The ORF3a is a large accessory protein in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which plays an important role in virulence and viral replication; especially in inflammasome activation and apoptosis. However,, the existing cryo-EM structure of SARS-CoV-2 ORF3a is incomplete, . making it challenging to understand its structural and functional features. The aim of this study is to investigate the dynamic behaviors of the full-sequence homology model of ORF3a and compare it with the cryo-EM structure using microsecond molecular dynamics simulations. The previous studies indicated that the unresolved residues of the cryo-EM structure are not only involved in the pathogenesis of the SARS-CoV-2 but also exhibit a significant antigenicity. The dynamics scenario of homology model revealed higher RMSD, Rg, and SASA values with stable pattern when compared to the cryo-EM structure. Moreover, the RMSF analysis demonstrated higher fluctuations at specific positions (1-43, 97-110, 172-180, 219-243) in the model structure, whereas the cryo-EM structure displayed lower overall drift (except 1-43) in comparison to the model structure.Secondary structural features indicated that a significant unfolding in the transmembrane domains and β-strand at positions 166 to 172, affecting the stability and compactness of the cryo-EM structure , whereas the model exhibited noticeable unfolding in transmembrane domains and small-coiled regions in the N-terminal. , The results from molecular docking and steered molecular dynamics investigations showed the model structure had a greater number of non-bonding interactions, leading to enhanced stability when compared to the cryo-EM structure. Consequently, higher forces were necessary for unbinding of the baricitinib and ruxolitinib inhibitors from the model structure.. Our findings can help better understanding of the significance of unresolved residues at the molecular level. Additionally, this information can guide researchers for experimental endeavors aimed at completing the full-sequence structure of the ORF3a.Communicated by Ramaswamy H. Sarma.

Comparative modelling studies of fruit bromelain using molecular dynamics simulation

Fruit bromelain is a cysteine protease accumulated in pineapple fruits. This proteolytic enzyme has received high demand for industrial and therapeutic applications. In this study, fruit bromelain sequences QIM61759, QIM61760 and QIM61761 were retrieved from the National Center for Biotechnology Information (NCBI) Genbank Database. The tertiary structure of fruit bromelain QIM61759, QIM61760 and QIM61761 was generated by using MODELLER. The result revealed that the local stereochemical quality of the generated models was improved by using multiple templates during modelling process. Moreover, by comparing with the available papain model, structural analysis provides an insight on how pro-peptide functions as a scaffold in fruit bromelain folding and contributing to inactivation of mature protein. The structural analysis also disclosed the similarities and differences between these models. Lastly, thermal stability of fruit bromelain was studied. Molecular dynamics simulation of fruit bromelain structures at several selected temperatures demonstrated how fruit bromelain responds to elevation of temperature.

Improved scFv Anti-LOX-1 Binding Activity by Fusion with LOX-1-Binding Peptides

The oxidized low-density lipoprotein receptor-1 (LOX-1) targeted single-chain variable fragment (scFvs) is a promising molecule for the targeted delivery of imaging and therapeutic molecules of atherosclerotic diseases; however, its applications are limited by the inherent low antigen affinity. In this study, the three-dimensional (3D) model of the anti-LOX-1 scFv was constructed and its docking with the LOX-1 protein was developed. To improve the LOX-1-binding activity, the anti-LOX-1 scFv was designed to fuse with one of three LOX-1-binding heptapeptides, LTPATAI, FQTPPQL, and LSIPPKA, at its N-terminus and C-terminus and in the linker region, which have different LOX-1-binding interfaces with the anti-LOX-1 scFv analyzed by an array of computational approaches. These scFv/peptide fusions were constructed, successfully expressed in Brevibacillus choshinensis hosts, and purified by a two-step column purification process. The antigen binding activity, structural characteristics, thermal stability, and stability in serum of these fusion proteins were examined. Results showed that the scFv with N-terminal fusing peptides proteins demonstrated increased LOX-1-binding activity without decrease in stability. These findings will help increase the application efficacy of LOX-1 targeting scFv in LOX-1-based therapy.

In silico approaches and chemical space of anti-P-type ATPase compounds for discovering new antituberculous drugs

Tuberculosis (TB) is one of the most important public health problems around the world. The emergence of multi-drug-resistant (MDR) and extensively drug-resistant (XDR) Mycobacterium tuberculosis strains has driven the finding of alternative anti-TB targets. In this context, P-type ATPases are interesting therapeutic targets due to their key role in ion homeostasis across the plasma membrane and the mycobacterial survival inside macrophages. In this review, in silico and experimental strategies used for the rational design of new anti-TB drugs are presented; in addition, the chemical space distribution based on the structure and molecular properties of compounds with anti-TB and anti-P-type ATPase activity is discussed. The chemical space distribution compared to public compound libraries demonstrates that natural product libraries are a source of novel chemical scaffolds with potential anti-P-type ATPase activity. Furthermore, compounds that experimentally display anti-P-type ATPase activity belong to a chemical space of molecular properties comparable to that occupied by those approved for oral use, suggesting that these kinds of molecules have a good pharmacokinetic profile (drug-like) for evaluation as potential anti-TB drugs.

General overview on structure prediction of twilight-zone proteins

Protein structure prediction from amino acid sequence has been one of the most challenging aspects in computational structural biology despite significant progress in recent years showed by critical assessment of protein structure prediction (CASP) experiments. When experimentally determined structures are unavailable, the predictive structures may serve as starting points to study a protein. If the target protein consists of homologous region, high-resolution (typically <1.5 Å) model can be built via comparative modelling. However, when confronted with low sequence similarity of the target protein (also known as twilight-zone protein, sequence identity with available templates is less than 30%), the protein structure prediction has to be initiated from scratch. Traditionally, twilight-zone proteins can be predicted via threading or ab initio method. Based on the current trend, combination of different methods brings an improved success in the prediction of twilight-zone proteins. In this mini review, the methods, progresses and challenges for the prediction of twilight-zone proteins were discussed.

In silico determination and validation of baumannii acinetobactin utilization a structure and ligand binding site

Acinetobacter baumannii is a deadly nosocomial pathogen. Iron is an essential element for the pathogen. Under iron-restricted conditions, the bacterium expresses iron-regulated outer membrane proteins (IROMPs). Baumannii acinetobactin utilization (BauA) is the most important member of IROMPs in A. baumannii. Determination of its tertiary structure could help deduction of its functions and its interactions with ligands. The present study unveils BauA 3D structure via in silico approaches. Apart from ab initio, other rational methods such as homology modeling and threading were invoked to achieve the purpose. For homology modeling, BLAST was run on the sequence in order to find the best template. The template was then served to model the 3D structure. All the models built were evaluated qualitatively. The best model predicted by LOMETS was selected for analyses. Refinement of 3D structure as well as determination of its clefts and ligand binding sites was carried out on the structure. In contrast to the typical trimeric arrangement found in porins, BauA is monomeric. The barrel is formed by 22 antiparallel transmembrane β -strands. There are short periplasmic turns and longer surface-located loops. An N-terminal domain referred to either as the cork, the plug, or the hatch domain occludes the β -barrel.

Quantitative structure-activity relationships and docking studies of calcitonin gene-related peptide antagonists

Defining the role of calcitonin gene-related peptide in migraine pathogenesis could lead to the application of calcitonin gene-related peptide antagonists as novel migraine therapeutics. In this work, quantitative structure-activity relationship modeling of biological activities of a large range of calcitonin gene-related peptide antagonists was performed using a panel of physicochemical descriptors. The computational studies evaluated different variable selection techniques and demonstrated shuffling stepwise multiple linear regression to be superior over genetic algorithm-multiple linear regression. The linear quantitative structure-activity relationship model revealed better statistical parameters of cross-validation in comparison with the non-linear support vector regression technique. Implementing only five peptide descriptors into this linear quantitative structure-activity relationship model resulted in an extremely robust and highly predictive model with calibration, leave-one-out and leave-20-out validation R(2) of 0.9194, 0.9103, and 0.9214, respectively. We performed docking of the most potent calcitonin gene-related peptide antagonists with the calcitonin gene-related peptide receptor and demonstrated that peptide antagonists act by blocking access to the peptide-binding cleft. We also demonstrated the direct contact of residues 28-37 of the calcitonin gene-related peptide antagonists with the receptor. These results are in agreement with the conclusions drawn from the quantitative structure-activity relationship model, indicating that both electrostatic and steric factors should be taken into account when designing novel calcitonin gene-related peptide antagonists.