Mechanistic insights into G-protein activation via phosphorylation mediated non-canonical pathway.

Activation of heterotrimeric G-proteins (Gαßγ) downstream to receptor tyrosine kinases (RTKs) is a well-established crosstalk between the signaling pathways mediated by G-protein coupled receptors (GPCRs) and RTKs. While GPCR serves as a guanine exchange factor (GEF) in the canonical activation of Gα that facilitates the exchange of GDP for GTP, the mechanism through which RTK phosphorylations induce Gα activation remains unclear. Recent experimental studies revealed that the epidermal growth factor receptor (EGFR), a well-known RTK, phosphorylates the helical domain tyrosine residues Y154 and Y155 and accelerates the GDP release from the Gαi3, a subtype of Gα-protein. Using well-tempered metadynamics and extensive unbiased molecular dynamics simulations, we captured the GDP release event and identified the intermediates between bound and unbound states through Markov state models. In addition to weakened salt bridges at the domain interface, phosphorylations induced the unfolding of helix αF, which contributed to increased flexibility near the hinge region, facilitating a greater distance between domains in the phosphorylated Gαi3. Although the larger domain separation in the phosphorylated system provided an unobstructed path for the nucleotide, the accelerated release of GDP was attributed to increased fluctuations in several conserved regions like P-loop, switch 1, and switch 2. Overall, this study provides atomistic insights into the activation of G-proteins induced by RTK phosphorylations and identifies the specific structural motifs involved in the process. The knowledge gained from the study could establish a foundation for targeting non-canonical signaling pathways and developing therapeutic strategies against the ailments associated with dysregulated G-protein signaling.

Biased Signaling in Mutated Variants of ß2-Adrenergic Receptor: Insights from Molecular Dynamics Simulations.

The molecular basis of receptor bias in G protein-coupled receptors (GPCRs) caused by mutations that preferentially activate specific intracellular transducers over others remains poorly understood. Two experimentally identified biased variants of ß2-adrenergic receptors (ß2AR), a prototypical GPCR, are a triple mutant (T68F, Y132A, and Y219A) and a single mutant (Y219A); the former bias the receptor toward the ß-arrestin pathway by disfavoring G protein engagement, while the latter induces G protein signaling explicitly due to selection against GPCR kinases (GRKs) that phosphorylate the receptor as a prerequisite of ß-arrestin binding. Though rigorous characterizations have revealed functional implications of these mutations, the atomistic origin of the observed transducer selectivity is not clear. In this study, we investigated the allosteric mechanism of receptor bias in ß2AR using microseconds of all-atom Gaussian accelerated molecular dynamics (GaMD) simulations. Our observations reveal distinct rearrangements in transmembrane helices, intracellular loop 3, and critical residues R1313.50 and Y3267.53 in the conserved motifs D(E)RY and NPxxY for the mutant receptors, leading to their specific transducer interactions. Moreover, partial dissociation of G protein from the receptor core is observed in the simulations of the triple mutant in contrast to the single mutant and wild-type receptor. The reorganization of allosteric communications from the extracellular agonist BI-167107 to the intracellular receptor-transducer interfaces drives the conformational rearrangements responsible for receptor bias in the single and triple mutants. The molecular insights into receptor bias of ß2AR presented here could improve the understanding of biased signaling in GPCRs, potentially opening new avenues for designing novel therapeutics with fewer side-effects and superior efficacy.

Molecular Insights into Phosphorylation-Induced Allosteric Conformational Changes in a ß2-Adrenergic Receptor.

The large conformational flexibility of G protein-coupled receptors (GPCRs) has been a puzzle in structural and pharmacological studies for the past few decades. Apart from structural rearrangements induced by ligands, enzymatic phosphorylations by GPCR kinases (GRKs) at the carboxy-terminal tail (C-tail) of a GPCR also make conformational alterations to the transmembrane helices and facilitates the binding of one of its transducer proteins named ß-arrestin. The phosphorylation-induced conformational transition of the receptor that causes specific binding to ß-arrestin but prevents the association of other transducers such as G proteins lacks atomistic understanding and is elusive to experimental studies. Using microseconds of all-atom conventional and Gaussian accelerated molecular dynamics (GaMD) simulations, we investigate the allosteric mechanism of phosphorylation induced-conformational changes in ß2-adrenergic receptor, a well-characterized GPCR model system. Free energy profiles reveal that the phosphorylated receptor samples a new conformational state in addition to the canonical active state corroborating with recent nuclear magnetic resonance experimental findings. The new state has a smaller intracellular cavity that is likely to accommodate ß-arrestin better than G protein. Using contact map and inter-residue interaction energy calculations, we found the phosphorylated C-tail adheres to the cytosolic surface of the transmembrane domain of the receptor. Transfer entropy calculations show that the C-tail residues drive the correlated motions of TM residues, and the allosteric signal is relayed via several residues at the cytosolic surface. Our results also illustrate how the redistribution of inter-residue nonbonding interaction couples with the allosteric communication from the phosphorylated C-tail to the transmembrane. Atomistic insight into phosphorylation-induced ß-arrestin specific conformation is therapeutically important to design drugs with higher efficacy and fewer side effects. Our results, therefore, open novel opportunities to fine-tune ß-arrestin bias in GPCR signaling.

Genome-scale metabolic reconstruction and metabolic versatility of an obligate methanotroph Methylococcus capsulatus str. Bath.

The increase in greenhouse gases with high global warming potential such as methane is a matter of concern and requires multifaceted efforts to reduce its emission and increase its mitigation from the environment. Microbes such as methanotrophs can assist in methane mitigation. To understand the metabolic capabilities of methanotrophs, a complete genome-scale metabolic model (GSMM) of an obligate methanotroph, Methylococcus capsulatus str. Bath was reconstructed. The model contains 535 genes, 899 reactions and 865 metabolites and is named iMC535. The predictive potential of the model was validated using previously-reported experimental data. The model predicted the Entner-Duodoroff pathway to be essential for the growth of this bacterium, whereas the Embden-Meyerhof-Parnas pathway was found non-essential. The performance of the model was simulated on various carbon and nitrogen sources and found that M. capsulatus can grow on amino acids. The analysis of network topology of the model identified that six amino acids were in the top-ranked metabolic hubs. Using flux balance analysis, 29% of the metabolic genes were predicted to be essential, and 76 double knockout combinations involving 92 unique genes were predicted to be lethal. In conclusion, we have reconstructed a GSMM of a methanotroph Methylococcus capsulatus str. Bath. This is the first high quality GSMM of a Methylococcus strain which can serve as an important resource for further strain-specific models of the Methylococcus genus, as well as identifying the biotechnological potential of M. capsulatus Bath.

IL17eScan: A Tool for the Identification of Peptides Inducing IL-17 Response.

IL-17 cytokines are pro-inflammatory cytokines and are crucial in host defense against various microbes. Induction of these cytokines by microbial antigens has been investigated in the case of ischemic brain injury, gingivitis, candidiasis, autoimmune myocarditis, etc. In this study, we have investigated the ability of amino acid sequence of antigens to induce IL-17 response using machine-learning approaches. A total of 338 IL-17-inducing and 984 IL-17 non-inducing peptides were retrieved from Immune Epitope Database. 80% of the data were randomly selected as training dataset and rest 20% as validation dataset. To predict the IL-17-inducing ability of peptides/protein antigens, different sequence-based machine-learning models were developed. The performance of support vector machine (SVM) and random forest (RF) was compared with different parameters to predict IL-17-inducing epitopes (IIEs). The dipeptide composition-based SVM-model displayed an accuracy of 82.4% with Matthews correlation coefficient = 0.62 at polynomial (t = 1) kernel on 10-fold cross-validation and outperformed RF. Amino acid residues Leu, Ser, Arg, Asn, and Phe and dipeptides LL, SL, LK, IL, LI, NL, LR, FK, SF, and LE are abundant in IIEs. The present tool helps in the identification of IIEs using machine-learning approaches. The induction of IL-17 plays an important role in several inflammatory diseases, and identification of such epitopes would be of great help to the immunologists. It is freely available at http://metagenomics.iiserb.ac.in/IL17eScan/ and http://metabiosys.iiserb.ac.in/IL17eScan/.

Prediction of anti-inflammatory proteins/peptides: an insilico approach.

BACKGROUND: The current therapy for inflammatory and autoimmune disorders involves the use of nonspecific anti-inflammatory drugs and other immunosuppressant, which are often accompanied with potential side effects. As an alternative therapy, anti-inflammatory peptides are recently being exploited as anti-inflammatory agents for treatment of various inflammatory diseases such as Alzheimer's disease and rheumatoid arthritis. Thus, understanding the correlation between amino acid sequence and its potential anti-inflammatory property is of great importance for the discovery of novel and efficient anti-inflammatory peptide-based therapeutics. METHODS: In this study, we have developed a prediction tool for the classification of peptides as anti-inflammatory epitopes or non anti-inflammatory epitopes. The training was performed using experimentally validated epitopes obtained from Immune epitope database and analysis resource database. Different sequence-based features and their hybrids with motif information were employed for development of support vector machine-based machine learning models. Similarly, machine learning models were also constructed using random forest. RESULTS: The composition and terminal residue conservation analysis of peptides revealed the dominance of leucine, serine, tyrosine and arginine residues in anti-inflammatory epitopes as compared to non anti-inflammatory epitopes. Similarly, the anti-inflammatory epitopes specific motifs were found to be rich in hydrophobic and polar residues. The hybrid of tripeptide composition-based support vector machine model and motif yielded the best performance on 10-fold cross validation with an accuracy of 78.1% and MCC of 0.58. The same displayed an accuracy of 72% and MCC of 0.45 on validation dataset, rejecting any possibility of over-fitting. The tripeptide composition-based random forest model displayed an accuracy of 0.8 and MCC of 0.59 on 10-fold cross validation, however, the accuracy (0.68) and MCC (0.31) was lower as compared to support vector machine model on validation dataset. Thus, the support vector machine model is implemented as the default model and an additional option of using the random forest model is provided. CONCLUSION: The prediction models along with tools for epitope mapping and similarity search have been provided as a web server which is freely accessible at http://metagenomics.iiserb.ac.in/antiinflam/ .

ProInflam: a webserver for the prediction of proinflammatory antigenicity of peptides and proteins.

BACKGROUND: Proinflammatory immune response involves a complex series of molecular events leading to inflammatory reaction at a site, which enables host to combat plurality of infectious agents. It can be initiated by specific stimuli such as viral, bacterial, parasitic or allergenic antigens, or by non-specific stimuli such as LPS. On counter with such antigens, the complex interaction of antigen presenting cells, T cells and inflammatory mediators like IL1α, IL1ß, TNFα, IL12, IL18 and IL23 lead to proinflammatory immune response and further clearance of infection. In this study, we have tried to establish a relation between amino acid sequence of antigen and induction of proinflammatory response. RESULTS: A total of 729 experimentally-validated proinflammatory and 171 non-proinflammatory epitopes were obtained from IEDB database. The A, F, I, L and V amino acids and AF, FA, FF, PF, IV, IN dipeptides were observed as preferred residues in proinflammatory epitopes. Using the compositional and motif-based features of proinflammatory and non-proinflammatory epitopes, we have developed machine learning-based models for prediction of proinflammatory response of peptides. The hybrid of motifs and dipeptide-based features displayed best performance with MCC = 0.58 and an accuracy of 87.6 %. CONCLUSION: The amino acid sequence-based features of peptides were used to develop a machine learning-based prediction tool for the prediction of proinflammatory epitopes. This is a unique tool for the computational identification of proinflammatory peptide antigen/candidates and provides leads for experimental validations. The prediction model and tools for epitope mapping and similarity search are provided as a comprehensive web server which is freely available at http://metagenomics.iiserb.ac.in/proinflam/ and http://metabiosys.iiserb.ac.in/proinflam/ .