A viscoactive constitutive modeling framework with variational updates for the myocardium.

We present a constitutive modeling framework for contractile cardiac mechanics by formulating a single variational principle from which incremental stress-strain relations and kinetic rate equations for active contraction and relaxation can all be derived. The variational framework seamlessly incorporates the hyperelastic behavior of the relaxed and contracted tissue along with the rate - and length - dependent generation of contractile force. We describe a three-element, Hill-type model that unifies the active tension and active deformation approaches. As in the latter approach, we multiplicatively decompose the total deformation gradient into active and elastic parts, with the active deformation parametrizing the contractile Hill element. We adopt as internal variables the fiber, cross-fiber, and sheet normal stretch ratios. The kinetics of these internal variables are modeled via definition of a kinetic potential function derived from experimental force-velocity relations. Additionally, we account for dissipation during tissue deformation by adding a Newtonian viscous potential. To model the force activation, the kinetic equations are coupled with the calcium transient obtained from a cardiomyocyte electrophysiology model. We first analyze our model at the material point level using stress and strain versus time curves for different viscosity values. Subsequently, we couple our constitutive framework with the finite element method (FEM) and study the deformation of three-dimensional tissue slabs with varying cardiac myocyte orientation. Finally, we simulate the contraction and relaxation of an ellipsoidal left ventricular model and record common kinematic measures, such as ejection fraction, and myocardial tissue volume changes.

Statistical Physics of Viral Capsids with Broken Symmetry.

We present a model to understand quantitatively the role of symmetry breaking in assembly of macromolecular aggregates in general, and the protein shells of viruses in particular. A simple dodecahedral lattice model with a quadrupolar order parameter allows us to demonstrate how symmetry breaking may reduce the probability of assembly errors and, consequently, enhance assembly efficiency. We show that the ground state is characterized by large-scale cooperative zero-energy modes. In analogy with other models, this suggests a general physical principle: the tendency of biological molecules to generate symmetric structures competes with the tendency to break symmetry in order to achieve specific functional goals.

Elasticity theory of macromolecular aggregates.

We present a version of continuum elasticity theory applicable to aggregates of functional biomolecules at length scales comparable to that of the component molecules. Unlike classical elasticity theory, the stress and strain fields have mathematical discontinuities along the interfaces of the macromolecules, due to conformational incompatibility and large scale conformational transitions. The method is applied to the P-II to EI shape transition of the protein shell of the virus HK97. We show that protein residual stresses generated by incompatibility drive a "reverse buckling" transition from an icosahedral to a dodecahedral shape via a "critical" spherical shape, which can be identified as the P-II state.

Unlocking internal prestress from protein nanoshells.

The capsids of icosahedral viruses are closed shells assembled from a hexagonal lattice of proteins with fivefold angular defects located at the icosahedral vertices. Elasticity theory predicts that these disclinations are subject to an internal compressive prestress, which provides an explanation for the link between size and shape of capsids. Using a combination of experiment and elasticity theory we investigate the question of whether macromolecular assemblies are subject to residual prestress, due to basic geometric incompatibility of the subunits. Here we report the first direct experimental test of the theory: by controlled removal of protein pentamers from the icosahedral vertices, we measure the mechanical response of so-called "whiffle ball" capsids of herpes simplex virus, and demonstrate the signature of internal prestress locked into wild-type capsids during assembly.

Morphological phase diagram for lipid membrane domains with entropic tension.

Circular domains in phase-separated lipid vesicles with symmetric leaflet composition commonly exhibit three stable morphologies: flat, dimpled, and budded. However, stable dimples (i.e., partially budded domains) present a puzzle since simple elastic theories of domain shape predict that only flat and spherical budded domains are mechanically stable in the absence of spontaneous curvature. We argue that this inconsistency arises from the failure of the constant surface tension ensemble to properly account for the effect of entropic bending fluctuations. Formulating membrane elasticity within an entropic tension ensemble, wherein tension represents the free energy cost of extracting membrane area from thermal bending of the membrane, we calculate a morphological phase diagram that contains regions of mechanical stability for each of the flat, dimpled, and budded domain morphologies.

Squeezing protein shells: how continuum elastic models, molecular dynamics simulations, and experiments coalesce at the nanoscale.

The current rapid growth in the use of nanosized particles is fueled in part by our increased understanding of their physical properties and ability to manipulate them, which is essential for achieving optimal functionality. Here we report detailed quantitative measurements of the mechanical response of nanosized protein shells (viral capsids) to large-scale physical deformations and compare them with theoretical descriptions from continuum elastic modeling and molecular dynamics (MD). Specifically, we used nanoindentation by atomic force microscopy to investigate the complex elastic behavior of Hepatitis B virus capsids. These capsids are hollow, approximately 30 nm in diameter, and conform to icosahedral (5-3-2) symmetry. First we show that their indentation behavior, which is symmetry-axis-dependent, cannot be reproduced by a simple model based on Föppl-von Kármán thin-shell elasticity with the fivefold vertices acting as prestressed disclinations. However, we can properly describe the measured nonlinear elastic and orientation-dependent force response with a three-dimensional, topographically detailed, finite-element model. Next, we show that coarse-grained MD simulations also yield good agreement with our nanoindentation measurements, even without any fitting of force-field parameters in the MD model. This study demonstrates that the material properties of viral nanoparticles can be correctly described by both modeling approaches. At the same time, we show that even for large deformations, it suffices to approximate the mechanical behavior of nanosized viral shells with a continuum approach, and ignore specific molecular interactions. This experimental validation of continuum elastic theory provides an example of a situation in which rules of macroscopic physics can apply to nanoscale molecular assemblies.

Useful scars: Physics of the capsids of archaeal viruses.

We propose a physical model for the capsids of tailed archaeal viruses as viscoelastic membranes under tension. The fluidity is generated by thermal motion of scarlike structures that are an intrinsic feature of the ground state of large particle arrays covering surfaces with nonzero Gauss curvature. The tension is generated by a combination of the osmotic pressure of the enclosed genome and an extension force generated by filamentous structure formation that drives the formation of the tails. In continuum theory, the capsid has the shape of a surface of constant mean curvature: an unduloid. Particle arrays covering unduloids are shown to exhibit pronounced subdiffusive and diffusive single-particle transport at temperatures that are well below the melting temperature of defect-free particle arrays on a surface with zero Gauss curvature.

Regional segmentation of ventricular models to achieve repolarization dispersion in cardiac electrophysiology modeling.

The electrocardiogram (ECG) is one of the most significant outputs of a computational model of cardiac electrophysiology because it relates the numerical results to clinical data and is a universal tool for diagnosing heart diseases. One key features of the ECG is the T-wave, which is caused by longitudinal and transmural heterogeneity of the action potential duration (APD). Thus, in order to model a correct wave of repolarization, different cell properties resulting in different APDs must be assigned across the ventricular wall and longitudinally from apex to base. To achieve this requirement, a regional parametrization of the heart is necessary. We propose a robust approach to obtain the transmural and longitudinal segmentation in a general heart geometry without relying on ad hoc procedures. Our approach is based on auxiliary harmonic lifting analyses, already used in the literature to generate myocardial fiber orientations. Specifically, the solution of a sequence of Laplace boundary value problems allows parametrically controlled segmentation of both heart ventricles. The flexibility and simplicity of the proposed method is demonstrated through several representative examples, varying the locations and extents of the epicardial, midwall, and endocardial layers. Effects of the control parameters on the T-wave morphology are illustrated via computed ECGs.

Nanoindentation studies of full and empty viral capsids and the effects of capsid protein mutations on elasticity and strength.

The elastic properties of capsids of the cowpea chlorotic mottle virus have been examined at pH 4.8 by nanoindentation measurements with an atomic force microscope. Studies have been carried out on WT capsids, both empty and containing the RNA genome, and on full capsids of a salt-stable mutant and empty capsids of the subE mutant. Full capsids resisted indentation more than empty capsids, but all of the capsids were highly elastic. There was an initial reversible linear regime that persisted up to indentations varying between 20% and 30% of the diameter and applied forces of 0.6-1.0 nN; it was followed by a steep drop in force that is associated with irreversible deformation. A single point mutation in the capsid protein increased the capsid stiffness. The experiments are compared with calculations by finite element analysis of the deformation of a homogeneous elastic thick shell. These calculations capture the features of the reversible indentation region and allow Young's moduli and relative strengths to be estimated for the empty capsids.