RESUMEN
Machine learning plays an important and growing role in molecular simulation. The newest version of the OpenMM molecular dynamics toolkit introduces new features to support the use of machine learning potentials. Arbitrary PyTorch models can be added to a simulation and used to compute forces and energy. A higher-level interface allows users to easily model their molecules of interest with general purpose, pretrained potential functions. A collection of optimized CUDA kernels and custom PyTorch operations greatly improves the speed of simulations. We demonstrate these features in simulations of cyclin-dependent kinase 8 (CDK8) and the green fluorescent protein chromophore in water. Taken together, these features make it practical to use machine learning to improve the accuracy of simulations with only a modest increase in cost.
Asunto(s)
Simulación de Dinámica Molecular , Agua , Aprendizaje AutomáticoRESUMEN
Machine learning plays an important and growing role in molecular simulation. The newest version of the OpenMM molecular dynamics toolkit introduces new features to support the use of machine learning potentials. Arbitrary PyTorch models can be added to a simulation and used to compute forces and energy. A higher-level interface allows users to easily model their molecules of interest with general purpose, pretrained potential functions. A collection of optimized CUDA kernels and custom PyTorch operations greatly improves the speed of simulations. We demonstrate these features on simulations of cyclin-dependent kinase 8 (CDK8) and the green fluorescent protein (GFP) chromophore in water. Taken together, these features make it practical to use machine learning to improve the accuracy of simulations at only a modest increase in cost.
RESUMEN
Liquid-liquid phase separation (LLPS) is an important mechanism that helps explain the membraneless compartmentalization of the nucleus. Because chromatin compaction and LLPS are collective phenomena, linking their modulation to the physicochemical features of nucleosomes is challenging. Here, we develop an advanced multiscale chromatin model-integrating atomistic representations, a chemically-specific coarse-grained model, and a minimal model-to resolve individual nucleosomes within sub-Mb chromatin domains and phase-separated systems. To overcome the difficulty of sampling chromatin at high resolution, we devise a transferable enhanced-sampling Debye-length replica-exchange molecular dynamics approach. We find that nucleosome thermal fluctuations become significant at physiological salt concentrations and destabilize the 30-nm fiber. Our simulations show that nucleosome breathing favors stochastic folding of chromatin and promotes LLPS by simultaneously boosting the transient nature and heterogeneity of nucleosome-nucleosome contacts, and the effective nucleosome valency. Our work puts forward the intrinsic plasticity of nucleosomes as a key element in the liquid-like behavior of nucleosomes within chromatin, and the regulation of chromatin LLPS.
Asunto(s)
Cromatina/metabolismo , ADN/metabolismo , Histonas/metabolismo , Conformación de Ácido Nucleico , Nucleosomas/metabolismo , Algoritmos , Cromatina/química , Cromatina/genética , Simulación por Computador , ADN/química , ADN/genética , Histonas/química , Modelos Genéticos , Simulación de Dinámica Molecular , Nucleosomas/química , Nucleosomas/genéticaRESUMEN
Neutrophils release their intracellular content, DNA included, into the bloodstream to form neutrophil extracellular traps (NETs) that confine and kill circulating pathogens. The mechanosensitive adhesive blood protein, von Willebrand Factor (vWF), interacts with the extracellular DNA of NETs to potentially immobilize them during inflammatory and coagulatory conditions. Here, we elucidate the previously unknown molecular mechanism governing the DNA-vWF interaction by integrating atomistic, coarse-grained, and Brownian dynamics simulations, with thermophoresis, gel electrophoresis, fluorescence correlation spectroscopy (FCS), and microfluidic experiments. We demonstrate that, independently of its nucleotide sequence, double-stranded DNA binds to a specific helix of the vWF A1 domain, via three arginines. This interaction is attenuated by increasing the ionic strength. Our FCS and microfluidic measurements also highlight the key role shear-stress has in enabling this interaction. Our simulations attribute the previously-observed platelet-recruitment reduction and heparin-size modulation, upon establishment of DNA-vWF interactions, to indirect steric hindrance and partial overlap of the binding sites, respectively. Overall, we suggest electrostatics-guiding DNA to a specific protein binding site-as the main driving force defining DNA-vWF recognition. The molecular picture of a key shear-mediated DNA-protein interaction is provided here and it constitutes the basis for understanding NETs-mediated immune and hemostatic responses.
Asunto(s)
ADN/química , Simulación del Acoplamiento Molecular , Factor de von Willebrand/química , Sitios de Unión , ADN/metabolismo , Humanos , Simulación de Dinámica Molecular , Concentración Osmolar , Unión Proteica , Electricidad Estática , Factor de von Willebrand/metabolismoRESUMEN
Protein flexibility and disorder is emerging as a crucial modulator of chromatin structure. Histone tail disorder enables transient binding of different molecules to the nucleosomes, thereby promoting heterogeneous and dynamic internucleosome interactions and making possible recruitment of a wide-range of regulatory and remodeling proteins. On the basis of extensive multiscale modeling we reveal the importance of linker histone H1 protein disorder for chromatin hierarchical looping. Our multiscale approach bridges microsecond-long bias-exchange metadynamics molecular dynamics simulations of atomistic 211-bp nucleosomes with coarse-grained Monte Carlo simulations of 100-nucleosome systems. We show that the long C-terminal domain (CTD) of H1-a ubiquitous nucleosome-binding protein-remains disordered when bound to the nucleosome. Notably, such CTD disorder leads to an asymmetric and dynamical nucleosome conformation that promotes chromatin structural flexibility and establishes long-range hierarchical loops. Furthermore, the degree of condensation and flexibility of H1 can be fine-tuned, explaining chromosomal differences of interphase versus metaphase states that correspond to partial and hyperphosphorylated H1, respectively. This important role of H1 protein disorder in large-scale chromatin organization has a wide range of biological implications.