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

Cytomolecular identification of individual wheat-wheat chromosome arm associations in wheat-rye hybrids

Chromosome pairing in the meiotic metaphase I of wheatrye hybrids has been characterized by sequential genomic and fluorescent in situ hybridization allowing not only the discrimination of wheat and rye chromosomes, but also the identification of the individual wheat and rye chromosome arms involved in the chromosome associations. The majority of associations (93.8%) were observed between the wheat chromosomes. The largest number of wheat-wheat chromosome associations (53%) was detected between the A and D genomes, while the frequency of B-D and A-B associations was significantly lower (32 and 8%, respectively). Among the A-D chromosome associations, pairing between the 3AL and 3DL arms was observed with the highest frequency, while the most frequent of all the chromosome associations (0.113/ cell) was found to be the 3DS-3BS. Differences in the pairing frequency of the individual chromosome arms of wheat-rye hybrids have been discussed in relation to the homoeologous relationships between the constituent genomes of hexaploid wheat

Developmental perspectives on interpersonal affective touch

In the last decade, philosophy, neuroscience and psychology alike have paid increasing attention to the study of interpersonal affective touch, which refers to the emotional and motivational facets of tactile sensation. Some aspects of affective touch have been linked to a neurophysiologically specialised system, namely the C tactile (CT) system. While the role of this sys-tem for affiliation, social bonding and communication of emotions have been widely investigated, only recently researchers have started to focus on the potential role of interpersonal affective touch in acquiring awareness of the body as our own, i.e. as belonging to our psychological ‘self’. We review and discuss recent developmental and adult findings, pointing to the central role of interpersonal affective touch in body awareness and social cognition in health and disorders. We propose that interpersonal affective touch, as an interoceptive modality invested of a social nature, can uniquely contribute to the ongoing debate in philosophy about the primacy of the relational nature of the minimal self

Multi-messenger observations of a binary neutron star merger

On 2017 August 17 a binary neutron star coalescence candidate (later designated GW170817) with merger time 12:41:04 UTC was observed through gravitational waves by the Advanced LIGO and Advanced Virgo detectors. The Fermi Gamma-ray Burst Monitor independently detected a gamma-ray burst (GRB 170817A) with a time delay of ~1.7 s with respect to the merger time. From the gravitational-wave signal, the source was initially localized to a sky region of 31 deg2 at a luminosity distance of 40+8-8 Mpc and with component masses consistent with neutron stars. The component masses were later measured to be in the range 0.86 to 2.26 Mo. An extensive observing campaign was launched across the electromagnetic spectrum leading to the discovery of a bright optical transient (SSS17a, now with the IAU identification of AT 2017gfo) in NGC 4993 (at ~40 Mpc) less than 11 hours after the merger by the One- Meter, Two Hemisphere (1M2H) team using the 1 m Swope Telescope. The optical transient was independently detected by multiple teams within an hour. Subsequent observations targeted the object and its environment. Early ultraviolet observations revealed a blue transient that faded within 48 hours. Optical and infrared observations showed a redward evolution over ~10 days. Following early non-detections, X-ray and radio emission were discovered at the transient’s position ~9 and ~16 days, respectively, after the merger. Both the X-ray and radio emission likely arise from a physical process that is distinct from the one that generates the UV/optical/near-infrared emission. No ultra-high-energy gamma-rays and no neutrino candidates consistent with the source were found in follow-up searches. These observations support the hypothesis that GW170817 was produced by the merger of two neutron stars in NGC4993 followed by a short gamma-ray burst (GRB 170817A) and a kilonova/macronova powered by the radioactive decay of r-process nuclei synthesized in the ejecta

Racial Disparity in Police Stop and Searches in England and Wales

Data published by the United Kingdom's Ministry for Justice clearly shows that, compared to persons who were White, members of racial minorities in England, particularly Blacks, were far more likely to be stopped and searched by the police. The question is whether such racial disparity in stops and searches could be justified by racial disparities in offending? Or whether the disparity in stop and searches exceeded the disparity in offending? This paper proposes a method for measuring the amount of excess in racial disparity in police stop and searches. Using the most recently published Ministry of Justice data (for 2007/08) for Police Areas in England and Wales it concludes that while in several Areas there was no excess to racial disparity in police stop and searches, there was, on the basis of the methodology proposed in the paper, evidence of such excess in some Police Areas of England and Wales

Tinkering Evolution of Post-Transcriptional RNA Regulons: Puf3p in Fungi as an Example

Genome-wide studies of post-transcriptional mRNA regulation in model organisms indicate a “post-transcriptional RNA regulon” model, in which a set of functionally related genes is regulated by mRNA–binding RNAs or proteins. One well-studied post-transcriptional regulon by Puf3p functions in mitochondrial biogenesis in budding yeast. The evolution of the Puf3p regulon remains unclear because previous studies have shown functional divergence of Puf3p regulon targets among yeast, fruit fly, and humans. By analyzing evolutionary patterns of Puf3p and its targeted genes in forty-two sequenced fungi, we demonstrated that, although the Puf3p regulon is conserved among all of the studied fungi, the dedicated regulation of mitochondrial biogenesis by Puf3p emerged only in the Saccharomycotina clade. Moreover, the evolution of the Puf3p regulon was coupled with evolution of codon usage bias in down-regulating expression of genes that function in mitochondria in yeast species after genome duplication. Our results provide a scenario for how evolution like a tinker exploits pre-existing materials of a conserved post-transcriptional regulon to regulate gene expression for novel functional roles

Double-stranded RNA-activated protein kinase PKR of fishes and amphibians: Varying the number of double-stranded RNA binding domains and lineage-specific duplications

BackgroundDouble-stranded (ds) RNA, generated during viral infection, binds and activates the mammalian anti-viral protein kinase PKR, which phosphorylates the translation initiation factor eIF2alpha leading to the general inhibition of protein synthesis. Although PKR-like activity has been described in fish cells, the responsible enzymes eluded molecular characterization until the recent discovery of goldfish and zebrafish PKZ, which contain Z-DNA-binding domains instead of dsRNA-binding domains (dsRBDs). Fish and amphibian PKR genes have not been described so far.ResultsHere we report the cloning and identification of 13 PKR genes from 8 teleost fish and amphibian species, including zebrafish, demonstrating the coexistence of PKR and PKZ in this latter species. Analyses of their genomic organization revealed up to three tandemly arrayed PKR genes, which are arranged in head-to-tail orientation. At least five duplications occurred independently in fish and amphibian lineages. Phylogenetic analyses reveal that the kinase domains of fish PKR genes are more closely related to those of fish PKZ than to the PKR kinase domains of other vertebrate species. The duplication leading to fish PKR and PKZ genes occurred early during teleost fish evolution after the divergence of the tetrapod lineage. While two dsRBDs are found in mammalian and amphibian PKR, one, two or three dsRBDs are present in fish PKR. In zebrafish, both PKR and PKZ were strongly upregulated after immunostimulation with some tissue-specific expression differences. Using genetic and biochemical assays we demonstrate that both zebrafish PKR and PKZ can phosphorylate eIF2alpha in yeast.ConclusionConsidering the important role for PKR in host defense against viruses, the independent duplication and fixation of PKR genes in different lineages probably provided selective advantages by leading to the recognition of an extended spectrum of viral nucleic acid structures, including both dsRNA and Z-DNA/RNA, and perhaps by altering sensitivity to viral PKR inhibitors. Further implications of our findings for the evolution of the PKR family and for studying PKR/PKZ interactions with viral gene products and their roles in viral infections are discussed

Pre-microRNA and Mature microRNA in Human Mitochondria

Chantier qualité GAInternational audienceBACKGROUND: Because of the central functions of the mitochondria in providing metabolic energy and initiating apoptosis on one hand and the role that microRNA (miRNA) play in gene expression, we hypothesized that some miRNA could be present in the mitochondria for post-transcriptomic regulation by RNA interference. We intend to identify miRNA localized in the mitochondria isolated from human skeletal primary muscular cells. METHODOLOGY/PRINCIPAL FINDINGS: To investigate the potential origin of mitochondrial miRNA, we in-silico searched for microRNA candidates in the mtDNA. Twenty five human pre-miRNA and 33 miRNA aligments (E-value35) for the smallest RNA input concentration and 204 miRNA for the maximum RNA input concentration. In silico analysis predicted 80 putative miRNA target sites in the mitochondrial genome (E-value<0.05). CONCLUSIONS/SIGNIFICANCE: The present study experimentally demonstrated for the first time the presence of pre-miRNA and miRNA in the human mitochondria isolated from skeletal muscular cells. A set of miRNA were significantly detected in mitochondria fraction. The origin of these pre-miRNA and miRNA should be further investigate to determine if they are imported from the cytosol and/or if they are partially processed in the mitochondria

OzDES Reverberation Mapping Programme: Mg II lags and R−L relation

The correlation between the broad line region radius and continuum luminosity (R-L relation) of active galactic nuclei (AGNs) is critical for single-epoch mass estimates of supermassive black holes (SMBHs). At z ∼ 1-2, where AGN activity peaks, the R-L relation is constrained by the reverberation mapping (RM) lags of the Mg II line. We present 25 Mg II lags from the Australian Dark Energy Survey RM project based on 6 yr of monitoring. We define quantitative criteria to select good lag measurements and verify their reliability with simulations based on both the damped random walk stochastic model and the rescaled, resampled versions of the observed light curves of local, well-measured AGN. Our sample significantly increases the number of Mg II lags and extends the R-L relation to higher redshifts and luminosities. The relative iron line strength RFe has little impact on the R-L relation. The best-fitting Mg II R-L relation has a slope α = 0.39 ± 0.08 with an intrinsic scatter σrl = 0.15+−000203. The slope is consistent with previous measurements and shallower than the H β R-L relation. The intrinsic scatter of the new R-L relation is substantially smaller than previous studies and comparable to the intrinsic scatter of the H β R-L relation. Our new R-L relation will enable more precise single-epoch mass estimates and SMBH demographic studies at cosmic noon