Generation of random numbers on graphics processors: forced indentation in silico of the bacteriophage HK97

Virtual indentation of the empty capsid of the minute virus of mice using a minimal coarse-grained model

A minimal coarse-grained model for T=1 viral capsids assembled from 20 protein rigid trimers has been designed by extending a previously proposed form of the interaction energy written as a sum of anisotropic pairwise interactions between the trimeric capsomers. The extension of the model has been performed to properly account for the coupling between two internal coordinates: the one that measures the intercapsomer distance and the other that gives the intercapsomer dihedral angle. The model has been able to fit with less than a 10% error the atomic force microscopy (AFM) indentation experimental data for the empty capsid of the minute virus of mice (MVM), providing in this way an admissible picture of the main mechanisms behind the capsid deformations. In this scenario, the bending of the intercapsomer dihedral angle is the angular internal coordinate that can support larger deformations away from its equilibrium values, determining important features of the AFM indentation experiments as the elastic constants along the three symmetry axes of the capsid and the critical indentations. From the value of one of the parameters of our model, we conclude that trimers in the MVM must be quite oblate tops, in excellent agreement with their known structure. The transition from the linear to the nonlinear regimes sampled in the indentation process appears to be an interesting topic for future research in physical virology.

Mechanical fatigue testing in silico: Dynamic evolution of material properties of nanoscale biological particles

Biological particles have evolved to possess mechanical characteristics necessary to carry out their functions. We developed a computational approach to "fatigue testing in silico", in which constant-amplitude cyclic loading is applied to a particle to explore its mechanobiology. We used this approach to describe dynamic evolution of nanomaterial properties and low-cycle fatigue in the thin spherical encapsulin shell, thick spherical Cowpea Chlorotic Mottle Virus (CCMV) capsid, and thick cylindrical microtubule (MT) fragment over 20 cycles of deformation. Changing structures and force-deformation curves enabled us to describe their damage-dependent biomechanics (strength, deformability, stiffness), thermodynamics (released and dissipated energies, enthalpy, and entropy) and material properties (toughness). Thick CCMV and MT particles experience material fatigue due to slow recovery and damage accumulation over 3-5 loading cycles; thin encapsulin shells show little fatigue due to rapid remodeling and limited damage. The results obtained challenge the existing paradigm: damage in biological particles is partially reversible owing to particle's partial recovery; fatigue crack may or may not grow with each loading cycle and may heal; and particles adapt to deformation amplitude and frequency to minimize the energy dissipated. Using crack size to quantitate damage is problematic as several cracks might form simultaneously in a particle. Dynamic evolution of strength, deformability, and stiffness, can be predicted by analyzing the cycle number (N) dependent damage, [Formula: see text] , where α is a power law and Nf is fatigue life. Fatigue testing in silico can now be used to explore damage-induced changes in the material properties of other biological particles. STATEMENT OF SIGNIFICANCE: Biological particles possess mechanical characteristics necessary to perform their functions. We developed "fatigue testing in silico" approach, which employes Langevin Dynamics simulations of constant-amplitude cyclic loading of nanoscale biological particles, to explore dynamic evolution of the mechanical, energetic, and material properties of the thin and thick spherical particles of encapsulin and Cowpea Chlorotic Mottle Virus, and the microtubule filament fragment. Our study of damage growth and fatigue development challenge the existing paradigm. Damage in biological particles is partially reversible as fatigue crack might heal with each loading cycle. Particles adapt to deformation amplitude and frequency to minimize energy dissipation. The evolution of strength, deformability, and stiffness, can be accurately predicted by analyzing the damage growth in particle structure.

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

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. 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 CellDynaMo package models, in 3D, any cellular subsystem where sufficient detail of the macromolecular players and the kinetics of relevant reactions are available. The package is based on the Stochastic Reaction-Diffusion-Dynamics model that combines the stochastic description of chemical kinetics, Brownian diffusion-based description of molecular transport, and Langevin dynamics-based representation of mechanical processes most pertinent to the system. We apply the model to test the Search-and-Capture mechanism of mitotic spindle assembly. We find that there is an optimal rate of microtubule-kinetochore detachments which maximizes the accuracy of chromosome connections, that chromosome arms improve the attachment accuracy by slowing down chromosome movements, that Aurora A kinase and kinetochore deformations have small positive effects on the accuracy, and that thermal fluctuations of the microtubules increase the rates of kinetochore capture and also improve the accuracy.

Single-molecule conductance of double-stranded RNA oligonucleotides

RNA oligonucleotides are crucial for a range of biological functions and in many biotechnological applications. Herein, we measured, for the first time, the conductance of individual double-stranded (ds)RNA molecules and compared it with the conductance of single DNA : RNA hybrids. The average conductance values are similar for both biomolecules, but the distribution of conductance values shows an order of magnitude higher variability for dsRNA, indicating higher molecular flexibility of dsRNA. Microsecond Molecular Dynamics simulations explain this difference and provide structural insights into the higher stability of DNA : RNA duplex with atomic level of detail. The rotations of 2'-OH groups of the ribose rings and the bases in RNA strands destabilize the duplex structure by weakening base stacking interactions, affecting charge transport, and making single-molecule conductance of dsRNA more variable (dynamic disorder). The results demonstrate that a powerful combination of state-of-the-art biomolecular electronics techniques and computational approaches can provide valuable insights into biomolecules' biophysics with unprecedented spatial resolution.

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-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 latticeU long , and the energies of longitudinal and lateral interactions required to complete MT to full cylinderU long add andU lat add . At high 250 μM tubulin concentration, the most important characteristics are L ,κ long , number of hydrolyzed αβ-tubulin dimersn hyd and number of lateral interactions per helical pitchn lat in MT lattice, energy of lateral interactions in MT latticeU lat , and energy of longitudinal interactions in MT tipu long . These results allow greater insights into what brings about kinetic state stability and the transitions between states involved in MT dynamic instability behavior.

Strength, deformability and toughness of uncrosslinked fibrin fibers from theoretical reconstruction of stress-strain curves

Structural mechanisms underlying the mechanical properties of fibrin fibers are elusive. We combined tensile testing of uncrosslinked fibrin polymers in vitro and in silico to explore their material properties. The experimental stress (σ) - strain (ε) curves for fibrin fibers are characterized by elastic deformations with a weaker elastic response for ε<160% due to unraveling of αC tethers and straightening of fibrin protofibrils, and a stronger response for ε>160% owing to unfolding of the coiled coils and γ nodules in fibrin monomers. Fiber rupture for strains ε>212% is due to dissociation of the knob-hole bonds and rupture of D:D interfaces. We developed the Fluctuating Bilinear Spring model to interpret the σ-ε profiles in terms of the free energy for protofibril alignment ΔG0 = 10.1-11.5 kBT, Young's moduli for protofibril alignment Yu = 1.9-3.2 MPa and stretching Ya = 5.7-9.7 MPa, strain scale ε˜≈ 12-40% for fiber rupture, and protofibril cooperativity m= 3.6-8. We applied the model to characterize the fiber strength σcr≈ 12-13 MPa, deformability εcr≈ 222%, and rupture toughness U≈ 9 MJ/m3, and to resolve thermodynamic state functions, 96.9 GJ/mol entropy change for protofibril alignment (at room temperature) and 113.6 GJ/mol enthalpy change for protofibril stretching, which add up to 210.5 GJ/mol free-energy change. Fiber elongation is associated with protofibril dehydration and sliding mechanism to create an ordered protofibril array. Fibrin fibers behave like a hydrogel; protofibril dehydration and water expulsion account for ∼94-98% of the total free-energy changes for fiber elongation and rupture. STATEMENT OF SIGNIFICANCE: Structural mechanisms underlying the mechanical properties of fibrin fibers, major components of blood clots and obstructive thrombi, are elusive. We performed tensile testing of uncrosslinked fibrin polymers in vitro and in silico to explore their material properties. Fluctuating Bilinear Spring theory was developed to interpret the stress-strain profiles in terms of the energy for protofibril alignment, elastic moduli for protofibril alignment and stretching, and strain scale for fiber rupture, and to probe the limits of fiber strength, extensibility and toughness. Fibrin fibers behave like a hydrogel. Fiber elongation is defined by the protofibril dehydration and sliding. Structural rearrangements in water matrix control fiber elasticity. These results contribute to fundamental understanding of blood clot breakage that underlies thrombotic embolization.

Minimal Design Principles for Icosahedral Virus Capsids

The geometrical structures of single- and multiple-shell icosahedral virus capsids are reproduced as the targets that minimize the cost corresponding to relatively simple design functions. Capsid subunits are first identified as building blocks at a given coarse-grained scale and then represented in these functions as point particles located on an appropriate number of concentric spherical surfaces. Minimal design cost is assigned to optimal spherical packings of the particles. The cost functions are inspired by the packings favored for the Thomson problem, which minimize the electrostatic potential energy between identical charged particles. In some cases, icosahedral symmetry constraints are incorporated as external fields acting on the particles. The simplest cost functions can be obtained by separating particles in disjoint nonequivalent sets with distinct interactions, or by introducing interacting holes (the absence of particles). These functions can be adapted to reproduce any capsid structure found in real viruses. Structures absent in Nature require significantly more complex designs. Measures of information content and complexity are assigned to both the cost functions and the capsid geometries. In terms of these measures, icosahedral structures and the corresponding cost functions are the simplest solutions.

The influence of random number generation in dissipative particle dynamics simulations using a cryptographic hash function

The tiny encryption algorithm (TEA) is widely used when performing dissipative particle dynamics (DPD) calculations in parallel, usually on distributed memory systems. In this research, we reduced the computational cost of the TEA hash function and investigated the influence of the quality of the random numbers generated on the results of DPD calculations. It has already been established that the randomness, or quality, of the random numbers depend on the number of processes from internal functions such as SHIFT, XOR and ADD, which are commonly referred to as “rounds”. Surprisingly, if we choose seed numbers from high entropy sources, with a minimum number of rounds, the quality of the random numbers generated is sufficient to successfully perform accurate DPD simulations. Although it is well known that using a minimal number of rounds is insufficient for generating high-quality random numbers, the combination of selecting good seed numbers and the robustness of DPD simulations means that we can reduce the random number generation cost without reducing the accuracy of the simulation results.

A minimal coarse-grained model for the low-frequency normal mode analysis of icosahedral viral capsids

The main goal of this work is the design of a coarse-grained theoretical model of minimal resolution for the study of the physical properties of icosahedral virus capsids within the linear-response regime. In this model the capsid is represented as an interacting many-body system whose composing elements are capsid subunits (capsomers), which are treated as three-dimensional rigid bodies. The total interaction potential energy is written as a sum of pairwise capsomer-capsomer interactions. Based on previous work [Gomez Llorente et al., Soft Matter, 2014, 10, 3560], a minimal and complete anisotropic binary interaction that includes a full Hessian matrix of independent force constants is proposed. In this interaction model, capsomers have rotational symmetry around an axis of order n > 2. The full coarse-grained model is applied to analyse the low-frequency normal-mode spectrum of icosahedral T = 1 capsids. The model performance is evaluated by fitting its predicted spectrum to the full-atom results for the Satellite Tobacco Necrosis Virus (STNV) capsid [Dykeman and Sankey, Phys. Rev. Lett., 2008, 100, 028101]. Two capsomer choices that are compatible with the capsid icosahedral symmetry are checked, namely pentamers (n = 5) and trimers (n = 3). Both subunit types provide fair fits, from which the magnitude of the coarse-grained force constants for a real virus is obtained. The model is able to uncover latent instabilities whose analysis is fully consistent with the current knowledge about the STNV capsid, which does not self-assemble in the absence of RNA and is thermally unstable. The straightforward generalisability of the model beyond the linear regime and its completeness make it a promising tool to theoretically interpret many experimental data such as those provided by the atomic force microscopy or even to better understand processes far from equilibrium such as the capsid self-assembly.

Kinetics of empty viral capsid assembly in a minimal model

The efficient construction of a protective protein shell or capsid is one of the most crucial steps in the replication cycle of a virus. The formation of the simplest capsid typically proceeds by the spontaneous assembly of identical building blocks. This process can also be achieved in vitro even in the absence of genetic material, thus opening the door to the production of artificial viral cages for a myriad of applications. In this work, we analyze the efficiency and the kinetic peculiarities of this self-assembly process using Brownian Dynamics simulations. We use a minimal model that considers identical assembly units and is able to reproduce successfully the correct final architecture of spherical capsids. The selection of a specific size and structure is achieved by changing a single parameter that imposes an angular anisotropy on the interaction. We analyze how the geometrical constraints of the interaction affect the efficiency of the assembly. We find that the optimal conditions for an efficient assembly from a kinetic point of view strongly depart from the lowest capsid energy corresponding to the minimum of the potential energy landscape. Our work illustrates the important differences between the equilibrium and dynamic characteristics of viral self-assembly, and provides important insights on how to design specific interactions for a successful assembly of artificial viral cages.

Viral nanomechanics with a virtual atomic force microscope

One of the most important components of a virus is the protein shell or capsid that encloses its genetic material. The main role of the capsid is to protect the viral genome against external aggressions, facilitating its safe and efficient encapsulation and delivery. As a consequence, viral capsids have developed astonishing mechanical properties that are crucial for viral function. These remarkable properties have started to be unveiled in single-virus nanoindentation experiments, and are opening the door to the use of viral-derived artificial nanocages for promising bio- and nano-technological applications. However, the interpretation of nanoindentation experiments is often difficult, requiring the support of theoretical and simulation analysis. Here we present a 'Virtual AFM' (VAFM), a Brownian Dynamics simulation of a coarse-grained model of virus aimed to mimic the standard setup of atomic force microscopy (AFM) nanoindentation experiments. Despite the heavy level of coarse-graining, these simulations provide valuable information which is not accessible in experiments. Rather than focusing on a specific virus, the VAFM will be used to analyze how the mechanical response and breaking of viruses depend on different parameters controlling the effective interactions between capsid's structural units. In particular, we will discuss the influence of adsorption, the tip radius, and the rigidity and shape of the shell on its mechanical response.

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

TensorCalculator: exploring the evolution of mechanical stress in the CCMV capsid

A new computational methodology for the accurate numerical calculation of the Cauchy stress tensor, stress invariants, principal stress components, von Mises and Tresca tensors is developed. The methodology is based on the atomic stress approach which permits the calculation of stress tensors, widely used in continuum mechanics modeling of materials properties, using the output from the MD simulations of discrete atomic and [Formula: see text]-based coarse-grained structural models of biological particles. The methodology mapped into the software package TensorCalculator was successfully applied to the empty cowpea chlorotic mottle virus (CCMV) shell to explore the evolution of mechanical stress in this mechanically-tested specific example of a soft virus capsid. We found an inhomogeneous stress distribution in various portions of the CCMV structure and stress transfer from one portion of the virus structure to another, which also points to the importance of entropic effects, often ignored in finite element analysis and elastic network modeling. We formulate a criterion for elastic deformation using the first principal stress components. Furthermore, we show that von Mises and Tresca stress tensors can be used to predict the onset of a viral capsid's mechanical failure, which leads to total structural collapse. TensorCalculator can be used to study stress evolution and dynamics of defects in viral capsids and other large-size protein assemblies.

Dynamics of microtubules: highlights of recent computational and experimental investigations

Microtubules are found in most eukaryotic cells, with homologs in eubacteria and archea, and they have functional roles in mitosis, cell motility, intracellular transport, and the maintenance of cell shape. Numerous efforts have been expended over the last two decades to characterize the interactions between microtubules and the wide variety of microtubule associated proteins that control their dynamic behavior in cells resulting in microtubules being assembled and disassembled where and when they are required by the cell. We present the main findings regarding microtubule polymerization and depolymerization and review recent work about the molecular motors that modulate microtubule dynamics by inducing either microtubule depolymerization or severing. We also discuss the main experimental and computational approaches used to quantify the thermodynamics and mechanics of microtubule filaments.

SOP-GPU: influence of solvent-induced hydrodynamic interactions on dynamic structural transitions in protein assemblies

Hydrodynamic interactions (HI) are incorporated into Langevin dynamics of the Cα -based protein model using the Truncated Expansion approximation (TEA) to the Rotne-Prager-Yamakawa diffusion tensor. Computational performance of the obtained GPU realization demonstrates the model's capability for describing protein systems of varying complexity (10(2) -10(5) residues), including biological particles (filaments, virus shells). Comparison of numerical accuracy of the TEA versus exact description of HI reveals similar results for the kinetics and thermodynamics of protein unfolding. The HI speed up and couple biomolecular transitions through cross-communication among protein domains, which result in more collective displacements of structure elements governed by more deterministic (less variable) dynamics. The force-extension/deformation spectra from nanomanipulations in silico exhibit sharper force signals that match well the experimental profiles. Hence, biomolecular simulations without HI overestimate the role of tension/stress fluctuations. Our findings establish the importance of incorporating implicit water-mediated many-body effects into theoretical modeling of dynamic processes involving biomolecules. © 2016 Wiley Periodicals, Inc.

Fluctuating Nonlinear Spring Model of Mechanical Deformation of Biological Particles

The mechanical properties of virus capsids correlate with local conformational dynamics in the capsid structure. They also reflect the required stability needed to withstand high internal pressures generated upon genome loading and contribute to the success of important events in viral infectivity, such as capsid maturation, genome uncoating and receptor binding. The mechanical properties of biological nanoparticles are often determined from monitoring their dynamic deformations in Atomic Force Microscopy nanoindentation experiments; but a comprehensive theory describing the full range of observed deformation behaviors has not previously been described. We present a new theory for modeling dynamic deformations of biological nanoparticles, which considers the non-linear Hertzian deformation, resulting from an indenter-particle physical contact, and the bending of curved elements (beams) modeling the particle structure. The beams’ deformation beyond the critical point triggers a dynamic transition of the particle to the collapsed state. This extreme event is accompanied by a catastrophic force drop as observed in the experimental or simulated force (F)-deformation (X) spectra. The theory interprets fine features of the spectra, including the nonlinear components of the FX-curves, in terms of the Young’s moduli for Hertzian and bending deformations, and the structural damage dependent beams’ survival probability, in terms of the maximum strength and the cooperativity parameter. The theory is exemplified by successfully describing the deformation dynamics of natural nanoparticles through comparing theoretical curves with experimental force-deformation spectra for several virus particles. This approach provides a comprehensive description of the dynamic structural transitions in biological and artificial nanoparticles, which is essential for their optimal use in nanotechnology and nanomedicine applications.

Dynamic force experiments, which have become available to explore the physical properties of biological assemblies, oftentimes reveal results that are difficult to understand without theoretical framework. We employed a multiscale modeling approach—a combination of Molecular Dynamics simulations of atomic structures with Langevin simulations of coarse-grained models of virus shells—to characterize the degrees of freedom defining the deformation and structural collapse of biological particles tested mechanically. This enabled us to develop an analytical model that provides meaningful interpretation of force-deformation spectra available from single-particle nanoindentation experiments. The Fluctuating Nonlinear Spring (FNS) model of uniaxial particle’s deformation captures essential features of the force-deformation spectra as observed in nanomanipulations in vitro and in silico: initial non-linearity, then a subsequent force decrease transition due to structural collapse. Our theory uniquely combines the elements of continuum mechanics with the statistics of extremes, enabling one to gather mechanical and statistical characteristics of nanoparticles, which determine the Hertzian deformation of the particle’s protein layer, and bending deformation and structural damage to the particle structure. We have demonstrated how the FNS theory can accurately model the deformation of several viral shells, showing promising model applications for describing a variety of natural and synthetic nanoparticles.

Tubulin bond energies and microtubule biomechanics determined from nanoindentation in silico

Microtubules, the primary components of the chromosome segregation machinery, are stabilized by longitudinal and lateral noncovalent bonds between the tubulin subunits. However, the thermodynamics of these bonds and the microtubule physicochemical properties are poorly understood. Here, we explore the biomechanics of microtubule polymers using multiscale computational modeling and nanoindentations in silico of a contiguous microtubule fragment. A close match between the simulated and experimental force–deformation spectra enabled us to correlate the microtubule biomechanics with dynamic structural transitions at the nanoscale. Our mechanical testing revealed that the compressed MT behaves as a system of rigid elements interconnected through a network of lateral and longitudinal elastic bonds. The initial regime of continuous elastic deformation of the microtubule is followed by the transition regime, during which the microtubule lattice undergoes discrete structural changes, which include first the reversible dissociation of lateral bonds followed by irreversible dissociation of the longitudinal bonds. We have determined the free energies of dissociation of the lateral (6.9 ± 0.4 kcal/mol) and longitudinal (14.9 ± 1.5 kcal/mol) tubulin–tubulin bonds. These values in conjunction with the large flexural rigidity of tubulin protofilaments obtained (18,000–26,000 pN·nm²) support the idea that the disassembling microtubule is capable of generating a large mechanical force to move chromosomes during cell division. Our computational modeling offers a comprehensive quantitative platform to link molecular tubulin characteristics with the physiological behavior of microtubules. The developed in silico nanoindentation method provides a powerful tool for the exploration of biomechanical properties of other cytoskeletal and multiprotein assemblies.

Structural transitions and energy landscape for Cowpea Chlorotic Mottle Virus capsid mechanics from nanomanipulation in vitro and in silico

Physical properties of capsids of plant and animal viruses are important factors in capsid self-assembly, survival of viruses in the extracellular environment, and their cell infectivity. Combined AFM experiments and computational modeling on subsecond timescales of the indentation nanomechanics of Cowpea Chlorotic Mottle Virus capsid show that the capsid's physical properties are dynamic and local characteristics of the structure, which change with the depth of indentation and depend on the magnitude and geometry of mechanical input. Under large deformations, the Cowpea Chlorotic Mottle Virus capsid transitions to the collapsed state without substantial local structural alterations. The enthalpy change in this deformation state ΔHind = 11.5-12.8 MJ/mol is mostly due to large-amplitude out-of-plane excitations, which contribute to the capsid bending; the entropy change TΔSind = 5.1-5.8 MJ/mol is due to coherent in-plane rearrangements of protein chains, which mediate the capsid stiffening. Direct coupling of these modes defines the extent of (ir)reversibility of capsid indentation dynamics correlated with its (in)elastic mechanical response to the compressive force. This emerging picture illuminates how unique physico-chemical properties of protein nanoshells help define their structure and morphology, and determine their viruses' biological function.

Order statistics inference for describing topological coupling and mechanical symmetry breaking in multidomain proteins

Cooperativity is a hallmark of proteins, many of which show a modular architecture comprising discrete structural domains. Detecting and describing dynamic couplings between structural regions is difficult in view of the many-body nature of protein-protein interactions. By utilizing the GPU-based computational acceleration, we carried out simulations of the protein forced unfolding for the dimer WW - WW of the all-β-sheet WW domains used as a model multidomain protein. We found that while the physically non-interacting identical protein domains (WW) show nearly symmetric mechanical properties at low tension, reflected, e.g., in the similarity of their distributions of unfolding times, these properties become distinctly different when tension is increased. Moreover, the uncorrelated unfolding transitions at a low pulling force become increasingly more correlated (dependent) at higher forces. Hence, the applied force not only breaks "the mechanical symmetry" but also couples the physically non-interacting protein domains forming a multi-domain protein. We call this effect "the topological coupling." We developed a new theory, inspired by order statistics, to characterize protein-protein interactions in multi-domain proteins. The method utilizes the squared-Gaussian model, but it can also be used in conjunction with other parametric models for the distribution of unfolding times. The formalism can be taken to the single-molecule experimental lab to probe mechanical cooperativity and domain communication in multi-domain proteins.

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. The critical extension of the α-helices is 0.25 nm, and the energy difference between the α-helices and β-sheets is 4.9 kcal/mol per helical pitch. The soft α-to-β phase transition in coiled coils might be a universal mechanism underlying mechanical properties of filamentous α-helical proteins.

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. This process is controlled by many cellular factors so that growth and shrinkage periods are correlated with the life cycle of a cell. Resolving the molecular basis for the action of these factors is of paramount importance for understanding the diverse functions of MTs. We employed a multiscale modeling approach to study the force-induced MT depolymerization by analyzing the mechanical response of a MT protofilament to external forces. We carried out self-organized polymer (SOP) model based simulations accelerated on Graphics Processing Units (GPUs). This approach enabled us to follow the mechanical behavior of the molecule on experimental time scales using experimental force loads. We resolved the structural details and determined the physical parameters that characterize the stretching and bending modes of motion of a MT protofilament. The central result is that the severing action of proteins, such as katanin and kinesin, can be understood in terms of their mechanical coupling to a protofilament. For example, the extraction of tubulin dimers from MT caps by katanin can be achieved by pushing the protofilament toward the axis of the MT cylinder, while the removal of large protofilaments curved into ″ram's horn″ structures by kinesin is the result of the outward bending of the protofilament. We showed that, at the molecular level, these types of deformations are due to the anisotropic, but homogeneous, micromechanical properties of MT protofilaments.

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. We have found that mechanical unraveling of fibrin(ogen) is determined by the combined molecular transitions that couple stepwise unfolding of the γ chain nodules and reversible extension-contraction of the α-helical coiled-coil connectors. These findings provide important characteristics of the fibrin(ogen) nanomechanics necessary to understand the molecular origins of fibrin viscoelasticity at the fiber and whole clot levels.