Artem Zhmurov

We introduce a Stochastic Reaction-Diffusion-Dynamics Model (SRDDM) for simulations of cellular mechanochemical processes with high spatial and temporal resolution. The SRDDM is mapped into the CellDynaMo package, which couples the spatially inhomogeneous reaction-diffusion master equation to account for biochemical reactions and molecular transport within the Langevin Dynamics (LD) framework to describe dynamic mechanical processes. This computational infrastructure allows the simulation of hours of molecular machine dynamics in reasonable wall-clock time. We apply SRDDM to test performance of the Search-and-Capture of mitotic spindle assembly by simulating, in three spatial dimensions, dynamic instability of elastic microtubules anchored in two centrosomes, movement and deformations of geometrically realistic centromeres with flexible kinetochores and chromosome arms. Furthermore, the SRDDM describes the mechanics and kinetics of Ndc80 linkers mediating transient attachments of microtubules to the chromosomal kinetochores. The rates of these attachments and detachments depend upon phosphorylation states of the Ndc80 linkers, which are regulated in the model by explicitly accounting for the reactions of Aurora A and B kinase enzymes undergoing restricted diffusion. We find that there is an optimal rate of microtubule-kinetochore detachments which maximizes the accuracy of the chromosome connections, that adding chromosome arms to kinetochores improve the accuracy by slowing down chromosome movements, that Aurora A and kinetochore deformations have a small positive effect on the attachment accuracy, and that thermal fluctuations of the microtubules increase the rates of kinetochore capture and also improve the accuracy of spindle assembly.

The blood clotting protein fibrin has extraordinary elastomeric properties due to its hierarchical structure. SAXS is combined with computational molecular modeling, providing insight in fibrin elasticity and guidelines for designing new polymers.

Fibrin is the major extracellular component of blood clots and a proteinaceous hydrogel used as a versatile biomaterial. Fibrin forms branched networks of polymeric fibers, built of laterally associated double-stranded protofibrils. This multiscale hierarchical structure is crucial for the extraordinary mechanical resilience of blood clots. Yet, the structural basis of clot mechanical properties remains largely unclear due, in part, to the unresolved molecular packing structure of fibrin fibers. Here we quantitatively assess the packing structure of fibrin fibers by combining Small Angle X-ray Scattering (SAXS) measurements of fibrin networks reconstituted under a wide range of conditions with computational molecular modeling of fibrin oligomers. The number, positions, and intensities of the Bragg peaks observed in the SAXS experiments were reproduced computationally based on the all-atom molecular structure of reconstructed fibrin protofibrils. Specifically, the model correctly predicts the intensities of the reflections of the 22.5 nm axial repeat, corresponding to the half-staggered longitudinal arrangement of fibrin molecules. In addition, the SAXS measurements showed that protofibrils within fibrin fibers have a partially ordered lateral arrangement with a characteristic transverse repeat distance of 13 nm, irrespective of the fiber thickness. These findings provide fundamental insights into the molecular structure of fibrin clots that underlies their biological and physical properties.

Given the need for modern researchers to produce open, reproducible scientific output, the lack of standards and best practices for sharing data and workflows used to produce and analyze molecular dynamics (MD) simulations have become an important issue in the field. There are now multiple well-established packages to perform molecular dynamics simulations, often highly tuned for exploiting specific classes of hardware, and each with strong communities surrounding them, but with very limited interoperability/transferability options. Thus, the choice of the software package often dictates the workflow for both simulation production and analysis. The level of detail in documenting the workflows and analysis code varies greatly in published work, hindering reproducibility of the reported results and the ability for other researchers to build on these studies. An increasing number of researchers are motivated to make their data available, but many challenges remain in order to effectively share and reuse simulation data. To discuss these and other issues related to best practices in the field in general, we organized a workshop in November 2018 ( https://bioexcel.eu/events/workshop-on-sharing-data-from-molecular-simulations/ ). Here, we present a brief overview of this workshop and topics discussed. We hope this effort will spark further conversation in the MD community to pave the way towards more open, interoperable and reproducible outputs coming from research studies using MD simulations.

We carried out dynamic force manipulations in silico on a variety of coiled-coil protein fragments from myosin, chemotaxis receptor, vimentin, fibrin, and phenylalanine zippers that vary in size and topology of their α-helical packing. When stretched along the superhelical axis, all superhelices show elastic, plastic, and inelastic elongation regimes and undergo a dynamic transition from the α-helices to the β-sheets, which marks the onset of plastic deformation. Using the Abeyaratne-Knowles formulation of phase transitions, we developed a new theoretical methodology to model mechanical and kinetic properties of protein coiled-coils under mechanical nonequilibrium conditions and to map out their energy landscapes. The theory was successfully validated by comparing the simulated and theoretical force-strain spectra. We derived the scaling laws for the elastic force and the force for α-to-β transition, which can be used to understand natural proteins’ properties as well as to rationally design novel biomaterials of required mechanical strength with desired balance between stiffness and plasticity.

This chapter demonstrates the utility of the approach to nanoindentation in silico by exploring the biomechanical properties of two specific examples of virus shell—the Cowpea chlorotic mottle virus (CCMV) capsid and the virus-like particle encapsulin. Ubiquitous in nature, viruses are biological infectious nanostructures that, in most instances, consist of just a nucleic acid packaged within a protein shell. Viruses themselves contain no metabolic machinery. Viruses appear to be ubiquitous, existing as specific infectious agents for all species of bacteria, archaebacteria, animals, and plants. Viruses that infect cells have biomedical and societal importance due to their roles as infectious agents. These viruses, plus the native cellular virus-like nanocompartment bioreactors, have potential applications in a number of areas. Recently discovered bacterial nanocompartments, called encapsulins, belong to a family of virus-like bioreactor structures. There exist a number of theoretical methods, which have been used to describe the dynamics of virus capsids.

We studied the hydrodynamic behavior of fibrinogen, a blood plasma protein involved in blood clotting, in a broad 0.3–60 mg/mL range of concentration and 5–42 °C temperature using pulsed-field gradient ¹H NMR-diffusometry. Arrhenius plots revealed the activation energy for fibrinogen diffusion E_d = 21.3 kJ/mol at 1.4 mg/mL and 28.4 kJ/mol at 38 mg/mL. We found a dramatic slowdown in fibrinogen self-diffusion with concentration beginning at 1.7–3.4 mg/mL, which deviated from the standard hard-particle behavior, suggesting a remarkable intermolecular entanglement. This concentration dependence was observed regardless of the absence or presence of the GPRP peptide (inhibitor of fibrin polymerization), and also in samples free of fibrin oligomers. By contrast, diffusivity of fibrinogen variant I-9 with truncated C-terminal portions of the Aα chains was much less concentration-dependent, indicating the importance of intermolecular linkages formed by the αC regions. Theoretical models combined with all-atom molecular dynamics simulations revealed partially bent fibrinogen solution conformations that interpolate between a flexible chain and a rigid rod observed in the crystal. The results obtained illuminate the important role of the αC regions in modulating the fibrinogen molecular shape through formation of weak intermolecular linkages that control the bulk properties of fibrinogen solutions.