Ab Initio Quantum Mechanics Analysis of Imidazole C-H…O Water Hydrogen Bonding and a Molecular Mechanics Forcefield Correction

Abstract

While it is well established that classical hydrogen bonds play an important role in enzyme structure, function and dynamics, the role of weaker, but ‘activated’ C-H donor hydrogen bonds is poorly understood. The most important such case involves histidine which often plays a direct role in enzyme catalysis and possesses the most acidic C-H donor group of the standard amino acids. In the present study, we obtained optimized geometries and hydrogen bond interaction energies for C-H…O hydrogen bonded complexes between methane, ethylene, benzene, acetylene, and imidazole with water at the MP2-FC/6-31++G(2d,2p) and MP2-FC/aug-cc-pVDZ//MP2-FC/6-31++G(2d,2p) levels of theory. A strong linear relationship is obtained between the stability of the various hydrogen bonded complexes and both separation distances for H…0 and C—O. In general, these calculations indicate that C-H…0 interactions can be classified as hydrogen bonding interactions, albeit significantly weaker than the classical hydrogen bonds, but significantly stronger than just van der Waals interactions. For instance, while the electronic energy of stabilization at the MP2-FC/aug-cc-pVDZ//MP2-FC/6-31++G(2d,2p) level of theory of a water C-H…O water hydrogen bond is 4.36 kcal/mol more stable than the methane C-H…O water interaction, the water-water hydrogen bond is only 2.06 kcal/mol more stable than the imidazole C_e?H…O water hydrogen bond. Neglecting this latter hydrogen bonding interaction is obviously unacceptable. We next compare the potential energy surfaces for the imidazole C_e?H…O water and imidazole N_d?H…O hydrogen bonded complexes computed at the MP2/6-31++G(2d,2p) level of theory with the potential energy surface computed using the AMBER molecular mechanics program and forcefields. While the Weiner et al and Cornell et al AMBER forcefields reasonably account for the imidazole N-H…O water interaction, these forcefields do not adequately account for the imidazole C_e?H…O water hydrogen bond. A forcefield modification is offered that results in excellent agreement between the ab initio and molecular mechanics geometry and energy for this C-H…O hydrogen bonded complex.     

Parasites as prey in aquatic food webs: implications for predator infection and parasite transmission

While the recent inclusion of parasites into food‐web studies has highlighted the role of parasites as consumers, there is accumulating evidence that parasites can also serve as prey for predators. Here we investigated empirical patterns of predation on parasites and their relationships with parasite transmission in eight topological food webs representing marine and freshwater ecosystems. Within each food web, we examined links in the typical predator–prey sub web as well as the predator–parasite sub web, i.e. the quadrant of the food web indicating which predators eat parasites. Most predator– parasite links represented ‘concomitant predation’ (consumption and death of a parasite along with the prey/host; 58–72%), followed by ‘trophic transmission’ (predator feeds on infected prey and becomes infected; 8–32%) and predation on free‐living parasite life‐cycle stages (4–30%). Parasite life‐cycle stages had, on average, between 4.2 and 14.2 predators. Among the food webs, as predator richness increased, the number of links exploited by trophically transmitted parasites increased at about the same rate as did the number of links where these stages serve as prey. On the whole, our analyses suggest that predation on parasites has important consequences for both predators and parasites, and food web structure. Because our analysis is solely based on topological webs, determining the strength of these interactions is a promising avenue for future research.     

New microbial fossils from ∼1.3 billion‐year‐old rocks of Eastern California

New types of microbial fossils and new occurrences of fossils previously reported only from the Beck Spring Dolomite of the Pahrump Group are now recognized from each of the three formations of the Pahrump Group (Crystal Spring Formation, Beck Spring Dolomite, and Kingston Peak Formation) approximately 1.3 X 10° years old. Comprising perhaps eight or nine distinctive forms, these fossils are characteristically preserved as faint ghostlike structures whose low‐contrast outlines are clearly revealed only when illuminated by a xenon lamp and recorded on high‐contrast film. They represent a distinctive, previously overlooked or neglected type of preservation that has significantly extended the known distribution of microbial fossils in the Pahrump. They include the oldest occurrence known to us of filaments designatable as Girvanella and apparently the first from rocks of pre‐Phanerozoic age. Similar fossils were also found, using the same techniques, in the Chuar Group of the Grand Canyon and in the Uluntui Suite of middle Riphean age in eastern Siberia. Although time correlation of pre‐Phanerozoic rocks based on similar microbial assemblages would be premature, similarity between such assemblages in all formations of the Pahrump Group and with that of the Uluntui Suite is consistent with the inferred unity and middle Riphean age of the Pahrump Group. In addition to the Girvanella we find two smaller types of filaments, two kinds of simple spheroids, and three composite forms (two spheroids and one stalked cluster) that attain diameters up to 80 μm and are probably eucaryotic.     

Comparison of the Genome Sequences of “Candidatus Portiera aleyrodidarum” Primary Endosymbionts of the Whitefly Bemisia tabaci B and Q Biotypes

“Candidatus Portiera aleyrodidarum” is the primary endosymbiont of whiteflies. We report two complete genome sequences of this bacterium from the worldwide invasive B and Q biotypes of the whitefly Bemisia tabaci. Differences in the two genome sequences may add insights into the complex differences in the biology of both biotypes.     

Proton transport facilitating water-oxidation: the role of second sphere ligands surrounding the catalytic metal cluster

The ability of PSII to extract electrons from water, with molecular oxygen as a by-product, is a remarkable biochemical and evolutionary innovation. From an evolutionary perspective, the invention of PSII approximately 2.7 Ga led to the accelerated accumulation of biomass in the biosphere and the accumulation of oxygen in the atmosphere, a combination that allowed for the evolution of a much more complex and extensive biosphere than would otherwise have been possible. From the biochemical and enzymatic perspective, PSII is remarkable because of the thermodynamic and kinetic obstacles that needed to have been overcome to oxidize water as the ultimate photosynthetic electron donor. This article focuses on how proton release is an integral part of how these kinetic and thermodynamic obstacles have been overcome: the sequential removal of protons from the active site of H₂O-oxidation facilitates the multistep oxidation of the substrate water at the Mn₄CaO _x , the catalytic heart of the H₂O-oxidation reaction. As noted previously, the facilitated deprotonation of the Mn₄CaO _x cluster exerts a redox-leveling function preventing the accumulation of excess positive charge on the cluster, which might otherwise hinder the already energetically difficult oxidation of water. Using recent results, including the characteristics of site-directed mutants, the role of the second sphere of amino acid ligands and the associated network of water molecules surrounding the Mn₄CaO _x is discussed in relation to proton transport in other systems. In addition to the redox-leveling function, a trapping function is assigned to the proton release step occurring immediately prior to the dioxygen chemistry. This trapping appears to involve a yet-to-be clarified gating mechanism that facilitates to coordinated release of a proton from the neighborhood of the active site thereby insuring that the backward charge-recombination reaction does not out-compete the forward reaction of dioxygen chemistry during this final step of H₂O-oxidation.     

The antitumor activity of hydrophobin SC3, a fungal protein

The use of mushroom extracts has been common practice in traditional medicine for centuries, including the treatment of cancer. Proteins called hydrophobins are very abundant in mushrooms. Here, it was examined whether they have antitumor activity. Hydrophobin SC3 of Schizophyllum commune was injected daily intraperitoneally starting 1 day after tumor induction in two tumor mouse models (sarcoma and melanoma). SC3 reduced the size and weight of the melanoma significantly, but the sarcoma seemed not affected. However, microscopic analysis of the tumors 12 days after induction revealed a strong antitumor effect of SC3 on both tumors. The mitotic activity of the tumor decreased 1.6- (melanoma) to 2.3-fold (sarcoma), while the vital mass decreased 2.3- (melanoma) to 4.3-fold (sarcoma) compared to the control. Treatment did not cause any signs of toxicity. Behavior, animal growth, and weight of organs were similar to animals injected with vehicle, and no histological abnormalities were found in the organs. In vitro cell culture studies revealed no direct cytotoxic effect of SC3 towards sarcoma cells, while cytotoxic activity was observed towards melanoma cells at a high SC3 concentration. Daily treatment with SC3 did not result in detectable levels of anti-SC3 antibodies in the plasma. Instead, a cellular immune response was observed. Incubation of spleen cells with SC3 resulted in a 1.5- to 2.5-fold increase in interleukin-10 and TNF-α mRNA levels. In conclusion, the nontoxic fungal hydrophobin SC3 showed tumor-suppressive activity possibly via immunomodulation and may be of benefit as adjuvant in combination with chemotherapy and radiation.     

Molecular Dynamics Simulations of Homo-oligomeric Bundles Embedded Within a Lipid Bilayer

Using molecular dynamics simulations, we studied the structure, interhelix interactions, and dynamics of transmembrane proteins. Specifically, we investigated homooligomeric helical bundle systems consisting of synthetic α-helices with either the sequence Ac-(LSLLLSL)₃-NH₂ (LS2) or Ac-(LSSLLSL)₃-NH₂ (LS3). The LS2 and LS3 helical peptides are designed to have amphipathic characteristics that form ion channels in membrane. We simulated bundles containing one to six peptides that were embedded in palmitoyl-oleoyl-phosphatidylcholine (POPC) lipid bilayer and placed between two lamellae of water. We aim to provide a fundamental understanding of how amphipathic helical peptides interact with each other and their dynamical behaviors in different homooligomeric states. To understand structural properties, we examined the helix lengths, tilt angles of individual helices and the entire bundle, interhelix distances, interhelix cross-angles, helix hydrophobic-to-hydrophilic vector projections, and the average number of interhelix hydrophilic (serine–serine) contacts lining the pore of the transmembrane channel. To analyze dynamical properties, we calculated the rotational autocorrelation function of each helix and the cross-correlation of the rotational velocity between adjacent helices. The observed structural and dynamical characteristics show that higher order bundles containing four to six peptides are composed of multiple lower order bundles of one to three peptides. For example, the LS2 channel was found to be stable in a tetrameric bundle composed of a “dimer of dimers.” In addition, we observed that there is a minimum of two strong hydrophilic contacts between a pair of adjacent helices in the dimer to tetramer systems and only one strong hydrophilic interhelix contact in helix pairs of the pentamer and hexamer systems. We believe these results are general and can be applied to more complex ion channels, providing insight into ion channel stability and assembly.     

Detecting and Characterizing Genomic Signatures of Positive Selection in Global Populations

Natural selection is a significant force that shapes the architecture of the human genome and introduces diversity across global populations. The question of whether advantageous mutations have arisen in the human genome as a result of single or multiple mutation events remains unanswered except for the fact that there exist a handful of genes such as those that confer lactase persistence, affect skin pigmentation, or cause sickle cell anemia. We have developed a long-range-haplotype method for identifying genomic signatures of positive selection to complement existing methods, such as the integrated haplotype score (iHS) or cross-population extended haplotype homozygosity (XP-EHH), for locating signals across the entire allele frequency spectrum. Our method also locates the founder haplotypes that carry the advantageous variants and infers their corresponding population frequencies. This presents an opportunity to systematically interrogate the whole human genome whether a selection signal shared across different populations is the consequence of a single mutation process followed subsequently by gene flow between populations or of convergent evolution due to the occurrence of multiple independent mutation events either at the same variant or within the same gene. The application of our method to data from 14 populations across the world revealed that positive-selection events tend to cluster in populations of the same ancestry. Comparing the founder haplotypes for events that are present across different populations revealed that convergent evolution is a rare occurrence and that the majority of shared signals stem from the same evolutionary event.     

Deep Whole-Genome Sequencing of 100 Southeast Asian Malays

Whole-genome sequencing across multiple samples in a population provides an unprecedented opportunity for comprehensively characterizing the polymorphic variants in the population. Although the 1000 Genomes Project (1KGP) has offered brief insights into the value of population-level sequencing, the low coverage has compromised the ability to confidently detect rare and low-frequency variants. In addition, the composition of populations in the 1KGP is not complete, despite the fact that the study design has been extended to more than 2,500 samples from more than 20 population groups. The Malays are one of the Austronesian groups predominantly present in Southeast Asia and Oceania, and the Singapore Sequencing Malay Project (SSMP) aims to perform deep whole-genome sequencing of 100 healthy Malays. By sequencing at a minimum of 30× coverage, we have illustrated the higher sensitivity at detecting low-frequency and rare variants and the ability to investigate the presence of hotspots of functional mutations. Compared to the low-pass sequencing in the 1KGP, the deeper coverage allows more functional variants to be identified for each person. A comparison of the fidelity of genotype imputation of Malays indicated that a population-specific reference panel, such as the SSMP, outperforms a cosmopolitan panel with larger number of individuals for common SNPs. For lower-frequency (<5%) markers, a larger number of individuals might have to be whole-genome sequenced so that the accuracy currently afforded by the 1KGP can be achieved. The SSMP data are expected to be the benchmark for evaluating the value of deep population-level sequencing versus low-pass sequencing, especially in populations that are poorly represented in population-genetics studies.