Strain Gradient Programming in 3D Fibrous Hydrogels to Direct Graded Cell Alignment.

Biological tissues experience various stretch gradients which act as mechanical signaling from the extracellular environment to cells. These mechanical stimuli are sensed by cells, triggering essential signaling cascades regulating cell migration, differentiation, and tissue remodeling. In most previous studies, a simple, uniform stretch to 2D elastic substrates has been applied to analyze the response of living cells. However, induction of nonuniform strains in controlled gradients, particularly in biomimetic 3D hydrogels, has proven challenging. In this study, 3D fibrin hydrogels of manipulated geometry are stretched by a silicone carrier to impose programmable strain gradients along different chosen axes. The resulting strain gradients are analyzed and compared to finite element simulations. Experimentally, the programmed strain gradients result in similar gradient patterns in fiber alignment within the gels. Additionally, temporal changes in the orientation of fibroblast cells embedded in the stretched fibrin gels correlate to the strain and fiber alignment gradients. The experimental and simulation data demonstrate the ability to custom-design mechanical gradients in 3D biological hydrogels and to control cell alignment patterns. It provides a new technology for mechanobiology and tissue engineering studies.

Does the Extracellular Matrix Support Cell-Cell Communication by Elastic Wave Packets?

The extracellular matrix (ECM) is a fibrous network supporting biological cells and provides them a medium for interaction. Cells modify the ECM by applying traction forces, and these forces can propagate to long ranges and establish a mechanism of mechanical communication between neighboring cells. Previous studies have mainly focused on analysis of static force transmission across the ECM. In this study, we explore the plausibility of dynamic mechanical interaction, expressed as vibrations or abrupt fluctuations, giving rise to elastic waves propagating along ECM fibers. We use a numerical mass-spring model to simulate the longitudinal and transversal waves propagating along a single ECM fiber and across a 2D random fiber network. The elastic waves are induced by an active contracting cell (signaler) and received by a passive neighboring cell (receiver). We show that dynamic wave propagation may amplify the signal at the receiver end and support up to an order of magnitude stronger mechanical cues and longer-ranged communication relative to static transmission. Also, we report an optimal impulse duration corresponding to the most effective transmission, as well as extreme fast impulses, in which the waves are encaged around the active cell and do not reach the neighboring cell, possibly due to the Anderson localization effect. Finally, we also demonstrate that extracellular fluid viscosity reduces, but still allows, dynamic propagation along embedded ECM fibers. Our results motivate future biological experiments in mechanobiology to investigate, on the one hand, the mechanosensitivity of cells to dynamic forces traveling and guided by the ECM and, on the other hand, the impact of ECM architecture and remodeling on dynamic force transmission and its spectral filtering, dispersion, and decay.

FITC-Dextran Release from Cell-Embedded Fibrin Hydrogels.

Fibrin hydrogel is a central biological material in tissue engineering and drug delivery applications. As such, fibrin is typically combined with cells and biomolecules targeted to the regenerated tissue. Previous studies have analyzed the release of different molecules from fibrin hydrogels; however, the effect of embedded cells on the release profile has yet to be quantitatively explored. This study focused on the release of Fluorescein isothiocyanate (FITC)-dextran (FD) 250 kDa from fibrin hydrogels, populated with different concentrations of fibroblast or endothelial cells, during a 48-h observation period. The addition of cells to fibrin gels decreased the overall release by a small percentage (by 7-15% for fibroblasts and 6-8% for endothelial cells) relative to acellular gels. The release profile was shown to be modulated by various cellular activities, including gel degradation and physical obstruction to diffusion. Cell-generated forces and matrix deformation (i.e., densification and fiber alignment) were not found to significantly influence the release profiles. This knowledge is expected to improve fibrin integration in tissue engineering and drug delivery applications by enabling predictions and ways to modulate the release profiles of various biomolecules.

Controlled Strain of 3D Hydrogels under Live Microscopy Imaging.

External forces are an important factor in tissue formation, development, and maintenance. The effects of these forces are often studied using specialized in vitro stretching methods. Various available systems use 2D substrate-based stretchers, while the accessibility of 3D techniques to strain soft hydrogels, is more restricted. Here, we describe a method that allows external stretching of soft hydrogels from their circumference, using an elastic silicone strip as the sample carrier. The stretching system utilized in this protocol is constructed from 3D-printed parts and low-cost electronics, making it simple and easy to replicate in other labs. The experimental process begins with polymerizing thick (>100 µm) soft fibrin hydrogels (Elastic Modulus of ~100 Pa) in a cut-out at the center of a silicone strip. Silicone-gel constructs are then attached to the printed-stretching device and placed on the confocal microscope stage. Under live microscopy the stretching device is activated, and the gels are imaged at various stretch magnitudes. Image processing is then used to quantify the resulting gel deformations, demonstrating relatively homogenous strains and fiber alignment throughout the gel's 3D thickness (Z-axis). Advantages of this method include the ability to strain extremely soft hydrogels in 3D while executing in situ microscopy, and the freedom to manipulate the geometry and size of the sample according to the user's needs. Additionally, with proper adaptation, this method can be used to stretch other types of hydrogels (e.g., collagen, polyacrylamide or polyethylene glycol) and can allow for analysis of cells and tissue response to external forces under more biomimetic 3D conditions.

Motion magnification analysis of microscopy videos of biological cells.

It is well recognized that isolated cardiac muscle cells beat in a periodic manner. Recently, evidence indicates that other, non-muscle cells, also perform periodic motions that are either imperceptible under conventional lab microscope lens or practically not easily amenable for analysis of oscillation amplitude, frequency, phase of movement and its direction. Here, we create a real-time video analysis tool to visually magnify and explore sub-micron rhythmic movements performed by biological cells and the induced movements in their surroundings. Using this tool, we suggest that fibroblast cells perform small fluctuating movements with a dominant frequency that is dependent on their surrounding substrate and its stiffness.

Straining 3D Hydrogels with Uniform Z-Axis Strains While Enabling Live Microscopy Imaging.

External forces play an important role in the development and regulation of many tissues. Such effects are often studied using specialized stretchers-standardized commercial and novel laboratory-designed. While designs for 2D stretchers are abundant, the range of available 3D stretcher designs is more limited, especially when live imaging is required. This work presents a novel method and a stretching device that allow straining of 3D hydrogels from their circumference, using a punctured elastic silicone strip as the sample carrier. The system was primarily constructed from 3D-printed parts and low-cost electronics, rendering it simple and cost-efficient to reproduce in other labs. To demonstrate the system functionality, > 100 µm thick soft fibrin gels (< 1 KPa) were stretched, while performing live confocal imaging. The subsequent strains and fiber alignment were analyzed and found to be relatively homogenous throughout the gel's thickness (Z axis). The uniform Z-response enabled by our approach was found to be in contrast to a previously reported approach that utilizes an underlying elastic substrate to convey strain to a 3D thick sample. This work advances the ability to study the role of external forces on biological processes under more physiological 3D conditions, and can contribute to the field of tissue engineering.

Force chains in cell-cell mechanical communication.

Force chains (FCs) are a key determinant of the micromechanical properties and behaviour of heterogeneous materials, such as granular systems. However, less is known about FCs in fibrous materials, such as the networks composing the extracellular matrix (ECM) of biological systems. Using a finite-element computational model, we simulated the contraction of a single cell and two nearby cells embedded in two-dimensional fibrous elastic networks and analysed the tensile FCs that developed in the ECM. The role of ECM nonlinear elasticity on FC formation was evaluated by considering linear and nonlinear, i.e. exhibiting 'buckling' and/or 'strain-stiffening', stress-strain curves. The effect of the degree of cell contraction and network coordination value was assessed. We found that nonlinear elasticity of the ECM fibres influenced the structure of the FCs, facilitating the transition towards more distinct chains that were less branched and more radially oriented than the chains formed in linear elastic networks. When two neighbouring cells contract, a larger number of FCs bridged between the cells in nonlinear networks, and these chains had a larger effective rigidity than the chains that did not reach a neighbouring cell. These results suggest that FCs function as a route for mechanical communication between distant cells and highlight the contribution of ECM fibre nonlinear elasticity to the formation of FCs.

Local Ordering at Mobile Sites in Proteins: Combining Perspectives from NMR Relaxation and Molecular Dynamics.

We report on progress toward improving NMR relaxation analysis in proteins in terms of the slowly relaxing local structure (SRLS) approach by developing a method that combines SRLS with molecular dynamics (MD) simulations. 15N-H bonds from the Rho GTPase binding domain of plexin-B1 are used as test case. We focus on the locally restricting/ordering potential of mean force (POMF), u(θ,φ), at the N-H site (θ and φ specify the orientation of the N-H bond in the protein). In SRLS, u(θ,φ) is expanded in the basis set of the real linear combinations of the Wigner rotation matrix elements with M = 0, D L,| K|(θ,φ). Because of limited data sensitivity, only the lowest ( L = 2) terms are preserved; this potential function is denoted by u(SRLS). In MD, the force-field-parametrized POMF is the potential, u(MD), defined in terms of the probability distribution, Peq(MD) ∝ exp(- u(MD)). Peq(MD), and subsequently u(MD), can be derived from the MD trajectory as histograms. One might contemplate utilizing u(MD) instead of u(SRLS); however, histograms cannot be used in SRLS analyses. Here, we approximate u(θ,φ) in terms of linear combinations of the D L,| K| functions with L = 1-4 and appropriate symmetry, denoted by u(DLK), and optimize the latter (via Peq) against u(MD). This yields for every N-H bond an analytical ordering potential, u(DLK-BEST), which exceeds u(SRLS) considerably in accuracy. u(DLK-BEST) can be used fixed in SRLS data fitting, thereby enabling the determination of additional parameters. This yields a substantially improved picture of structural dynamics, which is a significant benefit. The primary achievement of this work is to have employed for the first time MD data to derive a suitable (in terms of composition and symmetry) approximation to the SRLS POMF.

Nonlinear Elasticity of the ECM Fibers Facilitates Efficient Intercellular Communication.

Biological cells embedded in fibrous matrices have been observed to form intercellular bands of dense and aligned fibers through which they mechanically interact over long distances. Such matrix-mediated cellular interactions have been shown to regulate various biological processes. This study aimed to explore the effects of elastic nonlinearity of the fibers contained in the extracellular matrix (ECM) on the transmission of mechanical loads between contracting cells. Based on our biological experiments, we developed a finite-element model of two contracting cells embedded within a fibrous network. The individual fibers were modeled as showing linear elasticity, compression microbuckling, tension stiffening, or both of the latter two. Fiber compression buckling resulted in smaller loads in the ECM, which were primarily directed toward the neighboring cell. These loads decreased with increasing cell-to-cell distance; when cells were >9 cell diameters apart, no such intercellular interaction was observed. Tension stiffening further contributed to directing the loads toward the neighboring cell, though to a smaller extent. The contraction of two neighboring cells resulted in mutual attraction forces, which were considerably increased by tension stiffening and decayed with increasing cell-to-cell distances. Nonlinear elasticity contributed also to the onset of force polarity on the cell boundaries, manifested by larger contractile forces pointing toward the neighboring cell. The density and alignment of the fibers within the intercellular band were greater when fibers buckled under compression, with tension stiffening further contributing to this structural remodeling. Although previous studies have established the role of the ECM nonlinear mechanical behavior in increasing the range of force transmission, our model demonstrates the contribution of nonlinear elasticity of biological gels to directional and efficient mechanical signal transfer between distant cells, and rehighlights the importance of using fibrous gels in experimental settings for facilitating intercellular communication. VIDEO ABSTRACT.

Conformational Entropy from NMR Relaxation in Proteins: The SRLS Perspective.

Conformational entropy changes associated with bond-vector motions in proteins contribute to the free energy of ligand-binding. To derive such contributions, we apply the slowly relaxing local structure (SRLS) approach to NMR relaxation from 15N-H bonds or C-CDH2 moieties of several proteins in free and ligand-bound form. The spatial restraints on probe motion, which determine the extent of local order, are expressed in SRLS by a well-defined potential, u(θ). The latter yields the orientational probability density, Peq = exp(-u(θ)), and hence the related conformational entropy, S = -∫Peq(θ) ln[Peq(θ)] sin θ dθ (S is "entropy" in units of kBT, and θ represents the bond-vector orientation in the protein). SRLS is applied to 4-oxalocrotonate tautomerase (4-OT), the acyl-coenzyme A binding protein (ACBP), the C-terminal SH2 domain of phospholipase Cγ1 (PLCγ1C SH2), the construct dihydrofolate reductase-E:folate (DHFR-E:folate), and their complexes with appropriate ligands, to determine ΔS. Eglin C and its V18A and V34A mutants are also studied. Finally, SRLS is applied to the structurally homologous proteins TNfn3 and FNfn10 to characterize within its scope the unusual "dynamics" of the TNfn3 core. Upon ligand-binding, the backbones of 4-OT, ACBP, and PLCγ1C SH2 show limited, increased, and decreased order, respectively; the cores of DHFR-E:folate and PLCγ1C SH2 become more ordered. The V18A (V34A) mutation increases (decreases) the order within the eglin C core. The core of TNfn3 is less ordered structurally and more mobile kinetically. Secondary structure versus loops, surface-binding versus core insertion, and ligand size emerged as being important in rationalizing ΔS. The consistent and general tool developed herein is expected to provide further insights in future work.

An SRLS Study of 2H Methyl-Moiety Relaxation and Related Conformational Entropy in Free and Peptide-Bound PLCγ1C SH2.

The two-body (protein and probe) coupled-rotator slowly relaxing local structure (SRLS) approach for NMR relaxation in proteins is extended to derive conformational entropy, S. This version of SRLS is applied to deuterium relaxation from the C-CDH2 bonds of free and peptide-bound PLCγ1C SH2. Local C-CDH2 motion is described by a correlation time for local diffusion, τ2, and a Maier-Saupe potential, u. On average, τ2, which largely fulfills τ2 ≪ τ1 (τ1 - correlation time for global tumbling), is 270 ± 41 ps and u is 2 ± 0.1 kBT. The PLCγ1C SH2 data were analyzed previously with the model-free (MF) method. SRLS is a generalization of MF, assumed so far to yield the latter for τ2 ≪ τ1 and simple local geometry. Despite these conditions being fulfilled, we find here that τ2 and u differ substantially from their MF counterparts. This is shown to stem from MF (a) disregarding mode-coupling of the first type (see below) and (b) parametrizing the methyl-moiety-related spectral density function (SDF). Our main interest lies in ΔS, the conformational entropy difference between the peptide-bound and free PLCγ1C SH2 forms. We find that ΔS is rendered inaccurate in MF because factors a and b above impair the accuracy of Saxis, the parameter on which the calculation of ΔS is based. Conformational entropy was obtained previously using various simple system-specific models. SRLS is unique in obtaining this important thermodynamic quantity based on a general physically well-defined local potential. It is also unique in its ability to extract the information inherent in 2H relaxation parameters from methyl moieties in protein with accuracy commensurate with data sensitivity.

Local Ordering at Mobile Sites in Proteins from Nuclear Magnetic Resonance Relaxation: The Role of Site Symmetry.

Restricted motions in proteins (e.g., N-H bond dynamics) are studied effectively with NMR. By analogy with restricted motions in liquid crystals (LC), the local ordering has in the past been primarily represented by potentials comprising the L = 2, |K| = 0, 2 spherical harmonics. However, probes dissolved in LCs experience nonpolar ordering, often referred to as alignment, while protein-anchored probes experience polar ordering, often referred to as orientation. In this study we investigate the role of local (site) symmetry in the context of the polarity of the local ordering. We find that potentials comprising the L = 1, |K| = 0, 1 spherical harmonics represent adequately polar ordering. It is useful to characterize potential symmetry in terms of the irreducible representations of D2h point group, which is already implicit in the definition of the rotational diffusion tensor. Thus, the relevant rhombic L = 1 potentials have B1u and B3u symmetry whereas the relevant rhombic L = 2 potentials have Ag symmetry. A comprehensive scheme where local potentials and corresponding probability density functions (PDFs) are represented in Cartesian and spherical coordinates clarifies how they are affected by polar and nonpolar ordering. The Cartesian coordinates are chosen so that the principal axis of polar axial PDF is pointing along the z-axis, whereas the principal axis of the nonpolar axial PDF is pointing along ±z. Two-term axial potentials with 1 ≤ L ≤ 3 exhibit substantial diversity; they are expected to be useful in NMR-relaxation-data-fitting. It is shown how potential coefficients are reflected in the experimental order parameters. The comprehensive scheme representing local potentials and PDFs is exemplified for the L = 2 case using experimental data from (15)N-labeled plexin-B1 and thioredoxin, (2)H-, and (13)C-labeled benzenehexa-n-alkanoates, and nitroxide-labeled T4 lysozyme. Future prospects for improved ordering analysis based on combined atomistic and mesoscopic approaches are delineated.

Polar Versus Non-polar Local Ordering at Mobile Sites in Proteins: Slowly Relaxing Local Structure Analysis of (15)N Relaxation in the Third Immunoglobulin-Binding Domain of Streptococcal Protein G.

We developed recently the slowly relaxing local structure (SRLS) approach for studying restricted motions in proteins by NMR. The spatial restrictions have been described by potentials comprising the traditional L = 2, K = 0, 2 spherical harmonics. However, the latter are associated with non-polar ordering whereas protein-anchored probes experience polar ordering, described by odd-L spherical harmonics. Here we extend the SRLS potential to include the L = 1, K = 0, 1 spherical harmonics and analyze (15)N-(1)H relaxation from the third immunoglobulin-binding domain of streptococcal protein G (GB3) with the polar L = 1 potential (coefficients c0(1) and c1(1)) or the non-polar L = 2 potential (coefficients c0(2) and c2(2)). Strong potentials, with ⟨c0(1)⟩ ∼ 60 for L = 1 and ⟨c0(2)⟩ ∼ 20 for L = 2 (in units of kBT), are detected. In the α-helix of GB3 the coefficients of the rhombic terms are c1(1) ∼ c2(2) ∼ 0; in the preceding (following) chain segment they are ⟨c1(1)⟩ ∼ 6 for L = 1 and ⟨c2(2)⟩ ∼ 14 for L = 2 (⟨c1(1)⟩ ∼ 3 for L = 1 and ⟨c2(2)⟩ ∼ 7 for L = 2). The local diffusion rate, D2, lies in the 5 × 10(9)-1 × 10(11) s(-1) range; it is generally larger for L = 1. The main ordering axis deviates moderately from the N-H bond. Corresponding L = 1 and L = 2 potentials and probability density functions are illustrated for residues A26 of the α-helix, Y3 of the ß1-strand, and L12 of the ß1/ß2 loop; they differ considerably. Polar/orientational ordering is shown to be associated with GB3 binding to its cognate Fab fragment. The polarity of the local ordering is clearly an important factor.

Structural Dynamics of the Potassium Channel Blocker ShK: SRLS Analysis of (15)N Relaxation.

The 35-residue ShK peptide binds with high affinity to voltage-gated potassium channels. The dynamics of the binding surface was studied recently with (microsecond to millisecond) (15)N relaxation dispersion and (picosecond to nanosecond) (15)N spin relaxation of the N-H bonds. Relaxation dispersion revealed microsecond conformational-exchange-mediated exposure of the functionally important Y23 side chain to the peptide surface. The spin relaxation parameters acquired at 14.1 and 16.45 T have been subjected to model-free (MF) analysis, which yielded a squared generalized order parameter, S(2), of approximately 0.85 for virtually all of the N-H bonds. Only a "rigid backbone" evaluation could be inferred. We ascribe this limited information to the simplicity of MF in the context of challenging data. To improve the analysis, we apply the slowly relaxing local structure (SRLS) approach, which is a generalization of MF. SRLS describes N-H bond dynamics in ShK in terms of a local potential, u, ranging from 10 to 18.5 kBT, and a local diffusion rate, D2, ranging from 4.2 × 10(8) to 2.4 × 10(10) s(-1). This analysis shows that u is outstandingly strong for Y23 and relatively weak for K22, whereas D2 is slow for Y23 and fast for K22. These observations are relevant functionally because of the key role of the K22-Y23 dyad in ShK binding to potassium channels. The disulfide-bond network exhibits a medium-strength potential and an alternating wave-like D2 pattern. This is indicative of moderate structural restraints and motional plasticity, in support of, although not directly correlated with, the microsecond binding-related conformational exchange process detected previously. Thus, new information on functionally important residues in ShK and its overall conformational stability emerged from the SRLS analysis, as compared with the previous MF-based estimate of backbone dynamics as backbone rigidity.

An extended CCR5 ECL2 peptide forms a helix that binds HIV-1 gp120 through non-specific hydrophobic interactions.

UNLABELLED: C-C chemokine receptor 5 (CCR5) serves as a co-receptor for HIV-1. The CCR5 N-terminal segment, the second extracellular loop (ECL2) and the transmembrane helices have been implicated in binding the envelope glycoprotein gp120. Peptides corresponding to the sequence of the putative ECL2 as well as peptides containing extracellular loops 1 and 3 (ECL1 and ECL3) were found to inhibit HIV-1 infection. The aromatic residues in the C-terminal half of an ECL2 peptide were shown to interact with gp120. In the present study, we found that, in aqueous buffer, the segment Q188-Q194 in an elongated ECL2 peptide (R168-K197) forms an amphiphilic helix, which corresponds to the beginning of the fifth transmembrane helix in the crystal structure of CCR5. Two-dimensional saturation transfer difference NMR spectroscopy and dynamic filtering studies revealed involvement of Y187, F189, W190 and F193 of the helical segment in the interaction with gp120. The crystal structure of CCR5 shows that the aromatic side chains of F189, W190 and F193 point away from the binding pocket and interact with the membrane or with an adjacent CCR5 molecule, and therefore could not interact with gp120 in the intact CCR5 receptor. We conclude that these three aromatic residues of ECL2 peptides interact with gp120 through hydrophobic interactions that are not representative of the interactions of the intact CCR5 receptor. The HIV-1 inhibition by ECL2 peptides, as well as by ECL1 and ECL3 peptides and peptides corresponding to ECL2 of CXCR4, which serves as an alternative HIV-1 co-receptor, suggests that there is a hydrophobic surface in the envelope spike that could be a target for HIV-1 entry inhibitors. DATABASE: The structures and NMR data of ECL2S (Q186-T195) were deposited under Protein Data Bank ID 2mzx and BioMagResBank ID 25505.

Drosophila Raf's N terminus contains a novel conserved region and can contribute to torso RTK signaling.

Drosophila Raf (DRaf) contains an extended N terminus, in addition to three conserved regions (CR1-CR3); however, the function(s) of this N-terminal segment remains elusive. In this article, a novel region within Draf's N terminus that is conserved in BRaf proteins of vertebrates was identified and termed conserved region N-terminal (CRN). We show that the N-terminal segment can play a positive role(s) in the Torso receptor tyrosine kinase pathway in vivo, and its contribution to signaling appears to be dependent on the activity of Torso receptor, suggesting this N-terminal segment can function in signal transmission. Circular dichroism analysis indicates that DRaf's N terminus (amino acids 1-117) including CRN (amino acids 19-77) is folded in vitro and has a high content of helical secondary structure as predicted by proteomics tools. In yeast two-hybrid assays, stronger interactions between DRaf's Ras binding domain (RBD) and the small GTPase Ras1, as well as Rap1, were observed when CRN and RBD sequences were linked. Together, our studies suggest that DRaf's extended N terminus may assist in its association with the upstream activators (Ras1 and Rap1) through a CRN-mediated mechanism(s) in vivo.

Regulation of energy liberation during steady sarcomere shortening.

Energy liberation rate (E) during steady muscle shortening is a monotonic increasing or biphasic function of the shortening velocity (V). The study examines three plausible hypotheses for explaining the biphasic E-V relationship (EVR): 1) the cross-bridge (XB) turnover rate from non-force-generating (weak) to force-generating (strong) conformation decreases as V increases; 2) XB kinetics is determined by the number of strong XBs (XB-XB cooperativity); and 3) the affinity of troponin for calcium is modulated by the number of strong XBs (XB-Ca cooperativity). The relative role of the various energy-regulating mechanisms is not well defined. The hypotheses were tested by coupling calcium kinetics with XB cycling. All three hypotheses yield identical steady-state characteristics: 1) hyperbolic force-velocity relationship; 2) quasi-linear stiffness-force relationship; and 3) biphasic EVR, where E declines at high V due to decrease in the number of cycling XBs or in the weak-to-strong transition rate. The hypotheses differ in the ability to describe the existence of both monotonic and biphasic EVRs and in the effect of intracellular free calcium concentration ([Ca2+]i) on the EVR peak. Monotonic and biphasic EVRs with a shift in EVR peak to higher velocity at higher [Ca2+]i are obtained only by XB-Ca cooperativity. XB-XB cooperativity provides only biphasic EVRs. A direct effect of V on XB kinetics predicts that EVR peak is obtained at the same velocity independently of [Ca2+]i. The study predicts that measuring the dependence of the EVR on [Ca2+]i allows us to test the hypotheses and to identify the dominant energy-regulating mechanism. The established XB-XB and XB-Ca mechanisms provide alternative explanations to the various reported EVRs.

Is peptide bond cis/trans isomerization a key stage in the chemo-mechanical cycle of motor proteins?

Motor proteins such as myosin and kinesin are responsible for actively directed movement in vivo. The physicochemical mechanism underlying their function is still obscure. A novel and unifying model concerning the motors driving mechanism is suggested here. This model resides within the framework of the well-studied "swinging lever-arm" hypothesis, stating that cis/trans peptide bond isomerization (CTI) is a key stage in the chemo-mechanical coupling within actomyosin--the complex of the motor (myosin) and its specific track (actin). CTI is suggested to propel myosin's lever-arm swing. The model addresses on the submolecular level a broad spectrum of actomyosin's functional characteristics, such as kinetics, energetics, force exertion, stepping, and directionality. The model may be tested first with relative ease in kinesin--a smaller motor that could be specifically modified with unnatural amino acids using bacterial expression. Suggested modifications may be used for labeling and functional decoupling.