Multiscale impact of nucleotides and cations on the conformational equilibrium, elasticity and rheology of actin filaments and crosslinked networks

Cholate Conjugated Cationic Polymers for Regulation of Actin Dynamics

Cytoskeletal movement is a compulsory necessity for proper cell functioning and is largely controlled by actin filament dynamics. The actin dynamics can be finetuned by various natural and artificial materials...

Biomechanics of cells and subcellular components: A comprehensive review of computational models and applications

Cells are a fundamental structural, functional and biological unit for all living organisms. Up till now, considerable efforts have been made to study the responses of single cells and subcellular components to an external load, and understand the biophysics underlying cell rheology, mechanotransduction and cell functions using experimental and in silico approaches. In the last decade, computational simulation has become increasingly attractive due to its critical role in interpreting experimental data, analysing complex cellular/subcellular structures, facilitating diagnostic designs and therapeutic techniques, and developing biomimetic materials. Despite the significant progress, developing comprehensive and accurate models of living cells remains a grand challenge in the 21st century. To understand current state of the art, this review summarises and classifies the vast array of computational biomechanical models for cells. The article covers the cellular components at multi-spatial levels, that is, protein polymers, subcellular components, whole cells and the systems with scale beyond a cell. In addition to the comprehensive review of the topic, this article also provides new insights into the future prospects of developing integrated, active and high-fidelity cell models that are multiscale, multi-physics and multi-disciplinary in nature. This review will be beneficial for the researchers in modelling the biomechanics of subcellular components, cells and multiple cell systems and understanding the cell functions and biological processes from the perspective of cell mechanics.

Regulation of Actin Bundle Mechanics and Structure by Intracellular Environmental Factors

The mechanical and structural properties of actin cytoskeleton drive various cellular processes, including structural support of the plasma membrane and cellular motility. Actin monomers assemble into double-stranded helical filaments as well as higher-ordered structures such as bundles and networks. Cells incorporate macromolecular crowding, cation interactions, and actin-crosslinking proteins to regulate the organization of actin bundles. Although the roles of each of these factors in actin bundling have been well-known individually, how combined factors contribute to actin bundle assembly, organization, and mechanics is not fully understood. Here, we describe recent studies that have investigated the mechanisms of how intracellular environmental factors influence actin bundling. This review highlights the effects of macromolecular crowding, cation interactions, and actin-crosslinking proteins on actin bundle organization, structure, and mechanics. Understanding these mechanisms is important in determining in vivo actin biophysics and providing insights into cell physiology.

Dissociation Mechanisms of G-actin Subunits Govern Deformation Response of Actin Filament

Actin molecules are essential structural components of the cellular cytoskeleton. Here, we report a comprehensive analysis of F-actin's deformation behavior and highlight underlying mechanisms using steered molecular dynamics simulations (SMD). The investigation of F-actin was done under tension, compression, bending, and torsion. We report that the dissociation pattern of conformational locks at intrastrand and interstrand G-actin interfaces regulates the deformation response of F-actin. The conformational locks at the G-actin interfaces are portrayed by a spheroidal joint, interlocking serrated plates' analogy. Further, the SMD simulation approach was utilized to evaluate Young's modulus, flexural rigidity, persistent length, and torsional rigidity of F-actin, and the values obtained were found to be consistent with available experimental data. The evaluation of the mechanical properties of actin and the insight into the fundamental mechanisms contributing to its resilience described here are necessary for developing accurate models of eukaryotic cells and for assessing cellular viability and mobility.

Effects of Nucleotide and End-Dependent Actin Conformations on Polymerization

The regulation of actin is key for controlled cellular function. Filaments are regulated by actin-binding proteins, but the nucleotide state of actin is also an important factor. From extended molecular dynamics simulations, we find that both nucleotide states of the actin monomer have significantly less twist than their crystal structures and that the ATP monomer is flatter than the ADP form. We also find that the filament's pointed end is flatter than the remainder of the filament and has a conformation distinct from G-actin, meaning that incoming monomers would need to undergo isomerization that would weaken the affinity and slow polymerization. Conversely, the barbed end of the filament takes on a conformation nearly identical to the ATP monomer, enhancing ATP G-actin's ability to polymerize as compared with ADP G-actin. The thermodynamic penalty imposed by differences in isomerization for the ATP and ADP growth at the barbed end exactly matches experimental results.

When Stiffness Matters: Mechanosensing in Heart Development and Disease

During embryonic morphogenesis, the heart undergoes a complex series of cellular phenotypic maturations (e.g., transition of myocytes from proliferative to quiescent or maturation of the contractile apparatus), and this involves stiffening of the extracellular matrix (ECM) acting in concert with morphogenetic signals. The maladaptive remodeling of the myocardium, one of the processes involved in determination of heart failure, also involves mechanical cues, with a progressive stiffening of the tissue that produces cellular mechanical damage, inflammation, and ultimately myocardial fibrosis. The assessment of the biomechanical dependence of the molecular machinery (in myocardial and non-myocardial cells) is therefore essential to contextualize the maturation of the cardiac tissue at early stages and understand its pathologic evolution in aging. Because systems to perform multiscale modeling of cellular and tissue mechanics have been developed, it appears particularly novel to design integrated mechano-molecular models of heart development and disease to be tested in ex vivo reconstituted cells/tissue-mimicking conditions. In the present contribution, we will discuss the latest implication of mechanosensing in heart development and pathology, describe the most recent models of cell/tissue mechanics, and delineate novel strategies to target the consequences of heart failure with personalized approaches based on tissue engineering and induced pluripotent stem cell (iPSC) technologies.

Charge-dependent interactions of monomeric and filamentous actin with lipid bilayers

The cytoskeletal protein actin polymerizes into filaments that are essential for the mechanical stability of mammalian cells. In vitro experiments showed that direct interactions between actin filaments and lipid bilayers are possible and that the net charge of the bilayer as well as the presence of divalent ions in the buffer play an important role. In vivo, colocalization of actin filaments and divalent ions are suppressed, and cells rely on linker proteins to connect the plasma membrane to the actin network. Little is known, however, about why this is the case and what microscopic interactions are important. A deeper understanding is highly beneficial, first, to obtain understanding in the biological design of cells and, second, as a possible basis for the building of artificial cortices for the stabilization of synthetic cells. Here, we report the results of coarse-grained molecular dynamics simulations of monomeric and filamentous actin in the vicinity of differently charged lipid bilayers. We observe that charges on the lipid head groups strongly determine the ability of actin to adsorb to the bilayer. The inclusion of divalent ions leads to a reversal of the binding affinity. Our in silico results are validated experimentally by reconstitution assays with actin on lipid bilayer membranes and provide a molecular-level understanding of the actin-membrane interaction.

The Impact of Natural Compounds on S-Shaped Aβ42 Fibril: From Molecular Docking to Biophysical Characterization

The pursuit for effective strategies inhibiting the amyloidogenic process in neurodegenerative disorders, such as Alzheimer’s disease (AD), remains one of the main unsolved issues, and only a few drugs have demonstrated to delay the degeneration of the cognitive system. Moreover, most therapies induce severe side effects and are not effective at all stages of the illness. The need to find novel and reliable drugs appears therefore of primary importance. In this context, natural compounds have shown interesting beneficial effects on the onset and progression of neurodegenerative diseases, exhibiting a great inhibitory activity on the formation of amyloid aggregates and proving to be effective in many preclinical and clinical studies. However, their inhibitory mechanism is still unclear. In this work, ensemble docking and molecular dynamics simulations on S-shaped Aβ₄₂ fibrils have been carried out to evaluate the influence of several natural compounds on amyloid conformational behaviour. A deep understanding of the interaction mechanisms between natural compounds and Aβ aggregates may play a key role to pave the way for design, discovery and optimization strategies toward an efficient destabilization of toxic amyloid assemblies.

Molecular and Coarse-Grained Modeling to Characterize and Optimize Dendrimer-Based Nanocarriers for Short Interfering RNA Delivery

Dendrimer nanocarriers are unique hyper-branched polymers with biomolecule-like properties, representing a promising prospect as a nucleic acid delivery system. The design of effective dendrimer-based gene carriers requires considering several parameters, such as carrier morphology, size, molecular weight, surface chemistry, and flexibility/rigidity. In detail, the rational design of the dendrimer surface chemistry has been ascertained to play a crucial role on the efficiency of interaction with nucleic acids. Within this framework, advances in the field of organic chemistry have allowed us to design dendrimers with even small difference in the chemical structure of their surface terminals. In this study, we have selected two different cationic phosphorus dendrimers of generation 3 functionalized, respectively, with pyrrolidinium (DP) and morpholinium (DM) surface groups, which have demonstrated promising potential for short interfering RNA (siRNA) delivery. Despite DP and DM differing only for one atom in their chemical structure, in vitro and in vivo experiments have highlighted several differences between them in terms of siRNA complexation properties. In this context, we have employed coarse-grained molecular dynamics simulation techniques to shed light on the supramolecular characteristics of dendrimer-siRNA complexation, the so-called dendriplex formations. Our data provide important information on self-assembly dynamics driven by surface chemistry and competition mechanisms.

Aminoacid substitutions in the glycine zipper affect the conformational stability of amyloid beta fibrils

The aggregation of amyloid-beta peptides is associated with the pathogenesis of Alzheimer's disease. The hydrophobic core of the amyloid beta sequence contains a GxxxG repeated motif, called glycine zipper, which involves crucial residues for assuring stability and promoting the process of fibril formation. Mutations in this motif lead to a completely different oligomerization pathway and rate of fibril formation. In this work, we have tested G33L and G37L residue substitutions by molecular dynamics simulations. We found that both protein mutations may lead to remarkable changes in the fibril conformational stability. Results suggest the disruption of the glycine zipper as a possible strategy to reduce the aggregation propensity of amyloid beta peptides. On the basis of our data, further investigations may consider this key region as a binding site to design/discover novel effective inhibitors.Communicated by Ramaswamy H. Sarma.

The Role of Structural Polymorphism in Driving the Mechanical Performance of the Alzheimer's Beta Amyloid Fibrils

Alzheimer's Disease (AD) is related with the abnormal aggregation of amyloid β-peptides Aβ₁₋₄₀ and Aβ₁₋₄₂, the latter having a polymorphic character which gives rise to U- or S-shaped fibrils. Elucidating the role played by the nanoscale-material architecture on the amyloid fibril stability is a crucial breakthrough to better understand the pathological nature of amyloid structures and to support the rational design of bio-inspired materials. The computational study here presented highlights the superior mechanical behavior of the S-architecture, characterized by a Young's modulus markedly higher than the U-shaped architecture. The S-architecture showed a higher mechanical resistance to the enforced deformation along the fibril axis, consequence of a better interchain hydrogen bonds' distribution. In conclusion, this study, focusing the attention on the pivotal multiscale relationship between molecular phenomena and material properties, suggests the S-shaped Aβ₁₋₄₂ species as a target of election in computational screen/design/optimization of effective aggregation modulators.

Coarse-Grained Simulation of Full-Length Integrin Activation

Integrin conformational dynamics are critical to their receptor and signaling functions in many cellular processes, including spreading, adhesion, and migration. However, assessing integrin conformations is both experimentally and computationally challenging because of limitations in resolution and dynamic sampling. Thus, structural changes that underlie transitions between conformations are largely unknown. Here, focusing on integrin αvβ₃, we developed a modified form of the coarse-grained heterogeneous elastic network model (hENM), which allows sampling conformations at the onset of activation by formally separating local fluctuations from global motions. Both local fluctuations and global motions are extracted from all-atom molecular dynamics simulations of the full-length αvβ₃ bent integrin conformer, but whereas the former are incorporated in the hENM as effective harmonic interactions between groups of residues, the latter emerge by systematically identifying and treating weak interactions between long-distance domains with flexible and anharmonic connections. The new hENM model allows integrins and single-point mutant integrins to explore various conformational states, including the initiation of separation between α- and β-subunit cytoplasmic regions, headpiece extension, and legs opening.

Multiscale simulation of actin filaments and actin-associated proteins

Actin is an important cytoskeletal protein that serves as a building block to form filament networks that span across the cell. These networks are orchestrated by a myriad of other cytoskeletal entities including the unbranched filament-forming protein formin and branched network-forming protein complex Arp2/3. Computational models have been able to provide insights into many important structural transitions that are involved in forming these networks, and into the nature of interactions essential for actin filament formation and for regulating the behavior of actin-associated proteins. In this review, we summarize a subset of such models that focus on the atomistic features and those that can integrate atomistic features into a larger picture in a multiscale fashion.

Insights into the Cooperative Nature of ATP Hydrolysis in Actin Filaments

Actin filaments continually assemble and disassemble within a cell. Assembled filaments "age" as a bound nucleotide ATP within each actin subunit quickly hydrolyzes followed by a slower release of the phosphate Pi, leaving behind a bound ADP. This subtle change in nucleotide state of actin subunits affects filament rigidity as well as its interactions with binding partners. We present here a systematic multiscale ultra-coarse-graining approach that provides a computationally efficient way to simulate a long actin filament undergoing ATP hydrolysis and phosphate-release reactions while systematically taking into account available atomistic details. The slower conformational changes and their dependence on the chemical reactions are simulated with the ultra-coarse-graining model by assigning internal states to the coarse-grained sites. Each state is represented by a unique potential surface of a local heterogeneous elastic network. Internal states undergo stochastic transitions that are coupled to conformations of the underlying molecular system. The model reproduces mechanical properties of the filament and allows us to study whether conformational fluctuations in actin subunits produce cooperative filament aging. We find that the nucleotide states of neighboring subunits modulate the reaction kinetics, implying cooperativity in ATP hydrolysis and Pi release. We further systematically coarse grain the system into a Markov state model that incorporates assembly and disassembly, facilitating a direct comparison with previously published models. We find that cooperativity in ATP hydrolysis and Pi release significantly affects the filament growth dynamics only near the critical G-actin concentration, whereas far from it, both cooperative and random mechanisms show similar growth dynamics. In contrast, filament composition in terms of the bound nucleotide distribution varies significantly at all monomer concentrations studied. These results provide new insights, to our knowledge, into the cooperative nature of ATP hydrolysis and Pi release and the implications it has for actin filament properties, providing novel predictions for future experimental studies.

Crosslinking of fibrous hydrogels

In contrast to most synthetic hydrogels, biological gels are made of fibrous networks. This architecture gives rise to unique properties, like low concentration, high porosity gels with a high mechanical responsiveness as a result of strain-stiffening. Here, we used a synthetic polymer model system, based on polyisocyanides, that we crosslinked selectively inside the bundles. This approach allows us to lock in the fibrous network present at the crosslinking conditions. At minimum crosslink densities, we are able to freeze in the architecture, as well as the associated mechanical properties. Rheology and X-ray scattering experiments show that we able to accurately tailor network mechanics, not by changing the gel composition or architecture, but rather by tuning its (thermal) history. Selective crosslinking is a crucial step in making biomimetic networks with a controlled architecture.

Unlike synthetic hydrogels, biological gels are made of fibrous networks which give rise to unique properties, such as high porosity and mechanical responsiveness. Here the authors use polyisocyanide-based gels and selectively crosslink inside the bundles to lock the fibrous network and thus control the architecture and the mechanics.

Mechanical differences between ATP and ADP actin states: A molecular dynamics study

This paper aims to give a comprehensive atomistic modeling of the nanomechanical behavior of actin monomer. Actin is a ubiquitous and essential component of cytoskeleton which forms many different cellular structures. Despite for several years great effort has been devoted to the investigation of mechanical properties of the actin filament, studies on the nanomechanical behavior of actin monomer are still lacking. These scales are, however, important for a complete understanding of the role of actin as an important component in the cytoskeleton structure. Based on the accuracy of atomistic modeling methods such as molecular dynamics simulations, steered molecular dynamics method is performed to assess tension of monomeric G-actin molecule under different types of mechanical loading including axial and lateral. As a result, stress-strain curves are obtained in aqueous solution, with either ATP or ADP bound in the nucleotide binding pocket. The obtained results yield evaluation of the tensile stiffness of a single actin monomer in lateral and normal direction. In order to compare the behavior of ATP and ADP G-actins, the number of hydrogen bonds and nonbonded interactions between the nucleotide and the protein are analyzed. Moreover, The effect of virtual spring of steered molecular dynamics on the mechanical behavior of actin monomer is investigated. The results reveal increasing the virtual spring constant leads to convergence of the stiffness. Moreover, in this paper, a generalized model is proposed to extend the obtained results for the monomeric G-actin scale to the actin filament. Our modeling estimated a persistence length of actin filament 15.41 µm, close to experimental measurements. Moreover, In this paper, the breaking force actin-actin bond is evaluated using steered molecular dynamics simulation. By applying a tensile force, actin-actin bond ruptured at 4197.5 pN.

Nanomechanics of actin filament: A molecular dynamics simulation

Actin is known as the most abundant essentially protein in eukaryotic cells. Actin plays a crucial role in many cellular processes involving mechanical forces such as cell motility, adhesion, muscle contraction, and intracellular transport. However, little is known about the mechanical properties of this protein when subjected to mechanical forces in cellular processes. In this article, a series of large-scale molecular dynamics simulations are carried out to elucidate nanomechanical behavior such as elastic and viscoelastic properties of a single actin filament. Here, we used two individual methods namely, all-atoms and coarse-grained molecular dynamics, to evaluate elastic properties of a single actin filament. In the other word, based on Brownian motions of the filament and using the principle of the equipartition theorem, in aqueous solution, tensile stiffness, torsional rigidity, and bending rigidity of the single actin filament are studied. The results revealed that increasing the sampling window time leads to convergence of obtained mechanical properties to the experimental values. Moreover, in order to investigate viscoelastic properties of a single actin filament, constant force steered molecular dynamics method is used to apply different external tensile loads and perform five individual creep tests on the molecule. The strain-time response of the filament for each creep test is obtained. Based on the Kelvin-Voigt model, the results reveal that a single actin filament shows a nonlinear viscoelastic behavior, with a Young's modulus of 2.85 GPa, a viscosity of 4.06 GPa.ns, and a relaxation time in the range of 1.42 ns which were measured here for the first time at the single filament level. The findings of this article suggest that molecular dynamics simulations could also be a useful tool for investigating the mechanical behavior of bio-nanomaterials.

Probing eukaryotic cell mechanics via mesoscopic simulations

Cell mechanics has proven to be important in many biological processes. Although there is a number of experimental techniques which allow us to study mechanical properties of cell, there is still a lack of understanding of the role each sub-cellular component plays during cell deformations. We present a new mesoscopic particle-based eukaryotic cell model which explicitly describes cell membrane, nucleus and cytoskeleton. We employ Dissipative Particle Dynamics (DPD) method that provides us with the unified framework for modeling of a cell and its interactions in the flow. Data from micropipette aspiration experiments were used to define model parameters. The model was validated using data from microfluidic experiments. The validated model was then applied to study the impact of the sub-cellular components on the cell viscoelastic response in micropipette aspiration and microfluidic experiments.

Thermodynamic and kinetic stability of the Josephin Domain closed arrangement: evidences from replica exchange molecular dynamics

BACKGROUND

Molecular phenomena driving pathological aggregation in neurodegenerative diseases are not completely understood yet. Peculiar is the case of Spinocerebellar Ataxia 3 (SCA3) where the conformational properties of the AT-3 N-terminal region, also called Josephin Domain (JD), play a key role in the first step of aggregation, having the JD an amyloidogenic propensity itself. For this reason, unraveling the intimate relationship between JD structural features and aggregation tendency may lead to a step forward in understanding the pathology and rationally design a cure. In this connection, computational modeling has demonstrated to be helpful in exploring the protein molecular dynamics and mechanism of action.

RESULTS

Conformational dynamics of the JD is here finely investigated by replica exchange molecular dynamics simulations able to sample the microsecond time scale and to provide both a thermodynamic and kinetic description of the protein conformational changes. Accessible structural conformations of the JD have been identified in: open, intermediate and closed like arrangement. Data indicated the closed JD arrangement as the most likely protein arrangement. The protein transition from closed toward intermediate/open states was characterized by a rate constant higher than 700 ns. This result also explains the inability of classical molecular dynamics to explore transitions from closed to open JD configuration on a time scale of hundreds of nanoseconds.

CONCLUSION

This work provides the first kinetic estimation of the JD transition pathway from open-like to closed-like arrangement and vice-versa, indicating the closed-like arrangement as the most likely configuration for a JD in water environment. More widely, the importance of our results is also underscored considering that the ability to provide a kinetic description of the protein conformational changes is a scientific challenge for both experimental and theoretical approaches to date.

REVIEWERS

This article was reviewed by Oliviero Carugo, Bojan Zagrovic.

Surface coupling strength of gold nanoparticles affects cytotoxicity towards neurons

Weakly bound gold nanoparticles reveal awful toxicity towards neurons.

Morphological Transformation and Force Generation of Active Cytoskeletal Networks

Cells assemble numerous types of actomyosin bundles that generate contractile forces for biological processes, such as cytokinesis and cell migration. One example of contractile bundles is a transverse arc that forms via actomyosin-driven condensation of actin filaments in the lamellipodia of migrating cells and exerts significant forces on the surrounding environments. Structural reorganization of a network into a bundle facilitated by actomyosin contractility is a physiologically relevant and biophysically interesting process. Nevertheless, it remains elusive how actin filaments are reoriented, buckled, and bundled as well as undergo tension buildup during the structural reorganization. In this study, using an agent-based computational model, we demonstrated how the interplay between the density of myosin motors and cross-linking proteins and the rigidity, initial orientation, and turnover of actin filaments regulates the morphological transformation of a cross-linked actomyosin network into a bundle and the buildup of tension occurring during the transformation.

Cancer cells become less deformable and more invasive with activation of β-adrenergic signaling

Invasion by cancer cells is a crucial step in metastasis. An oversimplified view in the literature is that cancer cells become more deformable as they become more invasive. β-adrenergic receptor (βAR) signaling drives invasion and metastasis, but the effects on cell deformability are not known. Here, we show that activation of β-adrenergic signaling by βAR agonists reduces the deformability of highly metastatic human breast cancer cells, and that these stiffer cells are more invasive in vitro We find that βAR activation also reduces the deformability of ovarian, prostate, melanoma and leukemia cells. Mechanistically, we show that βAR-mediated cell stiffening depends on the actin cytoskeleton and myosin II activity. These changes in cell deformability can be prevented by pharmacological β-blockade or genetic knockout of the β₂-adrenergic receptor. Our results identify a β₂-adrenergic-Ca²⁺-actin axis as a new regulator of cell deformability, and suggest that the relationship between cell mechanical properties and invasion might be dependent on context.

Conformational fluctuations of the AXH monomer of Ataxin-1

In this paper, we report the results of molecular dynamics simulations of AXH monomer of Ataxin-1. The AXH domain plays a crucial role in Ataxin-1 aggregation, which accompanies the initiation and progression of Spinocerebellar ataxia type 1. Our simulations involving both classical and replica exchange molecular dynamics, followed by principal component analysis of the trajectories obtained, reveal substantial conformational fluctuations of the protein structure, especially in the N-terminal region. We show that these fluctuations can be generated by thermal noise since the free energy barriers between conformations are small enough for thermally stimulated transitions. In agreement with the previous experimental findings, our results can be considered as a basis for a future design of ataxin aggregation inhibitors that will require several key conformations identified in the present study as molecular targets for ligand binding.

Multiscale mechanobiology: computational models for integrating molecules to multicellular systems

Mechanical signals exist throughout the biological landscape. Across all scales, these signals, in the form of force, stiffness, and deformations, are generated and processed, resulting in an active mechanobiological circuit that controls many fundamental aspects of life, from protein unfolding and cytoskeletal remodeling to collective cell motions. The multiple scales and complex feedback involved present a challenge for fully understanding the nature of this circuit, particularly in development and disease in which it has been implicated. Computational models that accurately predict and are based on experimental data enable a means to integrate basic principles and explore fine details of mechanosensing and mechanotransduction in and across all levels of biological systems. Here we review recent advances in these models along with supporting and emerging experimental findings.