Nonlinear Finite Element Analysis of Nanoindentation of Viral Capsids

Melissa M. Gibbons,William S. Klug
DOI: https://doi.org/10.1103/PhysRevE.75.031901
2006-08-17
Abstract:Recent Atomic Force Microscope (AFM) nanoindentation experiments measuring mechanical response of the protein shells of viruses have provided a quantitative description of their strength and elasticity. To better understand and interpret these measurements, and to elucidate the underlying mechanisms, this paper adopts a course-grained modeling approach within the framework of three-dimensional nonlinear continuum elasticity. Homogeneous, isotropic, elastic, thick shell models are proposed for two capsids: the spherical Cowpea Chlorotic Mottle Virus (CCMV), and the ellipsocylindrical bacteriophage $\phi 29$. As analyzed by the finite element method, these models enable parametric characterization of the effects of AFM tip geometry, capsid dimensions, and capsid constitutive descriptions. The generally nonlinear force response of capsids to indentation is shown to be insensitive to constitutive details, and greatly influenced by geometry. Nonlinear stiffening and softening of the force response is dependent on the AFM tip dimensions and shell thickness. Fits of the models capture the roughly linear behavior observed in experimental measurements and result in estimates of Young's moduli of $\approx$280--360 MPa for CCMV and $\approx$4.5 GPa for $\phi 29$.
Biological Physics,Computational Physics
What problem does this paper attempt to address?
### What problems does this paper attempt to solve? This paper aims to study the mechanical responses of virus capsids in atomic force microscope (AFM) nano - indentation experiments through the nonlinear finite element analysis method. Specifically, the author attempts to solve the following key problems: 1. **Interpret experimental phenomena**: - Recent AFM experiments have shown that virus capsids still exhibit toughness and high elasticity under significant deformation, and their force - indentation responses are linear. This linear response is contrary to the phenomenon in traditional solid mechanics where large deformations are usually accompanied by nonlinear force responses. Therefore, the author hopes to explain why a linear response is observed in this case. 2. **Establish a more accurate model**: - Previous studies mostly used the linearized small - strain thin - shell theory to explain AFM experimental results. However, these models assume that the strain is small and are applicable to thin - shell structures where the thickness is much smaller than the overall structure size. For some thick - shell structures, such as the CCMV virus capsid, these assumptions may not hold. Therefore, the author attempts to abandon the linear assumption and use the three - dimensional nonlinear continuous elasticity theory to more comprehensively describe the mechanical behavior of the capsid. 3. **Evaluate the influence of geometric and material properties**: - Study the influence of different geometric parameters (such as AFM probe shape, capsid size and thickness) and material constitutive relations (such as linear and nonlinear hyperelastic models) on the force - indentation response curve. This helps to understand which factors are most critical to the mechanical response of the capsid and provides guidance for experimental design. 4. **Verify and compare experimental results**: - Through numerical simulation, verify the validity of existing experimental data and compare the simulation results with experimental measurements. This can not only verify the accuracy of the model but also provide theoretical support for further experiments. ### Formula summary To better understand the above problems, the following formulas are used in the paper to describe the material constitutive relations and mechanical responses: - **Strain energy function of the St. Venant - Kirchhoff model**: \[ W=\frac{\lambda_0}{2}(\text{tr}E)^2+\mu_0(\text{tr}E^2) \] where \(\lambda_0\) and \(\mu_0\) are the initial Lame constants, and \(E\) is the Green - strain tensor. - **Isovolume strain energy function of the Neo - Hookean model**: \[ W_{\text{iso}}=\frac{\mu_0}{2}(\bar{I}_1 - 3) \] where \(\mu_0\) is the initial shear modulus, and \(\bar{I}_1\) is the first invariant of the right Cauchy - Green tensor. - **Isovolume strain energy function of the Mooney - Rivlin model**: \[ W_{\text{iso}}=C_{10}(\bar{I}_1 - 3)+C_{01}(\bar{I}_2 - 3) \] where \(C_{10}\) and \(C_{01}\) are material constants, and \(\bar{I}_1\) and \(\bar{I}_2\) are the first and second invariants of the right Cauchy - Green tensor respectively. These formulas are used to describe the stress - strain relations of different material models and are solved by the finite element method to simulate the mechanical responses of the capsid under nano - indentation. ### Conclusion Through the above research, the author hopes to reveal the mechanical behavior of virus capsids under large deformations and provide valuable references for future experimental and theoretical research.