Antifouling activities of proteinase K and α-amylase enzymes: Laboratory bioassays and in silico analysis

The application of enzymes as antifoulants is one of the environment-friendly strategies in biofouling management. In this study, antifouling activities of commercially available proteinase K and α-amylase enzymes were evaluated using barnacle larva and biofilm-forming bacteria as test organisms. The enzymes were also tested against barnacle cement protein through in silico analysis. The results showed that both enzymes inhibited the attachment of bacteria and settlement of barnacle larvae on the test surface. The lowest minimum inhibitory concentration of 0.312 mg ml-1 was exhibited by proteinase K against biofilm-forming bacteria. The calculated LC50 values for proteinase K and α-amylase against the barnacle nauplii were 91.8 and 230.96 mg ml-1 respectively. While α-amylase showed higher antibiofilm activity, proteinase K exhibited higher anti-larval settlement activity. Similarly, in silico analysis of the enzymes revealed promising anti-settlement activity, as the enzymes showed good binding scores with barnacle cement protein. Overall, the results suggested that the enzymes proteinase K and α-amylase could be used in antifouling coatings to reduce the settlement of biofouling on artificial materials in the marine environment.

Interfacial residues in protein-protein complexes are in the eyes of the beholder

Interactions between proteins are vital in almost all biological processes. The characterization of protein-protein interactions helps us understand the mechanistic basis of biological processes, thereby enabling the manipulation of proteins for biotechnological and clinical purposes. The interface residues of a protein-protein complex are assumed to have the following two properties: (a) they always interact with a residue of a partner protein, which forms the basis for distance-based interface residue identification methods, and (b) they are solvent-exposed in the isolated form of the protein and become buried in the complex form, which forms the basis for Accessible Surface Area (ASA)-based methods. The study interrogates this popular assumption by recognizing interface residues in protein-protein complexes through these two methods. The results show that a few residues are identified uniquely by each method, and the extent of conservation, propensities, and their contribution to the stability of protein-protein interaction varies substantially between these residues. The case study analyses showed that interface residues, unique to distance, participate in crucial interactions that hold the proteins together, whereas the interface residues unique to the ASA method have a potential role in the recognition, dynamics, and specificity of the complex and can also be a hotspot. Overall, the study recommends applying both distance and ASA methods so that some interface residues missed by either method but crucial to the stability, recognition, dynamics, and function of protein-protein complexes are identified in a complementary manner.

'Computational studies on coumestrol-ArlR interaction to target ArlRS signaling cascade involved in MRSA virulence'

Two component signaling system ArlRS (Autolysis-related locus) regulates adhesion, biofilm formation and virulence in methicillin resistant Staphylococcus aureus. It consists of a histidine kinase ArlS and response regulator ArlR. ArlR is composed of a N-terminal receiver domain and DNA-binding effector domain at C-terminal. ArlR receiver domain dimerizes upon signal recognition and activates DNA binding by effector domain and subsequent virulence expression. In silico simulation and structural data suggest that coumestrol, a phytochemical found in Pueraria montana, forges a strong intermolecular interaction with residues involved in dimer formation and destabilizes ArlR dimerization, an essential conformational switch required for downstream effector domain to bind to virulent loci. Structural and energy profiles of simulated ArlR-coumestrol complexes suggest lower affinity between ArlR monomers due to structural rigidity at the dimer interface hindering the conformational rearrangements relevant for dimer formation. These analyses could be an attractive strategy to develop therapeutics and potent leads molecules response regulators of two component systems in which are involved in MRSA virulence as well as other drug-resistant pathogens.Communicated by Ramaswamy H. Sarma.

Leucine Motifs Stabilize Residual Helical Structure in Disordered Proteins

Many examples are known of regions of intrinsically disordered proteins that fold into α-helices upon binding to their targets. These helical binding motifs (HBMs) can be partially helical also in the unbound state, and this so-called residual structure can affect binding affinity and kinetics. To investigate the underlying mechanisms governing the formation of residual helical structure, we assembled a dataset of experimental helix contents of 65 peptides containing HBM that fold-upon-binding. The average residual helicity is 17% and increases to 60% upon target binding. The helix contents of residual and target-bound structures do not correlate, however the relative location of helix elements in both states shows a strong overlap. Compared to the general disordered regions, HBMs are enriched in amino acids with high helix preference and these residues are typically involved in target binding, explaining the overlap in helix positions. In particular, we find that leucine residues and leucine motifs in HBMs are the major contributors to helix stabilization and target-binding. For the two model peptides, we show that substitution of leucine motifs to other hydrophobic residues (valine or isoleucine) leads to reduction of residual helicity, supporting the role of leucine as helix stabilizer. From the three hydrophobic residues only leucine can efficiently stabilize residual helical structure. We suggest that the high occurrence of leucine motifs and a general preference for leucine at binding interfaces in HBMs can be explained by its unique ability to stabilize helical elements.

Distinct or Overlapping Areas of Mitochondrial Thioredoxin 2 May Be Used for Its Covalent and Strong Non-Covalent Interactions with Protein Ligands

In silico approaches were employed to examine the characteristics of interactions between human mitochondrial thioredoxin 2 (HsTrx2) and its 38 previously identified mitochondrial protein ligands. All interactions appeared driven mainly by electrostatic forces. The statistically significant residues of HsTrx2 for interactions were characterized as "contact hot spots". Since these were identical/adjacent to putative thermodynamic hot spots, an energy network approach identified their neighbors to highlight possible contact interfaces. Three distinct areas for binding emerged: (i) one around the active site for covalent interactions, (ii) another antipodal to the active site for strong non-covalent interactions, and (iii) a third area involved in both kinds of interactions. The contact interfaces of HsTrx2 were projected as respective interfaces for Escherichia coli Trx1 (EcoTrx1), 2, and HsTrx1. Comparison of the interfaces and contact hot spots of HsTrx2 to the contact residues of EcoTx1 and HsTrx1 from existing crystal complexes with protein ligands supported the hypothesis, except for a part of the cleft/groove adjacent to Trp30 preceding the active site. The outcomes of this study raise the possibility for the rational design of selective inhibitors for the interactions of HsTrx2 with specific protein ligands without affecting the entirety of the functions of the Trx system.

Modeling and monitoring the effects of three partly conserved Ile residues in the dimerization domain of a Mip-like virulence factor from Escherichia coli

FKBP22, an Escherichia coli-made peptidyl-prolyl cis-trans isomerase, has shown considerable homology with Mip-like virulence factors. While the C-terminal domain of this enzyme is used for executing catalytic function and binding inhibitor, the N-terminal domain is employed for its dimerization. To precisely determine the underlying factors of FKBP22 dimerization, its structural model, developed using a suitable template, was carefully inspected. The data show that the dimeric FKBP22, like dimeric Mip proteins, has a V-like shape. Further, it dimerizes using 40 amino acid residues including Ile 9, Ile 17, Ile 42, and Ile 65. All of the above Ile residues except Ile 9 are partly conserved in the Mip-like proteins. To confirm the roles of the partly conserved Ile residues, three FKBP22 mutants, constructed by substituting them with an Ala residue, were studied as well. The results together indicate that Ile 65 has little role in maintaining the dimeric state or enzymatic activity of FKBP22. Conversely, both Ile 17 and Ile 42 are essential for preserving the structure, enzymatic activity, and dimerization ability of FKBP22. Ile 42 in particular looks more essential to FKBP22. However, none of these two Ile residues is required for binding the cognate inhibitor. Additional computational studies also indicated the change of V-shape and the dimeric state of FKBP22 due to the Ala substitution at position 42. The ways Ile 17 and Ile 42 protect the structure, function, and dimerization of FKBP22 have been discussed at length.Communicated by Ramaswamy H. Sarma.

Computational designing of the ligands of Protein L affinity chromatography based on molecular docking and molecular dynamics simulations

Protein L is a multidomain protein from Peptostreptococcus magnus with binding affinity to kappa light chain of human immunoglobulin (Ig) which is used for the purification of antibody fragments by affinity chromatography. The advances in protein engineering and computational biology approaches lead to the development of engineered affinity ligands with improved properties including binding affinity. In this study, molecular dynamics simulations (MDs) and Osprey software were used to design single B domains of the Protein L with higher affinity to antibody fragments. The modified B domains were then polymerized to ligand with six B domains by homology modeling methods. The results showed that single B domain mutants of MB1 (Thr865Trp) and MB2 (Thr847Met-Thr865Trp) had higher binding affinity to Fab compared to the wild single B domain. Also, MDs and molecular docking results showed that the polymerized Proteins L including the wild and mutated six B domains (6B0, 6B1, and 6B2) were stable during MDs and the two mutants of 6B1 and 6B2 showed higher binding affinity to Fab relative to the wild type.Communicated by Ramaswamy H. Sarma.

Comparative analysis of permanent and transient domain-domain interactions in multi-domain proteins

Protein domains are structural, functional, and evolutionary units. These domains bring out the diversity of functionality by means of interactions with other co-existing domains and provide stability. Hence, it is important to study intra-protein inter-domain interactions from the perspective of types of interactions. Domains within a chain could interact over short timeframes or permanently, rather like protein-protein interactions (PPIs). However, no systematic study has been carried out between two classes, namely permanent and transient domain-domain interactions. In this work, we studied 263 two-domain proteins, belonging to either of these classes and their interfaces on the basis of several factors, such as interface area and details of interactions (number, strength, and types of interactions). We also characterized them based on residue conservation at the interface, correlation of residue motions across domains, its involvement in repeat formation, and their involvement in particular molecular processes. Finally, we could analyze the interactions arising from domains in two-domain monomeric proteins, and we observed significant differences between these two classes of domain interactions and a few similarities. This study will help to obtain a better understanding of structure-function and folding principles of multi-domain proteins.

In silico characterization, molecular docking, and dynamic simulation of a novel fungal cell-death suppressing effector, MoRlpA as potential cathepsin B-like cysteine protease inhibitor during rice blast infection

The blast fungus Magnaporthe oryzae is one of the most notorious pathogens affecting rice production worldwide. The cereal killer employs a special class of small secreted proteins called effectors to manipulate and perturb the host metabolism. In turn, the host plants trigger effector-triggered immunity (ETI) via localized cell death and hypersensitive response (HR). We have identified and characterized a novel secreted effector MoRlpA from M. oryzae by extensive in silico methods. The localization studies suggested that it is exclusively secreted in the host apoplasts. Interestingly, MoRlpA interacts with a protease, cathepsin B from rice with highest affinity. The 3D structural models of both the proteins were generated. Cathepsin B-like cysteine proteases are usually involved in programmed cell death (PCD) and autophagy in plants which lead to generation of HR upon infection. Our results suggest that MoRlpA interacts with rice cathepsin B-like cysteine protease and demolish the host counter-attack by suppressing cell death and HR during an active blast infection. This was further validated by molecular docking and molecular dynamic simulation analyses. The important residues involved in the rice-blast pathogen interactions were deciphered. Overall, this research highlights stable interactions between MoRlpA-OsCathB during rice blast pathogenesis and providing an insight into how this novel RlpA protease inhibitor-cum-effector modulates the host's apoplast to invade the host tissues and establish a successful infection. Thus, this research will help to develop potential fungicide to block the binding region of MoRlpA target so that the cryptic pathogen would be recognized by the host. HIGHLIGHTSFor the first time, a novel secreted effector protein, MoRlpA has been identified and characterised from M. oryzae in silicoMoRlpA contains a rare lipoprotein A-like DPBB domain which is often an enzymatic domain in other systemsMoRlpA as an apoplastic effector interacts with the rice protease OsCathB to suppress the cell death and hypersensitive response during rice blast infectionThe three-dimensional structures of both the MoRlpA and OsCathB proteins were predictedMoRlpA-OsCathB interactions were analysed by molecular docking and molecular dynamic simulation studiesCommunicated by Ramaswamy H. Sarma.

Immunoinformatics-guided design of a multi-valent vaccine against Rotavirus and Norovirus (ChRNV22)

Rotavirus (RV) and Norovirus (NV) are the main viral etiologic agents of acute gastroenteritis (AG), a serious pediatric condition associated with significant death rates and long-term complications. Anti-RV vaccination has been proved efficient in the reduction of severe AG worldwide, however, the available vaccines are all attenuated and have suboptimal efficiencies in developing countries, where AG leads to substantial disease burden. On the other hand, no NV vaccine has been licensed so far. Therefore, we used immunoinformatics tools to develop a multi-epitope vaccine (ChRNV22) to prevent severe AG by RV and NV. Epitopes were predicted against 17 prevalent genotypes of four structural proteins (NV's VP1, RV's VP4, VP6 and VP7), and then assembled in a chimeric protein, with two small adjuvant sequences (tetanus toxin P2 epitope and a conserved sequence of RV's enterotoxin, NSP4). Simulations of the immune response and interactions with immune receptors indicated the immunogenic properties of ChRNV22, including a Th1-biased response. In silico search for putative host-homologous, allergenic and toxic regions also indicated the vaccine safety. In summary, we developed a multi-epitope vaccine against different NV and RV genotypes that seems promising for the prevention of severe AG, which will be further assessed by in vivo tests.

Structural and dynamic insights into α-synuclein dimer conformations

Parkinson disease is associated with the aggregation of the protein α-synuclein. While α-synuclein can exist in multiple oligomeric states, the dimer has been a subject of extensive debates. Here, using an array of biophysical approaches, we demonstrate that α-synuclein in vitro exhibits primarily a monomer-dimer equilibrium in nanomolar concentrations and up to a few micromolars. We then use spatial information from hetero-isotopic cross-linking mass spectrometry experiments as restrains in discrete molecular dynamics simulations to obtain the ensemble structure of dimeric species. Out of eight structural sub-populations of dimers, we identify one that is compact, stable, abundant, and exhibits partially exposed β-sheet structures. This compact dimer is the only one where the hydroxyls of tyrosine 39 are in proximity that may promote dityrosine covalent linkage upon hydroxyl radicalization, which is implicated in α-synuclein amyloid fibrils. We propose that this α-synuclein dimer features etiological relevance to Parkinson disease.

Biophysical evaluation of the oligomerization and conformational properties of the N-terminal domain of TDP-43

TDP-43 is an RNA-binding protein that presents four domains comprising an N-terminal region, two RNA recognition motifs and a C-terminal region. The N-terminal domain (NTD) has a relevant role in the oligomerization and splicing activity of TDP-43. In this work, we have expressed, purified and biophysically characterized the region that includes residues 1 to 102 that contains the nuclear localization signal (residues 80-102, NLS). Furthermore, we have evaluated the oligomerization equilibrium for this protein fragment. Also, we have determined changes in the tertiary structure and its stability in a broad range of pH values by means of different spectroscopic methods. Additionally, we compared this fragment with the one that lacks the NLS employing experimental and computational methods. Finally, we evaluated the motion of dimeric forms to get insights into the conformational flexibility of this TDP-43 module in an oligomeric state. Our results suggest that this domain has a conformational plasticity in the vicinity of the single tryptophan of this domain (Trp68), which is enhanced by the presence of the nuclear localization signal. All these results help to understand the molecular features of the NTD of TDP-43.

EGFR and p38MAPK Contribute to the Apoptotic Effect of the Recombinant Lectin from Tepary Bean (Phaseolus acutifolius) in Colon Cancer Cells

Previous works showed that a Tepary bean lectin fraction (TBLF) induced apoptosis on colon cancer cells and inhibited early colonic tumorigenesis. One Tepary bean (TB) lectin was expressed in Pichia pastoris (rTBL-1), exhibiting similarities to one native lectin, where its molecular structure and in silico recognition of cancer-type N-glycoconjugates were confirmed. This work aimed to determine whether rTBL-1 retained its bioactive properties and if its apoptotic effect was related to EGFR pathways by studying its cytotoxic effect on colon cancer cells. Similar apoptotic effects of rTBL-1 with respect to TBLF were observed for cleaved PARP-1 and caspase 3, and cell cycle G0/G1 arrest and decreased S phase were observed for both treatments. Apoptosis induction on SW-480 cells was confirmed by testing HA2X, p53 phosphorylation, nuclear fragmentation, and apoptotic bodies. rTBL-1 increased EGFR phosphorylation but also its degradation by the lysosomal route. Phospho-p38 increased in a concentration- and time-dependent manner, matching apoptotic markers, and STAT1 showed activation after rTBL-1 treatment. The results show that part of the rTBL-1 mechanism of action is related to p38 MAPK signaling. Future work will focus further on the target molecules of this recombinant lectin against colon cancer.

The CDR3 region as the major driver of TREM-1 interaction with its ligands, an in silico characterization

The triggering receptor expressed on myeloid cells-1 (TREM-1) is a pattern recognition receptor heavily investigated in infectious and non-infectious diseases. Because of its role in amplifying inflammation, TREM-1 has been explored as a diagnostic/prognostic biomarker. Further, as the receptor has been implicated in the pathophysiology of several diseases, therapies aiming at modulating its activity represent a promising strategy to constrain uncontrolled inflammatory or infectious diseases. Despite this, several aspects concerning its interaction with ligands and activation process, remain unclear. Although many molecules have been suggested as TREM-1 ligands, only five have been confirmed to interact with the receptor: actin, eCIRP, HMGB1, Hsp70 and PGLYRP1. However, the domains involved in the interaction between the receptor and these proteins are not clarified yet. Therefore, here we used in silico approaches to investigate the putative binding domains in the receptor, using hot spots analysis, molecular docking and molecular dynamics simulations between TREM-1 and its five known ligands. Our results indicated the complementarity-determining regions (CDRs) of the receptor as the main mediators of antigen recognition, especially the CDR3 loop. We believe that our study could be used as structural basis for the elucidation of TREM-1's recognition process, and may be useful for prospective in silico and biological investigations exploring the receptor in different contexts.

Human TMPRSS2 non-catalytic ectodomain and SARS-CoV-2 S2' subunit interaction mediated SARS-CoV-2 endocytosis: a model proposal with virtual screening for potential drug molecules to inhibit this interaction

This study proposes a novel model for integration of SARS-CoV-2 into host cell via endocytosis as a possible alternative to the prevailing direct fusion model. It is known that the SARS-CoV-2 spike protein undergoes proteolytic cleavage at S1-S2 cleavage site and the cleaved S2 domain is primed by the activated serine protease domain (SPD) of humanTMPRSS2 to become S2'. The activated SPD of TMPRSS2 is formed after it is cleaved by autocatalysis from the membrane bound non-catalytic ectodomain (hNECD) comprising of LDLRA CLASS-I repeat and a SRCR domain. It is known that the SRCR domains as well as LDLRA repeat harboring proteins mediate endocytosis of viruses and certain ligands. Based on this, we put forward a hypothesis that the exposed hNECD binds to the S2' as both are at an interaction proximity soon after S2 is processed by the SPD and this interaction may lead to the endocytosis of virus. Based on this hypothesis we have modelled the hNECD structure, followed by docking studies with the known 3D structure of S2'. The interaction interface of hNECD with S2' was further used for virtual screening of FDA-approved drug molecules and Indian medicinal plant-based compounds. We also mapped the known mutations of concern and mutations of interest on interaction interface of S2' and found that none of the known mutations map onto the interaction interface. This indicates that targeting the interaction between the hNECD of TMPRSS2 and S2' may serve as an attractive therapeutic target.Communicated by Ramaswamy H. Sarma.

Anurans against SARS-CoV-2: A review of the potential antiviral action of anurans cutaneous peptides

At the end of 2019, in China, clinical signs and symptoms of unknown etiology have been reported in several patients whose sample sequencing revealed pneumonia caused by the SARS-CoV-2 virus. COVID-19 is a disease triggered by this virus, and in 2020, the World Health Organization declared it a pandemic. Since then, efforts have been made to find effective therapeutic agents against this disease. Identifying novel natural antiviral drugs can be an alternative to treatment. For this reason, antimicrobial peptides secreted by anurans' skin have gained attention for showing a promissory antiviral effect. Hence, this review aimed to elucidate how and which peptides secreted by anurans' skin can be considered therapeutic agents to treat or prevent human viral infectious diseases. Through a literature review, we attempted to identify potential antiviral frogs' peptides to combat COVID-19. As a result, the Magainin-1 and -2 peptides, from the Magainin family, the Dermaseptin-S9, from the Dermaseptin family, and Caerin 1.6 and 1.10, from the Caerin family, are molecules that already showed antiviral effects against SARS-CoV-2 in silico. In addition to these peptides, this review suggests that future studies should use other families that already have antiviral action against other viruses, such as Brevinins, Maculatins, Esculentins, Temporins, and Urumins. To apply these peptides as therapeutic agents, experimental studies with peptides already tested in silico and new studies with other families not tested yet should be considered.

Design, synthesis and evaluation of 4,7-disubstituted 8-methoxyquinazoline derivatives as potential cytotoxic agents targeting β-catenin/TCF4 signaling pathway

•

Novel 4,7-disubstituted 8-methoxyquinazoline based derivativesdisrupting β-catenin/TCF4 interactions were synthesized.

•

Cytotoxic potencies of these derivatives were found to be comparable with imatinib mesylate against constitutively activated β-catenin/TCF4 signaling cancer cells (IC₅₀ ranged from 5.64 ± 0.68 to 23.18 ±0.45 μM against HCT116 and HepG2 cells).

•

Compound (18B) was found to be the most potent compound.

•

Compound (18B) induced apoptosis and inhibited cell migration against HCT116 and HepG2 cells.

•

Compound (18B) significantly downregulated β-catenin/TCF4 signaling pathway (reduced TOP-luciferase activity), β-catenin and TCF4 protein expression, and mRNA levels of c-MYC and Cyclin D1 (Wnt target genes) against HCT116 cells and showed cytotoxicity against primary human gallbladder cancer cells (IC₅₀ 8.50±1.44).

Overactivation of Wnt/β-catenin signaling by accumulated β-catenin in the nucleus has been shown to play a crucial role in the etiology of cancer. Interaction of β-catenin with Transcription factor 4 (TCF4) is a key step for the activation of Wnt genes in response to upstream signals of the Wnt/β-catenin pathway. Hence, down regulation of Wnt/β-catenin signaling or targeting downstream events by selective β-catenin/TCF4 protein–protein interaction inhibitors could be a potential therapeutic strategy against such cancers. In this study structure-based drug design approach was followed to design novel 4,7-disubstituted 8-methoxyquinazoline-based derivatives which could act as potential cytotoxic agents inhibiting the β-catenin/TCF4 protein–protein interactions. Fifteen compounds possessing 4,7-disubstituted 8-methoxyquinazoline scaffold were synthesized. Cytotoxic potential of the synthesised derivatives were determined against constitutively activated β-catenin/TCF4 signaling pathway cancer cells (HCT116 and HepG2) using the sulforhodamine B assay. The most potent compound (18B) was selected for detailed biological evaluation. Cell morphology, Hoechst 33342 and Annexin V/PI staining were used to detect apoptosis, while inhibition of cell migration was assessed by in vitro wound healing assay against HCT116 and HepG2 cells. Effect on β-catenin/TCF mediated transcriptional activity was assessed by TOPFlash/FOPFlash assay, TCF4 and β-catenin protein expression by immunocytofluorescence, and Wnt target genes (like c-MYC and Cyclin D1) mRNA levels by RT-PCR against HCT116 cells. Cytotoxic potency of the most potential compound (18B) against primary human gallbladder cancer cells was also evaluated. The derivatives showed interactions with active site residues of β-catenin and were capable of hindering the TCF4 binding, thereby disrupting β-catenin/TCF4 interactions. Cytotoxic potencies (IC₅₀) of these derivatives ranged from 5.64 ± 0.68 to 23.18 ± 0.45 μM against HCT116 and HepG2 cells respectively. Compound (18B), the most potent compound among the series, induced apoptosis and inhibited cell migration against HCT116 and HepG2 cells. Mechanistic studies indicated that compound (18B) downregulated β-catenin/TCF4 signaling pathway, β-catenin and TCF4 protein expression, and mRNA levels of c-MYC andCyclin D1 in HCT116 cells and showed cytotoxicity against primary human gallbladder cancer cells with IC₅₀ value of 8.50 ± 1.44 μM. Thus, novel 4,7-disubstituted 8-methoxyquinazoline derivatives were identified as potential cytotoxic agents with potencies comparable to that of imatinib mesylate. Compound (18B) represents a promising lead molecule as anticancer agent against colon, hepatocellular and gallbladder cancers targeting β-catenin/TCF4 signaling pathway.

In silico investigations of heparin binding to SARS-CoV-2 variants with a focus at the RBD/ACE2 interface

The increased infectivity and transmissibility of SARS-CoV-2 new variants were contributed largely by increase binding of receptor binding domain (RBD) domain of the Spike (S) protein to its cellular receptor ACE2 (Angiotensin-Converting Enzyme 2). Several studies have indicated that heparin and its derivatives interact to SARS-CoV-2 S-RBD and inhibits the binding of ACE2 which blocks the viral invasion. However, it is largely unclear how these SARS-CoV-2 variants affects ACE2 binding in the presence of heparin. Herein, using the molecular docking and interaction energy analysis, we showed that N501Y, L452R-E484Q, and E484K mutations bind strongly with heparin in the range of - 7.4 to - 8.0 kcal/mol. The triple mutations, K417N-E484K-N501Y, and K417T-E484K-N501Y displayed weaker binding affinity to heparin (-6.6 kcal/mol). Further, we showed that most of the RBD mutations increased the binding affinity of ACE2 in the absence of heparin, with the maximum increase observed for N501Y (-13.7 kcal/mol). Also, in the presence of heparin, ACE2 binds strongly to the mutant RBD as compared to WT RBD. The strong RBD/ACE2 interaction was observed in case of triple variants (-11.3 kcal/mol) whereas, N501Y showed weakest binding of RBD/ACE2 in the presence of heparin (-9.2 kcal/mol). The strong binding of ACE2 to RBD-heparin complex in these variants will leads to strong inhibition of their entry into host cells.

Sequence Divergence and Functional Specializations of the Ancient Spliceosomal SF3b: Implications in Flexibility and Adaptations of the Multi-Protein Complex

Multi-protein assemblies are complex molecular systems that perform highly sophisticated biochemical functions in an orchestrated manner. They are subject to changes that are governed by the evolution of individual components. We performed a comparative analysis of the ancient and functionally conserved spliceosomal SF3b complex, to recognize molecular signatures that contribute to sequence divergence and functional specializations. For this, we recognized homologous sequences of individual SF3b proteins distributed across 10 supergroups of eukaryotes and identified all seven protein components of the complex in 578 eukaryotic species. Using sequence and structural analysis, we establish that proteins occurring on the surface of the SF3b complex harbor more sequence variation than the proteins that lie in the core. Further, we show through protein interface conservation patterns that the extent of conservation varies considerably between interacting partners. When we analyze phylogenetic distributions of individual components of the complex, we find that protein partners that are known to form independent subcomplexes are observed to share similar profiles, reaffirming the link between differential conservation of interface regions and their inter-dependence. When we extend our analysis to individual protein components of the complex, we find taxa-specific variability in molecular signatures of the proteins. These trends are discussed in the context of proline-rich motifs of SF3b4, functional and drug binding sites of SF3b1. Further, we report key protein-protein interactions between SF3b1 and SF3b6 whose presence is observed to be lineage-specific across eukaryotes. Together, our studies show the association of protein location within the complex and subcomplex formation patterns with the sequence conservation of SF3b proteins. In addition, our study underscores evolutionarily flexible elements that appear to confer adaptive features in individual components of the multi-protein SF3b complexes and may contribute to its functional adaptability.

SARS-CoV-2 pan-variant inhibitory peptides deter S1-ACE2 interaction and neutralize delta and omicron pseudoviruses

Approved neutralizing antibodies that target the prototype Spike are losing their potency against the emerging variants of concern (VOCs) of SARS-CoV-2, particularly Omicron. Although SARS-CoV-2 is continuously adapting the host environment, emerging variants recognize the same ACE2 receptor for cell entry. Protein and peptide decoys derived from ACE2 or Spike proteins may hold the pan-variant inhibitory potential. Here, we deployed interactive structure- and pharmacophore-based approaches to design short and stable peptides –Coronavirus Spike Neutralizing Peptides (CSNPs)– capable of neutralizing all SARS-CoV-2 VOCs. After in silico structural stability investigation and free energies perturbation of the isolated and target-bound peptides, nine candidate peptides were evaluated for the biophysical interaction through SPR assay. CSNP1, CSNP2, and Pep1 dose-dependently bind the S1 domain of the prototype Spike, whereas CSNP4 binds both S1 and ACE2. After safety and immunocytochemistry evaluation, peptides were probed for their pan-variant inhibitory effects. CSNP1, CSNP2, and CSNP4 inhibited all VOCs dose-dependently, whereas Pep1 had a moderate effect. CSNP2 and CSNP4 could neutralize the wild-type pseudovirus up to 80 % when treated at 0.5 µM. Furthermore, CSNP4 synergize the neutralization effect of monoclonal antibody and CSNP1 in Delta variant pseudovirus assay as they target different regions on the RBD. Thus, we suggest that CSNPs are SARS-CoV-2 pan-variant inhibitory candidates for COVID-19 therapy, which may pave the way for combating the emerging immune-escaping variants. We also propose that CSNP1/2-CSNP4 peptide cocktail or CSNP1/4 mAbs cocktail with no overlapping epitopes could be effective therapeutic strategies against COVID-19.

An integrated computational pipeline for designing high-affinity nanobodies with expanded genetic codes

Protein engineering and design principles employing the 20 standard amino acids have been extensively used to achieve stable protein scaffolds and deliver their specific activities. Although this confers some advantages, it often restricts the sequence, chemical space, and ultimately the functional diversity of proteins. Moreover, although site-specific incorporation of non-natural amino acids (nnAAs) has been proven to be a valuable strategy in protein engineering and therapeutics development, its utility in the affinity-maturation of nanobodies is not fully explored. Besides, current experimental methods do not routinely employ nnAAs due to their enormous library size and infinite combinations. To address this, we have developed an integrated computational pipeline employing structure-based protein design methodologies, molecular dynamics simulations and free energy calculations, for the binding affinity prediction of an nnAA-incorporated nanobody toward its target and selection of potent binders. We show that by incorporating halogenated tyrosines, the affinity of 9G8 nanobody can be improved toward epidermal growth factor receptor (EGFR), a crucial cancer target. Surface plasmon resonance (SPR) assays showed that the binding of several 3-chloro-l-tyrosine (3MY)-incorporated nanobodies were improved up to 6-fold into a picomolar range, and the computationally estimated binding affinities shared a Pearson's r of 0.87 with SPR results. The improved affinity was found to be due to enhanced van der Waals interactions of key 3MY-proximate nanobody residues with EGFR, and an overall increase in the nanobody's structural stability. In conclusion, we show that our method can facilitate screening large libraries and predict potent site-specific nnAA-incorporated nanobody binders against crucial disease-targets.

Unraveling the Thermodynamics of Ultra-tight Binding of Intrinsically Disordered Proteins

Protein interactions mediated by the intrinsically disordered proteins (IDPs) are generally associated with lower affinities compared to those between globular proteins. Here, we characterize the association between the intrinsically disordered HigA2 antitoxin and its globular target HigB2 toxin from Vibrio cholerae using competition ITC experiments. We demonstrate that this interaction reaches one of the highest affinities reported for IDP-target systems (K_D = 3 pM) and can be entirely attributed to a short, 20-residue-long interaction motif that folds into α-helix upon binding. We perform an experimentally based decomposition of the IDP-target association parameters into folding and binding contributions, which allows a direct comparison of the binding contribution with those from globular ultra-high affinity binders. We find that the HigA2-HigB2 interface is energy optimized to a similar extent as the interfaces of globular ultra-high affinity complexes, such as barnase-barstar. Evaluation of other ultra-high affinity IDP-target systems shows that a strategy based on entropy optimization can also achieve comparably high, picomolar affinities. Taken together, these examples show how IDP-target interactions achieve picomolar affinities either through enthalpy optimization (HigA2-HigB2), resembling the ultra-high affinity binding of globular proteins, or via bound-state fuzziness and entropy optimization (CcdA-CcdB, histone H1-prothymosin α).

Discovery of a Natural Product with Potent Efficacy Against SARS-CoV-2 by Drug Screening

The novel coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has spread worldwide for almost 2 years. It starts from viral adherence to host cells through an interaction between spike glycoprotein 1 (S1) containing a receptor-binding domain (RBD) and human angiotensin-converting enzyme-2 (ACE2). One of the useful strategies to prevent SARS-CoV-2 infection is to inhibit the attachment of RBD to ACE2. Therefore, the current work proposed potent peptides against SARS-CoV-2 infection by carrying out MM-PBSA calculation based on the binding of 52 antiviral peptides (AVPs) to RBD. Considering the binding free energies of AVPs to RBD, cyanovirin-N (CV-N) showed the strongest RBD binding affinity among 52 AVPs. Upon structural analysis of RBD complex with CV-N, it was observed that 12 of the 13 key residues of RBD binding to ACE2 were hijacked by CV-N. CV-N bound to RBD at a smaller affinity of 14.9 nM than that of ACE2 and inhibited the recruitment of S1 to human alveolar epithelial cells. Further analysis revealed that CV-N suppressed SARS-CoV-2 S pseudovirion infection with a half-maximal inhibitory concentration (IC₅₀) of 18.52 μg/mL. This study demonstrated a drug screening for AVPs against SARS-CoV-2 and discovered a peptide with inspiring antiviral properties, which provided a promising strategy for the COVID-19 therapeutic approach.

HD5 and LL-37 Inhibit SARS-CoV and SARS-CoV-2 Binding to Human ACE2 by Molecular Simulation

The coronavirus (COVID-19) pandemic is still spreading all over the world. As reported, angiotensin-converting enzyme-2 (ACE2) is a receptor of SARS-CoV-2 spike protein that initializes viral entry into host cells. Previously, the human defensin 5 (HD5) has been experimentally confirmed to be functional against the SARS-CoV-2. The present study proposes a human cathelicidin known as LL37 that strongly binds to the carboxypeptidase domain of human ACE2 compared to HD5. Therefore, LL37 bears a great potential to be tested as an anti-SARS-CoVD-2 peptide. We investigated the molecular interactions formed between the LL37 and ACE2 as well as HD5 and ACE2 tailed by their thermodynamic stability. The MM-PBSA and free energy landscape analysis outcomes confirmed its possible inhibitory effect against the SARS-CoV-2. The results obtained here could help propose a promising therapeutic strategy against the havoc caused by SARS-CoV-2 infections.

Structural Design and Analysis of the RHOA-ARHGEF1 Binding Mode: Challenges and Applications for Protein-Protein Interface Prediction

The interaction between two proteins may involve local movements, such as small side-chains re-positioning or more global allosteric movements, such as domain rearrangement. We studied how one can build a precise and detailed protein-protein interface using existing protein-protein docking methods, and how it can be possible to enhance the initial structures using molecular dynamics simulations and data-driven human inspection. We present how this strategy was applied to the modeling of RHOA-ARHGEF1 interaction using similar complexes of RHOA bound to other members of the Rho guanine nucleotide exchange factor family for comparative assessment. In parallel, a more crude approach based on structural superimposition and molecular replacement was also assessed. Both models were then successfully refined using molecular dynamics simulations leading to protein structures where the major data from scientific literature could be recovered. We expect that the detailed strategy used in this work will prove useful for other protein-protein interface design. The RHOA-ARHGEF1 interface modeled here will be extremely useful for the design of inhibitors targeting this protein-protein interaction (PPI).

Computational designing of a peptide that potentially blocks the entry of SARS-CoV, SARS-CoV-2 and MERS-CoV

Last decade has witnessed three major pandemics caused by SARS-CoV, SARS-CoV-2 and MERS-CoV that belong to Coronavirus family. Currently, there are no effective therapies available for corona virus infections. Since the three viruses belong to the same family and share many common features, we can theoretically design a drug that can be effective on all the three of them. In this study, using computational approach, we designed a peptide (Peptide 7) that can bind to the Receptor Binding Domain (RBD) of SARS-CoV, SARS-CoV-2 and MERS-CoV thereby preventing the entry of the viruses into the host cell. The peptide inhibitor was designed as a consensus peptide from three different peptides that might individually bind to the RBD of the three viruses. Docking studies and molecular dynamic simulations using Peptide 7 has shown that it binds with higher affinity than the native receptors of the RBD and forms a stable complex thereby preventing further viral-receptor interaction and inhibiting their cellular entry. This effective binding is observed for the three RBDs, despite the Peptide 7 interactions being slightly different. Hence; this peptide inhibitor can be used as a potential candidate for the development of peptide based anti-viral therapy against Corona viruses.

Identification of Critical Residues in the Carboxy Terminus of the Dopamine Transporter Involved in the G Protein βγ-Induced Dopamine Efflux

The dopamine transporter (DAT) plays a crucial role in the regulation of brain dopamine (DA) homeostasis through the re-uptake of DA back into the presynaptic terminal. In addition to re-uptake, DAT is also able to release DA through a process referred to as DAT-mediated DA efflux. This is the mechanism by which potent and highly addictive psychostimulants, such as amphetamine (AMPH) and its analogues, increase extracellular DA levels in motivational and reward areas of the brain. Recently, we discovered that G protein βγ subunits (Gβγ) binds to the DAT, and that activation of Gβγ results in DAT-mediated efflux - a similar mechanism as AMPH. Previously, we have shown that Gβγ binds directly to a stretch of 15 residues within the intracellular carboxy terminus of DAT (residues 582-596). Additionally, a TAT peptide containing residues 582 to 596 of DAT was able to block the Gβγ-induced DA efflux through DAT. Here, we use a combination of computational biology, mutagenesis, biochemical, and functional assays to identify the amino acid residues within the 582-596 sequence of the DAT carboxy terminus involved in the DAT-Gβγ interaction and Gβγ-induced DA efflux. Our in-silico protein-protein docking analysis predicted the importance of F587 and R588 residues in a network of interactions with residues in Gβγ. In addition, we observed that mutating R588 to alanine residue resulted in a mutant DAT which exhibited attenuated DA efflux induced by Gβγ activation. We demonstrate that R588, and to a lesser extent F5837, located within the carboxy terminus of DAT play a critical role in the DAT-Gβγ physical interaction and promotion of DA efflux. These results identify a potential new pharmacological target for the treatment of neuropsychiatric conditions in which DAT functionality is implicated including ADHD and substance use disorder.

An in-silico approach to identify the potential hot spots in SARS-CoV-2 spike RBD to block the interaction with ACE2 receptor

A novel acute viral pneumonia induced by SARS-CoV-2 exploded at the end of 2019, causing a severe medical and economic crisis. For developing specific pharmacotherapy against SARS-CoV-2, an in silico virtual screening was developed for the available in-house molecules. The conserved domain analysis was performed to identify the highly conserved and exposed amino acid regions in the SARS-CoV-2-S RBD sites. The Protein-Protein interaction analyses demonstrated the higher affinity between the SARS-CoV-2-S and ACE2 due to varieties of significant interactions between them. The computational alanine scanning mutation study has recognized the highly stabilized amino acids in the SARS-CoV-2-S RBD/ACE2 complex. The cumulative sequence investigations have inferred that Lys417, Phe486, Asn487, Tyr489, and Gln493 are perhaps the iconic target amino acids to develop a drug molecule or vaccine against SARS-CoV-2 infection. Most of the selected compounds include luteolin, zhebeirine, 3-dehydroverticine, embelin, andrographolide, ophiopogonin D, crocin-1, sprengerinin A, B, C, peimine, etc. were exhibited distinguish drug actions through the strong hydrogen bonding with the hot spots of the RBD. Besides, the 100 ns molecular dynamics simulation and free energy binding analysis showed the significant efficacy of luteolin to inhibit the infection of SARS-CoV-2.Communicated by Ramaswamy H. Sarma.

Unraveling the mechanism of arbidol binding and inhibition of SARS-CoV-2: Insights from atomistic simulations

The COVID-19 pandemic has spread rapidly and poses an unprecedented threat to the global economy and human health. Broad-spectrum antivirals are currently being administered to treat severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). China's prevention and treatment guidelines suggest the use of an anti-influenza drug, arbidol, for the clinical treatment of COVID-19. Reports indicate that arbidol could neutralize SARS-CoV-2. Monotherapy with arbidol is superior to lopinavir-ritonavir or favipiravir for treating COVID-19. In SARS-CoV-2 infection, arbidol acts by interfering with viral binding to host cells. However, the detailed mechanism by which arbidol induces the inhibition of SARS-CoV-2 is not known. Here, we present atomistic insights into the mechanism underlying membrane fusion inhibition of SARS-CoV-2 by arbidol. Molecular dynamics (MD) simulation-based analyses demonstrate that arbidol binds and stabilizes at the receptor-binding domain (RBD)/ACE2 interface with a high affinity. It forms stronger intermolecular interactions with the RBD than ACE2. Analyses of the detailed decomposition of energy components and binding affinities revealed a substantial increase in the affinity between the RBD and ACE2 in the arbidol-bound RBD/ACE2 complex, suggesting that arbidol generates favorable interactions between them. Based on our MD simulation results, we propose that the binding of arbidol induces structural rigidity in the viral glycoprotein, thus restricting the conformational rearrangements associated with membrane fusion and virus entry. Furthermore, key residues of the RBD and ACE2 that interact with arbidol were identified, opening the door for developing therapeutic strategies and higher-efficacy arbidol derivatives or lead drug candidates.

Classification and prediction of protein-protein interaction interface using machine learning algorithm

Structural insight of the protein-protein interaction (PPI) interface can provide knowledge about the kinetics, thermodynamics and molecular functions of the complex while elucidating its role in diseases and further enabling it as a potential therapeutic target. However, owing to experimental lag in solving protein-protein complex structures, three-dimensional (3D) knowledge of the PPI interfaces can be gained via computational approaches like molecular docking and post-docking analyses. Despite development of numerous docking tools and techniques, success in identification of native like interfaces based on docking score functions is limited. Hence, we employed an in-depth investigation of the structural features of the interface that might successfully delineate native complexes from non-native ones. We identify interface properties, which show statistically significant difference between native and non-native interfaces belonging to homo and hetero, protein-protein complexes. Utilizing these properties, a support vector machine (SVM) based classification scheme has been implemented to differentiate native and non-native like complexes generated using docking decoys. Benchmarking and comparative analyses suggest very good performance of our SVM classifiers. Further, protein interactions, which are proven via experimental findings but not resolved structurally, were subjected to this approach where 3D-models of the complexes were generated and most likely interfaces were predicted. A web server called Protein Complex Prediction by Interface Properties (PCPIP) is developed to predict whether interface of a given protein-protein dimer complex resembles known protein interfaces. The server is freely available at http://www.hpppi.iicb.res.in/pcpip/ .

A conserved arginine residue in a staphylococcal anti-sigma factor is required to preserve its kinase activity, structure, and stability

RsbW, σ^B, and RsbV, encoded by Staphylococcus aureus and related bacteria, act as an anti-sigma factor, an sigma factor, and an anti-anti-sigma factor, respectively. The interaction between RsbW and σ^B blocks the transcription initiation activity of the latter protein. RsbW also functions as a serine kinase and phosphorylates RsbV in the presence of ATP. Our modeling study indicates that the RsbW-RsbV complex is stabilized by twenty-four intermolecular non-covalent bonds. Of the bond-forming RsbW residues, Arg 23, and Glu 49 are conserved residues. To understand the roles of Arg 23 in RsbW, rRsbW[R23A], a recombinant S. aureus RsbW (rRsbW) harboring Arg to Ala change at position 23, was investigated using various probes. The results reveal that rRsbW[R23A], like rRsbW, exists as the dimers in the aqueous solution. However, rRsbW[R23A], unlike rRsbW, neither interacted with a chimeric RsbV (rRsbV) nor formed the phosphorylated rRsbV in the presence of ATP. Furthermore, the tertiary structure and hydrophobic surface area of rRsbW[R23A] matched little with those of rRsbW. Conversely, both rRsbW[R23A] and rRsbW showed interaction with a recombinant σ^B (rσ^B). rRsbW and rRsbW[R23A] were also unfolded via the formation of at least one intermediate in the presence of urea. However, the thermodynamic stability of rRsbW significantly differed from that of rRsbW[R23A]. Our molecular dynamics (MD) simulation study also reveals the substantial change of structure, dimension, and stability of RsbW due to the above mutation. The ways side chain of critical Arg 23 contributes to maintaining the tertiary structure, and stability of RsbW was elaborately discussed.Communicated by Ramaswamy H. Sarma.

Prediction of CD28-CD86 protein complex structure using different level of resolution approach

Immune system plays essential role in functioning of higher organisms. Its hyperactivity can lead to autoimmune diseases or even anaphylactic shock while hypoactivity leads to proneness to infections or even cancer. T-cells play crucial role in immunity mechanisms and their activation and inhibition is strictly controlled by the regulatory proteins, such as CD28 and CTLA-4. Activity of these proteins is controlled by a pair of ligands, named CD80 and CD86, which can non-covalently bound to their receptors. While structure of human CTLA-4-CD86 complex in known, there is still no available structure for the CD28-CD86 system. To obtain the reliable structure of CD28-CD86 complex we first validated our methodology on the CTLA-4-CD86 system. Then coarse-grained UNRES-dock molecular docking simulation was performed followed by all-atom molecular dynamics simulations. As a result, we obtained a complete CD28-CD86 complex structure on atomistic level, in which interaction interface is consistent with available data. We also determined the kinetic properties for CTLA4-CD86 and CD28-CD86 complexes with use of coarse-grained model and determined the key residues for complex formation with use of Robetta, PPCheck and HawkDock servers. Our results not only verify high accuracy of the UNRES-dock method, but also provide a highly reliable model of the CD28-CD86 complex, which can be used in further studies and drug design.

Experimentally based structural model of Yih1 provides insight into its function in controlling the key translational regulator Gcn2

Yeast impact homolog 1 (Yih1), or IMPACT in mammals, is part of a conserved regulatory module controlling the activity of General Control Nonderepressible 2 (Gcn2), a protein kinase that regulates protein synthesis. Yih1/IMPACT is implicated not only in many essential cellular processes, such as neuronal development, immune system regulation and the cell cycle, but also in cancer. Gcn2 must bind to Gcn1 in order to impair the initiation of protein translation. Yih1 hinders this key Gcn1-Gcn2 interaction by binding to Gcn1, thus preventing Gcn2-mediated inhibition of protein synthesis. Here, we solved the structures of the two domains of Saccharomyces cerevisiae Yih1 separately using Nuclear Magnetic Resonance and determined the relative positions of the two domains using a range of biophysical methods. Our findings support a compact structural model of Yih1 in which the residues required for Gcn1 binding are buried in the interface. This model strongly implies that Yih1 undergoes a large conformational rearrangement from a latent closed state to a primed open state to bind Gcn1. Our study provides structural insight into the interactions of Yih1 with partner molecules.

Interaction mechanism of Mycobacterium tuberculosis GroEL2 protein with macrophage Lectin-like, oxidized low-density lipoprotein receptor-1: An integrated computational and experimental study

BACKGROUND

Bacterial surface proteins act as potential adhesins or invasins. The GroEL is a signal peptide-free surface expressed protein that aids adhesion in Escherichia coli by binding to LOX-1 receptor of the host cells. Mycobacterium tuberculosis (Mtb) expresses GroEL2 protein, having high level sequence identity with E. coli GroEL. This study investigates the interaction mechanism of GroEL2 protein of Mtb with LOX-1 of macrophages using integrated computational and experimental approach.

METHODS

Mtb GroEL2 protein was purified as histidine tagged protein using Ni-NTA chromatography. Confocal and scanning electron microscopies were used to study the uptake of GroEL2 coated fluorescent latex beads through the LOX-1 receptor in RAW264.7 macrophage cell line. Docking studies were performed to understand the interaction between the GroEL2 and LOX-1 proteins. Polyinosinic acid (PIA) was used as a LOX-1 inhibitor in both in silico and in vitro experiments.

RESULTS

GroEL2 protein coating enhances uptake of latex beads into macrophages through LOX-1 receptor. LOX-1 inhibitor PIA decreased the uptake of GroEL2 coated latex beads. GroEL2 interacts with the key ligand binding regions of the LOX-1 receptor, such as the basic spine and the saddle hydrophobic patch. PIA molecule destabilized the LOX-1-GroEL2 docked complex.

CONCLUSION

Surface associated GroEL2 protein of Mtb is a potential ligand for macrophage LOX-1 receptor. Interaction between GroEL2 and LOX-1 receptor may be utilized by Mtb to gain its intracellular access.

GENERAL SIGNIFICANCE

Surface associated GroEL2 of Mtb may bind to the macrophage LOX-1 receptor, enabling the internalization of the bacteria and progression of the infection.

Computational Analysis of Non-covalent Interactions in Phycocyanin Subunit Interfaces

Protein-protein interactions are an important phenomenon in biological processes and functions. We used the manually curated non-redundant dataset of 118 phycocyanin interfaces to gain additional insight into this phenomenon using a robust inter-atomic non-covalent interaction analyzing tool PPCheck. Our observations indicate that there is a relatively high composition of hydrophobic residues at the interfaces. Most of the interface residues are clustered at the middle of the range which we call "standard-size" interfaces. Furthermore, the multiple interaction patterns founded in the present study indicate that more than half of the residues involved in these interactions participate in multiple and water-bridged hydrogen bonds. Thus, hydrogen bonds contribute maximally towards the stability of protein-protein complexes. The analysis shows that hydrogen bond energies contribute to about 88 % to the total energy and it also increases with interface size. Van der Waals (vdW) energy contributes to 9.3 %±1.7 % on average in these complexes. Moreover, there is about 1.9 %±1.5 % contribution by electrostatic energy. Nevertheless, the role by vdW and electrostatic energy could not be ignored in interface binding. Results show that the total binding energy is more for large phycocyanin interfaces. The normalized energy per residue was less than -16 kJ mol^-1 , while most of them have energy in the range from -6 to -14 kJ mol^-1 . The non-covalent interacting residues in these proteins were found to be highly conserved. Obtained results might contribute to the understanding of structural stability of this class of evolutionary essential proteins with increased practical application and future designs of novel protein-bioactive compound interactions.

Small protein-protein interfaces rich in electrostatic are often linked to regulatory function

Protein-protein interaction (PPI) is critical for several biological functions in living cells through the formation of an interface. Therefore, it is of interest to characterize protein-protein interfaces using an updated non-redundant structural dataset of 2557 homo (identical subunits) and 393 hetero (different subunits) dimer protein complexes determined by X-ray crystallography. We analyzed the interfaces using van der Waals (vdW), hydrogen bonding and electrostatic energies. Results show that on average homo and hetero interfaces are similar. Hence, we further grouped the 2950 interfaces based on percentage vdW to total energies into dominant (≥60%) and sub-dominant (<60%) vdW interfaces. Majority (92%) of interfaces have dominant vdW energy with large interface size (146 ± 87 (homo) and 137 ± 76 (hetero) residues) and interface area (1622 ± 1135 Ų (homo) and 1579 ± 1060 Ų (hetero)). However, a proportion (8%) of interfaces have sub-dominant vdW energy with small interface size (85 ± 46 (homo) and 88 ± 36 (hetero) residues) and interface area (823 ± 538 Ų (homo) and 881 ± 377 Ų (hetero)). It is found that large interfaces have two-fold more interface area and interface size than small interfaces with increasing hydrogen bonding energy to interface size. However, small interfaces have three-fold more electrostatics energy than large interfaces with increasing electrostatics to interface size. Thus, 8% of complexes having small interfaces with limited interface area and sub-dominant vdW energy are rich in electrostatics. It is interesting to observe that complexes having small interfaces are often associated with regulatory function. Hence, the observed structural features with known molecular function provide insights for the better understanding of PPI.Communicated by Ramaswamy H. Sarma.

Modeling of Protein⁻Protein Interactions in Cytokinin Signal Transduction

The signaling of cytokinins (CKs), classical plant hormones, is based on the interaction of proteins that constitute the multistep phosphorelay system (MSP): catalytic receptors-sensor histidine kinases (HKs), phosphotransmitters (HPts), and transcription factors-response regulators (RRs). Any CK receptor was shown to interact in vivo with any of the studied HPts and vice versa. In addition, both of these proteins tend to form a homodimer or a heterodimeric complex with protein-paralog. Our study was aimed at explaining by molecular modeling the observed features of in planta protein-protein interactions, accompanying CK signaling. For this purpose, models of CK-signaling proteins' structure from Arabidopsis and potato were built. The modeled interaction interfaces were formed by rather conserved areas of protein surfaces, complementary in hydrophobicity and electrostatic potential. Hot spots amino acids, determining specificity and strength of the interaction, were identified. Virtual phosphorylation of conserved Asp or His residues affected this complementation, increasing (Asp-P in HK) or decreasing (His-P in HPt) the affinity of interacting proteins. The HK-HPt and HPt-HPt interfaces overlapped, sharing some of the hot spots. MSP proteins from Arabidopsis and potato exhibited similar properties. The structural features of the modeled protein complexes were consistent with the experimental data.

UvrB protein of Corynebacterium pseudotuberculosis complements the phenotype of knockout Escherichia coli and recognizes DNA damage caused by UV radiation but not 8-oxoguanine in vitro

In prokaryotic cells, the UvrB protein plays a central role in nucleotide excision repair, which is involved in the recognition of bulky DNA lesions generated by chemical or physical agents. The present investigation aimed to characterize the uvrB gene of Corynebacterium pseudotuberculosis (CpuvrB) and evaluate its involvement in the DNA repair system of this pathogenic organism. In computational analysis, the alignment of the UvrB protein sequences of Escherichia coli, Mycobacterium tuberculosis, Bacillus caldotenax and Corynebacterium pseudotuberculosis showed high similarity and the catalytic amino acid residues and functional domains are preserved. A CpUvrB model was constructed by comparative modeling and presented structural similarity with the UvrB of E. coli. Moreover, in molecular docking analysis CpUvrB showed favorable interaction with EcUvrA and revealed a preserved ATP incorporation site. Heterologous functional complementation assays using E. coli uvrB-deficient cells exposed to UV irradiation showed that the CpUvrB protein contributed to an increased survival rate in relation to those in the absence of CpUvrB. Damaged oligonucleotides containing thymine dimer or 8-oxoguanine lesion were synthesized and incubated with CpUvrB protein, which was able to recognize and excise UV irradiation damage but not 8-oxoguanine. These results suggest that CpUvrB is involved in repairing lesions derived from UV light and encodes a protein orthologous to EcUvrB.

Computational Resources for Predicting Protein-Protein Interactions

Proteins are the essential building blocks and functional components of a cell. They account for the vital functions of an organism. Proteins interact with each other and form protein interaction networks. These protein interactions play a major role in all the biological processes and pathways. The previous methods of predicting protein interactions were experimental which focused on a small set of proteins or a particular protein. However, these experimental approaches are low-throughput as they are time-consuming and require a significant amount of human effort. This led to the development of computational techniques that uses high-throughput experimental data for analyzing protein-protein interactions. The main purpose of this review is to provide an overview on the computational advancements and tools for the prediction of protein interactions. The major databases for the deposition of these interactions are also described. The advantages, as well as the specific limitations of these tools, are highlighted which will shed light on the computational aspects that can help the biologist and researchers in their research.

Protein-protein interfaces are vdW dominant with selective H-bonds and (or) electrostatics towards broad functional specificity

Several catalysis, cellular regulation, immune function, cell wall assembly, transport, signaling and inhibition occur through Protein- Protein Interactions (PPI). This is possible with the formation of specific yet stable protein-protein interfaces. Therefore, it is of interest to understand its molecular principles using structural data in relation to known function. Several interface features have been documented using known X-ray structures of protein complexes since 1975. This has improved our understanding of the interface using structural features such as interface area, binding energy, hydrophobicity, relative hydrophobicity, salt bridges and hydrogen bonds. The strength of binding between two proteins is dependent on interface size (number of residues at the interface) and thus its corresponding interface area. It is known that large interfaces have high binding energy (sum of (van der Waals) vdW, H-bonds, electrostatics). However, the selective role played by each of these energy components and more especially that of vdW is not explicitly known. Therefore, it is important to document their individual role in known protein-protein structural complexes. It is of interest to relate interface size with vdW, H-bonds and electrostatic interactions at the interfaces of protein structural complexes with known function using statistical and multiple linear regression analysis methods to identify the prominent force. We used the manually curated non-redundant dataset of 278 hetero-dimeric protein structural complexes grouped using known functions by Sowmya et al. (2015) to gain additional insight to this phenomenon using a robust inter-atomic non-covalent interaction analyzing tool PPCheck (Anshul and Sowdhamini, 2015). This dataset consists of obligatory (enzymes, regulator, biological assembly), immune and nonobligatory (enzyme and regulator inhibitors) complexes. Results show that the total binding energy is more for large interfaces. However, this is not true for its individual energy factors. Analysis shows that vdW energies contribute to about 75% ± 11% on average among all complexes and it also increases with interface size (r2 ranging from 0.67 to 0.89 with p<0.01) at 95% confidence limit irrespective of molecular function. Thus, vdW is both dominant and proportional at the interface independent of molecular function. Nevertheless, H bond energy contributes to 15% ± 6.5% on average in these complexes. It also moderately increases with interface size (r2 ranging from 0.43 to 0.61 with p<0.01) only among obligatory and immune complexes. Moreover, there is about 11.3% ± 8.7% contribution by electrostatic energy. It increases with interface size specifically among non-obligatory regulator-inhibitors (r2 = 0.44). It is implied that both H-bonds and electrostatics are neither dominant nor proportional at the interface. Nonetheless, their presence cannot be ignored in binding. Therefore, H-bonds and (or) electrostatic energy having specific role for improved stability in complexes is implied. Thus, vdW is common at the interface stabilized further with selective H-bonds and (or) electrostatic interactions at an atomic level in almost all complexes. Comparison of this observation with residue level analysis of the interface is compelling. The role by H-bonds (14.83% ± 6.5% and r2 = 0.61 with p<0.01) among obligatory and electrostatic energy (8.8% ± 4.77% and r2 = 0.63 with p <0.01) among non-obligatory complexes within interfaces (class A) having more non-polar residues than surface is influencing our inference. However, interfaces (class B) having less non-polar residues than surface show 1.5 fold more electrostatic energy on average. The interpretation of the interface using inter-atomic (vdW, H-bonds, electrostatic) interactions combined with inter-residue predominance (class A and class B) in relation to known function is the key to reveal its molecular principles with new challenges.

HDAC inhibitors show differential epigenetic regulation and cell survival strategies on p53 mutant colon cancer cells

Besides inactivating tumour suppressor activity in cells, mutations in p53 confer significant oncogenic functions and promote metastasis and resistance to anticancer therapy. A variety of therapies involving genetic and epigenetic signalling events regulate tumorogenesis and progression in such cases. Pharmacological interventions with HDAC inhibitors have shown promise in therapy. This work explores the changes in efficacy of the four HDAC inhibitors SAHA, MS-275, valproic acid and sodium butyrate on a panel of colon cancer cell lines - HCT116 (p53 wt), HCT116 p53-/-, HT29 and SW480 (with mutations in p53). Clonogenic assays, gene profiling and epigenetic expression done on these cells point to p53 dependent differential activity of the 4 HDAC inhibitors which also elevate methylation levels in p53 mutant cell lines. In silico modelling establishes the alterations in interactions that lead to such differential activity of valproic acid, one of the inhibitors considered for the work. Molecular Dynamic simulations carried out on the valproic acid complex ensure stability of the complex. This work establishes a p53 dependent epigenetic signalling mechanism triggered by HDAC inhibition expanding the scope of HDAC inhibitors in adjuvant therapy for p53 mutant tumours.

Specificity and stability of transient protein-protein interactions

Remarkable features that are achieved in a protein-protein complex to precise levels are stability and specificity. Deviation from the normal levels of specificity and stability, which is often caused by mutations, could result in disease conditions. Chemical nature, 3-D arrangement and dynamics of interface residues code for both specificity and stability. This article reviews roles of interfacial residues in transient protein-protein complexes. It is proposed that aside from hotspot residues conferring stability to the complex, a small set of 'rigid' residues at the interface that maintain conformation between complexed and uncomplexed forms, play a major role in conferring specificity. Exceptionally, 'super hotspot' residues, which confer both stability and specificity, are attractive sites for interaction with small molecule inhibitors.

Comparative analyses of quaternary arrangements in homo-oligomeric proteins in superfamilies: Functional implications

A comprehensive analysis of the quaternary features of distantly related homo-oligomeric proteins is the focus of the current study. This study has been performed at the levels of quaternary state, symmetry, and quaternary structure. Quaternary state and quaternary structure refers to the number of subunits and spatial arrangements of subunits, respectively. Using a large dataset of available 3D structures of biologically relevant assemblies, we show that only 53% of the distantly related homo-oligomeric proteins have the same quaternary state. Considering these homologous homo-oligomers with the same quaternary state, conservation of quaternary structures is observed only in 38% of the pairs. In 36% of the pairs of distantly related homo-oligomers with different quaternary states the larger assembly in a pair shows high structural similarity with the entire quaternary structure of the related protein with lower quaternary state and it is referred as "Russian doll effect." The differences in quaternary state and structure have been suggested to contribute to the functional diversity. Detailed investigations show that even though the gross functions of many distantly related homo-oligomers are the same, finer level differences in molecular functions are manifested by differences in quaternary states and structures. Comparison of structures of biological assemblies in distantly and closely related homo-oligomeric proteins throughout the study differentiates the effects of sequence divergence on the quaternary structures and function. Knowledge inferred from this study can provide insights for improved protein structure classification and function prediction of homo-oligomers. Proteins 2016; 84:1190-1202. © 2016 Wiley Periodicals, Inc.

PIMA: Protein-Protein interactions in Macromolecular Assembly - a web server for its Analysis and Visualization

Protein-protein interactions are essential for the basic biological machinery of the cell. This is important for processes like protein synthesis, enzyme kinetics, molecular assembly and signal transduction. A high number of macromolecular structural complexes are known due to recent advances in structure determination techniques. Therefore, it is of interest to develop an interactive tool to objectively analyze large protein complexes. Hence, we describe the development and utility of a web enabled application named ‘Protein-Protein Interaction in Macro-molecular Assembly’ (PIMA) for the analysis of large protein assemblies. The intricate details of physical interactions amongst protein subunits in a large complex are presented as simple user preferred interactive network diagrams