Hemocyanin gene family evolution in spiders (Araneae), with implications for phylogenetic relationships and divergence times in the infraorder Mygalomorphae

Hemocyanins are multimeric copper-containing hemolymph proteins involved in oxygen binding and transport in all major arthropod lineages. Most arachnids have seven primary subunits (encoded by paralogous genes a–g), which combine to form a 24-mer (4 × 6) quaternary structure. Within some spider lineages, however, hemocyanin evolution has been a dynamic process with extensive paralog duplication and loss. We have obtained hemocyanin gene sequences from numerous representatives of the spider infraorders Mygalomorphae and Araneomorphae in order to infer the evolution of the hemocyanin gene family and estimate spider relationships using these conserved loci. Our hemocyanin gene tree is largely consistent with the previous hypotheses of paralog relationships based on immunological studies, but reveals some discrepancies in which paralog types have been lost or duplicated in specific spider lineages. Analyses of concatenated hemocyanin sequences resolved deep nodes in the spider phylogeny and recovered a number of clades that are supported by other molecular studies, particularly for mygalomorph taxa. The concatenated data set is also used to estimate dates of higher-level spider divergences and suggests that the diversification of extant mygalomorphs preceded that of extant araneomorphs. Spiders are diverse in behavior and respiratory morphology, and our results are beneficial for comparative analyses of spider respiration. Lastly, the conserved hemocyanin sequences allow for the inference of spider relationships and ancient divergence dates.     

The Genetic Algorithm and the Conformational Search of Polypeptides and Proteins

Abstract

The genetic algorithm is a technique of function optimization derived from the principles of evolutionary theory. We have adapted it to perform conformational search on polypeptides and proteins. The algorithm was first tested on several small polypeptides and the 46 amino acid protein crambin under the AMBER potential energy function. The probable global minimum conformations of the polypeptides were located 90% of the time and a non-native conformation of crambin was located that was 150kcal/mol lower in potential energy than the minimized crystal structure conformation. Next, we used a knowledge-based potential function to predict the structures of melittin, pancreatic polypeptide, and crambin. A 2.31 Å ΔRMS conformation of melittin and a 5.33 Å ΔRMS conformation of pancreatic polypeptide were located by genetic algorithm-based conformational search under the knowledge-based potential function. Although the ΔRMS of pancreatic polypeptide was somewhat high, most of the secondary structure was correct. The secondary structure of crambin was predicted correctly, but the potential failed to promote packing interactions. Finally, we tested the packing aspects of our potential function by attempting to predict the tertiary structure of cytochrome b ₅₆₂ given correct secondary structure as a constraint. The final predicted conformation of cytochrome b ₅₆₂ was an almost completely extended continuous helix which indicated that the knowledge-based potential was useless for tertiary structure prediction. This work serves as a warning against testing potential functions designed for tertiary structure prediction on small proteins.     

The design of a peptide sequence to inhibit HIV replication: a search algorithm combining Monte Carlo and self-consistent mean field techniques

We developed a search algorithm combining Monte Carlo (MC) and self-consistent mean field techniques to evolve a peptide sequence that has good binding capability to the anticodon stem and loop (ASL) of human lysine tRNA species, tRNA^Lys3, with the ultimate purpose of breaking the replication cycle of human immunodeficiency virus-1. The starting point is the 15-amino-acid sequence, RVTHHAFLGAHRTVG, found experimentally by Agris and co-workers to bind selectively to hypermodified tRNA^Lys3. The peptide backbone conformation is determined via atomistic simulation of the peptide-ASL^Lys3 complex and then held fixed throughout the search. The proportion of amino acids of various types (hydrophobic, polar, charged, etc.) is varied to mimic different peptide hydration properties. Three different sets of hydration properties were examined in the search algorithm to see how this affects evolution to the best-binding peptide sequences. Certain amino acids are commonly found at fixed sites for all three hydration states, some necessary for binding affinity and some necessary for binding specificity. Analysis of the binding structure and the various contributions to the binding energy shows that: 1) two hydrophilic residues (asparagine at site 11 and the cysteine at site 12) “recognize” the ASL^Lys3 due to the VDW energy, and thereby contribute to its binding specificity and 2) the positively charged arginines at sites 4 and 13 preferentially attract the negatively charged sugar rings and the phosphate linkages, and thereby contribute to the binding affinity.     

Methylerythritol phosphate pathway to isoprenoids: Kinetic modeling and in silico enzyme inhibitions in Plasmodium falciparum

The methylerythritol phosphate (MEP) pathway of Plasmodium falciparum (P. falciparum) has become an attractive target for anti-malarial drug discovery. This study describes a kinetic model of this pathway, its use in validating 1-deoxy-d-xylulose 5-phosphate reductoisomerase (DXR) as drug target from the systemic perspective, and additional target identification, using metabolic control analysis and in silico inhibition studies. In addition to DXR, 1-deoxy-d-xylulose 5-phosphate synthase (DXS) can be targeted because it is the first enzyme of the pathway and has the highest flux control coefficient followed by that of DXR. In silico inhibition of both enzymes caused large decrement in the pathway flux. An added advantage of targeting DXS is its influence on vitamin B1 and B6 biosynthesis. Two more potential targets, 2-C-methyl-d-erythritol 2,4-cyclodiphosphate synthase and 1-hydroxy-2-methyl-2-(E)-butenyl 4-diphosphate synthase, were also identified. Their inhibition caused large accumulation of their substrates causing instability of the system.     

The relationship between visual working memory and attention: retention of precise colour information in the absence of effects on perceptual selection

We examined the conditions under which a feature value in visual working memory (VWM) recruits visual attention to matching stimuli. Previous work has suggested that VWM supports two qualitatively different states of representation: an active state that interacts with perceptual selection and a passive (or accessory) state that does not. An alternative hypothesis is that VWM supports a single form of representation, with the precision of feature memory controlling whether or not the representation interacts with perceptual selection. The results of three experiments supported the dual-state hypothesis. We established conditions under which participants retained a relatively precise representation of a parcticular colour. If the colour was immediately task relevant, it reliably recruited attention to matching stimuli. However, if the colour was not immediately task relevant, it failed to interact with perceptual selection. Feature maintenance in VWM is not necessarily equivalent with feature-based attentional selection.     

Development of a pharmacophore model for the catecholamine release-inhibitory peptide catestatin: Virtual screening and functional testing identify novel small molecule therapeutics of hypertension

The endogenous catecholamine release-inhibitory peptide catestatin (CST) regulates events leading to hypertension and cardiovascular disease. Earlier we studied the structure of CST by NMR, molecular modeling, and amino acid scanning mutagenesis. That structure has now been exploited for elucidation of interface pharmacophores that mediate binding of CST to its target, with consequent secretory inhibition. Designed pharmacophore models allowed screening of 3D structural domains. Selected compounds were tested on both cultured catecholaminergic cells and an in vivo model of hypertension; in each case, the candidates showed substantial mimicry of native CST actions, with preserved or enhanced potency and specificity. The approach and compounds have thus enabled rational design of novel drug candidates for treatment of hypertension or autonomic dysfunction.     

InDiaMed: A Comprehensive Database of Indian Medicinal plants for Diabetes

According to International Diabetes Federation (IDF), India has 62.4 million people with diabetes and by 2030 it is predicted that the number will rise to 100 million. Studies claim that there are around 410 experimentally proven Indian medicinal plants which have anti-diabetic activity, of which the mechanism of action of 109 plants has been elucidated or reported. So, the need of the hour is to explore the claims of Indian medicinal flora and open up the facets of many Indian plants which are being examined for their beneficial role in diabetes. So, we created a database (InDiaMed) of Indian medicinal plants that captures their role in anti-diabetic activity. InDiaMed''s features include chemical, pharmacological, biochemical and geographical information of the medicinal plant, scientifically relevant information of the plant, and the coherent research done on it in the field of diabetes. The database also includes the list of poly-herbal formulations which are used for treatment of diabetes in India.

Availability

http://www.indiamed.info