Skip to main content

ORIGINAL RESEARCH article

Front. Mol. Biosci., 15 February 2021
Sec. Biological Modeling and Simulation
This article is part of the Research Topic Combining Simulations, Theory, and Experiments into Multiscale Models of Biological Events View all 23 articles

Coarse-Grained Simulation of Mechanical Properties of Single Microtubules With Micrometer Length

Jinyin ZhaJinyin Zha1Yuwei ZhangYuwei Zhang1Kelin Xia,Kelin Xia2,3Frauke Grter,,Frauke Gräter4,5,6Fei Xia,
Fei Xia1,7*
  • 1School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, China
  • 2Division of Mathematical Sciences, School of Physical and Mathematical Sciences, Nanyang Technological University, Singapore, Singapore
  • 3School of Biological Sciences, Nanyang Technological University, Singapore, Singapore
  • 4Interdisciplinary Centre for Scientific Computing (IWR), Heidelberg University, Heidelberg, Germany
  • 5Heidelberg Institute for Theoretical Studies (HITS), Schloβ-Wolfsbrunnenweg 35, Heidelberg, Germany
  • 6Max Planck School Matter to Life, Jahnstraβe 29, Heidelberg, Germany
  • 7Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, China

Microtubules are one of the most important components in the cytoskeleton and play a vital role in maintaining the shape and function of cells. Because single microtubules are some micrometers long, it is difficult to simulate such a large system using an all-atom model. In this work, we use the newly developed convolutional and K-means coarse-graining (CK-CG) method to establish an ultra-coarse-grained (UCG) model of a single microtubule, on the basis of the low electron microscopy density data of microtubules. We discuss the rationale of the micro-coarse-grained microtubule models of different resolutions and explore microtubule models up to 12-micron length. We use the devised microtubule model to quantify mechanical properties of microtubules of different lengths. Our model allows mesoscopic simulations of micrometer-level biomaterials and can be further used to study important biological processes related to microtubule function.

Introduction

Microtubule (Fratzl and Weinkamer, 2007; Conde and Caceres, 2009) is one of the primary components of the cytoskeleton which is responsible for the morphology of eukaryotic cells. The building blocks of microtubules are α- and β-tubulins which assemble into a hollow cylinder structure with grooves. The length of a single microtubule can reach 10 µm and they are associated with each other in a hierarchical organization such as in the axons (Lazarus et al., 2015; Soheilypour et al., 2015). A microtubule is generally composed of 13 protofilaments (Kollman et al., 2010) lying side by side along the longitudinal direction, where the position of each protofilament shifts slightly to each other. The mechanical properties (Franz et al., 2020) of microtubules play important roles in sustaining the cell shape, executing the mitotic division and related physiological functions. For instance, they serve as tracks for the motor proteins kinesin (Gittes et al., 1996) to transport the substances required for cells.

In order to uncover the microscopic mechanisms of microtubules for transporting substances in cells, much effort has been put into measuring the biomechanical properties (Dye et al., 1993; Gittes et al., 1993; Venier et al., 1994; Kurachi et al., 1995; Felgner et al., 1996; Vinckier, 1996; Wagner et al., 1999; Kis et al., 2002; Le Goff et al., 2002; VanBuren et al., 2002; de Pablo et al., 2003; Janson and Dogterom, 2004; Pampaloni et al., 2006; Schaap et al., 2006; Hawkins et al., 2013; Cross, 2019; Szatkowski et al., 2019) of microtubules associated with their function, such as Young’s modulus E, the flexural rigidity κ, and the persistence length lp. Previous studies have exploited a variety of advanced techniques including thermal fluctuation (Gittes et al., 1993), hydrodynamic flow (Venier et al., 1994), atomic force microscopy (Vinckier, 1996), and optical tweezers (Felgner et al., 1996) to measure these mechanical quantities. So far, reported Young’s modulus E, the flexural rigidity κ, and persistence length lp for microtubules with different lengths span a large range of 0.31–7.0 GPa (Kurachi et al., 1995; Vinckier, 1996; Kis et al., 2002; de Pablo et al., 2003; Schaap et al., 2006; Kis et al., 2008), 0.37–4.0 × 10−23 Nm2 (Gittes et al., 1993; Venier et al., 1994; Kurachi et al., 1995; Felgner et al., 1996), and 2.0–6.3 mm (Gittes et al., 1993; Venier et al., 1994; Janson and Dogterom, 2004; Pampaloni et al., 2006), respectively. Kis et al. (Kis et al., 2002) first reported that the observed phenomenon of the measured persistence length was length-dependent, caused by the anisotropic association of heterodimeric α- and β-tubulins. This finding was further demonstrated by Pampaloni et al. (Pampaloni et al., 2006) using the single-particle tracking method and they observed that the persistence length lp increased for short lengths varying from 2.6 to 47.5 μm. It approached a constant value of lp = 6.3 mm for lengths beyond 21 μm. All measured results unambiguously imply that microtubules do not follow the behavior of an ideal elastic rod and fail to be described with the worm-like chain model (Kroy and Frey, 1996; Taute et al., 2008) based on the isotropic hypothesis.

Theoretical modeling of microtubules (Kasas et al., 2004a; Kasas et al., 2004b; Molodtsov et al., 2005; VanBuren et al., 2005; Schaap et al., 2006; Deriu et al., 2007; Heuvel et al., 2007; Deriu et al., 2008; Dima and Joshi, 2008; Deriu et al., 2010; Sept and MacKintosh, 2010; Wells and Aksimentiev, 2010; Grafmuller and Voth, 2011; Ji and Feng, 2011a; b; Peter and Mofrad, 2012; Grafmuller et al., 2013; Theisen et al., 2013; Kononova et al., 2014; Lazarus et al., 2015; Soheilypour et al., 2015; Havelka et al., 2017; Stevens, 2017; Szatkowski et al., 2019; Tong and Voth, 2020) faces a huge challenge due to the micrometer length. All-atom molecular dynamics (MD) simulations have been primarily used to investigate the atomistic details of the interactions between the α- and β-tubulins (Deriu et al., 2007; Deriu et al., 2008; Sept and MacKintosh, 2010; Wells and Aksimentiev, 2010; Grafmuller and Voth, 2011; Grafmuller et al., 2013; Havelka et al., 2017) for short microtubules. Recent all-atom simulations (Tong and Voth, 2020) reached a length scale of a 1 μm microtubule including 3 protofilaments. Additionally, coarse-grained (CG) models (Dima and Joshi, 2008; Deriu et al., 2010; Ji and Feng, 2011a; Ji and Feng, 2011b; Theisen et al., 2013; Kononova et al., 2014; Lazarus et al., 2015; Soheilypour et al., 2015; Stevens, 2017; Szatkowski et al., 2019) have also been used to model longer microtubules. For instance, Deriu et al. (Deriu et al., 2010) constructed a 350 nm long microtubule to estimate its mechanical properties based on the elastic network model (Bahar et al., 1997; Atilgan et al., 2001; Xia and Lu, 2012; Xia et al., 2013a). Nevertheless, their model was still far shorter than the microtubules measured in experiments (Pampaloni et al., 2006). In some special CG models, each tubulin monomer of a microtubule was even represented by one CG bead (Ji and Feng, 2011a; Ji and Feng, 2011b; Lazarus et al., 2015) so that the whole microtubule model could reach the micrometer length. Alternatively, other theoretical methods based on continuous mechanics such as the finite element method (Kasas et al., 2004a; Kasas et al., 2004b; Schaap et al., 2006) were also used to evaluate the mechanical properties of macromaterials. However, these methods had an obvious disadvantage that they could not reflect the correct surface contour of microtubules, especially for the functional grooves serving as binding sites to interact with other proteins.

Ultra-coarse-grained (UCG) models (Dama et al., 2013; Davtyan et al., 2014; Dama et al., 2017; Zhang et al., 2017a; Jin et al., 2018; Jin and Voth, 2018) have been rapidly developed in recent years and applied to simulate the functions of biomolecular complexes such as the mechanical properties of F-actin (Katkar et al., 2018), the assembly of HIV capsids (Grime et al., 2016), and two-state conformational changes of proteins (Zhang et al., 2020). The bead resolution of UCG models is lower than the high-resolution models, with one UCG bead usually representing more than a few residues or even secondary structure elements. Previously, Voth and coworkers (Zhang et al., 2008; Zhang et al., 2009) developed the essential-dynamic coarse-graining (ED-CG) method to optimize the UCG representation (Sinitskiy et al., 2012) of proteins and later extended the framework to the density-based ED-CG method (Zhang and Voth, 2010). Xia and coworkers (Li et al., 2016a; Li et al., 2016b; Zhang et al., 2017b) developed a stepwise local iterative optimization-based fluctuation-maximization coarse-graining method to derive the optimal UCG representation (Wu et al., 2020) of large biomolecules. Further, Zhang et al. (Zhang et al., 2019) developed the convolutional and K-means coarse-graining (CK-CG) method to derive the UCG representation of biomolecules directly from the low-resolution cryo-EM density data (Ming et al., 2002; Esquivel-Rodriguez and Kihara, 2013), which allows the construction of models without any atomistic details needed further. Next, the UCG beads could be parameterized using the energy potential of the Multiscale Virtual Particle-based Anisotropic Network Model (MVP-ANM) (Xia et al., 2013b; Opron et al., 2014; Xia, 2017) for further functional studies. In the previous study (Zhang et al., 2019), we have validated the combined CK-CG and MVP-ANM methods in constructing UCG model for a correct estimation of the biomechanical properties of short segments F-actin and collagen fibrils. In this work, we attempt to verify the usefulness of UCG models in simulating the mechanical properties such as the length-dependent persistence length of long microtubules and expect to apply it to the simulation more complicated and hierarchically organized biomaterials (Zapp et al., 2020).

Computational and Simulation Details

Data Preprocessing and CK-CG Method

The density file EMD-4043(Grange et al., 2017) of microtubule used for model construction was downloaded from the cryo-EM data bank online. The original length of microtubule in the file EMD-4043 is only 32 nm, in a resolution of 18.5 Å. In order to construct long microtubules, we replicated the density data shown in Figure 1A along the longitudinal direction and extended the 32 nm length to the different lengths ranging from 1 to 12 μm. Then, we removed the density of solvents from the raw data with a technique described in the previous work (Zhang et al., 2019) and evaluated the normalized density µs (r) according to µs (r) = (µ(r) − µmin)/(µmax − µmin), where µ (r) denotes the cryo-EM density data. After the preprocessing on data, the remaining densities are positive and ready as the input for coarse-graining by using CK-CG method.

FIGURE 1
www.frontiersin.org

FIGURE 1. (A) The data reprocessing includes replicating the density, removing solvent, and normalizing density. (B) The process of coarse-graining cryo-EM data into 12 UCG beads using the CK-CG method. The density map is first divided into 9 blocks. The Hadamard product for each block is obtained by using a Gaussian kernel and the subsequent max-pooling generates the initial positions of 9 UCG beads. The positions of the remaining 3 UCG beads are determined by other maximal densities in the blocks. Then, the 12 positions serve as the centers for clustering to yield the final UCG beads.

Coarse-graining density data is actually to cluster the adjacent discretized density data into a group of single UCG particles. In the conventional K-means method, the positions of center points in initial guess are generated randomly, which renders the convergence quite slow. In the CK-CG method, we utilized an idea similar to the Convolutional Neutral Network model (Li et al., 2019) to generate the initial guess for clustering, including the steps of the domain partition, density convolution, and max-pooling, displayed in Figure 1B. In more detail, the cryo-EM density data is firstly divided into a series of 3D block matrixes with the voxel bijl, where the subscripts i, j, and l denote dimensions. In the next convolution operation, we used a 3D Gaussian kernel matrix with the voxel kijl to execute a Hadamard product based convolution operation and saved the matrix values hijl.

In the subsequent max-pooling, the position corresponding to the largest density value in each block was chosen as the initial position of UCG bead for clustering. It has to be mentioned that if the expected number of UCG beads N is larger than the number of divided blocks M, the excess number of (N-M) beads might be automatically assigned into the blocks with largest values. After generating the initial guess, the clustering using K-means was carried out until the convergence was achieved. More details about the procedure of CK-CG workflow are described in the previous work (Zhang et al., 2019).

Multiscale Virtual Particle-Based Anisotropic Network Model

UCG beads derived from CK-CG could be connected with harmonic springs and parameterized to be the MVP-ANM model developed by Xia et al. (Xia et al., 2013b; Opron et al., 2014; Xia, 2017). The energy function EMVP−ANM is a sum of harmonic interactions as shown in Eq. 1, where γij is a heterogeneous force constant, dij and dij0 denote the instantaneous and equilibrium distances of pairwise UCG beads i and j, and H(dcdij0) is a Heaviside function depending on a cutoff distance dc.

EMVPANM=ij>iγij(dijdij0)2H(dcdij0)(1)

The parameter γij contains both the contributions from the density and distance of pairwise UCG beads i and j, with the expression shown in Eq. 2:

γij=A(1+mi)(1+mj)edij2/R2(2)

The parameter A is a scaling factor to control the mechanical properties of whole model and can be determined by fitting the experimental Young’s modulus using force-clamped MD simulation (Isralewitz et al., 2001; Gräter and Grubmuller, 2007) as shown in Figure 2A. The term of (1+mi)(1+mj) is a function representing the mass contributions from the beads i and j. The distance contribution is defined as a Gaussian function in which R is a parameter to control the bond strength of pairwise beads i and j, and it is only taken into account if the equilibrium distance dij0 is less than the chosen cutoff dc. A minimal cutoff can be determined through MD simulations by judging whether the structures of microtubules remain stable all the time. Please see the details for determining a minimal dc described in Parameterization of UCG Model below. Additionally, the parameter γij decays significantly as dij>3R. Thus, we define the magnitude of R to be one-third of the minimal cutoff of dc.

FIGURE 2
www.frontiersin.org

FIGURE 2. (A) An illustration of force-clamped MD simulation for measuring Young’s modulus of microtubule. (B) The estimation of persistence length using a method of 3D linear regression based on thermal fluctuation.

Evaluation of Persistence Length From UCG MD Simulation

We performed the UCG MD simulation of microtubules at 300 K and calculated the persistence length according the definition in Eq. 3:

cosθ=eL/lp(3)

where cosθ means the average value of the cosine of θ (0θπ/2) and θ denotes the intersection angle of two tangent lines at the positions of the head and tail of microtubules. In order to calculate it, we recorded the Cartesian coordinates of all the beads every 12 ps in UCG MD simulation, and obtained a number of frames for further analysis. For the structure in each frame, we used a 3D liner regression method in Figure 2B to solve Eq. 3 and derived the persistence length from the calculated average value of cosθ. All the UCG MD simulations mentioned above were carried out with using the software LAMMPS (Plimpton, 1995) in the NVT ensembles at 300K. The Langevin dynamics simulation (Grønbech-Jensen, 2019) was performed under the nonperiodic condition with the damping factor being 10.0 ps. The time step of simulation is chosen to be 10.0 fs for the simulation of Young’s modulus. For the estimation of persistence length from a long time UCG simulation, since the lightest bead in the UCG model is around 6.0 kg/mol, nearly 6000 times greater than H atom, we chose a moderate value of 3.0 ps for the integration steps. The UCG MD simulation for a single microtubule with micrometer length could reach the time scale of 60 ms.

It has to be mentioned that the persistence length of microtubules was frequently estimated from the bending measurements in experiments (Kurachi et al., 1995; Felgner et al., 1996; Kis et al., 2008) and simulations (Kononova et al., 2014; Wu and Adnan, 2018). However, our UCG model is not suitable for being used in a direct simulation of bending deformation, since the MVP-ANM potential used in our model is harmonic and invalid in large deformation. A simulation of the complicated process such as the dissociation and association of tubulins from microtubule requires a more complex potential such as a Gō-like potential (Zhang et al., 2017a; Poma et al., 2017; Poma et al., 2018; Zhang et al., 2020).

Results and Discussions

Coarse-Graining Microtubules With Different Bead Resolutions and Lengths

Previously, we have constructed a UCG model for a short 160 nm microtubule (Zhang et al., 2019). In this work, we attempt to coarse-grain the microtubules with the lengths from 1 to 10 μm which are much longer than the previous one by 1 or 2 orders of magnitude. In order to demonstrate the efficiency of CK-CG over the conventional K-means methods, we test its time costs in coarse-graining the microtubules into different numbers of beads and lengths.

The first test given is the time cost for coarse-graining a 1 μm microtubule UCG model into different numbers of UCG beads. The test was done by running C++ language codes with one CPU on a CentOS platform. As shown in Figure 3A, the derived numbers of UCG beads for a fixed-length microtubule range from 7609 to 33480 beads per micrometer, corresponding to the averaged masses of UCG beads from 22.0 to 5.0 kg/mol. Figure 3A shows that the performances of the CK-CG and K-means methods are similar to each other in deriving the smaller numbers of UCG beads such as 7009, 9300, 11160, and 13950, whereas they exhibit distinct difference for large UCG numbers. The time cost of coarse-graining into 16740, 18600, 20925, and 27900 UCG beads by using CK-CG is almost half of that of K-means. For instance, coarse-graining into 27900 beads by CK-CG needs merely 149 minutes, while the K-means method costs 339 min. Also, it is seen that the time costs of both methods for deriving 23914 and 33480 UCG beads are comparable to each other. However, no case in which K-means was obviously superior to CK-CG was observed in our ten tests. Thus, the performance of CK-CG in coarse-graining a fixed-length microtubule into different beads outperforms the K-means method by 47.4 ± 44.4 hours with 95% confidence interval, from the statistical perspective. In this test, no simple correlation such as the linear relationship was found to exist between the time cost and the number of coarse-grained beads.

FIGURE 3
www.frontiersin.org

FIGURE 3. The time cost for building models with (A) different numbers of UCG beads and (B) different lengths. The blue squares and red dots represent the results from the K-means and CK-CG methods, respectively. The integers in plot (A) denote the exact time cost in minutes. In the plot (b), the set of data are fitted with two quadratic curves, respectively.

Next, we tested how the time cost changes with the increased lengths of microtubules under the same computational condition. The UCG models were built with the constraint of 7609 beads per micrometer, with the microtubule lengths changing from 1 to 10 μm. Figure 3B shows the comparison of time cost of using CK-CG and K-means methods for coarse-graining. It is surprisingly observed that the time cost exhibits a regular quadratic relationship with the increased length. It is found that the performances of two methods are nearly the same as the microtubules are short below 5 μm. Beyond that, the difference of time cost between them gradually becomes large with the increased length of microtubule. At the length of 10 μm, the CK-CG method saves 10 h more than K-means to optimize a total number of 76090 UCG beads. This comparison indicates that a good initial guess of center positions is important for accelerating the convergence of data clustering, especially for the biomolecules with the length at large scale.

Determination of a Minimal Model for Microtubule

After validating the efficiency of CK-CG for coarse-graining micrometer microtubules, we are ready to construct a minimal model of microtubule for its mechanical simulation. The so-call “minimal model” must meet two requirements for functional study: one is that the number of UCG beads for a given length is so minimal that it can reduce the simulation time to the maximal extent. Secondly, it is expected that this model is capable of maintaining the geometric feature of microtubule, especially for the functional grooves. According to the two requirements, we designed a minimal model with the possible minimal number of beads, which is illustrated in Figure 4.

FIGURE 4
www.frontiersin.org

FIGURE 4. (A) An illustration of the minimal number of UCG beads required for maintaining the grooves of microtubule. (B) The transverse and longitude views of the built models. (C) Zoomed top and side views of the density map of 8 adjacent tubulins, labeled with the upper cases A-H. The different colors denote the density regions which are clustered into the UCG beads, with the corresponding UCG bead centers shown.

As shown in Figure 4A, the groove along the longitudinal direction of a microtubule is formed by two adjacent protofilaments. The simplest and minimal UCG model to characterize the feature of a groove is to represent each tubulin with two connected beads, with one shared by two neighboring tubulins. Since the average molar mass of tubulin is 55 kg/mol (Fratzl and Weinkamer, 2007), the estimated average molar mass of UCG beads is 22.5 kg/mol. To derive the positions of UCG beads, we set the average mass with the value of 22 kg/mol in the coarse-graining by CK-CG, which lead to a number of 7609 beads per 1 μm in Coarse-Graining Microtubules With Different Bead Resolutions and Lengths. Figure 4B shows the side and transverse views of the resulting microtubule model, where each UCG bead is virtualized with solid dots and the 13 protofilaments are in different colors.

Figure 4C shows the zoomed top and side views of a segment of density map of 8 tubulins, labeled with the upper cases A-H, where the density points belonging to the same UCG bead centers are in the same color. From the top and the side views, we can see clearly that the bottoms of most grooves are well represented by one UCG bead, such as the density regions connecting the adjacent “AB”, “DE”, or “BC” tubulins. Meanwhile, two adjacent tubulins such as “DG”, “BE”, or “EH” along a protofilament share the same beads. From the view of UCG beads, it can be judged that the UCG beads representing the groove bottoms of tubulins usually reside in a lower position, while the beads representing the tubulins along the protofilaments stay in a high position. Thus, the coarse-grained beads between two adjacent protofilaments maintain a geometric “V” shape along the contour surface of the microtubules, which is potentially important for functional studies including interacting proteins.

Parameterization of Ultra-Coarse-Grained Model

Using the potential energy of MVP-ANM model of Eq. 1 to describe the interaction between pairwise UCG beads requires determining the two parameters R and A in Eq. 2. As mentioned in Multiscale Virtual Particle-Based Anisotropic Network Model, the parameter R could be defined as one-third of minimal cutoff distance dc. In a word, we need to determine a minimal cutoff at first. As shown in Figure 5A, we constructed the UCG models with the cutoff set as the tentative values of 6.0, 6.5, and 7.0 nm, respectively. Meanwhile, the parameter A adopts the empirical values of 6.0 × 10−3, 8.0 × 10−3, and 1.0 × 10−2 kg−1·s−1. After the 60 ms MD simulation, we found that all the three microtubule models with the cutoff of 6.0 nm have collapsed and their structures could not keep stable anymore. Besides, the UCG models with the cutoffs of 6.5 and 7.0 nm are stable in the long time MD simulation and maintain their tube shape well. The simulation results also indicate that the stabilities of microtubules did not depend on the choices of magnitudes of the parameter A.

FIGURE 5
www.frontiersin.org

FIGURE 5. (A) The structures of constructed UCG models for microtubules at the initial simulation time and after 60 ms. The cutoff distances adopt the values of 6.0, 6.5, and 7.0 nm for three models, with the tested parameter A adopting 6 × 10−3, 8 × 10−3, and 1 × 10−2 kg−1 s−1, respectively. (B) The UCG bond distances below 4.4 nm are colored in gray, while those bonds in 4.4–6.6 nm are in pink. (C) The UCG bonds below 4.0 nm in gray and those in 4.0–4.4 nm are colored in orange.

To rationalize the choice of a cutoff of 6.5 nm, we analyze the number of harmonic bonds formed between pairwise UCG beads. In Eq. 2, the strength of harmonic bonds is weakened with an increasing distance R. When dc adopts the cutoff 6.5 nm, we find that numerous bonds that describe strong interactions are defined between UCG beads, shown as gray lines in Figure 5B (Ji and Feng, 2011a; Kononova et al., 2014). Within the distance between 4.4 and 6.6 nm (2R–3R), bonds describing nonbonded interactions are found (shown in pink), and they are roughly 40 to 150 times weaker than the bonds within 4.4 nm. In contrast, Figure 5C shows the formed bonds within 4.0 nm in gray and the bonds in 4.0–6.0 nm in orange, as the cutoff adopts 6.0 nm. It is obvious that the two adjacent protofilaments in the UCG model do not feature sufficient bonding to keep the whole cylinder structure stable. Thus, the reason of collapse with the cutoff 6.0 nm lies in the lack of sufficient pairwise interactions to maintain the microtubule’s stability.

Eventually, we chose the cutoff to be 6.5 nm for all models and the resulting R value is 2.2 nm. The scaling factor A is determined by fitting the longitudinal Young’s modulus of a 100 nm short microtubule. There are two reasons that we need to fit Young’s modulus to derive the parameter A. One is that the longitudinal Young’s modulus is an independent property and has no influence on the measurement of the persistence length according to Eq. 3. The second is that the longitudinal Young’s modulus is independent of the length L of microtubule and only intrinsically determined by the strength of force constant γij. Thus, the parameter A in Eq. 2 is tightly correlated with the longitudinal Young’s modulus.

Kis et al. (2002) reported that a microtubule had a longitude Young’s modulus of 100 MPa under 300 K. In our force-clamped MD simulation, one terminal of the 100 nm microtubule is fixed, while the other one is pulled with a constant force of 2.36 nN. In Figure 6A, different values of the parameter A are used for UCG models to measure the difference between the simulated displacements and the experimental result around 13.0 nm. As the parameter A adopts the value of 1.0 × 10−2 kg−1·s−1, the simulated displacement of microtubule reaches the experimental value well. So far, both the parameters R and A for the models of microtubules are determined and ready for the subsequent simulation of mechanical properties. In addition, we also performed force-clamped UCG MD simulations with the same parameters (A = 1.0 × 10−2 kg−1·s−1 and R = 2.2 nm) under different pulling forces ranging from 1.04 nN to 2.36 nN, as shown in Figure 6B. It can be observed that the equilibrated displacements have a good linear relationship with the applied forces. Based on the slope of the linear relationship, Young’s modulus is calculated to be 99.2 MPa, which reproduces Young’s modulus under different forces and verifies the validity of UCG model (Kis et al., 2002).

FIGURE 6
www.frontiersin.org

FIGURE 6. (A) The relationship of the terminal displacements of microtubule models with the different values used for the parameter A in the units of kg−1 s−1. The black dashed line represents the experimental displacement and the colorful curves denote the displacements of microtubule models changing with time under an external force of 2360 pN in force-clamp MD simulations. (B) The curves of the displacements of microtubule models under the different forces ranging from 1.04 to 2.36 nN, where the black hollow circles denote the average values of the end-to-end displacements sampled from 100–200 ns. The black line denotes the linear fitting of the circles.

Evaluation of Persistence Length

The persistence length of microtubule is evaluated by using MD simulation based on its thermal fluctuation. It should be noted that the experimental measurement of thermal fluctuation usually occurred at time scale of seconds. However, it is impossible to perform such a long time MD simulation in silicon, even by employing the UCG model. The previous simulation performed by Ding et al. (Ding and Xu, 2011) for evaluating persistence length from thermal fluctuation was 50 μs. In our UCG MD simulation, we used an integration step of 3 ps and the simulation time for microtubules with micrometer length could reach milliseconds. A number of MD trajectories were saved for further analysis. Figure 7A shows the recorded intersection angle of a 2 μm microtubule changing with time in MD simulation. It can be seen that the angle fluctuates between 0.0 and 17.0 degrees and the angles almost exhibit a Gaussian distribution in the inset. The data in Figure 7A suggest that our MD samplings are sufficient to evaluate persistence length.

FIGURE 7
www.frontiersin.org

FIGURE 7. (A) The angle fluctuation of a 2 μm microtubule model within 60 ms. The inset displays the histogram distribution of angles with a black fitted Gaussian curve. (B) The orange solid squares show the relationship of the simulated persistence length and flexural rigidity with respect to the contour length, fitted with a blue curve based on Eq. 4. The green hollow circles denote the experimental measured results (Pampaloni et al., 2006). The simulation time ranges from microseconds to milliseconds for the different microtubules.

Based on the sampled structures of microtubules with different lengths, we calculated the persistence lengths corresponding to different lengths according to Eq. 3. Figure 7B shows the change of calculated persistence length lp changing with the contour length, with the longest microtubule being 12 µm. The calculated flexural rigidity κ according to the formula κ = lpkBT (Pampaloni et al., 2006) is plotted at the right side of y-axis, where kB is the Boltzmann constant and T is the simulation temperature. The simulated data, denoted by the yellow squares, clearly verify that the persistence length increases nonlinearly with the length. During the simulated range of length from 0.5 to 12 µm, the variation tendency of lp is in a good agreement with the experimental measurements (Pampaloni et al., 2006). In order to evaluate the magnitude of the persistence length lp, we fitted the data of calculated lp with the microtubule length L according to Eq. 4 derived by Pampaloni et al. (Pampaloni et al., 2006). As shown in Eq. 4 below, the lp has a complicated functional dependence on the lp at infinity:

lp=lp(1+3EIGkAL2)1(4)

where E is the longitude Young’s modulus and adopts the experimental value of 100 MPa (Kis et al., 2002) here, A is the area of cross section with 25 nm2 for a microtubule composed of 13 protofilaments, I is the area moment of the cross section with 16670 nm4, k is a geometric factor adopting 0.72, and G is the transverse sheer modulus. Figure 7B shows the fitted blue curve based on Eq. 4, similar to the fitted one in the previous work (Pampaloni et al., 2006). The persistence length lp derived from our fitting is 7.0 mm, consistent with the experimental value of 6.3 ± 0.8 mm for long microtubules. They also calculated the term of (3EI/GkA)0.5 to be 21 μm (Pampaloni et al., 2006), which was considered as a threshold between short and long microtubules. We further compare our simulation results with the experimental data denoted by the hollow green circles below 12 μm (Pampaloni et al., 2006). It appears that our simulation results fit quiet well with the experimental measurements within the contour length of 10 μm. Besides, in the experimental measurement of microtubules up to 47.5 μm, the persistence length lp seems to approach a constant beyond 21 μm. At present, our data obtained from UCG MD simulation based on the microtubule models with the lengths of 0.5–12 μm generated good results for lp and lp agreeable with known experimental results. In addition, the calculated values of the flexural rigidity κ for the length of 2–12 μm are in the range of 0.89–11.5 × 10−24 Nm2, falling in the range of the reported experimental values 3.7–40.0 × 10−24 Nm2 (Gittes et al., 1993; Venier et al., 1994; Kurachi et al., 1995; Felgner et al., 1996). However, a simulation of microtubules beyond 21 μm still faces a severe challenge memory requirement for calculating the interactions between pairwise particles in the UCG models.

Conclusion

In this work, we introduce a systematic strategy to build a UCG model for microtubule with micrometer length from the cryo-EM density data. For the large microtubule with micrometer length, the CK-CG serves as an efficient coarse-graining method to derive its UCG representation. The comparison of time cost in coarse-graining the microtubules with the lengths from 1 to 10 μm demonstrates that CK-CG is superior to the conventional K-means in coarse-graining large systems. A reasonable minimal model was then determined for modeling the mechanical properties of microtubules. By means of using the MVP-ANM force fields, we parameterized the UCG beads and constructed the corresponding models for the study of mechanical properties.

We estimated the persistence length of microtubules with different lengths by using the force-clamp MD simulation based on thermal fluctuation. The obtained relationship between the persistence length and contour length accords well with the previous theoretical and experimental results, which demonstrates the validity of our method in modeling the mechanical properties of biomaterials. It is highlighted that our method has the advantage of being capable of constructing models directly from the low-resolution density data of cryo-EM, without any atomistic details required. Thus, it is particularly suitable for the functional study of macrobiomolecules for which only low-resolution cryo-EM data are available. Our model relies on low-resolution cryo-EM data, which is sufficient given our choice for an ultra-coarse-grained nature of our model. However, it remains to be examined if higher resolution cryo-EM structures can further improve models, even those on the very coarse scale such as the UCG model presented here. However, the current UCG model is still a sort of ENMs and it is limited to be valid in the study of biomechanical properties of biomaterials with small deformations under external forces. Future work will be invested into developing more physically meaningful energy potentials for UCG model and using our method to investigate longer microtubules and eventually the mechanical properties of a bundle of microtubules such as those present in axons.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author Contributions

FX designed the project. JZ performed the simulation and analyzed the data. All people discussed and wrote the paper together.

Funding

This work was supported by the National Natural Science Foundation of China (Grant Nos. 21773065 and 22073029), Nanyang Technological University Startup Grant M4081842, Singapore Ministry of Education Academic Research fund Tier 1 RG31/18 and RG109/19, Tier 2 MOE2018-T2-1-033. FG acknowledges support by the Klaus Tschira Foundation and by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy–2082/1–390761711.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors acknowledge the support of the NYU-ECNU Center for Computational Chemistry at NYU Shanghai. They also thank the ECNU Public Platform for Innovation (001) for providing computer time.

References

Atilgan, A. R., Durell, S. R., Jernigan, R. L., Demirel, M. C., Keskin, O., and Bahar, I. (2001). Anisotropy of fluctuation dynamics of proteins with an elastic network model. Biophys J. 80 (1), 505–515. doi:10.1016/S0006-3495(01)76033-X

PubMed Abstract | CrossRef Full Text | Google Scholar

Bahar, I., Atilgan, A. R., and Erman, B. (1997). Direct evaluation of thermal fluctuations in proteins using a single-parameter harmonic potential. Fold Des. 2 (3), 173–181. doi:10.1016/s1359-0278(97)00024-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Conde, C., and Cáceres, A. (2009). Microtubule assembly, organization and dynamics in axons and dendrites. Nat Rev Neurosci. 10 (5), 319–332. doi:10.1038/nrn2631

PubMed Abstract | CrossRef Full Text | Google Scholar

Cross, R. A. (2019). Microtubule lattice plasticity. Curr Opin Cell Biol. 56, 88–93. doi:10.1016/j.ceb.2018.10.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Dama, J. F., Jin, J., and Voth, G. A. (2017). The theory of ultra-coarse-graining. 3. Coarse-grained sites with rapid local equilibrium of internal states. J Chem Theory Comput. 13 (3), 1010–1022. doi:10.1021/acs.jctc.6b01081

PubMed Abstract | CrossRef Full Text | Google Scholar

Dama, J. F., Sinitskiy, A. V., McCullagh, M., Weare, J., Roux, B., Dinner, A. R., et al. (2013). The theory of ultra-coarse-graining. 1. General principles. J Chem Theory Comput. 9 (5), 2466–2480. doi:10.1021/ct4000444

PubMed Abstract | CrossRef Full Text | Google Scholar

Davtyan, A., Dama, J. F., Sinitskiy, A. V., and Voth, G. A. (2014). The theory of ultra-coarse-graining. 2. Numerical implementation. J Chem Theory Comput. 10 (12), 5265–5275. doi:10.1021/ct500834t

PubMed Abstract | CrossRef Full Text | Google Scholar

de Pablo, P. J., Schaap, I. A., MacKintosh, F. C., and Schmidt, C. F. (2003). Deformation and collapse of microtubules on the nanometer scale. Phys Rev Lett. 91 (9), 098101–098104. doi:10.1103/PhysRevLett.91.098101

PubMed Abstract | CrossRef Full Text | Google Scholar

Deriu, M. A., Enemark, S., Soncini, M., Montevecchi, F. M., and Redaelli, A. (2007). Tubulin: from atomistic structure to supramolecular mechanical properties. J Mater Sci. 42 (21), 8864–8872. doi:10.1007/s10853-007-1784-6

CrossRef Full Text | Google Scholar

Deriu, M.A., Soncin, M., and Redaelli, A. (2008). Mechanical model of the tubulin dimer based on molecular dynamics simulations. J Biomech Eng. 130 (4), 41008–41014. doi:10.1115/1.2913330

CrossRef Full Text | Google Scholar

Deriu, M. A., Soncini, M., Orsi, M., Patel, M., Essex, J. W., Montevecchi, F. M., et al. (2010). Anisotropic elastic network modeling of entire microtubules. Biophys J. 99 (7), 2190–2199. doi:10.1016/j.bpj.2010.06.070

PubMed Abstract | CrossRef Full Text | Google Scholar

Dima, R. I., and Joshi, H. (2008). Probing the origin of tubulin rigidity with molecular simulations. Proc Natl Acad Sci USA. 105 (41), 15743–15748. doi:10.1073/pnas.0806113105

PubMed Abstract | CrossRef Full Text | Google Scholar

Ding, Y., and Xu, Z. (2011). Mechanics of microtubules from a coarse-grained model. BioNanoSci. 1 (4), 173–182. doi:10.1007/s12668-011-0027-0

CrossRef Full Text | Google Scholar

Dye, R. B., Fink, S. P., and Williams, J. (1993). Taxol-induced flexibility of microtubules and its reversal by MAP-2 and Tau. J Biol Chem. 268 (10), 6847–6850. doi:10.1016/s0021-9258(18)53113-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Esquivel-Rodríguez, J., and Kihara, D. (2013). Computational methods for constructing protein structure models from 3D electron microscopy maps. J Struct Biol. 184 (1), 93–102. doi:10.1016/j.jsb.2013.06.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Felgner, H., Frank, R., and Schliwa, M. (1996). Flexural rigidity of microtubules measured with the use of optical tweezers. J Cell Sci. 109 (2), 509–516. doi:10.1083/jcb.132.4.741

PubMed Abstract | CrossRef Full Text | Google Scholar

Franz, F., Daday, C., and Gräter, F. (2020). Advances in molecular simulations of protein mechanical properties and function. Curr Opin Struct Biol. 61, 132–138. doi:10.1016/j.sbi.2019.12.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Fratzl, P., and Weinkamer, R. (2007). Nature’s hierarchical materials. Progress in Materials Science. 52 (8), 1263–1334. doi:10.1016/j.pmatsci.2007.06.001

CrossRef Full Text | Google Scholar

Gittes, F., Meyhöfer, E., Baek, S., and Howard, J. (1996). Directional loading of the kinesin motor molecule as it buckles a microtubule. Biophys J. 70 (1), 418–429. doi:10.1016/S0006-3495(96)79585-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Gittes, F., Mickey, B., Nettleton, J., and Howard, J. (1993). Flexural rigidity of microtubules and actin filaments measured from thermal fluctuations in shape. J Cell Biol. 120 (4), 923–934. doi:10.1083/jcb.120.4.923

PubMed Abstract | CrossRef Full Text | Google Scholar

Grafmüller, A., Noya, E. G., and Voth, G. A. (2013). Nucleotide-dependent lateral and longitudinal interactions in microtubules. J Mol Biol. 425 (12), 2232–2246. doi:10.1016/j.jmb.2013.03.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Grafmüller, A., and Voth, G. A. (2011). Intrinsic bending of microtubule protofilaments. Structure. 19 (3), 409–417. doi:10.1016/j.str.2010.12.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Grange, M., Vasishtan, D., and Grünewald, K. (2017). Cellular electron cryo tomography and in situ sub-volume averaging reveal the context of microtubule-based processes. J Struct Biol. 197 (2), 181–190. doi:10.1016/j.jsb.2016.06.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Gräter, F., and Grubmüller, H. (2007). Fluctuations of primary ubiquitin folding intermediates in a force clamp. Journal of Structural Biology. 157 (3), 557–569. doi:10.1016/j.jsb.2006.11.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Grime, J. M., Dama, J. F., Ganser-Pornillos, B. K., Woodward, C. L., Jensen, G. J., Yeager, M., et al. (2016). Coarse-grained simulation reveals key features of HIV-1 capsid self-assembly. Nat Commun. 7, 11568–11578. doi:10.1038/ncomms11568

PubMed Abstract | CrossRef Full Text | Google Scholar

Grønbech-Jensen, N. (2019). Complete set of stochastic verlet-type thermostats for correct langevin simulations. Molecular Physics. 118 (8), e1662506–e1662525. doi:10.1080/00268976.2019.1662506

CrossRef Full Text | Google Scholar

Havelka, D., Deriu, M. A., Cifra, M., and Kučera, O. (2017). Deformation pattern in vibrating microtubule: structural mechanics study based on an atomistic approach. Sci Rep. 7 (1), 4227–4237. doi:10.1038/s41598-017-04272-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Hawkins, T. L., Sept, D., Mogessie, B., Straube, A., and Ross, J. L. (2013). Mechanical properties of doubly stabilized microtubule filaments. Biophys J. 104 (7), 1517–1528. doi:10.1016/j.bpj.2013.02.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Heuvel, M. G. L. v. d., Bolhuis, S., and Dekker, C. (2007). Persistence length measurements from stochastic single-microtubule trajectories. Nano Lett. 7 (10), 3138–3144. doi:10.1021/nl071696y

PubMed Abstract | CrossRef Full Text | Google Scholar

Isralewitz, B., Gao, M., and Schulten, K. (2001). Steered molecular dynamics and mechanical functions of proteins. Curr Opin Struct Biol. 11 (2), 224–2230. doi:10.1016/S0959-440X(00)00194-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Janson, M. E., and Dogterom, M. (2004). A bending mode analysis for growing microtubules: evidence for a velocity-dependent rigidity. Biophys J. 87 (4), 2723–2736. doi:10.1529/biophysj.103.038877

PubMed Abstract | CrossRef Full Text | Google Scholar

Ji, X. Y., and Feng, X. Q. (2011a). Coarse-grained mechanochemical model for simulating the dynamic behavior of microtubules. Phys Rev E Stat Nonlin Soft Matter Phys. 84 (3), 031933–031941. doi:10.1103/PhysRevE.84.031933

PubMed Abstract | CrossRef Full Text | Google Scholar

Ji, X. Y., and Feng, X. Q. (2011b). Mechanochemical modeling of dynamic microtubule growth involving sheet-to-tube transition. PLoS One. 6 (12), e29049–e29059. doi:10.1371/journal.pone.0029049

PubMed Abstract | CrossRef Full Text | Google Scholar

Jin, J., Han, Y., and Voth, G. A. (2018). Ultra-coarse-grained liquid state models with implicit hydrogen bonding. J Chem Theory Comput. 14 (12), 6159–6174. doi:10.1021/acs.jctc.8b00812

PubMed Abstract | CrossRef Full Text | Google Scholar

Jin, J., and Voth, G. A. (2018). Ultra-coarse-grained models allow for an accurate and transferable treatment of interfacial systems. J Chem Theory Comput. 14 (4), 2180–2197. doi:10.1021/acs.jctc.7b01173

PubMed Abstract | CrossRef Full Text | Google Scholar

Kasas, S., Cibert, C., Kis, A., De Los Rios, P., Riederer, B. M., Forró, L., et al. (2004a). Oscillation modes of microtubules. Biol Cell. 96 (9), 697–700. doi:10.1016/j.biolcel.2004.09.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Kasas, S., Kis, A., Riederer, B. M., Forró, L., Dietler, G., and Catsicas, S. (2004b). Mechanical properties of microtubules explored using the finite elements method. Chemphyschem. 5 (2), 252–257. doi:10.1002/cphc.200300799

PubMed Abstract | CrossRef Full Text | Google Scholar

Katkar, H. H., Davtyan, A., Durumeric, A. E. P., Hocky, G. M., Schramm, A. C., De La Cruz, E. M., et al. (2018). Insights into the cooperative nature of ATP hydrolysis in actin filaments. Biophys J. 115 (8), 1589–1602. doi:10.1016/j.bpj.2018.08.034

PubMed Abstract | CrossRef Full Text | Google Scholar

Kis, A., Kasas, S., Babić, B., Kulik, A J., Benoît, W., Briggs, G. A., et al. (2002). Nanomechanics of microtubules. Phys Rev Lett. 89 (24), 248101–248104. doi:10.1103/PhysRevLett.89.248101

PubMed Abstract | CrossRef Full Text | Google Scholar

Kis, A., Kasas, S., Kulik, A. J., Catsicas, S., and Forró, L. (2008). Temperature-dependent elasticity of microtubules. Langmuir. 24 (12), 6176–6181. doi:10.1021/la800438q

PubMed Abstract | CrossRef Full Text | Google Scholar

Kollman, J. M., Polka, J. K., Zelter, A., Davis, T. N., and Agard, D. A. (2010). Microtubule nucleating gamma-TuSC assembles structures with 13-fold microtubule-like symmetry. Nature. 466 (7308), 879–882. doi:10.1038/nature09207

PubMed Abstract | CrossRef Full Text | Google Scholar

Kononova, O., Kholodov, Y., Theisen, K. E., Marx, K. A., Dima, R. I., Ataullakhanov, F. I., et al. (2014). Tubulin bond energies and microtubule biomechanics determined from nanoindentation in silico. J Am Chem Soc. 136 (49), 17036–17045. doi:10.1021/ja506385p

PubMed Abstract | CrossRef Full Text | Google Scholar

Kroy, K., and Frey, E. (1996). Force-extension relation and plateau modulus for wormlike chains. Phys Rev Lett. 77 (2), 306–309. doi:10.1103/PhysRevLett.77.306

PubMed Abstract | CrossRef Full Text | Google Scholar

Kurachi, M., Hoshi, M., and Tashiro, H. (1995). Buckling of a single microtubule by optical trapping forces: direct measurement of microtubule rigidity. Cell Motil Cytoskeleton. 30 (3), 221–228. doi:10.1002/cm.970300306

PubMed Abstract | CrossRef Full Text | Google Scholar

Lazarus, C., Soheilypour, M., and Mofrad, M. R. (2015). Torsional behavior of axonal microtubule bundles. Biophys J. 109 (2), 231–239. doi:10.1016/j.bpj.2015.06.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Le Goff, L., Hallatschek, O., Frey, E., and Amblard, F. (2002). Tracer studies on f-actin fluctuations. Phys Rev Lett. 89 (25), 258101–258104. doi:10.1103/PhysRevLett.89.258101

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, M., Zhang, J. Z., and Xia, F. (2016a). Constructing optimal coarse-grained sites of huge biomolecules by fluctuation maximization. J Chem Theory Comput. 12 (4), 2091–2100. doi:10.1021/acs.jctc.6b00016

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, M., Zhang, J. Z., and Xia, F. (2016b). A new algorithm for construction of coarse-grained sites of large biomolecules. J Comput Chem. 37 (9), 795–804. doi:10.1002/jcc.24265

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Z., Cai, X., Liu, Y., and Zhu, B. (2019). A novel gaussian-bernoulli based convolutional deep belief networks for image feature extraction. Neural Process Lett. 49 (4), 305–319. doi:10.1007/s11063-017-9751-y

CrossRef Full Text | Google Scholar

Ming, D., Kong, Y., Lambert, M. A., Huang, Z., and Ma, J. (2002). How to describe protein motion without amino acid sequence and atomic coordinates. Proc Natl Acad Sci USA. 99 (13), 8620–8625. doi:10.1073/pnas.082148899

PubMed Abstract | CrossRef Full Text | Google Scholar

Molodtsov, M. I., Ermakova, E. A., Shnol, E. E., Grishchuk, E. L., McIntosh, J. R., and Ataullakhanov, F. I. (2005). A Molecular-Mechanical Model of the Microtubule. Biophys J. 88 (5), 3167–3179. doi:10.1529/biophysj.104.051789

PubMed Abstract | CrossRef Full Text | Google Scholar

Opron, K., Xia, K., and Wei, G. W. (2014). Fast and anisotropic flexibility-rigidity index for protein flexibility and fluctuation analysis. J Chem Phys. 140 (23), 234105–234123. doi:10.1063/1.4882258

PubMed Abstract | CrossRef Full Text | Google Scholar

Pampaloni, F., Lattanzi, G., Jonás, A., Surrey, T., Frey, E., and Florin, E. L. (2006). Thermal fluctuations of grafted microtubules provide evidence of a length-dependent persistence length. Proc Natl Acad Sci USA. 103 (27), 10248–10253. doi:10.1073/pnas.0603931103

PubMed Abstract | CrossRef Full Text | Google Scholar

Peter, S. J., and Mofrad, M. R. (2012). Computational modeling of axonal microtubule bundles under tension. Biophys J. 102 (4), 749–757. doi:10.1016/j.bpj.2011.11.4024

PubMed Abstract | CrossRef Full Text | Google Scholar

Plimpton, S. (1995). Fast parallel algorithms for short-range molecular dynamics. Journal of Computational Physics. 117 (1), 1–19. doi:10.1006/jcph.1995.1039

CrossRef Full Text | Google Scholar

Poma, A. B., Chwastyk, M., and Cieplak, M. (2017). Elastic moduli of biological fibers in a coarse-grained model: crystalline cellulose and β-amyloids. Phys Chem Chem Phys. 19 (41), 28195–28206. doi:10.1039/c7cp05269c

PubMed Abstract | CrossRef Full Text | Google Scholar

Poma, A. B., Li, M. S., and Theodorakis, P. E. (2018). Generalization of the elastic network model for the study of large conformational changes in biomolecules. Phys Chem Chem Phys. 20 (25), 17020–17028. doi:10.1039/c8cp03086c

PubMed Abstract | CrossRef Full Text | Google Scholar

Schaap, I. A., Carrasco, C., de Pablo, P. J., MacKintosh, F. C., and Schmidt, C. F. (2006). Elastic Response, Buckling, and Instability of Microtubules under Radial Indentation. Biophys J. 91 (4), 1521–1531. doi:10.1529/biophysj.105.077826

PubMed Abstract | CrossRef Full Text | Google Scholar

Sept, D., and MacKintosh, F. C. (2010). Microtubule elasticity: connecting all-atom simulations with continuum mechanics. Phys Rev Lett. 104 (1), 018101–018104. doi:10.1103/PhysRevLett.104.018101

PubMed Abstract | CrossRef Full Text | Google Scholar

Sinitskiy, A. V., Saunders, M. G., and Voth, G. A. (2012). Optimal number of coarse-grained sites in different components of large biomolecular complexes. J Phys Chem B. 116 (29), 8363–8374. doi:10.1021/jp2108895

PubMed Abstract | CrossRef Full Text | Google Scholar

Soheilypour, M., Peyro, M., Peter, S. J., and Mofrad, M. R. (2015). Buckling behavior of individual and bundled microtubules. Biophys J. 108 (7), 1718–1726. doi:10.1016/j.bpj.2015.01.030

PubMed Abstract | CrossRef Full Text | Google Scholar

Stevens, M. J. (2017). The long persistence length of model tubules. J Chem Phys. 147 (4), 044902–044907. doi:10.1063/1.4994913

PubMed Abstract | CrossRef Full Text | Google Scholar

Szatkowski, L., Merz, D. R., Jiang, N., Ejikeme, I., Belonogov, L., Ross, J. L., et al. (2019). Mechanics of the microtubule seam interface probed by molecular simulations and in vitro severing experiments. J Phys Chem B. 123 (23), 4888–4900. doi:10.1021/acs.jpcb.9b03059

PubMed Abstract | CrossRef Full Text | Google Scholar

Taute, K. M., Pampaloni, F., Frey, E., and Florin, E. L. (2008). Microtubule dynamics depart from the wormlike chain model. Phys Rev Lett. 100 (2), 028102–028105. doi:10.1103/PhysRevLett.100.028102

PubMed Abstract | CrossRef Full Text | Google Scholar

Theisen, K. E., Desai, N. J., Volski, A. M., and Dima, R. I. (2013). Mechanics of severing for large microtubule complexes revealed by coarse-grained simulations. J Chem Phys. 139 (12), 121926–121937. doi:10.1063/1.4819817

PubMed Abstract | CrossRef Full Text | Google Scholar

Tong, D., and Voth, G. A. (2020). Microtubule simulations provide insight into the molecular mechanism underlying dynamic instability. Biophys J. 118 (12), 2938–2951. doi:10.1016/j.bpj.2020.04.028

PubMed Abstract | CrossRef Full Text | Google Scholar

VanBuren, V., Cassimeris, L., and Odde, D. J. (2005). Mechanochemical model of microtubule structure and self-assembly kinetics. Biophys J. 89 (5), 2911–2926. doi:10.1529/biophysj.105.060913

PubMed Abstract | CrossRef Full Text | Google Scholar

VanBuren, V., Odde, D. J., and Cassimeris, L. (2002). Estimates of lateral and longitudinal bond energies within the microtubule lattice. Proc Natl Acad Sci USA. 99 (8), 6035–6040. doi:10.1073/pnas.040639310110.1073/pnas.092504999

PubMed Abstract | CrossRef Full Text | Google Scholar

Venier, P., Maggs, A. C., Carlier, M. F., and Pantaloni, D. (1994). Analysis of microtubule rigidity using hydrodynamic flow and thermal fluctuations. J Biol Chem. 269 (18), 13353–13360. doi:10.1016/0009-3084(94)02321-210.1016/s0021-9258(17)36840-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Vinckier, A. (1996). Dynamical and mechanical study of immobilized microtubules with atomic force microscopy. J. Vac. Sci. Technol. B. 14 (2), 1427–1431. doi:10.1116/1.589113

CrossRef Full Text | Google Scholar

Wagner, O., Zinke, J., Dancker, P., Grill, W., and Bereiter-Hahn, J. (1999). Viscoelastic properties of f-actin, microtubules, f-actin/alpha-actinin, and f-actin/hexokinase determined in microliter volumes with a novel nondestructive method. Biophys J. 76 (5), 2784–2796. doi:10.1016/S0006-3495(99)77432-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Wells, D. B., and Aksimentiev, A. (2010). Mechanical properties of a complete microtubule revealed through molecular dynamics simulation. Biophys J. 99 (2), 629–637. doi:10.1016/j.bpj.2010.04.038

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, Y. T., and Adnan, A. (2018). Damage and failure of axonal microtubule under extreme high strain rate: an in-silico molecular dynamics study. Sci Rep. 8 (1), 12260. doi:10.1038/s41598-018-29804-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, Z., Zhang, Y., Zhang, J. Z., Xia, K., and Xia, F. (2020). Determining optimal coarse‐grained representation for biomolecules using internal cluster validation indexes. J Comput Chem. 41 (1), 14–20. doi:10.1002/jcc.26070

PubMed Abstract | CrossRef Full Text | Google Scholar

Xia, F., and Lu, L. (2012). Multiscale coarse-graining via normal mode analysis. J Chem Theory Comput. 8 (11), 4797–4806. doi:10.1021/ct3006387

PubMed Abstract | CrossRef Full Text | Google Scholar

Xia, F., Tong, D., and Lu, L. (2013a). Robust heterogeneous anisotropic elastic network model precisely reproduces the experimental b-factors of biomolecules. J Chem Theory Comput. 9 (8), 3704–3714. doi:10.1021/ct4002575

PubMed Abstract | CrossRef Full Text | Google Scholar

Xia, K., Opron, K., and Wei, G. W. (2013b). Multiscale multiphysics and multidomain models--flexibility and rigidity. J Chem Phys. 139 (19), 194109–194124. doi:10.1063/1.4830404

PubMed Abstract | CrossRef Full Text | Google Scholar

Xia, K. (2017). Multiscale virtual particle based elastic network model (MVP-ENM) for normal mode analysis of large-sized biomolecules. Phys Chem Chem Phys. 20 (1), 658–669. doi:10.1039/c7cp07177a

PubMed Abstract | CrossRef Full Text | Google Scholar

Zapp, C., Obarska-Kosinska, A., Rennekamp, B., Kurth, M., Hudson, D. M., Mercadante, D., et al. (2020). Mechanoradicals in tensed tendon collagen as a source of oxidative stress. Nat Commun. 11 (1), 2315–2322. doi:10.1038/s41467-020-15567-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z., Lu, L., Noid, W. G., Krishna, V., Pfaendtner, J., and Voth, G. A. (2008). A systematic methodology for defining coarse-grained sites in large biomolecules. Biophys J. 95 (11), 5073–5083. doi:10.1529/biophysj.108.139626

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z., Pfaendtner, J., Grafmüller, A., and Voth, G. A. (2009). Defining coarse-grained representations of large biomolecules and biomolecular complexes from elastic network models. Biophys J. 97 (8), 2327–2337. doi:10.1016/j.bpj.2009.08.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Cao, Z., and Xia, F. (2017a). Construction of ultra-coarse-grained model of protein with a Gō-like potential. Chemical Physics Letters. 681, 1–6. doi:10.1016/j.cplett.2017.05.039

CrossRef Full Text | Google Scholar

Zhang, Y., Cao, Z., Zhang, J. Z., and Xia, F. (2017b). Performance comparison of systematic methods for rigorous definition of coarse-grained sites of large biomolecules. J Chem Inf Model. 57 (2), 214–222. doi:10.1021/acs.jcim.6b00683

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Cao, Z., Zhang, J. Z., and Xia, F. (2020). Double-well ultra-coarse-grained model to describe protein conformational transitions. J Chem Theory Comput. 16 (10), 6678–6689. doi:10.1021/acs.jctc.0c00551

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Xia, K., Cao, Z., Gräter, F., and Xia, F. (2019). A new method for the construction of coarse-grained models of large biomolecules from low-resolution cryo-electron microscopy data. Phys Chem Chem Phys. 21 (19), 9720–9727. doi:10.1039/c9cp01370a

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z., and Voth, G. A. (2010). Coarse-grained representations of large biomolecular complexes from low-resolution structural data. J Chem Theory Comput. 6 (9), 2990–3002. doi:10.1021/ct100374a

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: microtubule, persistence length, ultra-coarse-grained model, mechanical property, convolutional and K-means coarse-graining

Citation: Zha J, Zhang Y, Xia K, Gräter F and Xia F (2021) Coarse-Grained Simulation of Mechanical Properties of Single Microtubules With Micrometer Length. Front. Mol. Biosci. 7:632122. doi: 10.3389/fmolb.2020.632122

Received: 22 November 2020; Accepted: 30 December 2020;
Published: 15 February 2021.

Edited by:

Fabio Trovato, Freie Universität Berlin, Germany

Reviewed by:

Gregory Voth, University of Chicago, United States
Adolfo Poma, Institute of Fundamental Technological Research, Polish Academy of Sciences, Poland

Copyright © 2021 Zha, Zhang, Xia, Gräter and Xia. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Fei Xia, fxia@chem.ecnu.edu.cn

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.