Inferring molecular function: contributions from functional linkages

In the current era of high-throughput sequencing and structure determination, functional annotation has become a bottleneck in biomedical science. Here, we show that automated inference of molecular function using functional linkages among genes increases the accuracy of functional assignments by >= 8% and enriches functional descriptions in >= 34% of top assignments. Furthermore, biochemical literature supports >80% of automated inferences for previously unannotated proteins. These results emphasize the benefit of incorporating functional linkages in protein annotation.

Molecular phylogeny and biogeography of langurs and leaf monkeys of South Asia (Primates: Colobinae)

The two recently proposed taxonomies of the langurs and leaf monkeys (Subfamily Colobinae) provide different implications to our understanding of the evolution of Nilgiri and purple-faced langurs. Groves (2001) [Groves, C.P., 2001. Primate Taxonomy. Smithsonian Institute Press, Washington], placed Nilgiri and purple-faced langurs in the genus Trachypithecus, thereby suggesting disjunct distribution of the genus Trachypithecus. [Brandon-Jones, D., Eudey, A.A., Geissmann, T., Groves, C.P., Melnick, D.J., Morales, J.C., Shekelle, M., Stewart, C.-B., 2003. Asian primate classification. Int. J. Primatol. 25, 97–162] placed these langurs in the genus Semnopithecus, which suggests convergence of morphological characters in Nilgiri and purple-faced langurs with Trachypithecus. To test these scenarios, we sequenced and analyzed the mitochondrial cytochrome b gene and two nuclear DNA-encoded genes, lysozyme and protamine P1, from a variety of colobine species. All three markers support the clustering of Nilgiri and purple-faced langurs with Hanuman langur (Semnopithecus), while leaf monkeys of Southeast Asian (Trachypithecus) form a distinct clade. The phylogenetic position of capped and golden leaf monkeys is still unresolved. It is likely that this species group might have evolved due to past hybridization between Semnopithecus and Trachypithecus clades.

Melting of an Anchored Bilayer: Molecular Dynamics Simulations of the Structural Transition in (CnH2n+1NH3)(2)PbI4 (n=12, 14, 16, 18)

The thermally driven Structural phase transition in the organic-inorganic hybrid perovskite (CnH2n+1NH3)(2)PbI4 has been investigated using molecular dynamics (MD) simulations. This system consists of positively charged alkyl-amine chains anchored to a rigid negatively charged PbI4 sheet with the chains organized as bilayers with a herringbone arrangement. Atomistic simulations were performed using ail isothermal-isobaric ensemble over a wide temperature range from 65 to 665 K for different alkyl chain lengths, n = 12, 14, 16, and 18. The simulations are able to reproduce the essential Features of the experimental observations of this system, including the existence of a transition, the linear variation of the transition temperature with alkyl chain length, and the expansion of the bilayer thickness at the transition. By use of the distance fluctuation Criteria, it is Shown that the transition is associated With a Melting of the alkyl chains of the anchored bilayer. Ail analysis of the conformation of the alkyl chains shows increased disorder in the form of gauche defects above due melting transition. Simulations also show that the melting transition is characterized by the complete disappearance of all-trans alkyl chains in the anchored bilayer, in agreement with experimental observations. A conformationally disordered chain has a larger effective cross-sectional area, and above due transition a uniformly tilted arrangement of the anchored chains call no longer be Sustained. At the melt the angular distribution of the orientation of the chains are 110 longer uniform; the chains are splayed allowing for increased space for individual chains of the anchored bilayer. This is reflected in a sharp rise in the ratio of the mean head-to-head to tail-to-tail distance of the chains of the bilayer at the transition resulting in in expansion of the bilayer thickness. The present MD simulations provide a simple explanation as to how changes in conformation of individual alkyl-chains gives rise to the observed increase in the interlayer lattice spacing of (CnH2n+1NH3)(2)PbI4 at the melting transition.

Structures and molecular-dynamics studies of three active-site mutants of bovine pancreatic phospholipase A(2)

Phospholipase A(2) hydrolyzes phospholipids at the sn-2 position to cleave the fatty-acid ester bond of L-glycerophospholipids. The catalytic dyad (Asp99 and His48) along with a nucleophilic water molecule is responsible for enzyme hydrolysis. Furthermore, the residue Asp49 in the calcium-binding loop is essential for controlling the binding of the calcium ion and the catalytic action of phospholipase A2. To elucidate the structural role of His48 and Asp49, the crystal structures of three active-site single mutants H48N, D49N and D49K have been determined at 1.9 angstrom resolution. Although the catalytically important calcium ion is present in the H48N mutant, the crystal structure shows that proton transfer is not possible from the catalytic water to the mutated residue. In the case of the Asp49 mutants, no calcium ion was found in the active site. However, the tertiary structures of the three active-site mutants are similar to that of the trigonal recombinant enzyme. Molecular-dynamics simulation studies provide a good explanation for the crystallographic results.

Two elementary proofs of the Wigner theorem on symmetry in quantum mechanics

In quantum theory, symmetry has to be defined necessarily in terms of the family of unit rays, the state space. The theorem of Wigner asserts that a symmetry so defined at the level of rays can always be lifted into a linear unitary or an antilinear antiunitary operator acting on the underlying Hilbert space. We present two proofs of this theorem which are both elementary and economical. Central to our proofs is the recognition that a given Wigner symmetry can, by post-multiplication by a unitary symmetry, be taken into either the identity or complex conjugation. Our analysis often focuses on the behaviour of certain two-dimensional subspaces of the Hilbert space under the action of a given Wigner symmetry, but the relevance of this behaviour to the larger picture of the whole Hilbert space is made transparent at every stage.

The effect of halogen atom on the molecular quadratic nonlinearity of half sandwich complexes in the presence of an acceptor

Half sandwich complexes of the type [CpM(CO)(n)X] {X=Cl, Br, I; If, M=Fe, Ru; n=2 and if M=Mo; n=3} and [CpNiPPh3X] {X=Cl, Br, I} have been synthesized and their second order molecular nonlinearity (beta) measured at 1064 nm in CHCl3 by the hyper-Rayleigh scattering technique. Iron complexes consistently display larger beta values than ruthenium complexes while nickel complexes have marginally larger beta values than iron complexes. In the presence of an acceptor ligand such as CO or PPh3, the role of the halogen atom is that of a pi donor. The better overlap of Cl orbitals with Fe and Ni metal centres make Cl a better pi donor than Br or I in the respective complexes. Consequently, M-pi interaction is stronger in Fe/Ni-Cl complexes. The value of beta decreases as one goes down the halogen group. For the complexes of 4d metal ions where the metal-ligand distance is larger, the influence of pi orbital overlap appears to be less important, resulting in moderate changes in beta as a function of halogen substitution. (C) 2006 Elsevier B.V. All rights reserved.

Three-dimensional structure of heat shock protein 90 from Plasmodium falciparum: molecular modelling approach to rational drug design against malaria

We have recently implicated heat shock protein 90 from Plasmodium falciparum (PfHsp90) as a potential drug target against malaria. Using inhibitors specific to the nucleotide binding domain of Hsp90, we have shown potent growth inhibitory effects on development of malarial parasite in human erythrocytes. To gain better understanding of the vital role played by PfHsp90 in parasite growth, we have modeled its three dimensional structure using recently described full length structure of yeast Hsp90. Sequence similarity found between PfHsp90 and yeast Hsp90 allowed us to model the core structure with high confidence. The superimposition of the predicted structure with that of the template yeast Hsp90 structure reveals an RMSD of 3.31 angstrom. The N-terminal and middle domains showed the least RMSD (1.76 angstrom) while the more divergent C-terminus showed a greater RMSD (2.84 angstrom) with respect to the template. The structure shows overall conservation of domains involved in nucleotide binding, ATPase activity, co-chaperone binding as well as inter-subunit interactions. Important co-chaperones known to modulate Hsp90 function in other eukaryotes are conserved in malarial parasite as well. An acidic stretch of amino acids found in the linker region, which is uniquely extended in PfHsp90 could not be modeled in this structure suggesting a flexible conformation. Our results provide a basis to compare the overall structure and functional pathways dependent on PfHsp90 in malarial parasite. Further analysis of differences found between human and parasite Hsp90 may make it possible to design inhibitors targeted specifically against malaria.

Effect of substitution on molecular conformation and packing features in a series of aryl substituted ethyl-6-methyl-4-phenyl-2-thioxo-1,2,3,4-tetrahydropyrimidine-5-carboxylates

The supramolecular structures of eight aryl protected ethyl-6-methyl-4-phenyl-2-thioxo-1,2,3,4 tetrahydropyrimidine-5-carboxyl ates were analyzed in order to understand the effect of variations in functional groups on molecular geometry, conformation and packing of molecules in the crystalline lattice. It is observed that the existence of a short intra-molecular C-H center dot center dot center dot pi interaction between the aromatic hydrogen of the aryl ring with the isolated double bond of the six-membered tetrahydropyrimidine ring is a key feature which imparts additional stability to the molecular conformation in the solid state. The compounds pack via the cooperative involvement of both N-H center dot center dot center dot S=C and N-H center dot center dot center dot O=C intermolecular dimers forming a sheet like structure. In addition, weak C-H center dot center dot center dot O and C-H center dot center dot center dot pi intermolecular interactions provide additional stability to the crystal packing.

Diffusion of Pure CH4 and Its Binary Mixture with CO2 in Faujasite NaY: A Combination of Neutron Scattering Experiments and Molecular Dynamics Simulations

The self-diffusion properties of pure CH4 and its binary mixture with CO2 within MY zeolite have been investigated by combining an experimental quasi-elastic neutron scattering (QENS) technique and classical Molecular dynamics simulations. The QENS measurements carried out at 200 K led to an unexpected self-diffusivity profile for Pure CH4 with the presence of a maximum for a loading of 32 CH4/unit cell, which was never observed before for the diffusion of apolar species in azeolite system With large windows. Molecular dynamics simulations were performed using two distinct microscopic models for representing the CH4/NaY interactions. Depending on the model, we are able to fairly reproduce either the magnitude or the profile of the self-diffusivity.Further analysis allowed LIS to provide some molecular insight into the diffusion mechanism in play. The QENS measurements report only a slight decrease of the self-diffusivity of CH4 in the presence of CO2 when the CO2 loading increases. Molecular dynamics simulations successfully capture this experimental trend and suggest a plausible microscopic diffusion mechanism in the case of this binary mixture.

Understanding Membranes through the Molecular Design of Lipids

Lipids are amphiphilic molecules that are composed of hydrophilic and hydrophobic regions. A typical membranous aggregate (vesicles, water-filled lipid nanospheres) is formed upon the self-organization of lipids in water from a diverse collection of amphiphiles producing a dynamic supramolecular structure that shows phase behavior and ordering as required for specific biological functions. The determination of various physical properties of lipid aggregates is the key to determining structure-function relationships. Over the years, we have designed and synthesized a wide variety of lipid molecular systems for the investigation of their membrane-forming properties and have used them for purposes such as gene delivery and enzyme activation. In this feature article, we focus on our work on various types of lipids including ion-paired amphiphiles, cholesterol-based lipids, aromatic lipids, macrocyclic lipids containing disulfide tethers; cationic dimeric lipids, and so forth. The emphasis is oil experimental design and bottom-line conclusions.

Purification, characterization and molecular cloning of a monocot mannose-binding lectin from Remusatia vivipara with nematicidal activity

A mannose-binding lectin (RVL) was purified from the tubers of Remusatia vivipara, a monocot plant by single-step affinity chromatography on asialofetuin-Sepharose 4B. RVL agglutinated only rabbit erythrocytes and was inhibited by mucin, asialomucin, asialofetuin and thyroglobulin. Lectin activity was stable up to 80A degrees C and under wide range of pH (2.0-9.3). SDS-PAGE and gel filtration results showed the lectin is a homotetramer of Mr 49.5 kDa, but MALDI analysis showed two distinct peaks corresponding to subunit mass of 12 kDa and 12.7 kDa. Also the N-terminal sequencing gave two different sequences indicating presence of two polypeptide chains. Cloning of RVL gene indicated posttranslational cleavage of RVL precursor into two mature polypeptides of 116 and 117 amino-acid residues. Dynamic light scattering (DLS) and gel filtration studies together confirmed the homogeneity of the purified lectin and supported RVL as a dimer with Mr 49.5 kDa derived from single polypeptide precursor of 233 amino acids. Purified RVL exerts potent nematicidal activity on Meloidogyne incognita, a root knot nematode. Fluorescent confocal microscopic studies demonstrated the binding of RVL to specific regions of the alimentary-tract and exhibited a potent toxic effect on M. incognita. RVL-mucin complex failed to interact with the gut confirming the receptor mediated lectin interaction. Very high mortality (88%) rate was observed at lectin concentration as low as 30 A mu g/ml, suggesting its potential application in the development of nematode resistant transgenic-crops.

Ion channel formation by zervamicin-IIB. A molecular modelling study.

Zervamicin-IIB (Zrv-IIB) is a 16 residue peptaibol which forms voltage-activated, multiple conductance level channels in planar lipid bilayers. A molecular model of Zrv-IIB channels is presented. The structure of monomeric Zrv-IIB is based upon the crystal structure of Zervamicin-Leu. The helical backbone is kinked by a hydroxyproline residue at position 10. Zrv-IIB channels are modelled as helix bundles of from 4 to 8 parallel helices surrounding a central pore. The monomers are packed with their C-terminal helical segments in close contact, and the bundles are stabilized by hydrogen bonds between glutamine 11 and hydroxyproline 10 of adjacent helices. Interaction energy profiles for movement of three different probes species (K+, Cl- and water) through the central pore are analyzed. The conformations of: (a) the sidechain of glutamine 3; (b) the hydroxyl group of hydroxyproline 10; and (c) the C-terminal hydroxyl group are "optimized" in order to maximize favourable interactions between the channel and the probes, resulting in favourable interaction energy profiles for all three. This suggests that conformational flexibility of polar sidechains enables the channel lining to mimic an aqueous environment.

Charge density based classification of intermolecular interactions in molecular crystals

Evaluation of intermolecular interactions in terms of both experimental and theoretical charge density analyses has produced a unified picture with which to classify strong and weak hydrogen bonds, along with van der Waals interactions, into three regions.

Molecular markers for discriminating Streptococcus pyogenes and S.dysgalactiae subspecies equisimilis

Given the increasing aetiological importance of Streptococcus dysgalactiae subspecies equisimilis in diseases which are primarily attributed to S. pyogenes, molecular markers are essential to distinguish these species and delineate their epidemiology more precisely. Many clinical microbiology laboratories rely on agglutination reactivity and biochemical tests to distinguish them. These methods have limitations which are particularly exacerbated when isolates with mixed properties are encountered. In order to provide additional distinguishing parameters that could be used to unequivocally discriminate these two common pathogens, we assess here three molecular targets: the speB gene, intergenic region upstream of the scpG gene (IRSG) and virPCR. Of these, the former two respectively gave positive and negative results for S. pyogenes, and negative and positive results for S. dysgalactiae subsp. equisimilis. Thus,a concerted use of these nucleic acid-based methods is particularly helpful in epidemiological surveillance to accurately assess the relative contribution of these species to streptococcal infections and diseases.