Variational Methods for Biomolecular Modeling

Structure, function and dynamics of many biomolecular systems can be characterized by the energetic variational principle and the corresponding systems of partial differential equations (PDEs). This principle allows us to focus on the identification of essential energetic components, the optimal parametrization of energies, and the efficient computational implementation of energy variation or minimization. Given the fact that complex biomolecular systems are structurally non-uniform and their interactions occur through contact interfaces, their free energies are associated with various interfaces as well, such as solute-solvent interface, molecular binding interface, lipid domain interface, and membrane surfaces. This fact motivates the inclusion of interface geometry, particular its curvatures, to the parametrization of free energies. Applications of such interface geometry based energetic variational principles are illustrated through three concrete topics: the multiscale modeling of biomolecular electrostatics and solvation that includes the curvature energy of the molecular surface, the formation of microdomains on lipid membrane due to the geometric and molecular mechanics at the lipid interface, and the mean curvature driven protein localization on membrane surfaces. By further implicitly representing the interface using a phase field function over the entire domain, one can simulate the dynamics of the interface and the corresponding energy variation by evolving the phase field function, achieving significant reduction of the number of degrees of freedom and computational complexity. Strategies for improving the efficiency of computational implementations and for extending applications to coarse-graining or multiscale molecular simulations are outlined.Comment: 36 page

Reproductive protein evolution in two cryptic species of marine chordate

<p>Abstract</p> <p>Background</p> <p>Reproductive character displacement (RCD) is a common and taxonomically widespread pattern. In marine broadcast spawning organisms, behavioral and mechanical isolation are absent and prezygotic barriers between species often operate only during the fertilization process. Such barriers are usually a consequence of differences in the way in which sperm and egg proteins interact, so RCD can be manifest as faster evolution of these proteins between species in sympatry than allopatry. Rapid evolution of these proteins often appears to be a consequence of positive (directional) selection. Here, we identify a set of candidate gamete recognition proteins (GRPs) in the ascidian <it>Ciona intestinalis </it>and showed that these GRPs evolve more rapidly than control proteins (those not involved in gamete recognition). Choosing a subset of these gamete recognition proteins that show evidence of positive selection (CIPRO37.40.1, CIPRO60.5.1, CIPRO100.7.1), we then directly test the RCD hypothesis by comparing divergence (omega) and polymorphism (McDonald-Kreitman, Tajima's D, Fu and Li's D and F, Fay and Wu's H) statistics in sympatric and allopatric populations of two distinct forms of <it>C. intestinalis </it>(Types A and B) between which there are strong post-zygotic barriers.</p> <p>Results</p> <p>Candidate gamete recognition proteins from two lineages of <it>C. intestinalis </it>(Type A and B) are evolving more rapidly than control proteins, consistent with patterns seen in insects and mammals. However, ω (d_N/d_S) is not significantly different between the sympatric and allopatric populations, and none of the polymorphism statistics show significant differences between sympatric and allopatric populations.</p> <p>Conclusions</p> <p>Enhanced prezygotic isolation in sympatry has become a well-known feature of gamete recognition proteins in marine broadcast spawners. But in most cases the evolutionary process or processes responsible for this pattern have not been identified. Although gamete recognition proteins in <it>C. intestinalis </it>do appear to evolve more rapidly, on average, than proteins with other functions, rates of evolution are not different in allopatric and sympatric populations of the two reproductively isolated forms. That sympatry is probably human-mediated, and therefore recent, may explain the absence of RCD.</p

Measurement of the Bottom-Strange Meson Mixing Phase in the Full CDF Data Set

We report a measurement of the bottom-strange meson mixing phase \beta_s using the time evolution of B0_s -> J/\psi (->\mu+\mu-) \phi (-> K+ K-) decays in which the quark-flavor content of the bottom-strange meson is identified at production. This measurement uses the full data set of proton-antiproton collisions at sqrt(s)= 1.96 TeV collected by the Collider Detector experiment at the Fermilab Tevatron, corresponding to 9.6 fb-1 of integrated luminosity. We report confidence regions in the two-dimensional space of \beta_s and the B0_s decay-width difference \Delta\Gamma_s, and measure \beta_s in [-\pi/2, -1.51] U [-0.06, 0.30] U [1.26, \pi/2] at the 68% confidence level, in agreement with the standard model expectation. Assuming the standard model value of \beta_s, we also determine \Delta\Gamma_s = 0.068 +- 0.026 (stat) +- 0.009 (syst) ps-1 and the mean B0_s lifetime, \tau_s = 1.528 +- 0.019 (stat) +- 0.009 (syst) ps, which are consistent and competitive with determinations by other experiments.Comment: 8 pages, 2 figures, Phys. Rev. Lett 109, 171802 (2012

Hybrids of the bHLH and bZIP Protein Motifs Display Different DNA-Binding Activities In Vivo vs. In Vitro

Minimalist hybrids comprising the DNA-binding domain of bHLH/PAS (basic-helix-loop-helix/Per-Arnt-Sim) protein Arnt fused to the leucine zipper (LZ) dimerization domain from bZIP (basic region-leucine zipper) protein C/EBP were designed to bind the E-box DNA site, CACGTG, targeted by bHLHZ (basic-helix-loop-helix-zipper) proteins Myc and Max, as well as the Arnt homodimer. The bHLHZ-like structure of ArntbHLH-C/EBP comprises the Arnt bHLH domain fused to the C/EBP LZ: i.e. swap of the 330 aa PAS domain for the 29 aa LZ. In the yeast one-hybrid assay (Y1H), transcriptional activation from the E-box was strong by ArntbHLH-C/EBP, and undetectable for the truncated ArntbHLH (PAS removed), as detected via readout from the HIS3 and lacZ reporters. In contrast, fluorescence anisotropy titrations showed affinities for the E-box with ArntbHLH-C/EBP and ArntbHLH comparable to other transcription factors (Kd 148.9 nM and 40.2 nM, respectively), but only under select conditions that maintained folded protein. Although in vivo yeast results and in vitro spectroscopic studies for ArntbHLH-C/EBP targeting the E-box correlate well, the same does not hold for ArntbHLH. As circular dichroism confirms that ArntbHLH-C/EBP is a much more strongly α-helical structure than ArntbHLH, we conclude that the nonfunctional ArntbHLH in the Y1H must be due to misfolding, leading to the false negative that this protein is incapable of targeting the E-box. Many experiments, including protein design and selections from large libraries, depend on protein domains remaining well-behaved in the nonnative experimental environment, especially small motifs like the bHLH (60–70 aa). Interestingly, a short helical LZ can serve as a folding- and/or solubility-enhancing tag, an important device given the focus of current research on exploration of vast networks of biomolecular interactions

Contribution of discourse and morphosyntax skills to reading comprehension in Chinese dyslexic and typically developing children

This study aimed at identifying important skills for reading comprehension in Chinese dyslexic children and their typically developing counterparts matched on age (CA controls) or reading level (RL controls). The children were assessed on Chinese reading comprehension, cognitive, and reading-related skills. Results showed that the dyslexic children performed significantly less well than the CA controls but similarly to RL controls in most measures. Results of multiple regression analyses showed that word-level reading-related skills like oral vocabulary and word semantics were found to be strong predictors of reading comprehension among typically developing junior graders and dyslexic readers of senior grades, whereas morphosyntax, a text-level skill, was most predictive for typically developing senior graders. It was concluded that discourse and morphosyntax skills are particularly important for reading comprehension in the non-inflectional and topic-prominent Chinese system

Effects of methylphenidate on attention in Wistar rats treated with the neurotoxin N-(2-chloroethyl)-N-ethyl-2-bromobenzylamine (DSP4)

The aim of this study was to assess the effects of the neurotoxin N-(2-chloroethyl)-N-ethyl-2-bromobenzylamine (DSP4) on attention in rats as measured using the 5-choice-serial-reaction-time task (5CSRTT) and to investigate whether methylphenidate has effects on DSP4-treated rats. Methylphenidate is a noradrenaline and dopamine reuptake inhibitor and commonly used in the pharmacological treatment of individuals with attention deficit/hyperactivity disorder (ADHD). Wistar rats were trained in the 5CSRTT and treated with one of three doses of DSP4 or saline. Following the DSP4 treatment rats were injected with three doses of methylphenidate or saline and again tested in the 5CSRTT. The treatment with DSP4 caused a significant decline of performance in the number of correct responses and a decrease in response accuracy. A reduction in activity could also be observed. Whether or not the cognitive impairments are due to attention deficits or changes in explorative behaviour or activity remains to be investigated. The treatment with methylphenidate had no beneficial effect on the rats’ performance regardless of the DSP4 treatment. In the group without DSP4 treatment, methylphenidate led to a reduction in response accuracy and bidirectional effects in regard to parameters related to attention. These findings support the role of noradrenaline in modulating attention and call for further investigations concerning the effects of methylphenidate on attentional processes in rats

Molecular evolution of a gene cluster of serine proteases expressed in the Anopheles gambiae female reproductive tract

<p>Abstract</p> <p>Background</p> <p>Genes involved in post-mating processes of multiple mating organisms are known to evolve rapidly due to coevolution driven by sexual conflict among male-female interacting proteins. In the malaria mosquito <it>Anopheles gambiae </it>- a monandrous species in which sexual conflict is expected to be absent or minimal - recent data strongly suggest that proteolytic enzymes specifically expressed in the female lower reproductive tissues are involved in the processing of male products transferred to females during mating. In order to better understand the role of selective forces underlying the evolution of proteins involved in post-mating responses, we analysed a cluster of genes encoding for three serine proteases that are down-regulated after mating, two of which specifically expressed in the atrium and one in the spermatheca of <it>A. gambiae </it>females.</p> <p>Results</p> <p>The analysis of polymorphisms and divergence of these female-expressed proteases in closely related species of the <it>A. gambiae </it>complex revealed a high level of replacement polymorphisms consistent with relaxed evolutionary constraints of duplicated genes, allowing to rapidly fix novel replacements to perform new or more specific functions. Adaptive evolution was detected in several codons of the 3 genes and hints of episodic selection were also found. In addition, the structural modelling of these proteases highlighted some important differences in their substrate specificity, and provided evidence that a number of sites evolving under selective pressures lie relatively close to the catalytic triad and/or on the edge of the specificity pocket, known to be involved in substrate recognition or binding. The observed patterns suggest that these proteases may interact with factors transferred by males during mating (e.g. substrates, inhibitors or pathogens) and that they may have differently evolved in independent <it>A. gambiae </it>lineages.</p> <p>Conclusions</p> <p>Our results - also examined in light of constraints in the application of selection-inference methods to the closely related species of the <it>A. gambiae </it>complex - reveal an unexpectedly intricate evolutionary scenario. Further experimental analyses are needed to investigate the biological functions of these genes in order to better interpret their molecular evolution and to assess whether they represent possible targets for limiting the fertility of <it>Anopheles </it>mosquitoes in malaria vector control strategies.</p

Identification of Spt5 Target Genes in Zebrafish Development Reveals Its Dual Activity In Vivo

Spt5 is a conserved essential protein that represses or stimulates transcription elongation in vitro. Immunolocalization studies on Drosophila polytene chromosomes suggest that Spt5 is associated with many loci throughout the genome. However, little is known about the prevalence and identity of Spt5 target genes in vivo during development. Here, we identify direct target genes of Spt5 using fogsk8 zebrafish mutant, which disrupts the foggy/spt5 gene. We identified that fogsk8 and their wildtype siblings differentially express less than 5% of genes examined. These genes participate in diverse biological processes from stress response to cell fate specification. Up-regulated genes exhibit shorter overall gene length compared to all genes examined. Through chromatin immunoprecipitation in zebrafish embryos, we identified a subset of developmentally critical genes that are bound by both Spt5 and RNA polymerase II. The protein occupancy patterns on these genes are characteristic of both repressive and stimulatory elongation regulation. Together our findings establish Spt5 as a dual regulator of transcription elongation in vivo and identify a small but diverse set of target genes critically dependent on Spt5 during development

Explicitly searching for useful inventions: dynamic relatedness and the costs of connecting versus synthesizing

Inventions combine technological features. When features are barely related, burdensomely broad knowledge is required to identify the situations that they share. When features are overly related, burdensomely broad knowledge is required to identify the situations that distinguish them. Thus, according to my first hypothesis, when features are moderately related, the costs of connecting and costs of synthesizing are cumulatively minimized, and the most useful inventions emerge. I also hypothesize that continued experimentation with a specific set of features is likely to lead to the discovery of decreasingly useful inventions; the earlier-identified connections reflect the more common consumer situations. Covering data from all industries, the empirical analysis provides broad support for the first hypothesis. Regressions to test the second hypothesis are inconclusive when examining industry types individually. Yet, this study represents an exploratory investigation, and future research should test refined hypotheses with more sophisticated data, such as that found in literature-based discovery research