Data deposition and annotation at the worldwide protein data bank

Structural Characterization of Peptide Antibodies

The role of proteins as very effective immunogens for the generation of antibodies is indisputable. Nevertheless, cases in which protein usage for antibody production is not feasible or convenient compelled the creation of a powerful alternative consisting of synthetic peptides. Synthetic peptides can be modified to obtain desired properties or conformation, tagged for purification, isotopically labeled for protein quantitation or conjugated to immunogens for antibody production. The antibodies that bind to these peptides represent an invaluable tool for biological research and discovery. To better understand the underlying mechanisms of antibody-antigen interaction, here, we present a pipeline developed by us to structurally classify immunoglobulin antigen binding sites and to infer key sequence residues and other variables that have a prominent role in each structural class.

The Interactome of Protein, DNA, and RNA

Proteins participate in many processes of the organism and are very important for maintaining the health of the organism. However, proteins cannot function independently in the body. They must interact with proteins, DNA, RNA, and other substances to perform biological functions and maintain the body's health. At present, there are many experimental methods and software tools that can detect and predict the interaction between proteins and other substances. There are also many databases that record the interaction between proteins and other substances. This article mainly describes protein-protein, protein-DNA, and protein-RNA interactions in detail by introducing some commonly used experimental methods, the software tools produced with the accumulation of experimental data and the rapid development of machine learning, and the related databases that record the relationship between proteins and some substances. By this review, we hope that through the analysis and summary of various aspects, it will be convenient for researchers to conduct further research on protein interactions.

Structural Bioinformatic Survey of Protein-Small Molecule Interfaces Delineates the Role of Glycine in Surface Pocket Formation

With a structural bioinformatic approach, we have explored amino acid compositions at PISA defined interfaces between small molecules and proteins that are contained in an optimized subset of 11,351 PDB files. The use of a series of restrictions, to prevent redundancy and biases from interactions between amino acids with charged side chains and ions, yielded a final data set of 45,230 protein-small molecule interfaces. We have compared occurrences of natural amino acids in surface exposed regions and binding sites for all the proteins of our data set. From our structural bioinformatic survey, the most relevant signal arose from the unexpected Gly abundance at enzyme catalytic sites. This finding suggested that Gly must have a fundamental role in stabilizing concave protein surface moieties. Subsequently, we have tried to predict the effect of in silico Gly mutations in hen egg white lysozyme to optimize those conditions that can reshape the protein surface with the appearance of new pockets. Replacing amino acids having bulky side chains with Gly in specific protein regions seems a feasible way for designing proteins with additional surface pockets, which can alter protein surface dynamics, therefore, representing controllable switches for protein activity.

Proteomic Network Analysis of Bronchoalveolar Lavage Fluid in Ex-Smokers to Discover Implicated Protein Targets and Novel Drug Treatments for Chronic Obstructive Pulmonary Disease

Bronchoalveolar lavage of the epithelial lining fluid (BALF) can sample the profound changes in the airway lumen milieu prevalent in chronic obstructive pulmonary disease (COPD). We compared the BALF proteome of ex-smokers with moderate COPD who are not in exacerbation status to non-smoking healthy control subjects and applied proteome-scale translational bioinformatics approaches to identify potential therapeutic protein targets and drugs that modulate these proteins for the treatment of COPD. Proteomic profiles of BALF were obtained from (1) never-smoker control subjects with normal lung function (n = 10) or (2) individuals with stable moderate (GOLD stage 2, FEV1 50−80% predicted, FEV1/FVC < 0.70) COPD who were ex-smokers for at least 1 year (n = 10). After identifying potential crucial hub proteins, drug−proteome interaction signatures were ranked by the computational analysis of novel drug opportunities (CANDO) platform for multiscale therapeutic discovery to identify potentially repurposable drugs. Subsequently, a literature-based knowledge graph was utilized to rank combinations of drugs that most likely ameliorate inflammatory processes. Proteomic network analysis demonstrated that 233 of the >1800 proteins identified in the BALF were significantly differentially expressed in COPD versus control. Functional annotation of the differentially expressed proteins was used to detail canonical pathways containing the differential expressed proteins. Topological network analysis demonstrated that four putative proteins act as central node proteins in COPD. The drugs with the most similar interaction signatures to approved COPD drugs were extracted with the CANDO platform. The drugs identified using CANDO were subsequently analyzed using a knowledge-based technique to determine an optimal two-drug combination that had the most appropriate effect on the central node proteins. Network analysis of the BALF proteome identified critical targets that have critical roles in modulating COPD pathogenesis, for which we identified several drugs that could be repurposed to treat COPD using a multiscale shotgun drug discovery approach.

A chemical interpretation of protein electron density maps in the worldwide protein data bank

High-quality three-dimensional structural data is of great value for the functional interpretation of biomacromolecules, especially proteins; however, structural quality varies greatly across the entries in the worldwide Protein Data Bank (wwPDB). Since 2008, the wwPDB has required the inclusion of structure factors with the deposition of x-ray crystallographic structures to support the independent evaluation of structures with respect to the underlying experimental data used to derive those structures. However, interpreting the discrepancies between the structural model and its underlying electron density data is difficult, since derived sigma-scaled electron density maps use arbitrary electron density units which are inconsistent between maps from different wwPDB entries. Therefore, we have developed a method that converts electron density values from sigma-scaled electron density maps into units of electrons. With this conversion, we have developed new methods that can evaluate specific regions of an x-ray crystallographic structure with respect to a physicochemical interpretation of its corresponding electron density map. We have systematically compared all deposited x-ray crystallographic protein models in the wwPDB with their underlying electron density maps, if available, and characterized the electron density in terms of expected numbers of electrons based on the structural model. The methods generated coherent evaluation metrics throughout all PDB entries with associated electron density data, which are consistent with visualization software that would normally be used for manual quality assessment. To our knowledge, this is the first attempt to derive units of electrons directly from electron density maps without the aid of the underlying structure factors. These new metrics are biochemically-informative and can be extremely useful for filtering out low-quality structural regions from inclusion into systematic analyses that span large numbers of PDB entries. Furthermore, these new metrics will improve the ability of non-crystallographers to evaluate regions of interest within PDB entries, since only the PDB structure and the associated electron density maps are needed. These new methods are available as a well-documented Python package on GitHub and the Python Package Index under a modified Clear BSD open source license.

PDB2INS: bridging the gap between small-molecule and macromolecular refinement

The open-source Python program PDB2INS is designed to prepare a .ins file for refinement with SHELXL [Sheldrick (2015). Acta Cryst. C71, 3-8], taking atom coordinates and other information from a Protein Data Bank (PDB)-format file. If PDB2INS is provided with a four-character PDB code, both the PDB file and the accompanying mmCIF-format reflection data file (if available) are accessed via the internet from the PDB public archive [Read et al. (2011). Structure, 19, 1395-1412] or optionally from the PDB_REDO server [Joosten, Long, Murshudov & Perrakis (2014). IUCrJ, 1, 213-220]. The SHELX-format .ins (refinement instructions and atomic coordinates) and .hkl (reflection data) files can then be generated without further user intervention, appropriate restraints etc. being added automatically. PDB2INS was tested on the 23 974 X-ray structures deposited in the PDB between 2008 and 2018 that included reflection data to 1.7 Å or better resolution in a recognizable format. After creating the two input files for SHELXL without user intervention, ten cycles of conjugate-gradient least-squares refinement were performed. For 96% of these structures PDB2INS and SHELXL completed successfully without error messages.

Rodrigues GO, M. Lopes CE, Queiroz-Junior CM, Lima CX, Costa VV, Teixeira MM, Bozza FA, Bozza PT, Boechat N, Tanuri A, Souza TML. Yellow fever virus is susceptible to sofosbuvir both in vitro and in vivo

Yellow fever virus (YFV) is a member of the Flaviviridae family. In Brazil, yellow fever (YF) cases have increased dramatically in sylvatic areas neighboring urban zones in the last few years. Because of the high lethality rates associated with infection and absence of any antiviral treatments, it is essential to identify therapeutic options to respond to YFV outbreaks. Repurposing of clinically approved drugs represents the fastest alternative to discover antivirals for public health emergencies. Other Flaviviruses, such as Zika (ZIKV) and dengue (DENV) viruses, are susceptible to sofosbuvir, a clinically approved drug against hepatitis C virus (HCV). Our data showed that sofosbuvir docks onto YFV RNA polymerase using conserved amino acid residues for nucleotide binding. This drug inhibited the replication of both vaccine and wild-type strains of YFV on human hepatoma cells, with EC₅₀ values around 5 μM. Sofosbuvir protected YFV-infected neonatal Swiss mice and adult type I interferon receptor knockout mice (A129^-/-) from mortality and weight loss. Because of its safety profile in humans and significant antiviral effects in vitro and in mice, Sofosbuvir may represent a novel therapeutic option for the treatment of YF. Key-words: Yellow fever virus; Yellow fever, antiviral; sofosbuvir

Yellow fever virus is transmitted by mosquitoes and its infection may be asymptomatic or lead to a wide clinical spectrum ranging from a mild febrile illness to a potentially lethal viral hemorrhagic fever characterized by liver damage. Although a yellow fever vaccine is available, low coverage allows 80,000–200,000 cases and 30,000–60,000 deaths annually worldwide. There are no specific therapy and treatment relies on supportive care, reinforcing an urgent need for antiviral repourposing. Here, we showed that sofosbuvir, clinically approved against hepatitis C, inhibits yellow fever virus replication in liver cell lines and animal models. In vitro, sofosbuvir inhibits viral RNA replication, decreases the number of infected cells and the production of infectious virus particles. These data is particularly relevante since the liver is the main target of yellow fever infection. Sofosbuvir also protected infected animals from mortality, weight loss and liver injury, especially prophylatically. Our pre-clinical results supports a second use of sofosbuvir against yellow fever.

The Methods Employed in Mass Spectrometric Analysis of Posttranslational Modifications (PTMs) and Protein-Protein Interactions (PPIs)

Mass Spectrometry (MS) has revolutionized the way we study biomolecules, especially proteins, their interactions and posttranslational modifications (PTM). As such MS has established itself as the leading tool for the analysis of PTMs mainly because this approach is highly sensitive, amenable to high throughput and is capable of assigning PTMs to specific sites in the amino acid sequence of proteins and peptides. Along with the advances in MS methodology there have been improvements in biochemical, genetic and cell biological approaches to mapping the interactome which are discussed with consideration for both the practical and technical considerations of these techniques. The interactome of a species is generally understood to represent the sum of all potential protein-protein interactions. There are still a number of barriers to the elucidation of the human interactome or any other species as physical contact between protein pairs that occur by selective molecular docking in a particular spatiotemporal biological context are not easily captured and measured.PTMs massively increase the complexity of organismal proteomes and play a role in almost all aspects of cell biology, allowing for fine-tuning of protein structure, function and localization. There are an estimated 300 PTMS with a predicted 5% of the eukaryotic genome coding for enzymes involved in protein modification, however we have not yet been able to reliably map PTM proteomes due to limitations in sample preparation, analytical techniques, data analysis, and the substoichiometric and transient nature of some PTMs. Improvements in proteomic and mass spectrometry methods, as well as sample preparation, have been exploited in a large number of proteome-wide surveys of PTMs in many different organisms. Here we focus on previously published global PTM proteome studies in the Apicomplexan parasites T. gondii and P. falciparum which offer numerous insights into the abundance and function of each of the studied PTM in the Apicomplexa. Integration of these datasets provide a more complete picture of the relative importance of PTM and crosstalk between them and how together PTM globally change the cellular biology of the Apicomplexan protozoa. A multitude of techniques used to investigate PTMs, mostly techniques in MS-based proteomics, are discussed for their ability to uncover relevant biological function.

Homology-based loop modeling yields more complete crystallographic protein structures

Inherent protein flexibility, poor or low-resolution diffraction data or poorly defined electron-density maps often inhibit the building of complete structural models during X-ray structure determination. However, recent advances in crystallographic refinement and model building often allow completion of previously missing parts. This paper presents algorithms that identify regions missing in a certain model but present in homologous structures in the Protein Data Bank (PDB), and 'graft' these regions of interest. These new regions are refined and validated in a fully automated procedure. Including these developments in the PDB-REDO pipeline has enabled the building of 24 962 missing loops in the PDB. The models and the automated procedures are publicly available through the PDB-REDO databank and webserver. More complete protein structure models enable a higher quality public archive but also a better understanding of protein function, better comparison between homologous structures and more complete data mining in structural bioinformatics projects.

Validation of Structures in the Protein Data Bank

The Worldwide PDB recently launched a deposition, biocuration, and validation tool: OneDep. At various stages of OneDep data processing, validation reports for three-dimensional structures of biological macromolecules are produced. These reports are based on recommendations of expert task forces representing crystallography, nuclear magnetic resonance, and cryoelectron microscopy communities. The reports provide useful metrics with which depositors can evaluate the quality of the experimental data, the structural model, and the fit between them. The validation module is also available as a stand-alone web server and as a programmatically accessible web service. A growing number of journals require the official wwPDB validation reports (produced at biocuration) to accompany manuscripts describing macromolecular structures. Upon public release of the structure, the validation report becomes part of the public PDB archive. Geometric quality scores for proteins in the PDB archive have improved over the past decade.

Combining asymmetric 13C-labeling and isotopic filter/edit NOESY: a novel strategy for rapid and logical RNA resonance assignment

Although ∼98% of the human genomic output is transcribed as non-protein coding RNA, <2% of the protein data bank structures comprise RNA. This huge structural disparity stems from combined difficulties of crystallizing RNA for X-ray crystallography along with extensive chemical shift overlap and broadened linewidths associated with NMR of RNA. While half of the deposited RNA structures in the PDB were solved by NMR methods, the usefulness of NMR is still limited by the high cost of sample preparation and challenges of resonance assignment. Here we propose a novel strategy for resonance assignment that combines new strategic 13C labeling technologies with filter/edit type NOESY experiments to greatly reduce spectral complexity and crowding. This new strategy allowed us to assign important non-exchangeable resonances of proton and carbon (1', 2', 2, 5, 6 and 8) nuclei using only one sample and <24 h of NMR instrument time for a 27 nt model RNA. The method was further extended to assigning a 6 nt bulge from a 61 nt viral RNA element justifying its use for a wide range RNA chemical shift resonance assignment problems.

The clinically approved antiviral drug sofosbuvir inhibits Zika virus replication

Zika virus (ZIKV) is a member of the Flaviviridae family, along with other agents of clinical significance such as dengue (DENV) and hepatitis C (HCV) viruses. Since ZIKV causes neurological disorders during fetal development and in adulthood, antiviral drugs are necessary. Sofosbuvir is clinically approved for use against HCV and targets the protein that is most conserved among the members of the Flaviviridae family, the viral RNA polymerase. Indeed, we found that sofosbuvir inhibits ZIKV RNA polymerase, targeting conserved amino acid residues. Sofosbuvir inhibited ZIKV replication in different cellular systems, such as hepatoma (Huh-7) cells, neuroblastoma (SH-Sy5y) cells, neural stem cells (NSC) and brain organoids. In addition to the direct inhibition of the viral RNA polymerase, we observed that sofosbuvir also induced an increase in A-to-G mutations in the viral genome. Together, our data highlight a potential secondary use of sofosbuvir, an anti-HCV drug, against ZIKV.

Inference of Functionally-Relevant N-acetyltransferase Residues Based on Statistical Correlations

Over evolutionary time, members of a superfamily of homologous proteins sharing a common structural core diverge into subgroups filling various functional niches. At the sequence level, such divergence appears as correlations that arise from residue patterns distinct to each subgroup. Such a superfamily may be viewed as a population of sequences corresponding to a complex, high-dimensional probability distribution. Here we model this distribution as hierarchical interrelated hidden Markov models (hiHMMs), which describe these sequence correlations implicitly. By characterizing such correlations one may hope to obtain information regarding functionally-relevant properties that have thus far evaded detection. To do so, we infer a hiHMM distribution from sequence data using Bayes’ theorem and Markov chain Monte Carlo (MCMC) sampling, which is widely recognized as the most effective approach for characterizing a complex, high dimensional distribution. Other routines then map correlated residue patterns to available structures with a view to hypothesis generation. When applied to N-acetyltransferases, this reveals sequence and structural features indicative of functionally important, yet generally unknown biochemical properties. Even for sets of proteins for which nothing is known beyond unannotated sequences and structures, this can lead to helpful insights. We describe, for example, a putative coenzyme-A-induced-fit substrate binding mechanism mediated by arginine residue switching between salt bridge and π-π stacking interactions. A suite of programs implementing this approach is available (psed.igs.umaryland.edu).

Protein sequence data, when gathered in great quantity, contain important but implicit biological information manifest as statistical correlations. Here we describe an approach to access this information by comprehensively modeling and characterizing the distribution of sequences belonging to a major protein superfamily. This approach takes as input a large set of unaligned sequences belonging to the superfamily. By applying the minimum description length principle, it seeks the statistical model that best explains the sequences while avoiding over-fitting the data. It concurrently aligns the sequences and, to model evolutionary divergence, partitions them into subgroups that are hierarchically-arranged based upon correlated residue patterns. Auxiliary routines create PyMOL scripts to visualize the locations of correlated residues within available structures. Because these correlations likely arise from structural and biochemical constraints, they can help elucidate protein properties important for functional specificity. Comparing and contrasting sequence and structural features in this way may therefore suggest, in the light of published studies, plausible biological hypotheses for experimental investigation. We illustrate this approach with N-acetyltransferases.

Homology modeling and molecular dynamics simulation of the HIF2α degradation-related HIF2α-VHL complex

BACKGROUND

Hypoxia-inducible factor 2 alpha (HIF2α), prolyl hydroxylase domain protein 2 (PHD2), and the von Hippel Lindau tumor suppressor protein (pVHL) are three principal proteins in the oxygen-sensing pathway. Under normoxic conditions, a conserved proline in HIF2α is hydroxylated by PHD2 in an oxygen-dependent manner, and then pVHL binds and promotes the degradation of HIF2α. However, the crystal structure of the HIF2α-pVHL complex has not yet been established, and this has limited research on the interaction between HIF and pVHL. Here, we constructed a structural model of a 23-residue HIF2α peptide (528-550)-pVHL-ElonginB-ElonginC complex by using homology modeling and molecular dynamics simulations. We also applied these methods to HIF2α mutants (HYP531PRO, F540L, A530 V, A530T, and G537R) to reveal structural defects that explain how these mutations weaken the interaction with pVHL.

METHODS

Homology modeling and molecular dynamics simulations were used to construct a three-dimensional (3D) structural model of the HIF2α-VHL complex. Subsequently, MolProbity, an active validation tool, was used to analyze the reliability of the model. Molecular mechanics energies combined with the generalized Born and surface area continuum solvation (MM-GBSA) and solvated interaction energy (SIE) methods were used to calculate the binding free energy between HIF2a and pVHL, and the stability of the simulation system was evaluated by using root mean square deviation (RMSD) analysis. We also determined the secondary structure of the system by using the definition of secondary structure of proteins (DSSP) algorithm. Finally, we investigated the structural significance of specific point mutations known to have clinical implications.

RESULTS

We established a reliable structural model of the HIF2α-pVHL complex, which is similar to the crystal structure of HIF1α in 1LQB. Furthermore, we compared the structural model of the HIF2α-pVHL complex and the HIF2α (HYP531P, F540L, A530V, A530T, and G537R)-pVHL mutants on the basis of RMSD, DSSP, binding free energy, and hydrogen bonding. The experimental data indicate that the stability of the structural model of the HIF2α-pVHL complex is higher than that of the mutants, consistently with clinical observations.

CONCLUSIONS

The structural model of the HIF2α-pVHL complex presented in this study enhances understanding of how HIF2α is captured by pVHL. Moreover, the important contact amino acids that we identified may be useful in the development of drugs to treat HIF2a-related diseases.

Toward High-Throughput Predictive Modeling of Protein Binding/Unbinding Kinetics

One of the unaddressed challenges in drug discovery is that drug potency determined in vitro is not a reliable indicator of drug activity in vivo. Accumulated evidence suggests that in vivo activity is more strongly correlated with the binding/unbinding kinetics than the equilibrium thermodynamics of protein-ligand interactions (PLIs). However, existing experimental and computational techniques are insufficient in studying the molecular details of kinetics processes of PLIs on a large scale. Consequently, we not only have limited mechanistic understanding of the kinetic processes but also lack a practical platform for high-throughput screening and optimization of drug leads on the basis of their kinetic properties. For the first time, we address this unmet need by integrating coarse-grained normal mode analysis with multitarget machine learning (MTML). To test our method, HIV-1 protease is used as a model system. We find that computational models based on the residue normal mode directionality displacement of PLIs can not only recapitulate the results from all-atom molecular dynamics simulations but also predict protein-ligand binding/unbinding kinetics accurately. When this is combined with energetic features, the accuracy of combined kon and koff prediction reaches 74.35%. Furthermore, our integrated model provides us with new insights into the molecular determinants of the kinetics of PLIs. We propose that the coherent coupling of conformational dynamics and thermodynamic interactions between the receptor and the ligand may play a critical role in determining the kinetic rate constants of PLIs. In conclusion, we demonstrate that residue normal mode directionality displacement can serve as a kinetic fingerprint to capture long-time-scale conformational dynamics of the binding/unbinding kinetics. When this is coupled with MTML, it is possible to screen and optimize compounds on the basis of their binding/unbinding kinetics in a high-throughput fashion. The further development of such computational tools will bridge one of the critical missing links between in vitro compound screening and in vivo drug activity.

Combined use of optical spectroscopy and computational methods to study the binding and the photoinduced conformational modification of proteins when NMR and X-ray structural determinations are not an option

The functions of proteins depend on their interactions with various ligands and these interactions are controlled by the structure of the polypeptides. If one can manipulate the structure of proteins, their functions can in principle be modulated. The issue of protein structure-function relationship is not only a central problem in biophysics, but is becoming clear that the ability to "artificially" modify the structure of proteins could be relevant in fields beyond the biomedical area to provide, for instance, light responses in proteins which would not possess such properties in their native state. This chapter presents an overview of the combination of optical electronic and vibrational spectroscopy with various computational methods to investigate the binding between photoactive ligands and proteins.

Structure and stability of recombinant bovine odorant-binding protein: III. Peculiarities of the wild type bOBP unfolding in crowded milieu

Contrary to the majority of the members of the lipocalin family, which are stable monomers with the specific OBP fold (a β-barrel consisting of a 8-stranded anti-parallel β-sheet followed by a short α-helical segment, a ninth β-strand, and a disordered C-terminal tail) and a conserved disulfide bond, bovine odorant-binding protein (bOBP) does not have such a disulfide bond and forms a domain-swapped dimer that involves crossing the α-helical region from each monomer over the β-barrel of the other monomer. Furthermore, although natural bOBP isolated from bovine tissues exists as a stable domain-swapped dimer, recombinant bOBP has decreased dimerization potential and therefore exists as a mixture of monomeric and dimeric variants. In this article, we investigated the effect model crowding agents of similar chemical nature but different molecular mass on conformational stability of the recombinant bOBP. These experiments were conducted in order to shed light on the potential influence of model crowded environment on the unfolding-refolding equilibrium. To this end, we looked at the influence of PEG-600, PEG-4000, and PEG-12000 in concentrations of 80, 150, and 300 mg/mL on the equilibrium unfolding and refolding transitions induced in the recombinant bOBP by guanidine hydrochloride. We are showing here that the effect of crowding agents on the structure and conformational stability of the recombinant bOBP depends on the size of the crowder, with the smaller crowding agents being more effective in the stabilization of the bOBP native dimeric state against the guanidine hydrochloride denaturing action. This effect of the crowding agents is concentration dependent, with the high concentrations of the agents being more effective.

Structure and stability of recombinant bovine odorant-binding protein: II. Unfolding of the monomeric forms

In a family of monomeric odorant-binding proteins (OBPs), bovine OBP (bOBP), that lacks conserved disulfide bond found in other OBPs, occupies unique niche because of its ability to form domain-swapped dimers. In this study, we analyzed conformational stabilities of the recombinant bOBP and its monomeric variants, the bOBP-Gly121+ mutant containing an additional glycine residue after the residue 121 of the bOBP, and the GCC-bOBP mutant obtained from the bOBP-Gly121+ form by introduction of the Trp64Cys/His155Cys double mutation to restore the canonical disulfide bond. We also analyzed the effect of the natural ligand binding on the conformational stabilities of these bOBP variants. Our data are consistent with the conclusion that the unfolding-refolding pathways of the recombinant bOBP and its mutant monomeric forms bOBP-Gly121+ and GCC-bOBP are similar and do not depend on the oligomeric status of the protein. This clearly shows that the information on the unfolding-refolding mechanism is encoded in the structure of the bOBP monomers. However, the process of the bOBP unfolding is significantly complicated by the formation of the domain-swapped dimer, and the rates of the unfolding-refolding reactions essentially depend on the conditions in which the protein is located.

Structure and stability of recombinant bovine odorant-binding protein: I. Design and analysis of monomeric mutants

Bovine odorant-binding protein (bOBP) differs from other lipocalins by lacking the conserved disulfide bond and for being able to form the domain-swapped dimers. To identify structural features responsible for the formation of the bOBP unique dimeric structure and to understand the role of the domain swapping on maintaining the native structure of the protein, structural properties of the recombinant wild type bOBP and its mutant that cannot dimerize via the domain swapping were analyzed. We also looked at the effect of the disulfide bond by designing a monomeric bOBPs with restored disulfide bond which is conserved in other lipocalins. Finally, to understand which features in the microenvironment of the bOBP tryptophan residues play a role in the defining peculiarities of the intrinsic fluorescence of this protein we designed and investigated single-tryptophan mutants of the monomeric bOBP. Our analysis revealed that the insertion of the glycine after the residue 121 of the bOBP prevents domain swapping and generates a stable monomeric protein bOBP-Gly121+. We also show that the restored disulfide bond in the GCC-bOBP mutant leads to the noticeable stabilization of the monomeric structure. Structural and functional analysis revealed that none of the amino acid substitutions introduced to the bOBP affected functional activity of the protein and that the ligand binding leads to the formation of a more compact and stable state of the recombinant bOBP and its mutant monomeric forms. Finally, analysis of the single-tryptophan mutants of the monomeric bOBP gave us a unique possibility to find peculiarities of the microenvironment of tryptophan residues which were not previously described.

Modulation of Erythrocyte Plasma Membrane Redox System Activity by Curcumin

Plasma membrane redox system (PMRS) is an electron transport chain system ubiquitously present throughout all cell types. It transfers electron from intracellular substrates to extracellular acceptors for regulation of redox status. Curcumin, isolated from Curcuma longa, has modulatory effects on cellular physiology due to its membrane interaction ability and antioxidant potential. The present study investigates the effect of curcumin on PMRS activity of erythrocytes isolated from Wistar rats in vitro and in vivo and validated through an in silico docking simulation study using Molegro Virtual Docker (MVD). Effects of curcumin were also evaluated on level of glutathione (GSH) and the oxidant potential of plasma measured in terms of plasma ferric equivalent oxidative potentials (PFEOP). Results show that curcumin significantly (p < 0.01) downregulated the PMRS activity in a dose-dependent manner. Molecular docking results suggest that curcumin interacts with amino acids at the active site cavity of cytochrome b 5 reductase, a key constituent of PMRS. Curcumin also increased the GSH level in erythrocytes and plasma while simultaneously decreasing the oxidant potential (PFEOP) of plasma. Altered PMRS activity and redox status are associated with the pathophysiology of several health complications including aging and diabetes; hence, the above finding may explain part of the role of curcumin in health beneficial effects.

Allosteric effects of chromophore interaction with dimeric near-infrared fluorescent proteins engineered from bacterial phytochromes

Fluorescent proteins (FPs) engineered from bacterial phytochromes attract attention as probes for in vivo imaging due to their near-infrared (NIR) spectra and use of available in mammalian cells biliverdin (BV) as chromophore. We studied spectral properties of the iRFP670, iRFP682 and iRFP713 proteins and their mutants having Cys residues able to bind BV either in both PAS (Cys15) and GAF (Cys256) domains, in one of these domains, or without these Cys residues. We show that the absorption and fluorescence spectra and the chromophore binding depend on the location of the Cys residues. Compared with NIR FPs in which BV covalently binds to Cys15 or those that incorporate BV noncovalently, the proteins with BV covalently bound to Cys256 have blue-shifted spectra and higher quantum yield. In dimeric NIR FPs without Cys15, the covalent binding of BV to Сys256 in one monomer allosterically inhibits the covalent binding of BV to the other monomer, whereas the presence of Cys15 allosterically promotes BV binding to Cys256 in both monomers. The NIR FPs with both Cys residues have the narrowest blue-shifted spectra and the highest quantum yield. Our analysis resulted in the iRFP713/Val256Cys protein with the highest brightness in mammalian cells among available NIR FPs.

Data Mining of Macromolecular Structures

The use of macromolecular structures is widespread for a variety of applications, from teaching protein structure principles all the way to ligand optimization in drug development. Applying data mining techniques on these experimentally determined structures requires a highly uniform, standardized structural data source. The Protein Data Bank (PDB) has evolved over the years toward becoming the standard resource for macromolecular structures. However, the process selecting the data most suitable for specific applications is still very much based on personal preferences and understanding of the experimental techniques used to obtain these models. In this chapter, we will first explain the challenges with data standardization, annotation, and uniformity in the PDB entries determined by X-ray crystallography. We then discuss the specific effect that crystallographic data quality and model optimization methods have on structural models and how validation tools can be used to make informed choices. We also discuss specific advantages of using the PDB_REDO databank as a resource for structural data. Finally, we will provide guidelines on how to select the most suitable protein structure models for detailed analysis and how to select a set of structure models suitable for data mining.

Optimal ligand descriptor for pocket recognition based on the Beta-shape

Structure-based virtual screening is one of the most important and common computational methods for the identification of predicted hit at the beginning of drug discovery. Pocket recognition and definition is frequently a prerequisite of structure-based virtual screening, reducing the search space of the predicted protein-ligand complex. In this paper, we present an optimal ligand shape descriptor for a pocket recognition algorithm based on the beta-shape, which is a derivative structure of the Voronoi diagram of atoms. We investigate six candidates for a shape descriptor for a ligand using statistical analysis: the minimum enclosing sphere, three measures from the principal component analysis of atoms, the van der Waals volume, and the beta-shape volume. Among them, the van der Waals volume of a ligand is the optimal shape descriptor for pocket recognition and best tunes the pocket recognition algorithm based on the beta-shape for efficient virtual screening. The performance of the proposed algorithm is verified by a benchmark test.

Structural Characterization of Peptide Antibodies

The role of proteins as very effective immunogens for the generation of antibodies is indisputable. Nevertheless, cases in which protein usage for antibody production is not feasible or convenient compelled the creation of a powerful alternative consisting of synthetic peptides. Synthetic peptides can be modified to obtain desired properties or conformation, tagged for purification, isotopically labeled for protein quantitation or conjugated to immunogens for antibody production. The antibodies that bind to these peptides represent an invaluable tool for biological research and discovery. To better understand the underlying mechanisms of antibody-antigen interaction here we present a pipeline developed by us to structurally classify immunoglobulin antigen binding sites and to infer key sequence residues and other variables that have a prominent role in each structural class.

Tryptophan residue of the D-galactose/D-glucose-binding protein from E. Coli localized in its active center does not contribute to the change in intrinsic fluorescence upon glucose binding

Changes of the characteristics of intrinsic tryptophan fluorescence of the wild type of D-galactose/D-glucose-binding protein from Escherichia coli (GGBPwt) induced by D-glucose binding were examined by the intrinsic UV-fluorescence of proteins, circular dyhroism in the near-UV region, and acrylamide-induced fluorescence quenching. The analysis of the different characteristics of GGBPwt and its mutant form GGBP-W183A together with the analysis of the microenvironment of tryptophan residues of GGBPwt revealed that Trp 183, which is directly involved in sugar binding, has the least influence on the provoked by D-glucose blue shift and increase in the intensity of protein intrinsic fluorescence in comparison with other tryptophan residues of GGBP.

Small molecule annotation for the Protein Data Bank

The Protein Data Bank (PDB) is the single global repository for three-dimensional structures of biological macromolecules and their complexes, and its more than 100 000 structures contain more than 20 000 distinct ligands or small molecules bound to proteins and nucleic acids. Information about these small molecules and their interactions with proteins and nucleic acids is crucial for our understanding of biochemical processes and vital for structure-based drug design. Small molecules present in a deposited structure may be attached to a polymer or may occur as a separate, non-covalently linked ligand. During curation of a newly deposited structure by wwPDB annotation staff, each molecule is cross-referenced to the PDB Chemical Component Dictionary (CCD). If the molecule is new to the PDB, a dictionary description is created for it. The information about all small molecule components found in the PDB is distributed via the ftp archive as an external reference file. Small molecule annotation in the PDB also includes information about ligand-binding sites and about covalent and other linkages between ligands and macromolecules. During the remediation of the peptide-like antibiotics and inhibitors present in the PDB archive in 2011, it became clear that additional annotation was required for consistent representation of these molecules, which are quite often composed of several sequential subcomponents including modified amino acids and other chemical groups. The connectivity information of the modified amino acids is necessary for correct representation of these biologically interesting molecules. The combined information is made available via a new resource called the Biologically Interesting molecules Reference Dictionary, which is complementary to the CCD and is now routinely used for annotation of peptide-like antibiotics and inhibitors.

A drunken search in crystallization space

The REMARK280 field of the Protein Data Bank is the richest open source of successful crystallization information. The REMARK280 field is optional and currently uncurated, so significant effort needs to be applied to extract reliable data. There are well over 15 000 crystallization conditions available commercially from 12 different vendors. After putting the PDB crystallization information and the commercial cocktail data into a consistent format, these data are used to extract information about the overlap between the two sets of crystallization conditions. An estimation is made as to which commercially available conditions are most appropriate for producing well diffracting crystals by looking at which commercial conditions are found unchanged (or almost unchanged) in the PDB. Further analyses include which commercial kits are the most appropriate for shotgun or more traditional approaches to crystallization screening. This analysis suggests that almost 40% of the crystallization conditions found currently in the PDB are identical or very similar to a commercial condition.

Structural analysis of protein-ligand interactions: the binding of endogenous compounds and of synthetic drugs

The large number of macromolecular structures deposited with the Protein Data Bank (PDB) describing complexes between proteins and either physiological compounds or synthetic drugs made it possible a systematic analysis of the interactions occurring between proteins and their ligands. In this work, the binding pockets of about 4000 PDB protein-ligand complexes were investigated and amino acid and interaction types were analyzed. The residues observed with lowest frequency in protein sequences, Trp, His, Met, Tyr, and Phe, turned out to be the most abundant in binding pockets. Significant differences between drug-like and physiological compounds were found. On average, physiological compounds establish with respect to drugs about twice as many hydrogen bonds with protein atoms, whereas drugs rely more on hydrophobic interactions to establish target selectivity. The large number of PDB structures describing homologous proteins in complex with the same ligand made it possible to analyze the conservation of binding pocket residues among homologous protein structures bound to the same ligand, showing that Gly, Glu, Arg, Asp, His, and Thr are more conserved than other amino acids. Also in the cases in which the same ligand is bound to unrelated proteins, the binding pockets showed significant conservation in the residue types. In this case, the probability of co-occurrence of the same amino acid type in the binding pockets could be up to thirteen times higher than that expected on a random basis. The trends identified in this study may provide an useful guideline in the process of drug design and lead optimization.

Actinous enigma or enigmatic actin: Folding, structure, and functions of the most abundant eukaryotic protein

Being the most abundant protein of the eukaryotic cell, actin continues to keep its secrets for more than 60 years. Everything about this protein, its structure, functions, and folding, is mysteriously counterintuitive, and this review represents an attempt to solve some of the riddles and conundrums commonly found in the field of actin research. In fact, actin is a promiscuous binder with a wide spectrum of biological activities. It can exist in at least three structural forms, globular, fibrillar, and inactive (G-, F-, and I-actin, respectively). G-actin represents a thermodynamically instable, quasi-stationary state, which is formed in vivo as a result of the energy-intensive, complex posttranslational folding events controlled and driven by cellular folding machinery. The G-actin structure is dependent on the ATP and Mg²⁺ binding (which in vitro is typically substituted by Ca²⁺) and protein is easily converted to the I-actin by the removal of metal ions and by action of various denaturing agents (pH, temperature, and chemical denaturants). I-actin cannot be converted back to the G-form. Foldable and “natively folded” forms of actin are always involved in interactions either with the specific protein partners, such as Hsp70 chaperone, prefoldin, and the CCT chaperonin during the actin folding in vivo or with Mg²⁺ and ATP as it takes place in the G-form. We emphasize that the solutions for the mysteries of actin multifunctionality, multistructurality, and trapped unfolding can be found in the quasi-stationary nature of this enigmatic protein, which clearly possesses many features attributed to both globular and intrinsically disordered proteins.

Structure-based characterization of multiprotein complexes

Multiprotein complexes govern virtually all cellular processes. Their 3D structures provide important clues to their biological roles, especially through structural correlations among protein molecules and complexes. The detection of such correlations generally requires comprehensive searches in databases of known protein structures by means of appropriate structure-matching techniques. Here, we present a high-speed structure search engine capable of instantly matching large protein oligomers against the complete and up-to-date database of biologically functional assemblies of protein molecules. We use this tool to reveal unseen structural correlations on the level of protein quaternary structure and demonstrate its general usefulness for efficiently exploring complex structural relationships among known protein assemblies.

Coverage of protein domain families with structural protein-protein interactions: current progress and future trends

Protein interactions have evolved into highly precise and regulated networks adding an immense layer of complexity to cellular systems. The most accurate atomistic description of protein binding sites can be obtained directly from structures of protein complexes. The availability of structurally characterized protein interfaces significantly improves our understanding of interactomes, and the progress in structural characterization of protein-protein interactions (PPIs) can be measured by calculating the structural coverage of protein domain families. We analyze the coverage of protein domain families (defined according to CDD and Pfam databases) by structures, structural protein-protein complexes and unique protein binding sites. Structural PPI coverage of currently available protein families is about 30% without any signs of saturation in coverage growth dynamics. Given the current growth rates of domain databases and structural PPI deposition, complete domain coverage with PPIs is not expected in the near future. As a result of this study we identify families without any protein-protein interaction evidence (listed on a supporting website http://www.ncbi.nlm.nih.gov/Structure/ibis/coverage/) and propose them as potential targets for structural studies with a focus on protein interactions.

A knot in the protein structure - probing the near-infrared fluorescent protein iRFP designed from a bacterial phytochrome

The possibility of engineering near-infrared fluorescent proteins and biosensors from bacterial phytochrome photoreceptors (BphPs) has led to substantial interest in this family of proteins. The near-infrared fluorescent proteins have allowed non-invasive bio-imaging of deep tissues and whole organs in living animals. BphPs and derived near-infrared fluorescent proteins contain a structural element, called a knot, in their polypeptide chains. The formation of knot structures in proteins was refuted for a long time. Here, we studied the denaturation and renaturation processes of the near-infrared fluorescent probe iRFP, engineered from RpBphP2, which utilizes a heme-derived tetrapyrrole compound biliverdin as a chromophore. iRFP contains a unique figure-of-eight knot. The denaturation and renaturation curves of the iRFP apoform coincided well, suggesting efficient refolding. However, the iRFP holoform exhibited irreversible unfolding and aggregation associated with the bound chromophore. The knot structure in the apoform did not prevent subsequent binding of biliverdin, resulting in the functional iRFP holoform. We suggest that the irreversibility of protein unfolding is caused by post-translational protein modifications, such as chromophore binding, rather than the presence of the knot. These results are essential for future design of BphP-based near-infrared probes, and add important features to our knowledge of protein folding.

The quaternary structure of the recombinant bovine odorant-binding protein is modulated by chemical denaturants

A large group of odorant-binding proteins (OBPs) has attracted great scientific interest as promising building blocks in constructing optical biosensors for dangerous substances, such as toxic and explosive molecules. Native tissue-extracted bovine OBP (bOBP) has a unique dimer folding pattern that involves crossing the α-helical domain in each monomer over the other monomer’s β-barrel. In contrast, recombinant bOBP maintaining the high level of stability inherent to native tissue bOBP is produced in a stable native-like state with a decreased tendency for dimerization and is a mixture of monomers and dimers in a buffered solution. This work is focused on the study of the quaternary structure and the folding-unfolding processes of the recombinant bOBP in the absence and in the presence of guanidine hydrochloride (GdnHCl). Our results show that the recombinant bOBP native dimer is only formed at elevated GdnHCl concentrations (1.5 M). This process requires re-organizing the protein structure by progressing through the formation of an intermediate state. The bOBP dimerization process appears to be irreversible and it occurs before the protein unfolds. Though the observed structural changes for recombinant bOBP at pre-denaturing GdnHCl concentrations show a local character and the overall protein structure is maintained, such changes should be considered where the protein is used as a sensitive element in a biosensor system.

Protective effect of theaflavin on erythrocytes subjected to in vitro oxidative stress

Antioxidant and free radical scavenging effect of black tea theaflavins has been shown in many epidemiological studies. In the present work we report the protective mechanism of tea theaflavins on biomarkers of oxidative stress, which are elevated during stress conditions. We hereby report the in vitro effect of theaflavins on erythrocyte malondialdehyde (MDA), intracellular reduced glutathione (GSH), and plasma membrane redox system (PMRS) of rats. The effect of theaflavin on PMRS has also been validated through an in silico docking simulation study using Molegro Virtual Docker (MVD). We report that theaflavins show significant protection to erythrocyte against oxidative stress induced by tert-butyl hydroperoxide (t-BHP). The findings suggest a possible protective role of theaflavins as antioxidant.

Targeting UDP-galactopyranose mutases from eukaryotic human pathogens

UDP-Galactopyranose mutase (UGM) is a unique flavin-dependent enzyme that catalyzes the conversion of UDP-galactopyranose(UDP-Galp) to UDP-galactofuranose (UDP-Galf). The product of this reaction is the precursor to Galf, a major component of the cell wall and of cell surface glycoproteins and glycolipids in many eukaryotic and prokaryotic human pathogens. The function of UGM is important in the virulence of fungi, parasites, and bacteria. Its role in virulence and its absence in humans suggest that UGM is an ideal drug target. Significant structural and mechanistic information has been accumulated on the prokaryotic UGMs; however, in the past few years the research interest has shifted to UGMs from eukaryotic human pathogens such as fungi and protozoan parasites. It has become clear that UGMs from prokaryotic and eukaryotic organisms have different structural and mechanistic features. The amino acid sequence identity between these two classes of enzymes is low, resulting in differences in oligomeric states, substrate binding, active site flexibility, and interaction with redox partners. However, the unique role of the flavin cofactor in catalysis is conserved among this enzyme family. In this review, recent findings on eukaryotic UGMs are discussed and presented in comparison with prokaryotic UGMs.

CLIPS-4D: a classifier that distinguishes structurally and functionally important residue-positions based on sequence and 3D data

MOTIVATION

The precise identification of functionally and structurally important residues of a protein is still an open problem, and state-of-the-art classifiers predict only one or at most two different categories.

RESULT

We have implemented the classifier CLIPS-4D, which predicts in a mutually exclusively manner a role in catalysis, ligand-binding or protein stability for each residue-position of a protein. Each prediction is assigned a P-value, which enables the statistical assessment and the selection of predictions with similar quality. CLIPS-4D requires as input a multiple sequence alignment and a 3D structure of one protein in PDB format. A comparison with existing methods confirmed state-of-the-art prediction quality, even though CLIPS-4D classifies more specifically than other methods. CLIPS-4D was implemented as a multiclass support vector machine, which exploits seven sequence-based and two structure-based features, each of which was shown to contribute to classification quality. The classification of ligand-binding sites profited most from the 3D features, which were the assessment of the solvent accessible surface area and the identification of surface pockets. In contrast, five additionally tested 3D features did not increase the classification performance achieved with evolutionary signals deduced from the multiple sequence alignment.

ValiDichro: a website for validating and quality control of protein circular dichroism spectra

Circular dichroism (CD) spectroscopy is widely used in structural biology as a technique for examining the structure, folding and conformational changes of proteins. A new server, ValiDichro, has been developed for checking the quality and validity of CD spectral data and metadata, both as an aid to data collection and processing and as a validation procedure for spectra to be included in publications. ValiDichro currently includes 25 tests for data completeness, consistency and quality. For each test that is done, not only is a validation report produced, but the user is also provided with suggestions for correcting or improving the data. The ValiDichro server is freely available at http://valispec.cryst.bbk.ac.uk/circularDichroism/ValiDichro/upload.html.

Beta-barrel scaffold of fluorescent proteins: folding, stability and role in chromophore formation

This review focuses on the current view of the interaction between the β-barrel scaffold of fluorescent proteins and their unique chromophore located in the internal helix. The chromophore originates from the polypeptide chain and its properties are influenced by the surrounding protein matrix of the β-barrel. On the other hand, it appears that a chromophore tightens the β-barrel scaffold and plays a crucial role in its stability. Furthermore, the presence of a mature chromophore causes hysteresis of protein unfolding and refolding. We survey studies measuring protein unfolding and refolding using traditional methods as well as new approaches, such as mechanical unfolding and reassembly of truncated fluorescent proteins. We also analyze models of fluorescent protein unfolding and refolding obtained through different approaches, and compare the results of protein folding in vitro to co-translational folding of a newly synthesized polypeptide chain.

PLI: a web-based tool for the comparison of protein-ligand interactions observed on PDB structures

MOTIVATION

A large fraction of the entries contained in the Protein Data Bank describe proteins in complex with low molecular weight molecules such as physiological compounds or synthetic drugs. In many cases, the same molecule is found in distinct protein-ligand complexes. There is an increasing interest in Medicinal Chemistry in comparing protein binding sites to get insight on interactions that modulate the binding specificity, as this structural information can be correlated with other experimental data of biochemical or physiological nature and may help in rational drug design.

RESULTS

The web service protein-ligand interaction presented here provides a tool to analyse and compare the binding pockets of homologous proteins in complex with a selected ligand. The information is deduced from protein-ligand complexes present in the Protein Data Bank and stored in the underlying database.

AVAILABILITY

Freely accessible at http://bioinformatics.istge.it/pli/.

Distinct effects of guanidine thiocyanate on the structure of superfolder GFP

Having a high folding efficiency and a low tendency to aggregate, the superfolder GFP (sfGFP) offers a unique opportunity to study the folding of proteins that have a β-barrel topology. Here, we studied the unfolding–refolding of sfGFP that was induced by guanidine thiocyanate (GTC), which is a stronger denaturing agent than GdnHCl or urea. Structural perturbations of sfGFP were studied by spectroscopic methods (absorbance, fluorescence, and circular dichroism), by acrylamide quenching of tryptophan and green chromophore fluorescence, and by size-exclusion chromatography. Low concentrations of GTC (up to 0.1 M) induce subtle changes in the sfGFP structure. The pronounced changes in the visible absorption spectrum of sfGFP which are accompanied by a dramatic decrease in tryptophan and green chromophore fluorescence was recorded in the range 0–0.7 M GNC. These alterations of sfGFP characteristics that erroneously can be mixed up with appearance of intermediate state in fact have pure spectroscopic but not structural nature. Higher concentrations of GTC (from 0.7 to 1.7 M), induce a disruption of the sfGFP structure, that is manifested in simultaneous changes of all of the detected parameters. Full recovery of native properties of denaturated sfGFP was observed after denaturant removal. The refolding of sfGFP passes through the accumulation of the off-pathway intermediate state, in which inner alpha-helix and hence green chromophore and Trp57 are still not tuned up to (properly integrated into) the already formed β-barrel scaffold of protein. Incorporation of the chromophore in the β-barrel in the pathway of refolding and restoration of its ability to fluoresce occur in a narrow range of GTC concentrations from 1.0 to 0.7 M, and a correct insertion of Trp 57 occurs at concentrations ranging from 0.7 to 0 M GTC. These two processes determine the hysteresis of protein unfolding and refolding.

Proteomics approaches for the analysis of enriched microbial subpopulations and visualization of complex functional information

Advances in the separation of microbial subpopulations and in proteomics technologies have paved the way for the global molecular characterization of microbial cells that share common functional characteristics. Quantitative characterization of the dynamics of microbial proteomes enables an unprecedented view of the adaptive responses of microbes to environmental stimuli or during interaction with other species or host cells. However, the intrinsic complexity of such data requires sophisticated visualization methods for the display, mining, interpretation and efficient exploitation of these data resources. In this review, we discuss how new approaches in data visualization such as streamgraphs, network graphs or Voronoi treemaps are being used in the field to provide new insights into the functional complexity of microbial cells, populations and multispecies consortia.

Plant polyphenols as electron donors for erythrocyte plasma membrane redox system: validation through in silico approach

Background

The plasma membrane redox system (PMRS) has extensively been studied in erythrocytes. The PMRS plays an important role in maintaining plasma redox balance and provides a protective mechanism against oxidative stress. Earlier it was proposed that only NADH or NADPH provided reducing equivalents to PMRS; however, now it is acknowledged that some polyphenols also have the ability to donate reducing equivalents to PMRS.

Methods

Two different docking simulation softwares, Molegro Virtual Docker and Glide were used to study the interaction of certain plant polyphenols viz. quercetin, epigallocatechin gallate, catechin epicatechin and resveratrol with human erythroyte NADH-cytochrome b5 reductase, which is a component of PMRS and together with the identification of minimum pharmacophoric feature using Pharmagist.

Results

The derived common minimum pharmacophoric features show the presence of minimum bioactive component in all the selected polyphenols. Our results confirm wet lab findings which show that these polyphenols have the ability to interact and donate protons to the Human NADH-cytochrome b5 reductase.

Conclusion

With the help of these comparative results of docking simulation and pharmacophoric features, novel potent molecules can be designed with higher efficacy for activation of the PMRS system.