PocketAlign A Novel Algorithm for Aligning Binding Sites in Protein Structures

A fundamental task in bioinformatics involves a transfer of knowledge from one protein molecule onto another by way of recognizing similarities. Such similarities are obtained at different levels, that of sequence, whole fold, or important substructures. Comparison of binding sites is important to understand functional similarities among the proteins and also to understand drug cross-reactivities. Current methods in literature have their own merits and demerits, warranting exploration of newer concepts and algorithms, especially for large-scale comparisons and for obtaining accurate residue-wise mappings. Here, we report the development of a new algorithm, PocketAlign, for obtaining structural superpositions of binding sites. The software is available as a web-service at http://proline.physicslisc.emetin/pocketalign/. The algorithm encodes shape descriptors in the form of geometric perspectives, supplemented by chemical group classification. The shape descriptor considers several perspectives with each residue as the focus and captures relative distribution of residues around it in a given site. Residue-wise pairings are computed by comparing the set of perspectives of the first site with that of the second, followed by a greedy approach that incrementally combines residue pairings into a mapping. The mappings in different frames are then evaluated by different metrics encoding the extent of alignment of individual geometric perspectives. Different initial seed alignments are computed, each subsequently extended by detecting consequential atomic alignments in a three-dimensional grid, and the best 500 stored in a database. Alignments are then ranked, and the top scoring alignments reported, which are then streamed into Pymol for visualization and analyses. The method is validated for accuracy and sensitivity and benchmarked against existing methods. An advantage of PocketAlign, as compared to some of the existing tools available for binding site comparison in literature, is that it explores different schemes for identifying an alignment thus has a better potential to capture similarities in ligand recognition abilities. PocketAlign, by finding a detailed alignment of a pair of sites, provides insights as to why two sites are similar and which set of residues and atoms contribute to the similarity.

CO adsorption on ionic Pt, Pd and Cu sites in Ce(1-x)M(x)O(2-delta) (M = Pt(2+), Pd(2+), Cu(2+))

Noble metal ion substituted CeO(2) in the form of Ce(0.98)M(0.02)O(2-delta) solid solution (where M = Pt, Pd, Cu) are the new generation catalysts with applications in three-way exhaust catalysis. While adsorption of CO on noble metals ions is well-known, adsorption of CO on noble metal ions has not been studied because creating exclusive ionic sites has been difficult. Using first-principles density functional theory (DFT) we have shown that CO gets adsorbed on the noble metal Pt(2+), Pd(2+), Cu(2+) ionic sites in the respective compounds, and the net energy of the overall system decreases. Adsorption of CO on metal ions is also confirmed by Fourier transform infrared spectroscopy (FTIR).

Auxin-Responsive OsMGH3, a Common Downstream Target of OsMADS1 and OsMADS6, Controls Rice Floret Fertility

GH3 proteins control auxin homeostasis by inactivating excess auxin as conjugates of amino acids and sugars and thereby controlling cellular bioactive auxin. Since auxin regulates many aspects of plant growth and development, regulated expression of these genes offers a mechanism to control various developmental processes. OsMGH3/OsGH3-8 is expressed abundantly in rice florets and is regulated by two related and redundant transcription factors, OsMADS1 and OsMADS6, but its contribution to flower development is not known. We functionally characterize OsMGH3 by overexpression and knock-down analysis and show a partial overlap in these phenotypes with that of mutants in OsMADS1 and OsMADS6. The overexpression of OsMGH3 during the vegetative phase affects the overall plant architecture, whereas its inflorescence-specific overexpression creates short panicles with reduced branching, resembling in part the effects of OsMADS1 overexpression. In contrast, the down-regulation of endogenous OsMGH3 caused phenotypes consistent with auxin overproduction or activated signaling, such as ectopic rooting from aerial nodes. Florets in OsMGH3 knock-down plants were affected in carpel development and pollen viability, both of which reduced fertility. Some of these floret phenotypes are similar to osmads6 mutants. Taken together, we provide evidence for the functional significance of auxin homeostasis and its transcriptional regulation during rice panicle branching and floret organ development.

Density matrix renormalization group algorithm for Bethe lattices of spin-1/2 or spin-1 sites with Heisenberg antiferromagnetic exchange

A density matrix renormalization group (DMRG) algorithm is presented for the Bethe lattice with connectivity Z = 3 and antiferromagnetic exchange between nearest-neighbor spins s = 1/2 or 1 sites in successive generations g. The algorithm is accurate for s = 1 sites. The ground states are magnetic with spin S(g) = 2(g)s, staggered magnetization that persists for large g > 20, and short-range spin correlation functions that decrease exponentially. A finite energy gap to S > S(g) leads to a magnetization plateau in the extended lattice. Closely similar DMRG results for s = 1/2 and 1 are interpreted in terms of an analytical three-site model.

C2 domain is responsible for targeting rice phosphoinositide specific phospholipase C

Phosphoinositide-specific phospholipase C (PLC) is involved in Ca2+ mediated signalling events that lead to altered cellular status. Using various sequence-analysis methods, we identified two conserved motifs in known PLC sequences. The identified motifs are located in the C2 domain of plant PLCs and are not found in any other protein. These motifs are specifically found in the Ca2+ binding loops and form adjoining beta strands. Further, we identified certain conserved residues that are highly distinct from corresponding residues of animal PLCs. The motifs reported here could be used to annotate plant-specific phospholipase C sequences. Furthermore, we demonstrated that the C2 domain alone is capable of targeting PLC to the membrane in response to a Ca2+ signal. We also showed that the binding event results from a change in the hydrophobicity of the C2 domain upon Ca2+ binding. Bioinformatic analyses revealed that all PLCs from Arabidopsis and rice lack a transmembrane domain, myristoylation and GPI-anchor protein modifications. Our bioinformatic study indicates that plant PLCs are located in the cytoplasm, the nucleus and the mitochondria. Our results suggest that there are no distinct isoforms of plant PLCs, as have been proposed to exist in the soluble and membrane associated fractions. The same isoform could potentially be present in both subcellular fractions, depending on the calcium level of the cytosol. Overall, these data suggest that the C2 domain of PLC plays a vital role in calcium signalling.

Homo- or Heterosynthon? A Crystallographic Study on a Series of New Cocrystals Derived from Pyrazinecarboxamide and Various Carboxylic Acids Equipped with Additional Hydrogen Bonding Sites

Four new cocrystals of a well studied active pharmaceutical ingredient (API), namely, pyrazinecarboxamide (PZA), with various monocarboxylic acids equipped with additional hydrogen bonding sites such as vanillic acid (VA), gallic acid (GA), 1-hydroxy-2-naphthoic acid (1HNA), and indole-2-carboxylic acid (I2CA) have been successfully prepared and characterized by FTIR, H-1 NMR, differential scanning calorimetry (DSC), and single crystal and powder X-ray diffraction (SXRD and PXRD, respectively) techniques. In the majority of the cases, preferential occurrence amide amide and acid acid homosynthons has been observed. Since the heterosynthon is energetically preferred to homosynthon, such unusual occurrence of homosynthon in these cocrystals is intriguing.

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.

Rice LHS1/OsMADS1 Controls Floret Meristem Specification by Coordinated Regulation of Transcription Factors and Hormone Signaling Pathways

SEPALLATA (SEP) MADS box transcription factors mediate floral development in association with other regulators. Mutants in five rice (Oryza sativa) SEP genes suggest both redundant and unique functions in panicle branching and floret development. LEAFY HULL STERILE1/OsMADS1, from a grass-specific subgroup of LOFSEP genes, is required for specifying a single floret on the spikelet meristem and for floret organ development, but its downstream mechanisms are unknown. Here, key pathways and directly modulated targets of OsMADS1 were deduced from expression analysis after its knockdown and induction in developing florets and by studying its chromatin occupancy at downstream genes. The negative regulation of OsMADS34, another LOFSEP gene, and activation of OsMADS55, a SHORT VEGETATIVE PHASE-like floret meristem identity gene, show its role in facilitating the spikelet-to-floret meristem transition. Direct regulation of other transcription factor genes like OsHB4 (a class III homeodomain Leu zipper member), OsBLH1 (a BEL1-like homeodomain member), OsKANADI2, OsKANADI4, and OsETTIN2 show its role in meristem maintenance, determinacy, and lateral organ development. We found that the OsMADS1 targets OsETTIN1 and OsETTIN2 redundantly ensure carpel differentiation. The multiple effects of OsMADS1 in promoting auxin transport, signaling, and auxin-dependent expression and its direct repression of three cytokinin A-type response regulators show its role in balancing meristem growth, lateral organ differentiation, and determinacy. Overall, we show that OsMADS1 integrates transcriptional and signaling pathways to promote rice floret specification and development.

Triarylborane-dipyrromethane conjugates bearing dual receptor sites: the synthesis and evaluation of the anion binding site preference

The synthesis and optical properties of four new triarylborane-dipyrromethane (TAB-DPM) conjugates (3a-d) containing dual binding sites (hydrogen bond donor and Lewis acid) have been reported. The new compounds exhibit a selective fluorogenic response towards the F-ion. The NMR titrations show that the anions bind to the TAB-DPM conjugates via the Lewis acidic triarylborane centre in preference to the hydrogen bond donor (dipyrromethane) units.

Phenotypic characters of rice landraces reveal independent lineages of short-grain aromatic indica rice

Rice landraces are lineages developed by farmers through artificial selection during the long-term domestication process. Despite huge potential for crop improvement, they are largely understudied in India. Here, we analyse a suite of phenotypic characters from large numbers of Indian landraces comprised of both aromatic and non-aromatic varieties. Our primary aim was to investigate the major determinants of diversity, the strength of segregation among aromatic and non-aromatic landraces as well as that within aromatic landraces. Using principal component analysis, we found that grain length, width and weight, panicle weight and leaf length have the most substantial contribution. Discriminant analysis can effectively distinguish the majority of aromatic from non-aromatic landraces. More interestingly, within aromatic landraces long-grain traditional Basmati and short-grain non-Basmati aromatics remain morphologically well differentiated. The present research emphasizes the general patterns of phenotypic diversity and finds out the most important characters. It also confirms the existence of very unique short-grain aromatic landraces, perhaps carrying signatures of independent origin of an additional aroma quantitative trait locus in the indica group, unlike introgression of specific alleles of the BADH2 gene from the japonica group as in Basmati. We presume that this parallel origin and evolution of aroma in short-grain indica landraces are linked to the long history of rice domestication that involved inheritance of several traits from Oryza nivara, in addition to O. rufipogon. We conclude with a note that the insights from the phenotypic analysis essentially comprise the first part, which will likely be validated with subsequent molecular analysis.

PocketMatch (version 2.0): A parallel algorithm for the detection of structural similarities between protein ligand binding-sites

Knowledge of protein-ligand interactions is essential to understand several biological processes and important for applications ranging from understanding protein function to drug discovery and protein engineering. Here, we describe an algorithm for the comparison of three-dimensional ligand-binding sites in protein structures. A previously described algorithm, PocketMatch (version 1.0) is optimised, expanded, and MPI-enabled for parallel execution. PocketMatch (version 2.0) rapidly quantifies binding-site similarity based on structural descriptors such as residue nature and interatomic distances. Atomic-scale alignments may also be obtained from amino acid residue pairings generated. It allows an end-user to compute database-wide, all-to-all comparisons in a matter of hours. The use of our algorithm on a sample dataset, performance-analysis, and annotated source code is also included.