Learning probabilistic models of hydrogen bond stability from molecular dynamics simulation trajectories

Hydrogen bonds (H-bonds) play a key role in both the formation and stabilization of protein structures. H-bonds involving atoms from residues that are close to each other in the main-chain sequence stabilize secondary structure elements. H-bonds between atoms from distant residues stabilize a protein’s tertiary structure. However, H-bonds greatly vary in stability. They form and break while a protein deforms. For instance, the transition of a protein from a nonfunctional to a functional state may require some H-bonds to break and others to form. The intrinsic strength of an individual H-bond has been studied from an energetic viewpoint, but energy alone may not be a very good predictor. Other local interactions may reinforce (or weaken) an H-bond. This paper describes inductive learning methods to train a protein-independent probabilistic model of H-bond stability from molecular dynamics (MD) simulation trajectories. The training data describes H-bond occurrences at successive times along these trajectories by the values of attributes called predictors. A trained model is constructed in the form of a regression tree in which each non-leaf node is a Boolean test (split) on a predictor. Each occurrence of an H-bond maps to a path in this tree from the root to a leaf node. Its predicted stability is associated with the leaf node. Experimental results demonstrate that such models can predict H-bond stability quite well. In particular, their performance is roughly 20 % better than that of models based on H-bond energy alone. In addition, they can accurately identify a large fraction of the least stable H-bonds in a give

Quantum interference of electromechanically stabilized emitters in nanophotonic devices

Photon-mediated coupling between distant matter qubits may enable secure communication over long distances, the implementation of distributed quantum computing schemes, and the exploration of new regimes of many-body quantum dynamics. Nanophotonic devices coupled to solid-state quantum emitters represent a promising approach towards realization of these goals, as they combine strong light-matter interaction and high photon collection efficiencies. However, the scalability of these approaches is limited by the frequency mismatch between solid-state emitters and the instability of their optical transitions. Here we present a nano-electromechanical platform for stabilization and tuning of optical transitions of silicon-vacancy (SiV) color centers in diamond nanophotonic devices by dynamically controlling their strain environments. This strain-based tuning scheme has sufficient range and bandwidth to alleviate the spectral mismatch between individual SiV centers. Using strain, we ensure overlap between color center optical transitions and observe an entangled superradiant state by measuring correlations of photons collected from the diamond waveguide. This platform for tuning spectrally stable color centers in nanophotonic waveguides and resonators constitutes an important step towards a scalable quantum network

Strain engineering of the silicon-vacancy center in diamond

We control the electronic structure of the silicon-vacancy (SiV) color-center in diamond by changing its static strain environment with a nano-electro-mechanical system. This allows deterministic and local tuning of SiV optical and spin transition frequencies over a wide range, an essential step towards multiqubit networks. In the process, we infer the strain Hamiltonian of the SiV revealing large strain susceptibilities of order 1 PHz/strain for the electronic orbital states. We identify regimes where the spin-orbit interaction results in a large strain susceptibility of order 100 THz/strain for spin transitions, and propose an experiment where the SiV spin is strongly coupled to a nanomechanical resonator

Several Distinct Polycomb Complexes Regulate and Co-Localize on the INK4a Tumor Suppressor Locus

Misexpression of Polycomb repressive complex 1 (PRC1) components in human cells profoundly influences the onset of cellular senescence by modulating transcription of the INK4a tumor suppressor gene. Using tandem affinity purification, we find that CBX7 and CBX8, two Polycomb (Pc) homologs that repress INK4a, both participate in PRC1-like complexes with at least two Posterior sex combs (Psc) proteins, MEL18 and BMI1. Each complex contains a single representative of the Pc and Psc families. In primary human fibroblasts, CBX7, CBX8, MEL18 and BMI1 are present at the INK4a locus and shRNA-mediated knockdown of any one of these components results in de-repression of INK4a and proliferative arrest. Sequential chromatin immunoprecipitation (ChIP) reveals that CBX7 and CBX8 bind simultaneously to the same region of chromatin and knockdown of one of the Pc or Psc proteins results in release of the other, suggesting that the binding of PRC1 complexes is interdependent. Our findings provide the first evidence that a single gene can be regulated by several distinct PRC1 complexes and raise important questions about their configuration and relative functions

RNAi-Mediated c-Rel Silencing Leads to Apoptosis of B Cell Tumor Cells and Suppresses Antigenic Immune Response In Vivo

c-Rel is a member of the Rel/NF-κB transcription factor family and is predominantly expressed in lymphoid and myeloid cells, playing a critical role in lymphocyte proliferation and survival. Persistent activation of the c-Rel signal transduction pathway is associated with allergies, inflammation, autoimmune diseases, and a variety of human malignancies. To explore the potential of targeting c-Rel as a therapeutic agent for these disorders, we designed a small interfering RNA (siRNA) to silence c-Rel expression in vitro and in vivo. C-Rel-siRNA expression via a retroviral vector in a B cell tumor cell line leads to growth arrest and apoptosis of the tumor cells. Silencing c-Rel in primary B cells in vitro compromises their proliferative and survival response to CD40 activation signals, similar to the impaired response of c-Rel knockout B cells. Most important, in vivo silencing of c-Rel results in significant impairment in T cell-mediated immune responses to antigenic stimulation. Our study thus validates the efficacy of c-Rel-siRNA, and suggests the development of siRNA-based therapy, as well as small molecular inhibitors for the treatment of B cell tumors as well as autoimmune diseases

Search for eccentric black hole coalescences during the third observing run of LIGO and Virgo

Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effects of eccentricity. Here, we present observational results for a waveform-independent search sensitive to eccentric black hole coalescences, covering the third observing run (O3) of the LIGO and Virgo detectors. We identified no new high-significance candidates beyond those that were already identified with searches focusing on quasi-circular binaries. We determine the sensitivity of our search to high-mass (total mass M>70 M⊙) binaries covering eccentricities up to 0.3 at 15 Hz orbital frequency, and use this to compare model predictions to search results. Assuming all detections are indeed quasi-circular, for our fiducial population model, we place an upper limit for the merger rate density of high-mass binaries with eccentricities 0<e≤0.3 at 0.33 Gpc−3 yr−1 at 90\% confidence level

Ultralight vector dark matter search using data from the KAGRA O3GK run

Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we present the result of a search for U(1)B−L gauge boson DM using the KAGRA data from auxiliary length channels during the first joint observation run together with GEO600. By applying our search pipeline, which takes into account the stochastic nature of ultralight DM, upper bounds on the coupling strength between the U(1)B−L gauge boson and ordinary matter are obtained for a range of DM masses. While our constraints are less stringent than those derived from previous experiments, this study demonstrates the applicability of our method to the lower-mass vector DM search, which is made difficult in this measurement by the short observation time compared to the auto-correlation time scale of DM