Summary In recent years even with remarkable scientific advancements and significant increase of global R&D spending, drugs are frequently withdrawn from markets. This is primarily due to their side-effects or toxicities. Drug molecules often interact with multiple targets, coined as polypharmacology, and the unintended drug-target interactions could cause side-effects. Polypharmacology remains to be one of the major challenges in drug development, and it opens novel avenues to rationally design next generation of more effective but less toxic therapeutic agents. This review outlines the latest progress and challenges in polypharmacology studies.

We report a structural study on the membrane binding of ovispirin using 2D IR line shape analysis, isotope labeling, and molecular dynamics simulations. Ovispirin is an antibiotic polypeptide that binds to the surfaces of membranes as an alpha-helix. By resolving individual backbone vibrational modes (amide I) using 1-(13)C=(18)O labeling, we measured the 2D IR line shapes for 15 of the 18 residues in this peptide. A comparison of the line shapes reveals an oscillation in the inhomogeneous line width that has a period equal to that of an alpha-helix (3.6 amino acids). The periodic trend is caused by the asymmetric environment of the membrane bilayer that exposes one face of the alpha-helix to much stronger environmental electrostatic forces than the other. We compare our experimental results to 2D IR line shapes calculated using the lowest free energy structure identified from molecular dynamics simulations. These simulations predict a periodic trend similar to the experiment and lead us to conclude that ovispirin lies in the membrane just below the headgroups, is tilted, and may be kinked. Besides providing insight into the antibiotic mechanism of ovispirin, our procedure provides an infrared method for studying peptide and protein structures that relies on the natural vibrational modes of the backbone. It is a complementary method to other techniques that utilize line shapes, such as fluorescence, NMR, and ESR spectroscopies, because it does not require mutations, the spectra can be quantitatively simulated using molecular dynamics, and the technique can be applied to difficult-to-study systems like ion channels, aggregated proteins, and kinetically evolving systems.

The geographic origin and time of dispersal of Austroasiatic (AA) speakers, presently settled in south and southeast Asia, remains disputed. Two rival hypotheses, both assuming a demic component to the language dispersal, have been proposed. The first of these places the origin of Austroasiatic speakers in southeast Asia with a later dispersal to south Asia during the Neolithic, whereas the second hypothesis advocates pre-Neolithic origins and dispersal of this language family from south Asia. To test the two alternative models, this study combines the analysis of uniparentally inherited markers with 610,000 common single nucleotide polymorphism loci from the nuclear genome. Indian AA speakers have high frequencies of Y chromosome haplogroup O2a; our results show that this haplogroup has significantly higher diversity and coalescent time (17-28 thousand years ago) in southeast Asia, strongly supporting the first of the two hypotheses. Nevertheless, the results of principal component and "structure-like" analyses on autosomal loci also show that the population history of AA speakers in India is more complex, being characterized by two ancestral components-one represented in the pattern of Y chromosomal and EDAR results and the other by mitochondrial DNA diversity and genomic structure. We propose that AA speakers in India today are derived from dispersal from southeast Asia, followed by extensive sex-specific admixture with local Indian populations.

Patients with type II diabetes exhibit fibrillar deposits of human amylin protein in the pancreas. It has been proposed that amylin oligomers arising along the aggregation or fibril-formation pathways are important in the genesis of the disease. In a step toward understanding these aggregation pathways, in this work we report the conformational preferences of human amylin monomer in solution using molecular simulations and infrared experiments. In particular, we identify a stable conformer that could play a key role in aggregation. We find that amylin adopts three stable conformations: one with an α-helical segment comprising residues 9-17 and a short antiparallel β-sheet comprising residues 24-28 and 31-35; one with an extended antiparallel β-hairpin with the turn region comprising residues 20-23; and one with no particular structure. Using detailed calculations, we determine the relative stability of these various conformations, finding that the β-hairpin conformation is the most stable, followed by the α-helical conformation, and then the unstructured coil. To test our predicted structure, we calculate its infrared spectrum in the amide I stretch regime, which is sensitive to secondary structure through vibrational couplings and linewidths, and compare it to experiment. We find that theoretically predicted spectra are in good agreement with the experimental line shapes presented herein. The implications of the monomer secondary structures on its aggregation pathway and on its interaction with cell membranes are discussed.

Ab initio (HF, MP2, and CCSD(T)) and DFT (B3LYP) calculations were done in modeling the cation (H(+), Li(+), Na(+), K(+), Ca(2+), Mg(2+), NH(4)(+), and NMe(4)(+)) interaction with aromatic side chain motifs of four amino acids (viz., phenylalanine, tyrosine, tryptophan and histidine). As the metal ion approaches the pi-framework of the model systems, they form strongly bound cation-pi complexes, where the metal ion is symmetrically disposed with respect to all ring atoms. In contrast, proton prefers to bind covalently to one of the ring carbons. The NH(4)(+) and NMe(4)(+) ions have shown N-H...pi interaction and C-H...pi interaction with the aromatic motifs. The interaction energies of N-H...pi and C-H...pi complexes are higher than hydrogen bonding interactions; thus, the orientation of aromatic side chains in protein is effected in the presence of ammonium ions. However, the regioselectivity of metal ion complexation is controlled by the affinity of the site of attack. In the imidazole unit of histidine the ring nitrogen has much higher metal ion (as well as proton) affinity as compared to the pi-face, facilitating the in-plane complexation of the metal ions. The interaction energies increase in the order of 1-M < 2-M < 3-M < 4-M < 5-M for all the metal ion considered. Similarly, the complexation energies with the model systems decrease in the following order: Mg(2+) > Ca(2+) > Li(+) > Na(+) > K(+) congruent with NH(4)(+) > NMe(4)(+). The variation of the bond lengths and the extent of charge transfer upon complexation correlate well with the computed interaction energies.