Riboswitch Detection Using Profile Hidden Markov Models

Background  

Riboswitches are a type of noncoding RNA that regulate gene expression by switching from one structural conformation to another on ligand binding. The various classes of riboswitches discovered so far are differentiated by the ligand, which on binding induces a conformational switch. Every class of riboswitch is characterized by an aptamer domain, which provides the site for ligand binding, and an expression platform that undergoes conformational change on ligand binding. The sequence and structure of the aptamer domain is highly conserved in riboswitches belonging to the same class. We propose a method for fast and accurate identification of riboswitches using profile Hidden Markov Models (pHMM). Our method exploits the high degree of sequence conservation that characterizes the aptamer domain.     

A global optimization algorithm for protein surface alignment

Background  

A relevant problem in drug design is the comparison and recognition of protein binding sites. Binding sites recognition is generally based on geometry often combined with physico-chemical properties of the site since the conformation, size and chemical composition of the protein surface are all relevant for the interaction with a specific ligand. Several matching strategies have been designed for the recognition of protein-ligand binding sites and of protein-protein interfaces but the problem cannot be considered solved.     

Fragment-based identification of determinants of conformational and spectroscopic change at the ricin active site

Background  

Ricin is a potent toxin and known bioterrorism threat with no available antidote. The ricin A-chain (RTA) acts enzymatically to cleave a specific adenine base from ribosomal RNA, thereby blocking translation. To understand better the relationship between ligand binding and RTA active site conformational change, we used a fragment-based approach to find a minimal set of bonding interactions able to induce rearrangements in critical side-chain positions.     

Elastic network model of allosteric regulation in protein kinase PDK1

Background  

Structural switches upon binding of phosphorylated moieties underpin many signalling networks. The ligand activation is a form of allosteric modulation of the protein, where the binding site is remote from the structural change in the protein. Recently this structural switch has been elegantly demonstrated with the crystallisation of the activated form of 3-phosphoinositide-dependent protein kinase-1 (PDK1). The purpose of the present work is to determine whether the allosteric coupling in PDK1 emerges at the level of a simple coarse grained model of protein dynamics.     

Crystallographic Study of Novel Transthyretin Ligands Exhibiting Negative-Cooperativity between Two Thyroxine Binding Sites

Background

Transthyretin (TTR) is a homotetrameric serum and cerebrospinal fluid protein that transports thyroxine (T4) and retinol by binding to retinol binding protein. Rate-limiting tetramer dissociation and rapid monomer misfolding and disassembly of TTR lead to amyloid fibril formation in different tissues causing various amyloid diseases. Based on the current understanding of the pathogenesis of TTR amyloidosis, it is considered that the inhibition of amyloid fibril formation by stabilization of TTR in native tetrameric form is a viable approach for the treatment of TTR amyloidosis.

Methodology and Principal Findings

We have examined interactions of the wtTTR with a series of compounds containing various substitutions at biphenyl ether skeleton and a novel compound, previously evaluated for binding and inhibiting tetramer dissociation, by x-ray crystallographic approach. High resolution crystal structures of five ligands in complex with wtTTR provided snapshots of negatively cooperative binding of ligands in two T4 binding sites besides characterizing their binding orientations, conformations, and interactions with binding site residues. In all complexes, the ligand has better fit and more potent interactions in first T4 site i.e. (AC site) than the second T4 site (BD site). Together, these results suggest that AC site is a preferred ligand binding site and retention of ordered water molecules between the dimer interfaces further stabilizes the tetramer by bridging a hydrogen bond interaction between Ser117 and its symmetric copy.

Conclusion

Novel biphenyl ether based compounds exhibit negative-cooperativity while binding to two T4 sites which suggests that binding of only single ligand molecule is sufficient to inhibit the TTR tetramer dissociation.     

Auto-FACE: An NMR Based Binding Site Mapping Program for Fast Chemical Exchange Protein-Ligand Systems

Background

Nuclear Magnetic Resonance (NMR) spectroscopy offers a variety of experiments to study protein-ligand interactions at atomic resolution. Among these experiments, N Heteronuclear Single Quantum Correlation (HSQC) experiment is simple, less time consuming and highly informative in mapping the binding site of the ligand. The interpretation of N HSQC becomes ambiguous when the chemical shift perturbations are caused by non-specific interactions like allosteric changes and local structural rearrangement. Under such cases, detailed chemical exchange analysis based on chemical shift perturbation will assist in locating the binding site accurately.

Methodology/Principal Findings

We have automated the mapping of binding sites for fast chemical exchange systems using information obtained from N HSQC spectra of protein serially titrated with ligand of increasing concentrations. The automated program Auto-FACE (Auto-FAst Chemical Exchange analyzer) determines the parameters, e.g. rate of change of perturbation, binding equilibrium constant and magnitude of chemical shift perturbation to map the binding site residues. Interestingly, the rate of change of perturbation at lower ligand concentration is highly sensitive in differentiating the binding site residues from the non-binding site residues. To validate this program, the interaction between the protein and the ligand BH3I-1 was studied. Residues in the hydrophobic BH3 binding groove of were easily identified to be crucial for interaction with BH3I-1 from other residues that also exhibited perturbation. The geometrically averaged equilibrium constant () calculated for the residues present at the identified binding site is consistent with the values obtained by other techniques like isothermal calorimetry and fluorescence polarization assays (). Adjacent to the primary site, an additional binding site was identified which had an affinity of 3.8 times weaker than the former one. Further NMR based model fitting for individual residues suggest single site model for residues present at these binding sites and two site model for residues present between these sites. This implies that chemical shift perturbation can represent the local binding event much more accurately than the global binding event.

Conclusion/Significance

Detail NMR chemical shift perturbation analysis enabled binding site residues to be distinguished from non-binding site residues for accurate mapping of interaction site in complex fast exchange system between small molecule and protein. The methodology is automated and implemented in a program called “Auto-FACE”, which also allowed quantitative information of each interaction site and elucidation of binding mechanism.     

Identification of ATP binding residues of a protein from its primary sequence

Background  

One of the major challenges in post-genomic era is to provide functional annotations for large number of proteins arising from genome sequencing projects. The function of many proteins depends on their interaction with small molecules or ligands. ATP is one such important ligand that plays critical role as a coenzyme in the functionality of many proteins. There is a need to develop method for identifying ATP interacting residues in a ATP binding proteins (ABPs), in order to understand mechanism of protein-ligands interaction.     

Exon 6 of human JAG1 encodes a conserved structural unit

Background  

Notch signaling drives developmental processes in all metazoans. The receptor binding region of the human Notch ligand Jagged-1 is made of a DSL (Delta/Serrate/Lag-2) domain and two atypical epidermal growth factor (EGF) repeats encoded by two exons, exon 5 and 6, which are out of phase with respect to the EGF domain boundaries.     

Mannosylerythritol lipid, a yeast extracellular glycolipid, shows high binding affinity towards human immunoglobulin G

Background  

There have been many attempts to develop new materials with stability and high affinity towards immunoglobulins. Some of glycolipids such as gangliosides exhibit a high affinity toward immunoglobulins. However, it is considerably difficult to develop these glycolipids into the practical separation ligand due to their limited amounts. We thus focused our attention on the feasible use of "mannosylerythritol lipid A", a yeast glycolipid biosurfactant, as an alternative ligand for immunoglobulins, and undertook the investigation on the binding between mannosylerythritol lipid A (MEL-A) and human immunoglobulin G (HIgG).     

An automated stochastic approach to the identification of the protein specificity determinants and functional subfamilies

Background  

Recent progress in sequencing and 3 D structure determination techniques stimulated development of approaches aimed at more precise annotation of proteins, that is, prediction of exact specificity to a ligand or, more broadly, to a binding partner of any kind.     

A novel chemogenomics analysis of G protein-coupled receptors (GPCRs) and their ligands: a potential strategy for receptor de-orphanization

Background  

G protein-coupled receptors (GPCRs) represent a family of well-characterized drug targets with significant therapeutic value. Phylogenetic classifications may help to understand the characteristics of individual GPCRs and their subtypes. Previous phylogenetic classifications were all based on the sequences of receptors, adding only minor information about the ligand binding properties of the receptors. In this work, we compare a sequence-based classification of receptors to a ligand-based classification of the same group of receptors, and evaluate the potential to use sequence relatedness as a predictor for ligand interactions thus aiding the quest for ligands of orphan receptors.     

Saturation Behavior: a general relationship described by a simple second-order differential equation

Background  

The numerous natural phenomena that exhibit saturation behavior, e.g., ligand binding and enzyme kinetics, have been approached, to date, via empirical and particular analyses. This paper presents a mechanism-free, and assumption-free, second-order differential equation, designed only to describe a typical relationship between the variables governing these phenomena. It develops a mathematical model for this relation, based solely on the analysis of the typical experimental data plot and its saturation characteristics. Its utility complements the traditional empirical approaches.     

Structural motifs recurring in different folds recognize the same ligand fragments

Background  

The structural analysis of protein ligand binding sites can provide information relevant for assigning functions to unknown proteins, to guide the drug discovery process and to infer relations among distant protein folds. Previous approaches to the comparative analysis of binding pockets have usually been focused either on the ligand or the protein component. Even though several useful observations have been made with these approaches they both have limitations. In the former case the analysis is restricted to binding pockets interacting with similar ligands, while in the latter it is difficult to systematically check whether the observed structural similarities have a functional significance.     

re-TAMD: exploring interactions between H3 peptide and YEATS domain using enhanced sampling

Background

Analysis of preferred binding regions of a ligand on a protein is important for detecting cryptic binding pockets and improving the ligand selectivity.

Result

The enhanced sampling approach TAMD has been adapted to allow a ligand to unbind from its native binding site and explore the protein surface. This so-called re-TAMD procedure was then used to explore the interaction between the N terminal peptide of histone H3 and the YEATS domain. Depending on the length of the peptide, several regions of the protein surface were explored. The peptide conformations sampled during the re-TAMD correspond to peptide free diffusion around the protein surface.

Conclusions

The re-TAMD approach permitted to get information on the relative influence of different regions of the N terminal peptide of H3 on the interaction between H3 and YEATS.