Gopichandran Sowmya

Linking structural features of protein complexes and biological function

Protein–protein interaction (PPI) establishes the central basis for complex cellular networks in a biological cell. Association of proteins with other proteins occurs at varying affinities, yet with a high degree of specificity. PPIs lead to diverse functionality such as catalysis, regulation, signaling, immunity, and inhibition, playing a crucial role in functional genomics. The molecular principle of such interactions is often elusive in nature. Therefore, a comprehensive analysis of known protein complexes from the Protein Data Bank (PDB) is essential for the characterization of structural interface features to determine structure–function relationship. Thus, we analyzed a nonredundant dataset of 278 heterodimer protein complexes, categorized into major functional classes, for distinguishing features. Interestingly, our analysis has identified five key features (interface area, interface polar residue abundance, hydrogen bonds, solvation free energy gain from interface formation, and binding energy) that are discriminatory among the functional classes using Kruskal‐Wallis rank sum test. Significant correlations between these PPI interface features amongst functional categories are also documented. Salt bridges correlate with interface area in regulator‐inhibitors (r = 0.75). These representative features have implications for the prediction of potential function of novel protein complexes. The results provide molecular insights for better understanding of PPIs and their relation to biological functions.

Insights from the structural analysis of protein heterodimer interfaces

Protein heterodimer complexes are often involved in catalysis, regulation, assembly, immunity and inhibition. This involves the formation of stable interfaces between the interacting partners. Hence, it is of interest to describe heterodimer interfaces using known structural complexes. We use a non-redundant dataset of 192 heterodimer complex structures from the protein databank (PDB) to identify interface residues and describe their interfaces using amino-acids residue property preference. Analysis of the dataset shows that the heterodimer interfaces are often abundant in polar residues. The analysis also shows the presence of two classes of interfaces in heterodimer complexes. The first class of interfaces (class A) with more polar residues than core but less than surface is known. These interfaces are more hydrophobic than surfaces, where protein-protein binding is largely hydrophobic. The second class of interfaces (class B) with more polar residues than core and surface is shown. These interfaces are more polar than surfaces, where binding is mainly polar. Thus, these findings provide insights to the understanding of protein-protein interactions.

Characterization of the interaction between heterodimeric αvβ6 integrin and urokinase plasminogen activator receptor (uPAR) using functional proteomics

Urokinase plasminogen activator receptor (uPAR) and the epithelial integrin αvβ6 are thought to individually play critical roles in cancer metastasis. These observations have been highlighted by the recent discovery (by proteomics) of an interaction between these two molecules, which are also both implicated in the epithelial-mesenchymal transition (EMT) that facilitates escape of cells from tissue barriers and is a common signature of cancer metastases. In this study, orthogonal in cellulo and in vitro functional proteomic approaches were used to better characterize the uPAR·αvβ6 interaction. Proximity ligation assays (PLA) confirmed the uPAR·αvβ6 interaction on OVCA429 (ovarian cancer line) and four different colon cancer cell lines including positive controls in cells with de novo β6 subunit expression. PLA studies were then validated using peptide arrays, which also identified potential physical sites of uPAR interaction with αvβ6, as well as verifying interactions with other known uPAR ligands (e.g., uPA, vitronectin) and individual integrin subunits (i.e., αv, β1, β3, and β6 alone). Our data suggest that interaction with uPAR requires expression of the complete αβ heterodimer (e.g., αvβ6), not individual subunits (i.e., αv, β1, β3, or β6). Finally, using in silico structural analyses in concert with these functional proteomics studies, we propose and demonstrate that the most likely unique sites of interaction between αvβ6 and uPAR are located in uPAR domains II and III.

Protein-protein interactions and prediction: a comprehensive overview.

Molecular function in cellular processes is governed by protein-protein interactions (PPIs) within biological networks. Selective yet specific association of these protein partners contributes to diverse functionality such as catalysis, regulation, assembly, immunity, and inhibition in a cell. Therefore, understanding the principles of protein-protein association has been of immense interest for several decades. We provide an overview of the experimental methods used to determine PPIs and the key databases archiving this information. Structural and functional information of existing protein complexes confers knowledge on the principles of PPI, based on which a classification scheme for PPIs is then introduced. Obtaining high-quality non-redundant datasets of protein complexes for interaction characterisation is an essential step towards deciphering their underlying binding principles. Analysis of physicochemical features and their documentation has enhanced our understanding of the molecular basis of protein-protein association. We describe the diverse datasets created/collected by various groups and their key findings inferring distinguishing features. The currently available interface databases and prediction servers have also been compiled.

A site for direct integrin αvβ6·uPAR interaction from structural modelling and docking

Integrin αvβ6 is an epithelially-restricted heterodimeric transmembrane glycoprotein, known to interact with the urokinase plasminogen activating receptor (uPAR), playing a critical role in cancer progression. While the X-ray crystallographic structures of segments of other integrin heterodimers are known, there is no structural information for the complete αvβ6 integrin to assess its direct interaction with uPAR. We have performed structural analysis of αvβ6·uPAR interactions using model data with docking simulations to pinpoint their interface, in accord with earlier reports of the β-propeller region of integrin α-chain interacting with uPAR. Interaction of αvβ6·uPAR was demonstrated by our previous study using immunoprecipitation coupled with proteomic analysis by mass spectrometry. Recently this interaction was validated with proximity ligation assays and peptide arrays. The data suggested that two potential peptide regions from domain II and one peptide region from domain III of uPAR, interact with αvβ6 integrin. Only the peptide region from domain III is consistent with the three-dimensional interaction site proposed in this study. The molecular basis of integrin αvβ6·uPAR binding using structural data is discussed for its implications as a potential therapeutic target in cancer management.

HIV-1 envelope accessible surface and polarity: clade, blood, and brain.

The human immunodeficiency virus type-1 (HIV-1) gp160 (gp120-gp41 complex) trimer envelope (ENV) protein is a potential vaccine candidate for HIV/AIDS. HIV-1 vaccine development has been problematic and charge polarity as well as sequence variation across clades may relate to the difficulties. Further obstacles are caused by sequence variation between blood and brain-derived sequences, since the brain is a separate compartment for HIV-1 infection. We utilize a threedimensional residue measure of solvent exposure, accessible surface area (ASA), which shows that major segments of gp120 and gp41 known structures are solvent exposed across clades. We demonstrate a large percent sequence polarity for solvent exposed residues in gp120 and gp41. The range of sequence polarity varies across clades, blood, and brain from different geographical locations. Regression analysis shows that blood and brain gp120 and gp41 percent sequence polarity range correlate with mean Shannon entropy. These results point to the use of protein modifications to enhance HIV-1 ENV vaccines across multiple clades, blood, and brain. It should be noted that we do not address the issue of protein glycosylation here; however, this is an important issue for vaccine design and development. Abbreviations HIV-1 - human immunodeficiency virus type 1, AIDS - acquired immunodeficiency syndrome, ENV - envelope, gp160 - 160,000d glycoprotein, gp120 - 120,000d glycoprotein, gp41 - 41,000d glycoprotein, LANL - Los Alamos National Laboratories, PDB - Protein Data Bank, HVTN - STEP HIV vaccine trial, AA - amino acids, MSA - multiple sequence alignment, ASA - accessible surface area, SNPs- single nucleotide polymorphisms, HAART - Highly Active Antiretroviral Therapy, CCR5 - C-C chemokine receptor type 5, CNS - central nervous system, HIVE - HIV encephalitis, P - polarity, NP - non-polarity, CTL - cytotoxic T lymphocyte, NIAID - National Institute of Allergy and Infectious Diseases.

Discrete structural features among interface residue-level classes

BackgroundProtein-protein interaction (PPI) is essential for molecular functions in biological cells. Investigation on protein interfaces of known complexes is an important step towards deciphering the driving forces of PPIs. Each PPI complex is specific, sensitive and selective to binding. Therefore, we have estimated the relative difference in percentage of polar residues between surface and the interface for each complex in a non-redundant heterodimer dataset of 278 complexes to understand the predominant forces driving binding.ResultsOur analysis showed ~60% of protein complexes with surface polarity greater than interface polarity (designated as class A). However, a considerable number of complexes (~40%) have interface polarity greater than surface polarity, (designated as class B), with a significantly different p-value of 1.66E-45 from class A. Comprehensive analyses of protein complexes show that interface features such as interface area, interface polarity abundance, solvation free energy gain upon interface formation, binding energy and the percentage of interface charged residue abundance distinguish among class A and class B complexes, while electrostatic visualization maps also help differentiate interface classes among complexes.ConclusionsClass A complexes are classical with abundant non-polar interactions at the interface; however class B complexes have abundant polar interactions at the interface, similar to protein surface characteristics. Five physicochemical interface features analyzed from the protein heterodimer dataset are discriminatory among the interface residue-level classes. These novel observations find application in developing residue-level models for protein-protein binding prediction, protein-protein docking studies and interface inhibitor design as drugs.

Interaction modes at protein hetero-dimer interfaces.

Hetero dimer (different monomers) interfaces are involved in catalysis and regulation through the formation of interface active sites. This is critical in cell and molecular biology events. The physical and chemical factors determining the formation of the interface active sites is often large in numbers. The combined role of interacting features is frequently combinatorial and additive in nature. Therefore, it is important to determine the physical and chemical features of such interactions. A number of such features have been documented in literature since 1975. However, the use of such interaction features in the prediction of interaction partners and sites given their sequences is still a challenge. In a non-redundant dataset of 156 hetero-dimer structures determined by X-ray crystallography, the interacting partners are often varying in size and thus, size variation between subunits is an important factor in determining the mode of interface formation. The size of protein subunits interacting are either small-small, largelarge, medium-medium, large-small, large-medium and small-medium. It should also be noted that the interface formed between subunits have physical interactions at N terminal (N), C terminal (C) and middle (M) region of the protein with reference to their sequences in one dimension. These features are believed to have application in the prediction of interaction partners and sites from sequences. However, the use of such features for interaction prediction from sequence is not currently clear.

Short Oligopeptide T-Cell Epitopes in HIV-1/AIDS Vaccine Development: Current Status, Design, Promises and Challenges

The development of a HIV-1/AIDS vaccine is a technological challenge over the last two decades due to HIV-1 viral diversity, sequence similarity between human and HIV proteins and HIV-1/human host molecular mimicry. The challenges and complexities associated with an ENV-gp120 subunit-based recombinant protein vaccine are discussed in detail for over a decade. Clinical trials using RV144 (Env-gp120, Gag and Pro) candidate vaccine are marginally promising yet unsuccessful. This is attributed to viral mutation, protein size, protein type (unstable membrane protein) and the native structure of glycosylated gp160 as a trimer complex spike on the envelope (ENV). The structural resolution of gp120 trimer spike pre-fusion is highly promising for recombinant vaccine development. However, the production of the native spike trimer complex is a biotechnological challenge for protein engineering, refolding and assembly. Hence, there is a need for the development of an alternative vaccine type for HIV-1/AIDS.

Short Peptide Vaccine Design and Development: Promises and Challenges

Vaccine development for viral diseases is a challenge where subunit vaccines are often ineffective. Therefore, the need for alternative solutions is crucial. Thus, short peptide vaccine candidates promise effective answers under such circumstances. Short peptide vaccine candidates are linear T-cell epitopes (antigenic determinants that are recognized by the immune system) that specifically function by binding human leukocyte antigen (HLA) alleles of different ethnicities (including Black, Caucasian, Oriental, Hispanic, Pacific Islander, American Indian, Australian aboriginal, and mixed ethnicities). The population-specific allele-level HLA sequence data in the public IMGT/HLA database contains approximately 12542 nomenclature defined class I (9437) and class II (3105) HLA alleles as of March 2015 present in several ethnic populations. The bottleneck in short peptide vaccine design and development is HLA polymorphism on the one hand and viral diversity on the other hand. Hence, a crucial step in its design and development is HLA allele-specific binding of short antigen peptides. This is usually combinatorial and computationally labor intensive. Mathematical models utilizing structure-defined pockets are currently available for class I and class II HLA-peptide-binding peptides. Frameworks have been developed to design protocols to identify the most feasible short peptide cocktails as vaccine candidates with superantigen properties among known HLA supertypes. This approach is a promising solution to develop new viral vaccines given the current advancement in T-cell immuno-informatics, yet challenging in terms of prediction efficiency and protocol development.