The origins of the evolutionary signal used to predict protein-protein interactions

Background: The correlation of genetic distances between pairs of protein sequence alignments has been used to infer protein-protein interactions. It has been suggested that these correlations are based on the signal of co-evolution between interacting proteins. However, although mutations in different proteins associated with maintaining an interaction clearly occur (particularly in binding interfaces and neighbourhoods), many other factors contribute to correlated rates of sequence evolution. Proteins in the same genome are usually linked by shared evolutionary history and so it would be expected that there would be topological similarities in their phylogenetic trees, whether they are interacting or not. For this reason the underlying species tree is often corrected for. Moreover processes such as expression level, are known to effect evolutionary rates. However, it has been argued that the correlated rates of evolution used to predict protein interaction explicitly includes shared evolutionary history; here we test this hypothesis. Results: In order to identify the evolutionary mechanisms giving rise to the correlations between interaction proteins, we use phylogenetic methods to distinguish similarities in tree topologies from similarities in genetic distances. We use a range of datasets of interacting and non-interacting proteins from Saccharomyces cerevisiae. We find that the signal of correlated evolution between interacting proteins is predominantly a result of shared evolutionary rates, rather than similarities in tree topology, independent of evolutionary divergence. Conclusions: Since interacting proteins do not have tree topologies that are more similar than the control group of non-interacting proteins, it is likely that coevolution does not contribute much to, if any, of the observed correlations.

Structural and molecular basis of interaction of HCV non-structural protein 5A with human casein kinase 1 alpha and PKR

Background: Interaction of non-structural protein 5A (NS5A) of Hepatitis C virus (HCV) with human kinases namely, casein kinase 1 alpha (ck1 alpha) and protein kinase R (PKR) have different functional implications such as regulation of viral replication and evasion of interferon induced immune response respectively. Understanding the structural and molecular basis of interactions of the viral protein with two different human kinases can be useful in developing strategies for treatment against HCV. Results: Serine 232 of NS5A is known to be phosphorylated by human ck1 alpha. A structural model of NS5A peptide containing phosphoacceptor residue Serine 232 bound to ck1 alpha has been generated using the known 3-D structures of kinase-peptide complexes. The substrate interacting residues in ck1 alpha has been identified from the model and these are found to be conserved well in the ck1 family. ck1 alpha - substrate peptide complex has also been used to understand the structural basis of association between ck1 alpha and its other viral stress induced substrate, tumour suppressor p53 transactivation domain which has a crystal structure available. Interaction of NS5A with another human kinase PKR is primarily genotype specific. NS5A from genotype 1b has been shown to interact and inhibit PKR whereas NS5A from genotype 2a/3a are unable to bind and inhibit PKR efficiently. This is one of the main reasons for the varied response to interferon therapy in HCV patients across different genotypes. Using PKR crystal structure, sequence alignment and evolutionary trace analysis some of the critical residues responsible for the interaction of NS5A 1b with PKR have been identified. Conclusions: The substrate interacting residues in ck1 alpha have been identified using the structural model of kinase substrate peptide. The PKR interacting NS5A 1b residues have also been predicted using PKR crystal structure, NS5A sequence analysis along with known experimental results. Functional significance and nature of interaction of interferon sensitivity determining region and variable region 3 of NS5A in different genotypes with PKR which was experimentally shown are also supported by the findings of evolutionary trace analysis. Designing inhibitors to prevent this interaction could enable the HCV genotype 1 infected patients respond well to interferon therapy.

Computational studies on histidine kinase protein BaeS to target multidrug-resistant Salmonella

Multidrug-resistant Salmonella serovars have been a recent concern in curing infectious diseases like typhoid. Salmonella BaeS and BaeR are the two-component system (TCS) that signal transduction proteins found to play an important role in its multidrug resistance. A canonical TCS comprises a histidine kinase (HK) and its cognate partner response regulator (RR). The general approaches for therapeutic targeting are either the catalytic ATP-binding domain or the dimerization domain HisKA (DHp) of the HK, and in some cases, the receiver or the regulatory domain of the RR proteins. Earlier efforts of identifying novel drugs targeting the signal transduction protein have not been quite successful, as it shares similar ATP-binding domain with the key house keeping gene products of the mammalian GHL family. However, targeting the dimerization domain of HisKA through which the signals are received from the RR can be a better approach. In this article, we show stepwise procedure to specifically identify the key interacting residues involved in the dimerization with the RR along with effective targeting by ligands screened from the public database. We have found a few inhibitors which target effectively the important residues for the dimerization activity. Our results suggest a plausible de novo design of better DHp domain inhibitors.

Current trends in modeling hostpathogen interactions

The rapid emergence of infectious diseases calls for immediate attention to determine practical solutions for intervention strategies. To this end, it becomes necessary to obtain a holistic view of the complex hostpathogen interactome. Advances in omics and related technology have resulted in massive generation of data for the interacting systems at unprecedented levels of detail. Systems-level studies with the aid of mathematical tools contribute to a deeper understanding of biological systems, where intuitive reasoning alone does not suffice. In this review, we discuss different aspects of hostpathogen interactions (HPIs) and the available data resources and tools used to study them. We discuss in detail models of HPIs at various levels of abstraction, along with their applications and limitations. We also enlist a few case studies, which incorporate different modeling approaches, providing significant insights into disease. (c) 2013 Wiley Periodicals, Inc.

Reconstructing the properties of dark energy using standard sirens

Future space-based gravity wave (GW) experiments such as the Big Bang Observatory (BBO), with their excellent projected, one sigma angular resolution, will measure the luminosity distance to a large number of GW sources to high precision, and the redshift of the single galaxies in the narrow solid angles towards the sources will provide the redshifts of the gravity wave sources. One sigma BBO beams contain the actual source in only 68% of the cases; the beams that do not contain the source may contain a spurious single galaxy, leading to misidentification. To increase the probability of the source falling within the beam, larger beams have to be considered, decreasing the chances of finding single galaxies in the beams. Saini et al. T.D. Saini, S.K. Sethi, and V. Sahni, Phys. Rev. D 81, 103009 (2010)] argued, largely analytically, that identifying even a small number of GW source galaxies furnishes a rough distance-redshift relation, which could be used to further resolve sources that have multiple objects in the angular beam. In this work we further develop this idea by introducing a self-calibrating iterative scheme which works in conjunction with Monte Carlo simulations to determine the luminosity distance to GW sources with progressively greater accuracy. This iterative scheme allows one to determine the equation of state of dark energy to within an accuracy of a few percent for a gravity wave experiment possessing a beam width an order of magnitude larger than BBO (and therefore having a far poorer angular resolution). This is achieved with no prior information about the nature of dark energy from other data sets such as type Ia supernovae, baryon acoustic oscillations, cosmic microwave background, etc. DOI:10.1103/PhysRevD.87.083001

Ratiometric, reversible, and parts per billion level detection of multiple toxic transition metal ions using a single probe in micellar media

We present the selective sensing of multiple transition metal ions in water using a synthetic single probe. The probe is made up of pyrene and pyridine as signaling and interacting moiety, respectively. The sensor showed different responses toward metal ions just by varying the medium of detection. In organic solvent (acetonitrile), the probe showed selective detection of Hg2+ ion. In water, the fluorescence quenching was observed with three metal ions, Cu2+, Hg2+, and Ni2+. Further, just by varying the surface charge on the micellar aggregates, the probe could detect and discriminate the above-mentioned three different toxic metal ions appropriately. In neutral micelles (Brij 58), the probe showed a selective interaction with Hg2+ ion as observed in acetonitrile medium. However, in anionic micellar medium (sodium dodecyl sulfate, SDS), the probe showed changes with both Cu2+ and Ni2+. under UV-vis absorption spectroscopy. The discrimination between these two ions was achieved by recording their emission spectra, where it showed selective quenching with Cu2+.

The evolution of complexity in social organization-A model using dominance-subordinate behavior in two social wasp species

Dominance and subordinate behaviors are important ingredients in the social organizations of group living animals. Behavioral observations on the two eusocial species Ropalidia marginata and Ropalidia cyathiformis suggest varying complexities in their social systems. The queen of R. cyathiformis is an aggressive individual who usually holds the top position in the dominance hierarchy although she does not necessarily show the maximum number of acts of dominance, while the R. marginata queen rarely shows aggression and usually does not hold the top position in the dominance hierarchy of her colony. In R. marginata, more workers are involved in dominance-subordinate interactions as compared to R. cyathiformis. These differences are reflected in the distribution of dominance-subordinate interactions among the hierarchically ranked individuals in both the species. The percentage of dominance interactions decreases gradually with hierarchical ranks in R. marginata while in R. cyathiformis it first increases and then decreases. We use an agent-based model to investigate the underlying mechanism that could give rise to the observed patterns for both the species. The model assumes, besides some non-interacting individuals, the interaction probabilities of the agents depend on their pre-differentiated winning abilities. Our simulations show that if the queen takes up a strategy of being involved in a moderate number of dominance interactions, one could get the pattern similar to R. cyathiformis, while taking up the strategy of very low interactions by the queen could lead to the pattern of R. marginata. We infer that both the species follow a common interaction pattern, while the differences in their social organization are due to the slight changes in queen as well as worker strategies. These changes in strategies are expected to accompany the evolution of more complex societies from simpler ones.

Chiral Mott insulator with staggered loop currents in the fully frustrated Bose-Hubbard model

Motivated by experiments on Josephson junction arrays in a magnetic field and ultracold interacting atoms in an optical lattice in the presence of a ``synthetic'' orbital magnetic field, we study the ``fully frustrated'' Bose-Hubbard model and quantum XY model with half a flux quantum per lattice plaquette. Using Monte Carlo simulations and the density matrix renormalization group method, we show that these kinetically frustrated boson models admit three phases at integer filling: a weakly interacting chiral superfluid phase with staggered loop currents which spontaneously break time-reversal symmetry, a conventional Mott insulator at strong coupling, and a remarkable ``chiral Mott insulator'' (CMI) with staggered loop currents sandwiched between them at intermediate correlation. We discuss how the CMI state may be viewed as an exciton condensate or a vortex supersolid, study a Jastrow variational wave function which captures its correlations, present results for the boson momentum distribution across the phase diagram, and consider various experimental implications of our phase diagram. Finally, we consider generalizations to a staggered flux Bose-Hubbard model and a two-dimensional (2D) version of the CMI in weakly coupled ladders.

Why are some galaxy discs extremely thin?

Some low-surface-brightness galaxies are known to have extremely thin stellar discs with the vertical-to-planar axes ratio 0.1 or less, often referred to as superthin galaxies. Although their existence is now known for over three decades, the physical origin of the superthin discs is still not understood. We model the vertical thickness of the stellar disc using our model of a two-component (gravitationally coupled stars and gas) disc embedded in a dark matter halo, for a bulgeless, superthin galaxy UGC 7321 which has a dense, compact halo, and is compare with a typical dwarf irregular galaxy Holmberg II which has a low-density, non-compact halo. We show that while the presence of gas does constrain the stellar disc thickness and hence its axial ratio, it is the compact dark matter halo which plays the decisive role in determining the mean distribution of stars in the vertical direction in low-luminosity bulgeless galaxies like UGC 7321, and causes the stellar disc to be superthin. Thus, the compactness of the dark matter halo significantly affects the disc structure and this could be important for the early evolution of galaxies.

Superparamagnetic behaviour and T-1, T-2 relaxivity of ZnFe2O4 nanoparticles for magnetic resonance imaging

In the present study, ZnFe2O4 nanoparticles were synthesized by the chemical co-precipitation followed by calcinations at 473 and 673K for 4h. Particle sizes obtained were 4 and 6nm for the calcination temperatures of 473 and 673K, respectively. To study the origin of system's low temperature spin dynamic behaviour, temperature dependence of susceptibility was investigated as a function of particle size and frequency. Slight increase in the grain size from 4nm at 473K to 6nm at 673K has led to a peak shift of temperature dependence of susceptibility measured at a constant frequency of 400Hz. Temperature dependence of at different frequencies also resulted in peak shift. Relaxation time dependence of peak temperature obeys a power law, which provides the fitting parameters within the range of superparamagnetic nature of the particles. Further, dependence of relaxation time and peak temperature obeys VogelFulcher law rather than NeelBrown equation demonstrating that the particles follow the behaviour of superparamagnetism of slightly interacting system. Spinlattice, T-1 and spinspin, T-2 relaxivity of proton of the water molecule in the presence of chitosan-coated superparamagnetic ZnFe2O4 nanoparticle yields the values of 0.002 and 0.360s(1)perppm.

Fermions in synthetic non-Abelian gauge potentials: rashbon condensates to novel Hamiltonians

Recent advances in the generation of synthetic gauge fields in cold atomic systems have stimulated interest in the physics of interacting bosons and fermions in them. In this paper, we discuss interacting two-component fermionic systems in uniform non-Abelian gauge fields that produce a spin-orbit interaction and uniform spin potentials. Two classes of gauge fields discussed include those that produce a Rashba spin-orbit interaction and the type of gauge fields (SM gauge fields) obtained in experiments by the Shanxi and MIT groups. For high symmetry Rashba gauge fields, a two-particle bound state exists even for a vanishingly small attractive interaction described by a scattering length. Upon increasing the strength of a Rashba gauge field, a finite density of weakly interacting fermions undergoes a crossover from a BCS like ground state to a BEC state of a new kind of boson called the rashbon whose properties are determined solely by the gauge field and not by the interaction between the fermions. The rashbon Bose-Einstein condensate (RBEC) is a quite intriguing state with the rashbon-rashbon interactions being independent of the fermion-fermion interactions (scattering length). Furthermore, we show that the RBEC has a transition temperature of the order of the Fermi temperature, suggesting routes to enhance the transition temperatures of weakly interacting superfluids by tuning the spin-orbit coupling. For the SM gauge fields, we show that in a regime of parameters, a pair of particles with finite centre-of-mass momentum is the most strongly bound. In other regimes of centre-of-mass momenta, there is no two-body bound state, but a resonance like feature appears in the scattering continuum. In the many-body setting, this results in flow enhanced pairing. Also, strongly interacting normal states utilizing the scattering resonance can be created opening the possibility of studying properties of helical Fermi liquids. This paper contains a general discussion of the physics of Feshbach resonance in a non-Abelian gauge field, where several novel features such as centre-of-mass-momentum-dependent effective interactions are shown. It is also shown that a uniform non-Abelian gauge field in conjunction with a spatial potential can be used to generate novel Hamiltonians; we discuss an explicit example of the generation of a monopole Hamiltonian.

Engineering Gold Nanoparticle Interaction by PAMAM Dendrimer

Bare faceted gold nanoparticles (AuNPs) have a tendency to aggregate through a preferred attachment of the 111] surfaces. We have used fully atomistic classical molecular dynamics simulations to obtain a quantitative estimate of this surface interaction using umbrella sampling (US) at various temperatures. To tune this surface interaction, we use polyamidoamine (PAMAM) dendrimer to coat the gold surface under various conditions. We observe a spontaneous adsorption of the protonated as well as nonprotonated PAMAM dendrimer on the AuNP surface. The adsorbed dendrimer on the nanoparticle surface strongly alters the interaction between the nanoparticles. We calculate the interaction between dendrimercoated AuNPs using US and show how the interaction between two bare faceted AuNPs can be tuned as a function of dendrimer concentration and charge (pH dependent) With appropriate choice of the dendrimer concentration and charge, two strongly interacting AuNPs can be made effectively noninteracting. Our simulation results demonstrate a strategy to tune the nanoparticle interaction, which can help in engineering self-assembly of such nanoparticles.

Heat and SDS insensitive NDK dimers are largely stabilised by hydrophobic interaction to form functional hexamer in Mycobacterium smegmatis

The primary structure and function of nucleoside diphosphate kinase (NDK), a substrate non-specific enzyme involved in the maintenance of nucleotide pools is also implicated to play pivotal roles in many other cellular processes. NDK is conserved from bacteria to human and forms a homotetramer or hexamer to exhibit its biological activity. However, the nature of the functional oligomeric form of the enzyme differs among different organisms. The functional form of NDKs from many bacterial systems, including that of the human pathogen, Mycobacterium tuberculosis (MtuNDK), is a hexamer, although some bacterial NDKs are tetrameric in nature. The present study addresses the oligomeric property of MsmNDK and how a dimer, the basic subunit of a functional hexamer, is stabilized by hydrogen bonds and hydrophobic interactions. Homology modeling was generated using the three-dimensional structure of MtuNDK as a template; the residues interacting at the monomer-monomer interface of MsmNDK were mapped. Using recombinant enzymes of wild type, catalytically inactive mutant, and monomer-monomer interactive mutants of MsmNDK, the stability of the dimer was verified under heat, SDS, low pH, and methanol. The predicted residues (Gln17, Ser24 and Glu27) were engaged in dimer formation, however the mutated proteins retained the ATPase and GTPase activity even after introducing single (MsmNDK- Q17A, MsmNDK-E27A, and MsmNDK-E27Q) and double (MsmNDK-E27A/Q17A) mutation. However, the monomer monomer interaction could be abolished using methanol, indicating the stabilization of the monomer-monomer interaction by hydrophobic interaction.

The effect of base roughness on the development of a dense granular flow down an inclined plane

The development of the flow of a granular material down an inclined plane starting from rest is studied as a function of the base roughness. In the simulations, the particles are rough frictional spheres interacting via the Hertz contact law. The rough base is made of a random configuration of fixed spheres with diameter different from the flowing particles, and the base roughness is decreased by decreasing the diameter of the base particles. The transition from an ordered to a disordered flowing state at a critical value of the base particle diameter, first reported by Kumaran and Maheshwari Phys. Fluids 24, 053302 (2012)] for particles with the linear contact model, is observed for the Hertzian contact model as well. The flow development for the ordered and disordered flows is very different. During the development of the disordered flow for the rougher base, there is shearing throughout the height. During the development of the ordered flow for the smoother base, there is a shear layer at the bottom and a plug region with no internal shearing above. In the shear layer, the particles are layered and hexagonally ordered in the plane parallel to the base, and the velocity profile is well approximated by Bagnold law. The flow develops in two phases. In the first phase, the thickness of the shear layer and the maximum velocity increase linearly in time till the shear front reaches the top. In the second phase, after the shear layer encompasses the entire flow, there is a much slower increase in the maximum velocity until the steady state is reached. (C) 2013 AIP Publishing LLC.

Stabilizing effect of electrostatic vs. aromatic interactions in diproline nucleated peptide beta-hairpins

The contribution of Tyr-His vs. Cys-His interacting pairs to the scaffold stability of (D)Pro-(L)Pro nucleated peptide beta-hairpins has been examined. We present direct evidence for the superiority of the Cys-His pairs, mediated by sulphur-imidazole interactions, as added stabilizing agents of the beta-hairpin scaffold.