Fibrin polymerization involves thrombin-mediated exposure of knobs on one monomer that bind to holes available on another, leading to the formation of fibers. In silico evidence has suggested that the classical A:a knob-hole interaction is enhanced by surrounding residues not directly involved in the binding pocket of hole a, via noncovalent interactions with knob A. We assessed the importance of extended knob-hole interactions by performing biochemical, biophysical, and in silico modeling studies on recombinant human fibrinogen variants with mutations at residues responsible for the extended interactions. Three single fibrinogen variants, γD297N, γE323Q, and γK356Q, and a triple variant γDEK (γD297N/γE323Q/γK356Q) were produced in a CHO (Chinese Hamster Ovary) cell expression system. Longitudinal protofibril growth probed by atomic force microscopy was disrupted for γD297N and enhanced for the γK356Q mutation. Initial polymerization rates were reduced for all variants in turbidimetric studies. Laser scanning confocal microscopy showed that γDEK and γE323Q produced denser clots, whereas γD297N and γK356Q were similar to wild type. Scanning electron microscopy and light scattering studies showed that fiber thickness and protofibril packing of the fibers were reduced for all variants. Clot viscoelastic analysis showed that only γDEK was more readily deformable. In silico modeling suggested that most variants displayed only slip-bond dissociation kinetics compared with biphasic catch-slip kinetics characteristics of wild type. These data provide new evidence for the role of extended interactions in supporting the classical knob-hole bonds involving catch-slip behavior in fibrin formation, clot structure, and clot mechanics.

Microtubules (MTs), a cellular structure element, exhibit dynamic instability and can switch stochastically from growth to shortening; but the factors that trigger these processes at the molecular level are not understood. We developed a 3D Microtubule Assembly and Disassembly DYnamics (MADDY) model, based upon a bead-per-monomer representation of the αβ-tubulin dimers forming an MT lattice, stabilized by the lateral and longitudinal interactions between tubulin subunits. The model was parameterized against the experimental rates of MT growth and shortening, and pushing forces on the Dam1 protein complex due to protofilaments splaying out. Using the MADDY model, we carried out GPU-accelerated Langevin simulations to access dynamic instability behavior. By applying Machine Learning techniques, we identified the MT characteristics that distinguish simultaneously all four kinetic states: growth, catastrophe, shortening, and rescue. At the cellular 25 μM tubulin concentration, the most important quantities are the MT length L, average longitudinal curvature κlong, MT tip width w, total energy of longitudinal interactions in MT lattice Ulong, and the energies of longitudinal and lateral interactions required to complete MT to full cylinder Ulongadd and Ulatadd. At high 250 μM tubulin concentration, the most important characteristics are L, κlong, number of hydrolyzed αβ-tubulin dimers nhyd and number of lateral interactions per helical pitch nlat in MT lattice, energy of lateral interactions in MT lattice Ulat, and energy of longitudinal interactions in MT tip ulong. These results allow greater insights into what brings about kinetic state stability and the transitions between states involved in MT dynamic instability behavior.

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 introduction of accelerator devices such as graphics processing units (GPUs) has had profound impact on molecular dynamics simulations and has enabled order-of-magnitude performance advances using commodity hardware. To fully reap these benefits, it has been necessary to reformulate some of the most fundamental algorithms, including the Verlet list, pair searching, and cutoffs. Here, we present the heterogeneous parallelization and acceleration design of molecular dynamics implemented in the GROMACS codebase over the last decade. The setup involves a general cluster-based approach to pair lists and non-bonded pair interactions that utilizes both GPU and central processing unit (CPU) single instruction, multiple data acceleration efficiently, including the ability to load-balance tasks between CPUs and GPUs. The algorithm work efficiency is tuned for each type of hardware, and to use accelerators more efficiently, we introduce dual pair lists with rolling pruning updates. Combined with new direct GPU-GPU communication and GPU integration, this enables excellent performance from single GPU simulations through strong scaling across multiple GPUs and efficient multi-node parallelization.

Fibrin is the major extracellular component of blood clots and a proteinaceous hydrogel used as a versatile biomaterial. Fibrin forms branched networks 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 of fibrin fibers. Here the packing structure of fibrin fibers is quantitatively assessed by combining Small Angle X-ray Scattering (SAXS) measurements of fibrin reconstituted under a wide range of conditions with computational molecular modeling of fibrin protofibrils. 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.

Thirteen tubulin protofilaments, made of αβ-tubulin heterodimers, interact laterally to produce cytoskeletal microtubules. Microtubules exhibit the striking property of dynamic instability, manifested in their intermittent growth and shrinkage at both ends. This behavior is key to many cellular processes, such as cell division, migration, maintenance of cell shape, etc. Although assembly and disassembly of microtubules is known to be linked to hydrolysis of a guanosine triphosphate molecule in the pocket of β-tubulin, detailed mechanistic understanding of corresponding conformational changes is still lacking. Here we take advantage of the recent generation of in-microtubule structures of tubulin to examine the properties of protofilaments, which serve as important microtubule assembly and disassembly intermediates. We find that initially straight tubulin protofilaments, relax to similar non-radially curved and slightly twisted conformations. Our analysis further suggests that guanosine triphosphate hydrolysis primarily affects the flexibility and conformation of the inter-dimer interface, without a strong impact on the shape or flexibility of αβ-heterodimer. Inter-dimer interfaces are significantly more flexible compared to intra-dimer interfaces. We argue that such a difference in flexibility could be key for distinct stability of the plus and minus microtubule ends. The higher flexibility of the inter-dimer interface may have implications for development of pulling force by curving tubulin protofilaments during microtubule disassembly, a process of major importance for chromosome motions in mitosis.

AbstractEven when modern computational platforms and parallel techniques are used, conventional all-atom simulations are limited both in terms of reachable timescale and number of atoms in the biomolecular system of interest. On the other hand, coarse-grained models, which allow to overcome this limitation, rely on proper and rigorous parametrization of the underlying force field. Here, we present a novel iterative approach for parametrization of coarse-grained models based on direct comparison of equilibrium simulations at all-atom and coarse-grained resolutions. In order to assess the accuracy of our method, we have built and parametrized an elastic network model (ENM) of the tubulin protofilament consisting of four monomers. For this system, our method shows good convergence and the parametrized ENM reproduces protein dynamics in a finer way when compared to ENMs parametrized using the conventional approach. The presented method can be extended to other coarse-grained models with a slight adjustment of the equations describing the iterative scheme.

The lifetimes of non-covalent A:a knob-hole bonds in fibrin probed with the optical trap-based force-clamp first increases ("catch bonds") and then decreases ("slip bonds") with increasing tensile force. Molecular modeling of "catch-to-slip" transition using the atomic structure of the A:a complex reveals that the movable flap serves as tension-dependent molecular switch. Flap dissociation from the regulatory B-domain in $\gamma$-nodule and translocation from the periphery to knob `A' triggers the hole `a' closure and interface remodeling, which results in the increased binding affinity and prolonged bond lifetimes. Fluctuating bottleneck theory is developed to understand the "catch-to-slip" transition in terms of the interface stiffness $\kappa =$ 15.7 pN nm $^{-1}$, interface size fluctuations 0.7-2.7 nm, knob `A' escape rate constant $k_0 =$ 0.11 nm$^2$ s$^{-1}$, and transition distance for dissociation $\sigma_y =$ 0.25 nm. Strengthening of the A:a knob-hole bonds under small tension might favor formation and reinforcement of nascent fibrin clots under hydrodynamic shear.

Fibrinogen is a blood plasma protein that polymerizes to form the fibrin clot on cleavage of its fibrinopeptides. This work provides quantitative characteristics of the molecular hydrodynamics of fibrinogen in a broad 0.3-60 mg/ml range of concentration and 5-42°C temperature obtained using pulsed-field gradient 1H NMR. Arrhenius plots revealed the activation energy for fibrinogen diffusion Ed = 21.3 kJ/mol at 1.4 mg/ml and 28.4 kJ/mol at 38 mg/ml. The diffusive motion of fibrinogen underwent a remarkable slowdown with concentrations beginning at 1.7-3.4 mg/ml, which deviated from the standard hard-particle behavior, suggesting concentration-dependent intermolecular entanglement. 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. The remarkable concentration dependence was observed regardless of the absence or presence of the GPRP peptide (inhibitor of fibrin polymerization) and in samples free of fibrin oligomers, confirming that the observed steep decrease in fibrinogen diffusivity was not due to potential contamination by oligomeric fibrin-fibrinogen complexes. Theoretical models combined with all-atom Molecular Dynamics simulations revealed that fibrinogen in solution has a bendable conformation that interpolates 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. The results emphasize the importance of fibrinogen self-assembly in modulating the hydrodynamic properties of semi-diluted and concentrated fibrinogen preparations, and point to the potential role of fibrinogen self-assembly in biological and clinical applications.

Fibrin formation and mechanical stability are essential in thrombosis and hemostasis. To reveal how mechanical load impacts fibrin, we carried out optical trap-based single-molecule forced unbinding experiments. The strength of noncovalent A:a knob-hole bond stabilizing fibrin polymers first increases with tensile force (catch bonds) and then decreases with force when the force exceeds a critical value (slip bonds). To provide the structural basis of catch-slip-bond behavior, we analyzed crystal structures and performed molecular modeling of A:a knob-hole complex. The movable flap (residues γ295 to γ305) containing the weak calcium-binding site γ2 serves as a tension sensor. Flap dissociation from the B domain in the γ-nodule and translocation to knob 'A' triggers hole 'a' closure, resulting in the increase of binding affinity and prolonged bond lifetimes. The discovery of biphasic kinetics of knob-hole bond rupture is quantitatively explained by using a theory, formulated in terms of structural transitions in the binding pocket between the low-affinity (slip) and high-affinity (catch) states. We provide a general framework to understand the mechanical response of protein pairs capable of tension-induced remodeling of their association interface. Strengthening of the A:a knob-hole bonds at 30- to 40-pN forces might favor formation of nascent fibrin clots subject to hydrodynamic shear in vivo.

The space-filling fibrin network is a major part of clots and thrombi formed in blood. Fibrin polymerization starts when fibrinogen, a plasma protein, is proteolytically converted to fibrin, which self-assembles to form double-stranded protofibrils. When reaching a critical length, these intermediate species aggregate laterally to transform into fibers arranged into branched fibrin network. We combined multiscale modeling in silico with atomic force microscopy (AFM) imaging to reconstruct complete atomic models of double-stranded fibrin protofibrils with γ-γ crosslinking, A:a and B:b knob-hole bonds, and αC regions-all important structural determinants not resolved crystallographically. Structures of fibrin oligomers and protofibrils containing up to 19 monomers were successfully validated by quantitative comparison with high-resolution AFM images. We characterized the protofibril twisting, bending, kinking, and reversibility of A:a knob-hole bonds, and calculated hydrodynamic parameters of fibrin oligomers. Atomic structures of protofibrils provide a basis to understand mechanisms of early stages of fibrin polymerization.

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.

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 1H NMR-diffusometry. Arrhenius plots revealed the activation energy for fibrinogen diffusion Ed = 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.

Antibodies to DNA play an important role in the pathogenesis of autoimmune diseases. The elucidation of structural mechanisms of both the antigen recognition and the interaction of anti-DNA antibodies with DNA will help to understand the role of DNA-containing immune complexes in various pathologies and can provide a basis for new treatment modalities. Moreover, the DNA-antibody complex is an analog of specific intracellular DNA-protein interactions. In this work, we used in silico molecular dynamic simulations of bimolecular complexes of the dsDNA segment containing the Fab fragment of an anti-DNA antibody to obtain the detailed thermodynamic and structural characteristics of dynamic intermolecular interactions. Using computationally modified crystal structure of the Fab-DNA complex (PDB ID: 3VW3), we studied the equilibrium molecular dynamics of the 64M-5 antibody Fab fragment associated with the dsDNA fragment containing the thymine dimer, the product of DNA photodamage. Amino acid residues that constitute paratopes and the complementary nucleotide epitopes for the Fab-DNA construct were identified. Stacking and electrostatic interactions were found to play the main role in mediating the most specific antibodydsDNA contacts, while hydrogen bonds were less significant. These findings may shed light on the formation and properties of pathogenic anti-DNA antibodies in autoimmune diseases, such as systemic lupus erythematosus associated with skin photosensitivity and DNA photodamage.

Fibrin is a filamentous network made in blood to stem bleeding; it forms when fibrinogen is converted into fibrin monomers that self-associate into oligomers and then to polymers. To gather structural insights into fibrin formation and properties, we combined high-resolution atomic force microscopy of fibrin(ogen) oligomers and molecular modeling of crystal structures of fibrin(ogen) and its fragments. We provided a structural basis for the intermolecular flexibility of single-stranded fibrin(ogen) oligomers and identified a hinge region at the D:D inter-monomer junction. Following computational reconstruction of the missing portions, we recreated the full-atomic structure of double-stranded fibrin oligomers that was validated by quantitative comparison with the experimental images. We characterized previously unknown intermolecular binding contacts at the D:D and D:E:D interfaces, which drive oligomerization and reinforce the intra-and inter-strand connections in fibrin besides the known knob-hole bonds. The atomic models provide valuable insights into the submolecular mechanisms of fibrin polymerization.

Fibrin is a transient extracellular matrix that forms upon vascular injury as a first step towards wound healing. Fibrin fibers can be extended several times their own length without rupture, though the origin of this extraordinary extensibility is still poorly understood. In this chapter, we study the origin of fibrin extensibility by using small-angle X-ray scattering (SAXS) in combination with tensile tests. SAXS measurements provide insight into both the axial and lateral molecular packing structure of the fibers. We show that the degree of axial and radial packing order is dependent on fibrin bundle size, which we relate to the mesh size of the network that puts an upper bound on the range of crystalline order. By comparing our SAXS measurements with full-atom simulations of protofibrils, we show that the second order reflection of the half-staggering axial repeat distance is suppressed due to the symmetric structure of the fibrinogen monomer. Upon stretching, however, several new peaks appear, including the second order reflection. At high strains of ∼100%, the peaks indicative of axial order disappear, showing that stretching causes forced molecular unfolding. To trace the molecular origin of these changes, we compare our SAXS results with predicted SAXS spectra based on full-atom simulations of stretched protofibrils. Taken together, our experiments and simulations directly show that protofibril elongation is mediated by 1) elongation of the alpha-helical coiled-coil regions and 2) γ-nodule unfolding, which starts at 30% strain. These results are important to understand the reduced extensibility in case of disease-related fibrinogen mutations and the influence of FXIIIa crosslinking on fibrin fiber extensibility.

Antibodies against double-stranded DNA play an important role in the pathogenesis of autoimmune diseases. Structural analysis of antibody-DNA complexes contributes to our understanding of the role of DNA-containing immune complexes in human pathologies and may help in designing novel treatments. In this paper, we study dynamics of the full-atomic structure of the molecular complexes formed by an antibody Fab fragment with double-stranded DNA, containing or not containing a thymine dimer. Molecular dynamics simulations are used in conjunction with the Principle Component Analysis technique. We found that removing a covalent bond from the thymine dimer results in changes of the structural dynamics of the light and heavy chains of Fab as well as in the DNA strands. A significant increase in mobility of the Fab light chain was observed throughout the entire simulation runs with a higher amplitude of fluctuations at the interface with DNA. Essential dynamics analysis of simulation trajectories of the antibody-dsDNA complexes shows that fluctuations in the low-frequency eigenvectors are localized at the ends of the DNA sequences, suggesting that these bending motions are important for the DNA-antibody interactions.