Learning generative models of molecular dynamics

We introduce three algorithms for learning generative models of molecular structures from molecular dynamics simulations. The first algorithm learns a Bayesian-optimal undirected probabilistic model over user-specified covariates (e.g., fluctuations, distances, angles, etc). L1 regularization is used to ensure sparse models and thus reduce the risk of over-fitting the data. The topology of the resulting model reveals important couplings between different parts of the protein, thus aiding in the analysis of molecular motions. The generative nature of the model makes it well-suited to making predictions about the global effects of local structural changes (e.g., the binding of an allosteric regulator). Additionally, the model can be used to sample new conformations. The second algorithm learns a time-varying graphical model where the topology and parameters change smoothly along the trajectory, revealing the conformational sub-states. The last algorithm learns a Markov Chain over undirected graphical models which can be used to study and simulate kinetics. We demonstrate our algorithms on multiple molecular dynamics trajectories.     

Rad50 Zinc Hook Is Important for the Mre11 Complex to Bind Chromosomal DNA Double-stranded Breaks and Initiate Various DNA Damage Responses

The Mre11-Rad50-Nbs1 (MRN) complex plays critical roles in checkpoint activation and double-stranded break (DSB) repair. The Rad50 zinc hook domain mediates zinc-dependent intercomplex associations of MRN, which is important for DNA tethering. Studies in yeast suggest that the Rad50 zinc hook domain is essential for MRN functions, but its role in mammalian cells is not clear. We demonstrated that the human Rad50 hook mutants are severely defective in various DNA damage responses including ATM (Ataxia telangiectasia mutated) activation, homologous recombination, sensitivity to IR, and activation of the ATR pathway. By using live cell imaging, we observed that the Rad50 hook mutants fail to be recruited to chromosomal DSBs, suggesting a novel mechanism underlying the severe defects observed for the Rad50 hook mutants. In vitro analysis showed that Zn(2+) promotes wild type but not the hook mutant of MR to bind double-stranded DNA. In vivo, the Rad50 hook mutants are defective in being recruited to chromosomal DSBs in both H2AX-proficient and -deficient cells, suggesting that the Rad50 hook mutants are impaired in direct binding to chromosomal DSB ends. We propose that the Rad50 zinc hook domain is important for the initial binding of MRN to DSBs, leading to ATM activation to phosphorylate H2AX, which recruits more MRN to the DSB-flanking chromosomal regions. Our studies reveal a critical role for the Rad50 zinc hook domain in establishing and maintaining MRN recruitment to chromosomal DSBs and suggest an important mechanism of how the Rad50 zinc hook domain contributes to DNA repair and checkpoint activation.     

Assessment of red blood cell aggregation with dextran by ultrasonic interferometry.

Aggregation of human red blood cells (RBCs) induced by dextrans of various molecular weight has been studied by using a new ultrasonic interferometry method. This method, based on A-mode echography, allowed for the measurement of the accumulation rate of particles on a solid plate which is related to their sedimentation rate (i.e., to their mean size). The initial aggregation process, the mean and the maximum sedimentation rate of aggregates and the packing of the sedimented RBCs have been investigated. Effects of hematocrit, molecular weight of dextrans and inhibition by dextran 40 on the RBC aggregation induced by dextran of higher molecular weight have been determined by analysing variations of the aggregate size. Results obtained confirm the aggregation effect of dextrans of molecular weights equal or higher than 70,000 dalton and disaggregation effect of dextran 40,000 dalton on aggregation by dextrans of higher molecular weight.