Stretching the α-helix: a direct measure of the hydrogen-bond energy of a single-peptide molecule

Peptide-Crosslinked, Highly Entangled Hydrogels with Excellent Mechanical Properties but Ultra-Low Solid Content

Low solid content is the ultimate reason for the brittleness and weakness of ordinary hydrogels. Here, hydrogels with ultra-low solid content but good mechanical properties are successfully synthesized using high monomer concentrations and low cross-linker/monomer ratios to obtain highly entangled structure and poly(l-lysine)-based cross-linker to introduce peptide chains. Compared with hydrogel cross-linked with N,N'-methylenebisacrylamide (BIS), the peptide-crosslinked one has a larger swelling degree in water, leading to fully swollen gel with ultra-low solid content (5.8%). However, it still exhibits excellent mechanical properties, including high stretchability (440%), high tensile strength (220 KPa), superb resilience (99%), high fracture toughness (2100 J m^-2 ), excellent fatigue resistance (720 J m^-2 ), low friction (0.0059), and high wear resistance. These properties are comparable to or even better than the BIS-crosslinked hydrogel, although the former has much lower solid content. The excellent mechanical properties of the peptide-crosslinked gel are attributed to its highly entangled structure and also to the introduction of a novel mechanism for energy dissipation, that is, energy dissipation via breakage of intramolecular hydrogen bonds stabilizing the helical structure of the peptide.

Structural role of osteocalcin and its modification in bone fracture

Osteocalcin (OC), an abundant non-collagenous protein in bone extracellular matrix, plays a vital role in both its biological and mechanical function. OC undergoes post-translational modification, such as glycation; however, it remains unknown whether glycation of OC affects bone's resistance to fracture. Here, for the first time, we demonstrate the formation of pentosidine, an advanced glycation end-product (AGE) cross-link on mouse OC analyzed by ultra-performance liquid chromatography. Next, we establish that the presence of OC in mouse bone matrix is associated with lower interlamellar separation (distance) and thicker bridges spanning the lamellae, both of which are critical for maintaining bone's structural integrity. Furthermore, to determine the impact of modification of OC by glycation on bone toughness, we glycated bone samples in vitro from wild-type (WT) and osteocalcin deficient (Oc-/-) mice, and compared the differences in total fluorescent AGEs and fracture toughness between the Oc -/- glycated and control mouse bones and the WT glycated and control mouse bones. We determined that glycation resulted in significantly higher AGEs in WT compared to Oc-/- mouse bones (delta-WT > delta-OC, p = 0.025). This observed change corresponded to a significant decrease in fracture toughness between WT and Oc-/- mice (delta-WT vs delta-OC, p = 0.018). Thus, we propose a molecular deformation and fracture mechanics model that corroborates our experimental findings and provides evidence to support a 37%-90% loss in energy dissipation of OC due to formation of pentosidine cross-link by glycation. We anticipate that our study will aid in elucidating the effects of a major non-collagenous bone matrix protein, osteocalcin, and its modifications on bone fragility and help identify potential therapeutic targets for maintaining skeletal health.

Peptide-enhanced tough, resilient and adhesive eutectogels for highly reliable strain/pressure sensing under extreme conditions

Natural gels and biomimetic hydrogel materials have been able to achieve outstanding integrated mechanical properties due to the gain of natural biological structures. However, nearly every natural biological structure relies on water as solvents or carriers, which limits the possibility in extreme conditions, such as sub-zero temperatures and long-term application. Here, peptide-enhanced eutectic gels were synthesized by introducing α-helical "molecular spring" structure into deep eutectic solvent. The gel takes full advantage of the α-helical structure, achieving high tensile/compression, good resilience, superior fracture toughness, excellent fatigue resistance and strong adhesion, while it also inherits the benefits of the deep eutectic solvent and solves the problems of solvent volatilization and freezing. This enables unprecedentedly long and stable sensing of human motion or mechanical movement. The electrical signal shows almost no drift even after 10,000 deformations for 29 hours or in the -20 °C to 80 °C temperature range.

Contribution of Hydrophobic Interactions to Protein Mechanical Stability

The role of hydrophobic and polar interactions in providing thermodynamic stability to folded proteins has been intensively studied, but the relative contribution of these interactions to the mechanical stability is less explored. We used steered molecular dynamics simulations with constant-velocity pulling to generate force-extension curves of selected protein domains and monitor hydrophobic surface unravelling upon extension. Hydrophobic contribution was found to vary between one fifth and one third of the total force while the rest of the contribution is attributed primarily to hydrogen bonds. Moreover, hydrophobic force peaks were shifted towards larger protein extensions with respect to the force peaks attributed to hydrogen bonds. The higher importance of hydrogen bonds compared to hydrophobic interactions in providing mechanical resistance is in contrast with the relative importance of the hydrophobic interactions in providing thermodynamic stability of proteins. The different contributions of these interactions to the mechanical stability are explained by the steeper free energy dependence of hydrogen bonds compared to hydrophobic interactions on the relative positions of interacting atoms. Comparative analyses for several protein domains revealed that the variation of hydrophobic forces is modest, while the contribution of hydrogen bonds to the force peaks becomes increasingly important for mechanically resistant protein domains.

Thermal fluctuations assist mechanical signal propagation in coiled-coil proteins

Recently, it has been shown that the long coiled-coil membrane tether protein early endosome antigen 1 (EEA1) switches from a rigid to a flexible conformation upon binding of a signaling protein to its free end. This flexibility switch represents a motorlike activity, allowing EEA1 to generate a force that moves vesicles closer to the membrane they will fuse with. It was hypothesized that the binding-induced signal could propagate along the coiled coil and lead to conformational changes through the localized domains of the protein chain that deviate from a perfect coiled-coil structure. To elucidate, if upon binding of a single protein the corresponding mechanical signal could propagate through the whole 200-nm-long chain, we propose a simplified description of the coiled coil as a one-dimensional Frenkel-Kontorova chain. Using numerical simulations, we find that an initial perturbation of the chain can propagate along its whole length in the presence of thermal fluctuations. This may enable the change of the configuration of the entire molecule and thereby affect its stiffness. Our work sheds light on intramolecular communication and force generation in long coiled-coil proteins.

Unfolding mechanism and free energy landscape of single, stable, alpha helices at low pull speeds

Single alpha helices (SAHs) stable in isolated form are often found in motor proteins where they bridge functional domains. Understanding the mechanical response of SAHs is thus critical to understand their function. The quasi-static force-extension relation of a small number of SAHs is known from single-molecule experiments. Unknown, or still controversial, are the molecular scale details behind those observations. We show that the deformation mechanism of SAHs pulled from the termini at pull speeds approaching the quasi-static limit differs from that of typical helices found in proteins, which are stable only when interacting with other protein domains. Using molecular dynamics simulations with atomistic resolution at low pull speeds previously inaccessible to simulation, we show that SAHs start unfolding from the termini at all pull speeds we investigated. Unfolding proceeds residue-by-residue and hydrogen bond breaking is not the main event determining the barrier to unfolding. We use the molecular simulation data to test the cooperative sticky chain model. This model yields excellent fits of the force-extension curves and quantifies the distance, xE = 0.13 nm, to the transition state, the natural frequency of bond vibration, ν0 = 0.82 ns-1, and the height, V0 = 2.9 kcal mol-1, of the free energy barrier associated with the deformation of single residues. Our results demonstrate that the sticky chain model could advantageously be used to analyze experimental force-extension curves of SAHs and other biopolymers.

Structural role of osteocalcin and osteopontin in energy dissipation in bone

Non-collagenous proteins are a vital component of bone matrix. Amongst them, osteocalcin (OC) and osteopontin (OPN) hold special significance due to their intimate interaction with the mineral and collagenous matrix in bone. Both proteins have been associated with microdamage and fracture, but their structural role in energy dissipation is unclear. This study used bone tissue from genetic deficient mice lacking OC and/or OPN and subjected them to a series of creep-fatigue-creep tests. To this end, whole tibiae were loaded in four-point bending to 70% stiffness loss which captured the three characteristic phases of fatigue associated with initiation, propagation, and coalescence of microdamage. Fatigue loading preceded and followed creep tests to determine creep and dampening parameters. Microdamage in the form of linear microcracks and diffuse damage were analyzed by histology. It was shown that OC and OPN were 'activated' following stiffness loss associated with fatigue damage where they facilitated creep and dampening parameters (i.e. increased energy dissipation). More specifically, post-fatigue creep rate and dampening were significantly greater in wild-types (WTs) than genetic deficient mice (p < 0.05). These results were supported by microdamage analysis which showed significant increase in creep-associated diffuse damage formation in WTs compared to genetic deficient groups (p < 0.05). Based on these findings, we propose that during local yield events, OC and OPN rely on ionic interactions of their charged side chains and on hydrogen bonding to dissipate energy in bone.

Quantitative Evaluation of Viral Protein Binding to Phosphoinositide Receptors and Pharmacological Inhibition

There is significant interest in developing analytical methods to characterize molecular recognition events between proteins and phosphoinositides, which are a medically important class of carbohydrate-functionalized lipids. Within this scope, one area of high priority involves quantitatively evaluating drug candidates that pharmacologically inhibit protein-phosphoinositide interactions. As full-length proteins are often difficult to produce, establishing methods to study these interactions with shorter, bioactive peptides would be advantageous. Herein, we report an atomic force microscopy (AFM)-based force spectroscopic approach to detect the specific interaction between an amphipathic, α-helical (AH) peptide derived from the hepatitis C virus NS5A protein and its biological target, the phosphatidylinositol (4,5)-bisphosphate [PI(4,5)P₂] phosphoinositide receptor. After optimization of the peptide tethering strategy and measurement parameters, the binding specificity of AH peptide for PI(4,5)P₂ receptors was comparatively evaluated across a panel of phosphoinositides and the influence of ionic strength on AH-PI(4,5)P₂ binding strength was tested. Importantly, these capabilities were translated into the development of a novel experimental methodology to determine the inhibitory activity of a small-molecule drug candidate acting against the AH-PI(4,5)P₂ interaction, and extracted kinetic parameters agree well with literature values obtained by conventional biochemical methods. Taken together, our findings provide a nanomechanical basis for explaining the high binding specificity of the NS5A AH to PI(4,5)P₂ receptors, in turn establishing an analytical framework to study phosphoinositide-binding viral peptides and proteins as well as a broadly applicable approach to evaluate candidate inhibitors of protein-phosphoinositide interactions.

Molecular basis of the mechanical hierarchy in myomesin dimers for sarcomere integrity

Myomesin is one of the most important structural molecules constructing the M-band in the force-generating unit of striated muscle, and a critical structural maintainer of the sarcomere. Using molecular dynamics simulations, we here dissect the mechanical properties of the structurally known building blocks of myomesin, namely α-helices, immunoglobulin (Ig) domains, and the dimer interface at myomesin's 13th Ig domain, covering the mechanically important C-terminal part of the molecule. We find the interdomain α-helices to be stabilized by the hydrophobic interface formed between the N-terminal half of these helices and adjacent Ig domains, and, interestingly, to show a rapid unfolding and refolding equilibrium especially under low axial forces up to ∼ 15 pN. These results support and yield atomic details for the notion of recent atomic-force microscopy experiments, namely, that the unique helices inserted between Ig domains in myomesin function as elastomers and force buffers. Our results also explain how the C-terminal dimer of two myomesin molecules is mechanically outperforming the helices and Ig domains in myomesin and elsewhere, explaining former experimental findings. This study provides a fresh view onto how myomesin integrates elastic helices, rigid immunoglobulin domains, and an extraordinarily resistant dimer into a molecular structure, to feature a mechanical hierarchy that represents a firm and yet extensible molecular anchor to guard the stability of the sarcomere.

The nanocrystal superlattice pressure cell: a novel approach to study molecular bundles under uniaxial compression

Ordered assemblies of inorganic nanocrystals coated with organic linkers present interesting scientific challenges in hard and soft matter physics. We demonstrate that a nanocrystal superlattice under compression serves as a nanoscopic pressure cell to enable studies of molecular linkers under uniaxial compression. We developed a method to uniaxially compress the bifunctional organic linker by attaching both ends of aliphatic chains to neighboring PbS nanocrystals in a superlattice. Pressurizing the nanocrystal superlattice in a diamond anvil cell thus results in compression of the molecular linkers along their chain direction. Small-angle and wide-angle X-ray scattering during the compression provide insights into the structure of the superlattice and nanocrystal cores under compression, respectively. We compare density functional theory calculations of the molecular linkers as basic Hookean springs to the experimental force-distance relationship. We determine the density of linkers on the nanocrystal surfaces. We demonstrate our method to probe the elastic force of single molecule as a function of chain length. The methodology introduced in this paper opens doors to investigate molecular interactions within organic molecules compressed within a nanocrystal superlattice.

Persistence length and stochastic fragmentation of supramolecular nanotubes under mechanical force

Cyclic peptide nanotubes (CPNs) exhibit impressive structural, mechanical and chemical properties in resemblance to beta-sheet proteins found in silks and amyloids, and show potential as supramolecular nanotubes that can be utilized to generate novel nanocomposites and nanoporous thin films. Quantifying the persistence length and thermomechanical fragmentation of CPNs is of great importance for establishing a theoretical basis of how to generate rectilinear nanostructures with controlled aspect ratio and rigidity. However, factors governing the elasticity and dynamical breaking of these supramolecular nanostructures remain to be fully understood. Here we present a statistical analysis of the Young's modulus and persistence length of CPNs using fully-atomistic molecular dynamic simulations in explicit solvent. We show that the measured properties exhibit a dependence on the magnitude of the shear force applied, and extrapolation to the quasi-static deformation case yields 0.46 μm for the persistence length and 7.8 GPa for the Young's modulus, in agreement with our experimental observations from TEM images. We establish a theoretical model for the spatial and temporal distribution of stochastic fracture, which we use to explain the simulation-based observations of spontaneous fragmentation under an applied shear force. Our methodology, blending theory, simulation and experiments provide a framework that can be utilized to investigate the mechanical behavior of self-assembling protein materials, paving the way for their design towards biological and industrial applications.

Combined Poisson and soft-particle DLVO analysis of the specific and nonspecific adhesion forces measured between L. monocytogenes grown at various temperatures and silicon nitride

Adhesion forces between pathogenic L. monocytogenes EGDe and silicon nitride (Si(3)N(4)) were measured using atomic force microscopy (AFM) under water and at room temperature for cells grown at five different temperatures (10, 20, 30, 37, and 40 °C). Adhesion forces were then decoupled into specific (hydrogen bonding) and nonspecific (electrostatic and Lifshitz-van der Waals) force components using Poisson statistical analysis. The strongest specific and nonspecific attraction forces were observed for cells grown at 30 °C, compared to those observed for cells grown at higher or lower temperatures, respectively. By combining the results of Poisson analysis with the results obtained through soft-particle Derjaguin-Landau-Verwey-Overbeek (DLVO) analysis, the contributions of the Lifshitz-van der Waals and electrostatic forces to the overall nonspecific interaction forces were determined. Our results showed that the Lifshitz-van der Waals attraction forces dominated the total nonspecific adhesion forces for all investigated thermal conditions. However, irrespective of the temperature of growth investigated, hydrogen bonding forces were always stronger than the nonspecific forces. Finally, by combining Poisson analysis with soft-particle analysis of DLVO forces, the closest separation distances where the irreversible bacterial adhesion takes place can be determined relatively easily. For all investigated thermal conditions, the closest separation distances were <1 nm.

Gas-liquid chromatography-obtained differences in the dissolution enthalpy between two positional isomers in a polar stationary phase: a measure of the inter- or intramolecular hydrogen bond energy?

Previous GLC work with several 2- and 4-substituted phenols and anilines, as well as with a pyrrolizidine alcohol, had determined the difference between the heats of dissolution of two positional isomers in a strong polar stationary liquid phase; one of these isomers forms an intramolecular hydrogen bond (intra-HB) and the other has no such bond for steric reasons. The energies of the intermolecular hydrogen bonds (inter-HBs), ΔH(inter-HB), formed by the 1,2- and 1,4-isomers with the molecules of a polar phase had been assumed approximately equal, so the difference between them could be ignored. The same assumption had been made for the energies of nonspecific interactions (NSIs), ΔH(NSI). It had been concluded that the found difference can be considered as an intra-HB energy (enthalpy), ΔH(intra-HB), when the energies (enthalpies) of inter-HBs formed by the 1,2- and 1,4-isomers under study with the molecules of a polar phase are much greater in absolute value than ΔH(intra-HB). And, conversely, when |ΔH(intra-HB)|>|ΔH(inter-HB)|, an inter-HB enthalpy will result. With the same assumptions, we here obtained an extended thermodynamic equation and corrected this above conclusion on the basis of a general consideration of the dissolution thermodynamics for two isomers of a molecule in a polar phase. Account was taken of the coefficients of isomer partitioning between the liquid and the gaseous phase at the experimental temperature. The conclusion made previously was adjusted for ΔH(NSI) and formulated as follows: The GLC method determines the intra-HB energy at |ΔH(intra-HB)|≤|ΔH(inter-HB)+ΔH(NSI)|. If |ΔH(intra-HB)|>|ΔH(inter-HB)+ΔH(NSI)|, the method yields the values of ΔH(inter-HB)+ΔH(NSI). This new conclusion was illustrated with virtual (numerical) experiments in which various ΔH(intra-HB), ΔH(inter-HB), and ΔH(NSI) values were postulated and results were obtained that would have been achieved by GLC if it had been done. Using a capillary column with the PEG 20M stationary phase, we measured the differences between the heats of dissolution for six pairs of isomers of phenolic compounds and for seven pairs of disubstituted benzene derivatives, which have similar structures but cannot form an intra-HB. The benzene derivatives served to make an approximate experimental estimate of the difference between the energies of NSIs of the isomers under study with a polar phase, ΔH(NSI).

Electrostatic fields near the active site of human aldose reductase: 2. New inhibitors and complications caused by hydrogen bonds

Vibrational Stark effect spectroscopy was used to measure electrostatic fields in the hydrophobic region of the active site of human aldose reductase (hALR2). A new nitrile-containing inhibitor was designed and synthesized, and the X-ray structure of its complex, along with cofactor NADP(+), with wild-type hALR2 was determined at 1.3 Å resolution. The nitrile is found to be in the proximity of T113, consistent with a hydrogen bond interaction. Two vibrational absorption peaks were observed at room temperature in the nitrile region when the inhibitor binds to wild-type hALR2, indicating that the nitrile probe experiences two different microenvironments, and these could be empirically separated into a hydrogen-bonded and non-hydrogen-bonded population by comparison with the T113A mutant, in which a hydrogen bond to the nitrile is not present. Classical molecular dynamics simulations based on the structure predict a double-peak distribution in protein electric fields projected along the nitrile probe. The interpretation of these two peaks as a hydrogen bond formation-dissociation process between the probe nitrile group and a nearby amino acid side chain is used to explain the observation of two IR bands, and the simulations were used to investigate the molecular details of this conformational change. Hydrogen bonding complicates the simplest analysis of vibrational frequency shifts as being due solely to electrostatic interactions through the vibrational Stark effect, and the consequences of this complication are discussed.

How protein materials balance strength, robustness, and adaptability

Proteins form the basis of a wide range of biological materials such as hair, skin, bone, spider silk, or cells, which play an important role in providing key functions to biological systems. The focus of this article is to discuss how protein materials are capable of balancing multiple, seemingly incompatible properties such as strength, robustness, and adaptability. To illustrate this, we review bottom-up materiomics studies focused on the mechanical behavior of protein materials at multiple scales, from nano to macro. We focus on alpha-helix based intermediate filament proteins as a model system to explain why the utilization of hierarchical structural features is vital to their ability to combine strength, robustness, and adaptability. Experimental studies demonstrating the activation of angiogenesis, the growth of new blood vessels, are presented as an example of how adaptability of structure in biological tissue is achieved through changes in gene expression that result in an altered material structure. We analyze the concepts in light of the universality and diversity of the structural makeup of protein materials and discuss the findings in the context of potential fundamental evolutionary principles that control their nanoscale structure. We conclude with a discussion of multiscale science in biology and de novo materials design.

Energy landscape, structure and rate effects on strength properties of alpha-helical proteins

The strength of protein domains is crucial to identify the mechanical role of protein domains in biological processes such as mechanotransduction, tissue mechanics and tissue remodeling. Whereas the concept of strength has been widely investigated for engineered materials, the strength of fundamental protein material building blocks and how it depends on structural parameters such as the chemical bonding, the protein filament length and the timescale of observation or deformation velocity remains poorly understood. Here we report a systematic analysis of the influence of key parameters that define the energy landscape of the strength properties of alpha-helical protein domains, including energy barriers, unfolding and refolding distances, the locations of folded and unfolded states, as well as variations of the length and pulling velocity of alpha-helical protein filaments. The analysis is facilitated by the development of a double-well mesoscale potential formulation, utilized here to carry out a systematic numerical analysis of the behavior of alpha-helices. We compare the results against widely used protein strength models based on the Bell model, one of the simplest models used to characterize the strength of protein filaments. We find that, whereas Bell-type models are a reasonable approximation to describe the rupture of alpha-helical protein domains for a certain range of pulling speeds and values of energy barriers, the model ceases to hold for very large energy barriers and for very small pulling speeds, in agreement with earlier findings. We conclude with an application of our mesoscale model to investigate the effect of the length of alpha-helices on their mechanical strength. We find a weakening effect as the length of alpha-helical proteins increases, followed by an asymptotic regime in which the strength remains constant. We compare strand lengths found in biological proteins with the scaling law of strength versus alpha-helix filament length. The mesoscale model reported here is generally applicable to other protein filaments that feature a serial array of domains that unfold under applied strain, where a similar length-dependent strength could be observed.

Free energies and forces in helix-coil transition of homopolypeptides under stretching

We show here that constant velocity steered molecular dynamics (SMD) simulations of alpha-helices in a vacuum present a well defined plateau in the force-extension relationship for homopolypeptides having more than (approximately) twenty residues. With the processes being far away from equilibrium, the energies strongly depend on the stretching velocity. Importantly, for a given velocity variation, the energy variation depends also on the helix sequence. Additionally, our observations show that homopolypeptides made of ten different amino acids (Ala, Cys, Gln, Ile, Leu, Met, Phe, Ser, Thr and Val) present a linear helix-coil transition.

Tensile mechanics of alanine-based helical polypeptide: force spectroscopy versus computer simulations

In nature, an alpha-helix is commonly used to build thermodynamically stable and mechanically rigid protein conformations. In view of growing interest in the mechanical rigidity of proteins, we measured the tensile profile of an alanine-based alpha-helical polypeptide on an atomic-force microscope to investigate the basic mechanics of helix extension with minimal interference from side-chain interactions. The peptide was extended to its maximum contour length with much less force than in reported cases of poly-L-Glu or poly-L-Lys, indicating that chain stiffness strongly depended on the physicochemical properties of side chains, such as their bulkiness. The low tensile-force extension originated presumably in locally unfolded parts because of spontaneous structural fluctuations. In 50% trifluoroethanol, the well-known helix-promoting agent, the rigidity of the sample polypeptide was markedly increased. Computer simulations of the peptide-stretching process showed that a majority of constituent residues underwent a transition from an alpha-helical to an extended conformation by overcoming an energy barrier around psi approximately 0 degrees on the Ramachandran plot. The observed lability of an isolated helix signified the biological importance of the lateral bundling of helices to maintain a rigid protein structure.

Deformation and failure of protein materials in physiologically extreme conditions and disease

Biological protein materials feature hierarchical structures that make up a diverse range of physiological materials. The analysis of protein materials is an emerging field that uses the relationships between biological structures, processes and properties to probe deformation and failure phenomena at the molecular and microscopic level. Here we discuss how advanced experimental, computational and theoretical methods can be used to assess structure-process-property relations and to monitor and predict mechanisms associated with failure of protein materials. Case studies are presented to examine failure phenomena in the progression of disease. From this materials science perspective, a de novo basis for understanding biological processes can be used to develop new approaches for treating medical disorders. We highlight opportunities to use knowledge gained from the integration of multiple scales with physical, biological and chemical concepts for potential applications in materials design and nanotechnology.

Nanomechanical strength mechanisms of hierarchical biological materials and tissues

Biological protein materials (BPMs), intriguing hierarchical structures formed by assembly of chemical building blocks, are crucial for critical functions of life. The structural details of BPMs are fascinating: They represent a combination of universally found motifs such as alpha-helices or beta-sheets with highly adapted protein structures such as cytoskeletal networks or spider silk nanocomposites. BPMs combine properties like strength and robustness, self-healing ability, adaptability, changeability, evolvability and others into multi-functional materials at a level unmatched in synthetic materials. The ability to achieve these properties depends critically on the particular traits of these materials, first and foremost their hierarchical architecture and seamless integration of material and structure, from nano to macro. Here, we provide a brief review of this field and outline new research directions, along with a review of recent research results in the development of structure-property relationships of biological protein materials exemplified in a study of vimentin intermediate filaments.

A multi-timescale strength model of alpha-helical protein domains

Here we report a constitutive model that characterizes the strength of an alpha-helical protein domain subjected to tensile deformation, covering more than ten orders of magnitude in timescales. The model elucidates multiple physical mechanisms of failure in dependence on the deformation rate, quantitatively linking atomistic simulation results with experimental strength measurements of alpha-helical protein domains. The model provides a description of the strength of alpha-helices based on fundamental physical parameters such as the H-bond energy and the polypeptide's persistence length, showing that strength is controlled by energetic, nonequilibrium processes at high rates and by thermodynamical, equilibrium processes at low rates. Our model provides a novel perspective on the strength of protein domains at ultra-slow pulling speeds relevant under physiologic and experimental conditions.

Importance of LPS structure on protein interactions with Pseudomonas aeruginosa

Atomic force microscopy (AFM) was used to quantify the adhesion forces between Pseudomonas aeruginosa PAO1 and AK1401, and a representative model protein, bovine serum albumin (BSA). The two bacteria strains differ in terms of the structure of their lipopolysaccharide (LPS) layers. While PAO1 is the wild-type expressing a complete LPS and two types of saccharide units in the O-antigen (A(+) B(+)), the mutant AK1401 expresses only a single unit of the A-band saccharide (A(+) B(-)). The mean adhesion force (F(adh)) between BSA and AK1401 was 1.12 nN, compared to 0.40 nN for F(adh) between BSA and PAO1. In order to better understand the fundamental forces that would control bacterial-protein interactions at equilibrium conditions, we calculated interfacial free energies using the van Oss-Chaudhury-Good (VCG) thermodynamic modeling approach. The hydrogen bond strength was also calculated using a Poisson statistical analysis. AK1401 has a higher ability to participate in hydrogen bonding with BSA than does PAO1, which may be because the short A-band and absence of B-band polymer allowed the core oligosaccharides and lipid A regions to be more exposed and to participate in hydrogen and chemical bonding. Interactions between PAO1 and BSA were weak due to the dominance of neutral and hydrophilic sugars of the A-band polymer. These results show that bacterial interactions with protein-coated surfaces will depend on the types of bonds that can form between bacterial surface macromolecules and the protein. We suggest that strategies to prevent bacterial colonization of biomaterials can focus on inhibiting these bonds.

Hierarchies, multiple energy barriers, and robustness govern the fracture mechanics of alpha-helical and beta-sheet protein domains

The fundamental fracture mechanisms of biological protein materials remain largely unknown, in part, because of a lack of understanding of how individual protein building blocks respond to mechanical load. For instance, it remains controversial whether the free energy landscape of the unfolding behavior of proteins consists of multiple, discrete transition states or the location of the transition state changes continuously with the pulling velocity. This lack in understanding has thus far prevented us from developing predictive strength models of protein materials. Here, we report direct atomistic simulation that over four orders of magnitude in time scales of the unfolding behavior of alpha-helical (AH) and beta-sheet (BS) domains, the key building blocks of hair, hoof, and wool as well as spider silk, amyloids, and titin. We find that two discrete transition states corresponding to two fracture mechanisms exist. Whereas the unfolding mechanism at fast pulling rates is sequential rupture of individual hydrogen bonds (HBs), unfolding at slow pulling rates proceeds by simultaneous rupture of several HBs. We derive the hierarchical Bell model, a theory that explicitly considers the hierarchical architecture of proteins, providing a rigorous structure-property relationship. We exemplify our model in a study of AHs, and show that 3-4 parallel HBs per turn are favorable in light of the protein's mechanical and thermodynamical stability, in agreement with experimental findings that AHs feature 3.6 HBs per turn. Our results provide evidence that the molecular structure of AHs maximizes its robustness at minimal use of building materials.

Direct observation of active protein folding using lock-in force spectroscopy

Direct observation of the folding of a single polypeptide chain can provide important information about the thermodynamic states populated along its folding pathway. In this study, we present a lock-in force-spectroscopy technique that improves resolution of atomic-force microscopy force spectroscopy to 400 fN. Using this technique we show that immunoglobulin domain 4 from Dictyostelium discoideum filamin (ddFLN4) refolds against forces of approximately 4 pN. Our data show folding of this domain proceeds directly from an extended state and no thermodynamically distinct collapsed state of the polypeptide before folding is populated. Folding of ddFLN4 under load proceeds via an intermediate state. Three-state folding allows ddFLN4 to fold against significantly larger forces than would be possible for a mere two-state folder. We present a general model for protein folding kinetics under load that can predict refolding forces based on chain-length and zero force refolding rate.

FTIR spectra of solid poly-l-lysine in the stretching NH mode range

Three bands at 3270 cm(-1), 3200 cm(-1) and 3030 cm(-1) are found in the IR stretching proton (nu(1)) mode spectral range in spectra of solid poly-l-lysine (PLL). Strong quantitative changes of these bands are observed in samples dried from water solutions with different pH. The narrow band at 3270 cm(-1), which is strong in the spectrum of PLL precipitated from pH=12 alkaline medium, is assigned to the nu(1) peptide proton mode of NH-CO (amide A) of the beta-sheet structure type. The band at 3200 cm(-1), which is intensified in PLL precipitated from pH=1 acidic medium, relates to the nu(1) peptide mode in the random coil structure. The band at 3030 cm(-1), whose peak intensity increases two-fold in going from alkaline to acidic medium, is assigned to the nu(1) modes of protonated NH(3)(+) side chain groups. The frequencies of all bands were used for estimating H-bond energy relying on an empirical correlation between this property and the red shift of the nu(1) band. The enthalpy of the secondary structure transition from beta-sheet to the random coil, which is observed in PLL at the change of pH from 11 to 1 amounts to 4.7 kJ mol(-1).

Specific and nonspecific interaction forces between Escherichia coli and silicon nitride, determined by poisson statistical analysis

The nature of the physical interactions between Escherichia coli JM109 and a model surface (silicon nitride) was investigated in water via atomic force microscopy (AFM). AFM force measurements on bacteria can represent the combined effects of van der Waals and electrostatic forces, hydrogen bonding, steric interactions, and perhaps ligand-receptor type bonds. It can be difficult to decouple these forces into their individual components since both specific (chemical or short-range forces such as hydrogen bonding) and nonspecific (long-range colloidal) forces may be present in the overall profiles. An analysis is presented based on the application of Poisson statistics to AFM adhesion data, to decouple the specific and nonspecific interactions. Comparisons with classical DLVO theory and a modified form of a van der Waals expression for rough surfaces were made in order to help explain the nature of the interactions. The only specific forces in the system were due to hydrogen bonding, which from the Poisson analysis were found to be -0.125 nN. The nonspecific forces of 0.155 nN represent an overall repulsive interaction. These nonspecific forces are comparable to the forces calculated from DLVO theory, in which electrostatic-double layer interactions are added to van der Waals attractions calculated at the distance of closest approach, as long as the van der Waals model for "rough" spherical surfaces is used. Calculated electrostatic-double layer and van der Waals interactions summed to 0.116 nN. In contrast, if the classic (i.e., smooth) sphere-sphere model was used to predict the van der Waals forces, the sum of electrostatic and van der Waals forces was -7.11 nN, which appears to be a large overprediction. The Poisson statistical analysis of adhesion forces may be very useful in applications of bacterial adhesion, because it represents an easy way to determine the magnitude of hydrogen bonding in a given system and it allows the fundamental forces to be easily broken into their components.

Coiled-coil nanomechanics and uncoiling and unfolding of the superhelix and alpha-helices of myosin

The nanomechanical properties of the coiled-coils of myosin are fundamentally important in understanding muscle assembly and contraction. Force spectra of single molecules of double-headed myosin, single-headed myosin, and coiled-coil tail fragments were acquired with an atomic force microscope and displayed characteristic triphasic force-distance responses to stretch: a rise phase (R) and a plateau phase (P) and an exponential phase (E). The R and P phases arise mainly from the stretching of the coiled-coils, with the hinge region being the main contributor to the rise phase at low force. Only the E phase was analyzable by the worm-like chain model of polymer elasticity. Restrained molecular mechanics simulations on an existing x-ray structure of scallop S2 yielded force spectra with either two or three phases, depending on the mode of stretch. It revealed that coiled-coil chains separate completely near the end of the P phase and the stretching of the unfolded chains gives rise to the E phase. Extensive conformational searching yielded a P phase force near 40 pN that agreed well with the experimental value. We suggest that the flexible and elastic S2 region, particularly the hinge region, may undergo force-induced unfolding and extend reversibly during actomyosin powerstroke.

Sacrificial bonds and hidden length: unraveling molecular mesostructures in tough materials

Sacrificial bonds and hidden length in structural molecules and composites have been found to greatly increase the fracture toughness of biomaterials by providing a reversible, molecular-scale energy-dissipation mechanism. This mechanism relies on the energy, of order 100 eV, needed to reduce entropy and increase enthalpy as molecular segments are stretched after being released by the breaking of weak bonds, called sacrificial bonds. This energy is relatively large compared to the energy needed to break the polymer backbone, of order a few eV. In many biological cases, the breaking of sacrificial bonds has been found to be reversible, thereby additionally providing a "self-healing" property to the material. Due to the nanoscopic nature of this mechanism, single molecule force spectroscopy using an atomic force microscope has been a useful tool to investigate this mechanism. Especially when investigating natural molecular constructs, force versus distance curves quickly become very complicated. In this work we propose various types of sacrificial bonds, their combination, and how they appear in single molecule force spectroscopy measurements. We find that by close analysis of the force spectroscopy curves, additional information can be obtained about the molecules and their bonds to the native constructs.

Frequency modulation atomic force microscopy reveals individual intermediates associated with each unfolded I27 titin domain

In this study, we apply a dynamic atomic force microscopy (AFM) technique, frequency modulation (FM) detection, to the mechanical unfolding of single titin I27 domains and make comparisons with measurements made using the AFM contact or static mode method. Static mode measurements revealed the well-known force transition occurring at 100-120 pN in the first unfolding peak, which was less clear, or more often absent, in the subsequent unfolding peaks. In contrast, some FM-AFM curves clearly resolved a force transition associated with each of the unfolding peaks irrespective of the number of observed unfolded domains. As expected for FM-AFM, the frequency shift response of the main unfolding peaks and their intermediates could only be detected when the oscillation amplitudes used were smaller than the interaction lengths being measured. It was also shown that the forces measured for the dynamical interaction of the FM-AFM technique were significantly lower than those measured using the static mode. This study highlights the potential for using dynamic AFM for investigating biological interactions, including protein unfolding and the detection of novel unfolding intermediates.

A new scheme for determining the intramolecular seven-membered ring N–H⋯OC hydrogen-bonding energies of glycine and alanine peptides

In this paper a new scheme was proposed to calculate the intramolecular hydrogen-bonding energies in peptides and was applied to calculate the intramolecular seven-membered ring N-H...O=C hydrogen-bonding energies of the glycine and alanine peptides. The density-functional theory B3LYP6-31G(d) and B3LYP6-311G(d,p) methods and the second-order Moller-Plesset perturbation theory MP26-31G(d) method were used to calculate the optimal geometries and frequencies of glycine and alanine peptides and related structures. MP26-311++G(d,p), MP26-311++G(3df,2p), and MP2/aug-cc-pVTZ methods were then used to evaluate the single-point energies. It was found that the B3LYP6-31G(d), MP26-31G(d), and B3LYP6-311G(d,p) methods yield almost similar structural parameters for the conformers of the glycine and alanine dipeptides. MP2/aug-cc-pVTZ predicts that the intramolecular seven-membered ring N-H...O=C hydrogen-bonding strength has a value of 5.54 kcal/mol in glycine dipeptide and 5.73 and 5.19 kcal/mol in alanine dipeptides, while the steric repulsive interactions of the seven-membered ring conformers are 4.13 kcal/mol in glycine dipeptide and 6.62 and 3.71 kcal/mol in alanine dipeptides. It was also found that MP26-311++G(3df,2p) gives as accurate intramolecular N-H...O=C hydrogen-bonding energies and steric repulsive interactions as the much more costly MP2/aug-cc-pVTZ does.

Nonlinear elasticity of an alpha -helical polypeptide

We study a minimal extension of the wormlike chain model to describe polypeptides having alpha-helical secondary structure. In this model the presence or absence of secondary structure enters as a scalar variable that controls the local chain bending modulus. Using this model we compute the extensional compliance of an alpha-helix under tensile stress, the bending compliance of the molecule under externally imposed torques, and the nonlinear interaction of such torques and forces on the molecule. We find that, due to coupling of the "internal" secondary structure variables to the conformational degrees of freedom of the polymer, the molecule has a highly nonlinear response to applied stress and force couples. In particular we demonstrate a sharp lengthening transition under applied force and a buckling transition under applied torque. Finally, we speculate that the inherent bistability of the molecule may underlie protein conformational change in vivo.

Effects of temperature and pH on the helicity of a peptide adsorbed to colloidal silica

The conformation of a cationic alpha-helical peptide (DDDDAAAARRRRR) adsorbed to anionic colloidal silica has been investigated by circular dichroism (CD) spectroscopy as a function of temperature and pH in order to examine how the structure of an adsorbed molecule responds to two simultaneous perturbations. Increased temperature destabilizes the helicity of the peptide in solution, while pH changes alter the substrate surface charge and the corresponding strength of the interaction with the peptide. Near neutral pH, the helicity of the adsorbed peptide, which is determined from the intensity of the CD signal at 222 nm, decreases with increasing temperature, similarly to the temperature-dependent behavior observed for the peptide in aqueous solution. By contrast, at basic pH and a strongly negative surface charge, the helicity of the adsorbed peptide increases with temperature. In order to elucidate the origin of the reversal of the temperature dependence of helicity, a statistical model for the conformation of the adsorbed peptide has been formulated based on the Lifson-Roig model for the helix-coil transition of the peptide in solution. The model provides insight into the trends in fractional helicity and reveals that the temperature dependence of the helicity of the adsorbed peptide results from a competition between the intramolecular interactions that promote helicity and the intermolecular interactions with the surface. The statistical model also enables estimation of the free energy contributions from specific aspects of the adsorption process. Through identification of a connection between the conformation of adsorbed peptide and the interactions of the peptide with the surface, this work suggests a route for the control of adsorbate conformation through peptide and surface engineering.

Mechanical force analysis of peptide interactions using atomic force microscopy

Some peptides have previously been reported to bind low molecular weight chemicals. One such peptide with the amino acid sequence His-Ala-Ser-Tyr-Ser was selectively screened from a phage library and bound to a cationic porphyrin, 5,10,15,20-tetrakis(N-methylpyridinium-4-yl)-21H,23H-porphine (TMpyP), with a binding constant of 10(5) M(-1) (J. Kawakami, T. Kitano, and N. Sugimoto, Chemical Communications, 1999, pp. 1765-1766). The proposed binding was due to pi-electron stacking from two aromatic amino acids of histidine and tyrosine. In this study, the weak interactions between TMpyP and the peptide were further investigated by force curve analysis using atomic force microscopy (AFM). The mechanical force required to unbind the peptide-porphyrin complex was measured by vertical movement of the AFM tip. Peptide self-assembled monolayers were formed on both a gold-coated mica substrate and a gold-coated AFM tip. The TMpyPs could bind between the two peptide layers when the peptide-immobilized AFM tip contacted the peptide-immobilized substrate in solution containing TMpyP. In the retracting process a force that ruptured the interaction between TMpyPs and peptides was observed. The unbinding force values correlated to the concentration of TMpyP. A detection limit of 100 ng/mL porphyrin was obtained for the force measurement, and was similar to surface plasmon resonance sensor detection limits. Furthermore, we calculated the product of the observed force and the length of the molecular elongation to determine the work required to unbind the complexes. The obtained values of unbinding work were in a reasonable range compared to the binding energy of porphyrin-peptide.

Redox-driven proton pumping by heme-copper oxidases

One of the key problems of molecular bioenergetics is the understanding of the function of redox-driven proton pumps on a molecular level. One such class of proton pumps are the heme-copper oxidases. These enzymes are integral membrane proteins in which proton translocation across the membrane is driven by electron transfer from a low-potential donor, such as, e.g. cytochrome c, to a high-potential acceptor, O(2). Proton pumping is associated with distinct exergonic reaction steps that involve gradual reduction of oxygen to water. During the process of O(2) reduction, unprotonated high pK(a) proton acceptors are created at the catalytic site. Initially, these proton acceptors become protonated as a result of intramolecular proton transfer from a residue(s) located in the membrane-spanning part of the enzyme, but removed from the catalytic site. This residue is then reprotonated from the bulk solution. In cytochrome c oxidase from Rhodobacter sphaeroides, the proton is initially transferred from a glutamate, E(I-286), which has an apparent pK(a) of 9.4. According to a recently published structure of the enzyme, the deprotonation of E(I-286) is likely to result in minor structural changes that propagate to protonatable groups on the proton output (positive) side of the protein. We propose that in this way, the free energy available from the O(2) reduction is conserved during the proton transfer. On the basis of the observation of these structural changes, a possible proton-pumping model is presented in this paper. Initially, the structural changes associated with deprotonation of E(I-286) result in the transfer of a proton to an acceptor for pumped protons from the input (negative) side of the membrane. After reprotonation of E(I-286) this acceptor releases a proton to the output side of the membrane.

Measurement of the length of the a helical section of a peptide directly using atomic force microscopy

Using atomic force microscopy (AFM), the length of the alpha-helix structure of poly-L-lysine was investigated by stretching the peptide directly, one molecule at a time. In the absence of urea, many rupturing points that seemed to be due to the breaking of some hydrogen bonds were observed in force-extension curves, while these points were never observed in the presence of 8 M urea. In the presence of 0.4 or 1.6 M urea, both force-extension curve types were observed. Total peptide elongation for each condition was calculated from force-extension curves reflecting the alpha-helix rupturing process. The experimental value of total elongation divided by the theoretical value of total alpha-helix elongation yields the alpha-helix content. This value was compatible with circular dichroism (CD) measurement results. This suggests that peptide conformation and content of the alpha-helix on a single molecule scale can be investigated by direct mechanical measurement using atomic force microscopy.

Oligomerization of light-harvesting I antenna peptides of Rhodospirillum rubrum

We investigated the oligomerization of the core light-harvesting complex (LH1) of Rhodospirillum rubrum from the separated alpha beta BChl(2) subunits (B820) and the oligomerization of the B820 subunit from its monomeric peptides. The full LH1 complex was reversibly associated from B820 subunits by either varying the temperature in the range 277-300 K or by varying the detergent concentration in the buffer from 0.36 to 0.52% n-octyl-beta-D-glucopyranoside. Temperature-induced transition measurements showed hysteresis: raising the temperature induced dissociation of B873 directly into B820 subunits whereas upon recooling an intermediate spectral form was observed with an absorption maximum located around 850 nm. This intermediate form was also observed in detergent-induced transitions. It is speculated that the B850 form is a small aggregate of B820, for instance a dimer. Additionally, during a temperature-mediated transition at low detergent concentration, a set of spectral forms with maxima slightly blue-shifted from 873 nm were observed, possibly due to opened rings with one or only a few alpha beta BChl(2) units missing. The temperature-induced transition of LH1 is discussed in terms of a simple assembly model. It is concluded that a moderately cooperative assembly explains the formation of small aggregates of B820 as well as of incomplete rings. Furthermore, the B820 subunits were reversibly dissociated into the monomeric B777 form by increasing either the temperature or the detergent concentration. Estimations of the enthalpy and entropy changes for the dimeric association reaction of B777 into B820 yielded an enthalpy change of -216 kJ mol(-1) and an entropy change of -0.59 kJ mol(-1)K(-1), at a detergent concentration of 0.8% n-octyl-beta-D-glucopyranoside.

Stretching a macromolecule in an atomic force microscope: statistical mechanical analysis

We formulate the proper statistical mechanics to describe the stretching of a macromolecule under a force provided by the cantilever of an Atomic Force Microscope. In the limit of a soft cantilever, the generalized ensemble of the coupled molecule-cantilever system reduces to the Gibbs ensemble for an isolated molecule subject to a constant force in which the extension is fluctuating. For a stiff cantilever, one obtains the Helmholtz ensemble for an isolated molecule held at a fixed extension with the force fluctuating. Numerical examples and predictions for experiments with cantilevers of differing stiffness are given for short and long chains of poly (ethylene glycol), based on parameter-free ab initio calculations.

Helix-coil transition in homopolypeptides under stretching

We consider the effect of an external applied force on the alpha-helix-coil transition of a single-stranded homopolypeptide chain. An annealed scenario is assumed, where the building amino acid monomers may interconvert between random-coiled and ordered alpha-helical configurations. By exact evaluation of the partition function of the freely jointed chain with helix-coil internal degrees of freedom in the thermodynamic limit, we obtain the result that the stress-strain characteristic has an asymmetrical sigmoid shape with a prominent pseudoplateau. Because of the one-dimensional nature of this system, fluctuations dominate over the mean-field approximation, which incorrectly predicts a second-order phase transition.

Stretching a macromolecule in an atomic force microscope: statistical mechanical analysis

We formulate the proper statistical mechanics to describe the stretching of a macromolecule under a force provided by the cantilever of an atomic force microscope. In the limit of a soft cantilever the generalized ensemble of the coupled molecule/cantilever system reduces to the Gibbs ensemble for an isolated molecule subject to a constant force in which the extension is fluctuating. For a stiff cantilever we obtain the Helmholtz ensemble for an isolated molecule held at a fixed extension with the force fluctuating. Numerical examples are given for poly (ethylene glycol) chains.

Advances in the characterization of supported lipid films with the atomic force microscope

During the past decade, the atomic force microscope (AFM) has become a key technique in biochemistry and biophysics to characterize supported lipid films, as testified by the continuous growth in the number of papers published in the field. The unique capabilities of AFM are: (i) capacity to probe, in real time and in aqueous environment, the surface structure of lipid films; (ii) ability to directly measure physical properties at high spatial resolution; (iii) possibility to modify the film structure and biophysical processes in a controlled way. Such experiments, published up to June 2000, are the focus of the present review. First, we provide a general introduction on the preparation and characterization of supported lipid films as well as on the principles of AFM. The section 'Structural properties' focuses on the various applications of AFM for characterizing the structure of supported lipid films: visualization of molecular structure, formation of structural defects, effect of external agents, formation of supported films, organization of phase-separated films (coexistence region, mixed films) and, finally, the use of supported lipid bilayers for anchoring biomolecules such as DNA, enzymes and crystalline protein arrays. The section 'Physical properties' introduces the principles of force measurements by AFM, interpretation of these measurements and their recent application to supported lipid films and related structures. Finally, we highlight the major achievements brought by the technique and some of the current limitations.

Spring mechanics of alpha-helical polypeptide

To design protein- and polymer-based micro-machineries, it is important to understand the mechanical properties of basic structural elements such as the alpha-helix of polypeptides. We employed the force measurement mode of an atomic force microscope (AFM) to investigate the spring mechanics of poly-L-glutamic acid (PGA) in its helical and randomly coiled states. After covalently anchoring the polypeptide between a silicon substrate and an AFM tip, the force required to stretch the polymer was measured. The results indicated that PGA in its helical conformation could be stretched almost fully with a continuous increase in the stretching force, suggesting that it can be used as a reliable coil-spring in the future design of spring-loaded molecular machineries.