Olga Kononova

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.

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

Assembly and Mechanical Properties of the Cargo-Free and Cargo-Loaded Bacterial Nanocompartment Encapsulin

Prokaryotes mostly lack membranous compartments that are typical of eukaryotic cells, but instead, they have various protein-based organelles. These include bacterial microcompartments like the carboxysome and the virus-like nanocompartment encapsulin. Encapsulins have an adaptable mechanism for enzyme packaging, which makes it an attractive platform to carry a foreign protein cargo. Here we investigate the assembly pathways and mechanical properties of the cargo-free and cargo-loaded nanocompartments, using a combination of native mass spectrometry, atomic force microscopy and multiscale computational molecular modeling. We show that encapsulin dimers assemble into rigid single-enzyme bacterial containers. Moreover, we demonstrate that cargo encapsulation has a mechanical impact on the shell. The structural similarity of encapsulins to virus capsids is reflected in their mechanical properties. With these robust mechanical properties encapsulins provide a suitable platform for the development of nanotechnological applications.

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

Background: Knob-hole interactions underlie formation and properties of fibrin polymer, the scaffold of blood clots and thrombi. Results: The structural mechanisms, dissociation kinetics, and thermodynamic parameters of the A:a and B:b knob-hole interactions have been determined. Conclusion: The knob-hole bonds are inherently variable and sensitive to pH and temperature. Significance: The emerging molecular picture offers mechanistic insights into fibrin polymerization. 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). There were similar structural changes in holes ‘a’ and ‘b’ during forced dissociation of the knob-hole complexes: elongation of loop I, stretching of the interior region, and translocation of the moveable flap. The disruption of the knob-hole interactions was not an “all-or-none” transition as it occurred through distinct two-step or single step pathways with or without intermediate states. The knob-hole bonds were stronger, tighter, and more brittle at pH 7 than at pH 5. The B:b knob-hole bonds were weaker, looser, and more compliant than the A:a knob-hole bonds at pH 7 but stronger, tighter, and less compliant at pH 5. Surprisingly, the knob-hole bonds were stronger, not weaker, at elevated temperature (T = 37 °C) compared with T = 25 °C due to the helix-to-coil transition in loop I that helps stabilize the bonds. These results provide detailed qualitative and quantitative characteristics underlying the most significant non-covalent interactions involved in fibrin polymerization.

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.

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 I2 U, with two intermediate states I1 and I2. BoNT/A loses almost all its secondary structure in 3.75 M urea (I1), yet it displays a native-like secondary structure in 5 M urea (I2). This agrees with the results of theoretical modeling, which helped to determine the molecular basis of unique behavior of BoNT/A in solution. Except for I2, all the states revert back to full enzymatic activity for SNAP-25 including the unfolded state U stable in 7 M urea. Our results stress the importance of structural flexibility in the toxin’s mechanism of survival and action, an unmatched evolutionary trait from billion-year-old bacteria, which also correlates with the long-lasting enzymatic activity of BoNT inside neuronal cells. BoNT/A provides a rich model to explore protein folding in relation to functional activity.

Mechanistic Basis for the Binding of RGD- and AGDV-Peptides to the Platelet Integrin αIIbβ3

Binding of soluble fibrinogen to the activated conformation of the integrin αIIbβ3 is required for platelet aggregation and is mediated exclusively by the C-terminal AGDV-containing dodecapeptide (γC-12) sequence of the fibrinogen γ chain. However, peptides containing the Arg-Gly-Asp (RGD) sequences located in two places in the fibrinogen Aα chain inhibit soluble fibrinogen binding to αIIbβ3 and make substantial contributions to αIIbβ3 binding when fibrinogen is immobilized and when it is converted to fibrin. Here, we employed optical trap-based nanomechanical measurements and computational molecular modeling to determine the kinetics, energetics, and structural details of cyclic RGDFK (cRGDFK) and γC-12 binding to αIIbβ3. Docking analysis revealed that NMR-determined solution structures of cRGDFK and γC-12 bind to both the open and closed αIIbβ3 conformers at the interface between the αIIb β-propeller domain and the β3 βI domain. The nanomechanical measurements revealed that cRGDFK binds to αIIbβ3 at least as tightly as γC-12. A subsequent analysis of molecular force profiles and the number of peptide-αIIbβ3 binding contacts revealed that both peptides form stable bimolecular complexes with αIIbβ3 that dissociate in the 60-120 pN range. The Gibbs free energy profiles of the αIIbβ3-peptide complexes revealed that the overall stability of the αIIbβ3-cRGDFK complex was comparable with that of the αIIbβ3-γC-12 complex. Thus, these results provide a mechanistic explanation for previous observations that RGD- and AGDV-containing peptides are both potent inhibitors of the αIIbβ3-fibrinogen interactions and are consistent with the observation that RGD motifs, in addition to AGDV, support interaction of αIIbβ3 with immobilized fibrinogen and fibrin.

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 models capability for describing protein systems of varying complexity (102–105 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.

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

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

Significance The fibrin network is the scaffold of blood clots that prevent bleeding. Fibrin formation is driven by noncovalent “knob-hole” bonds that determine the mechanical properties of clots. Counterintuitively, the strength of knob-hole bonds first increased (catch bond) and then decreased (slip bond) with tensile force, thereby displaying a rare dual catch–slip behavior. We resolved the structural mechanism of tension-activated strengthening of the knob-hole bond. Our theory quantitatively accounts for the biphasic catch–slip force response, which comprises an important aspect of fibrin structural bionanomechanics. The shear-enhanced strengthening of noncovalent knob-hole bonds might favor fibrin polymerization under blood flow and might prevent breakup of nascent blood clots in the circulation. 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.