CellDynaMo-stochastic reaction-diffusion-dynamics model: Application to search-and-capture process of mitotic spindle assembly.

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.

Fibrin protofibril packing and clot stability are enhanced by extended knob-hole interactions and catch-slip bonds.

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.

Microtubule assembly and disassembly dynamics model: Exploring dynamic instability and identifying features of Microtubules' Growth, Catastrophe, Shortening, and Rescue.

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--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.

Heterogeneous parallelization and acceleration of molecular dynamics simulations in GROMACS.

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.

Sharing Data from Molecular Simulations.

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 has 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, 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.

Mechanical properties of tubulin intra- and inter-dimer interfaces and their implications for microtubule dynamic instability.

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.

Regulatory element in fibrin triggers tension-activated transition from catch to slip bonds.

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).

Atomic Structural Models of Fibrin Oligomers.

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.

Dynamic Transition from α-Helices to β-Sheets in Polypeptide Coiled-Coil Motifs.

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.

Conformational Flexibility and Self-Association of Fibrinogen in Concentrated Solutions.

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 = 21.

Structural Basis of Interfacial Flexibility in Fibrin Oligomers.

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.

Botulinum neurotoxin: unique folding of enzyme domain of the most-poisonous poison.

Botulinum neurotoxin (BoNT), the most toxic substance known to mankind, is the first example of the fully active molten globule state. To understand its folding mechanism, we performed urea denaturation experiments and theoretical modeling using BoNT serotype A (BoNT/A). We found that the extent of BoNT/A denaturation from the native state (N) shows a nonmonotonic dependence on urea concentration indicating a unique multistep denaturation process, N → I1 [Formula: see text] I2 [Formula: see text] U, with two intermediate states I1 and I2.

Molecular mechanisms, thermodynamics, and dissociation kinetics of knob-hole interactions in fibrin.

Polymerization of fibrin, the primary structural protein of blood clots and thrombi, occurs through binding of knobs 'A' and 'B' in the central nodule of fibrin monomer to complementary holes 'a' and 'b' in the γ- and β-nodules, respectively, of another monomer. We characterized the A:a and B:b knob-hole interactions under varying solution conditions using molecular dynamics simulations of the structural models of fibrin(ogen) fragment D complexed with synthetic peptides GPRP (knob 'A' mimetic) and GHRP (knob 'B' mimetic). The strength of A:a and B:b knob-hole complexes was roughly equal, decreasing with pulling force; however, the dissociation kinetics were sensitive to variations in acidity (pH 5-7) and temperature (T = 25-37 °C).

Mechanical transition from α-helical coiled coils to β-sheets in fibrin(ogen).

We characterized the α-to-β transition in α-helical coiled-coil connectors of the human fibrin(ogen) molecule using biomolecular simulations of their forced elongation and theoretical modeling. The force (F)-extension (X) profiles show three distinct regimes: (1) the elastic regime, in which the coiled coils act as entropic springs (F < 100-125 pN; X < 7-8 nm); (2) the constant-force plastic regime, characterized by a force-plateau (F ≈ 150 pN; X ≈ 10-35 nm); and (3) the nonlinear regime (F > 175-200 pN; X > 40-50 nm). In the plastic regime, the three-stranded α-helices undergo a noncooperative phase transition to form parallel three-stranded β-sheets.

Comparative simulation studies of native and single-site mutant human beta-defensin-1 peptides.

Human defensins play important roles in a broad range of biological functions, such as microbial defense and immunity. Yet, little is known about their molecular properties, i.e.

Multiscale modeling of the nanomechanics of microtubule protofilaments.

Large-size biomolecular systems that spontaneously assemble, disassemble, and self-repair by controlled inputs play fundamental roles in biology. Microtubules (MTs), which play important roles in cell adhesion and cell division, are a prime example. MTs serve as ″tracks″ for molecular motors, and their biomechanical functions depend on dynamic instability-a stochastic switching between periods of rapid growing and shrinking.

Mechanism of fibrin(ogen) forced unfolding.

Fibrinogen, upon enzymatic conversion to monomeric fibrin, provides the building blocks for fibrin polymer, the scaffold of blood clots and thrombi. Little has been known about the force-induced unfolding of fibrin(ogen), even though it is the foundation for the mechanical and rheological properties of fibrin, which are essential for hemostasis. We determined mechanisms and mapped the free energy landscape of the elongation of fibrin(ogen) monomers and oligomers through combined experimental and theoretical studies of the nanomechanical properties of fibrin(ogen), using atomic force microscopy-based single-molecule unfolding and simulations in the experimentally relevant timescale.

Exploring the mechanical stability of the C2 domains in human synaptotagmin 1.

Human synaptotagmin 1 (Syt1) plays a crucial role in the bending of the membrane during neurotransmitter release at the synapse. Hence, resolving the structural details of Syt1 that underlie its biological function is fundamental for providing mechanistic insights into the nature of the synaptic response. We explored the unfolding micromechanics of Syt1 by analyzing the free energy landscape of the whole molecule and its C2A and C2B domains.