A common approach to inducing selective mechanochemical transformations relies on embedding the target molecules (called mechanophores) within elastomeric polymer networks. Mechanical properties of such elastomers can also be modulated through the mechanochemical response of the constituent polymer chains. The inherent randomness in the molecular structure of such materials leads to heterogeneity of the local forces exerted on individual mechanophores. Here we use coarse-grained simulations to study the force distributions within random elastomeric networks and show that those distributions are close to exponential regardless of the applied macroscopic load, entanglement effects, or network parameters. Exponential form of the distribution allows one to completely characterize the mechanophore kinetics in terms of the mean value of the force. At the same time, heterogeneity of the local force affects the kinetics qualitatively: While a narrow force distribution around the mean would lead to exponential kinetics, exponential force distribution results in highly nonexponential kinetics, with a fast kinetic phase involving highly loaded molecules, followed by a slow phase dominated by unloaded molecules.

Protein aggregation and amyloid formation are implicated in many diseases as well as in other biological phenomena. Recent studies have suggested that amyloid formation of tumor suppressor p53 can lead to loss of its physiological function, resulting in accelerated cancer progression. Design of cancer therapeutics, therefore, requires understanding of the mechanism of p53 aggregation. Here, we have employed atomistic simulations to characterize the aggregation process of the aggregation-prone (as suggested by experimental studies) p53 fragment (LTIITLE, 252–258) and to assess the efficiency of its I254R mutant as an aggregation suppressor. We show that the wild-type sequence attains stable β-sheet rich structure in the parallely arranged dimeric form, which dissociates in a sequential manner under mechanical force. The wild-type sequence further displays high aggregation propensity self-assembling into structures with parallel peptide arrangement. The I254R mutation destabilizes the dimer, changes the mechanical dissociation of the dimer to cooperative unfolding, reduces the aggregation propensity of the sequence, and alters the relative orientation of the peptides in the aggregate. Addition of the wild-type sequence, however, partially restores the aggregation propensity of the I254R mutant.

The notions of a reaction path and a reaction coordinate are central to chemistry as they provide low-dimensional descriptions of complex molecular processes. Here we discuss how to define, compute, and use the reaction paths for chemical transformations in molecules that are subjected to mechanical stress and thus driven toward regions of conformational space that are otherwise inaccessible both in computational studies and in reality. We further show that the circuitous nature of mechanochemical pathways often makes their one-dimensional description impossible and describe how multidimensional effects can be detected experimentally.

While the field of polymer mechanochemistry has traditionally focused on the use of mechanical forces to accelerate chemical processes, theoretical considerations predict an underexplored alternative: the suppression of reactivity through mechanical perturbation. Here, we use electronic structure calculations to analyze the mechanical reactivity of six mechanophores, or chemical functionalities that respond to mechanical stress in a controlled manner. Our computational results indicate that appropriately directed tensile forces could attenuate (as opposed to facilitate) mechanochemical phenomena. Accompanying experimental studies supported the theoretical predictions and demonstrated that relatively simple computational models may be used to design new classes of mechanically responsive materials. In addition, our computational studies and theoretical considerations revealed the prevalence of the anti-Hammond (as opposed to Hammond) effect (i.e., the increased structural dissimilarity between the reactant and transition state upon lowering of the reaction barrier) in the mechanical activation of polyatomic molecules.

The dynamics of surface-attached polymers play a key role in the operation of a number of biological sensors, yet its current understanding is rather limited. Here we use computer simulations to study the dynamics of a reaction between the free end of a polymer chain and a surface, to which its other end has been attached. We consider two limiting cases, the diffusion-controlled limit, where the reaction is accomplished whenever the free chain end diffuses to within a specified distance from the surface, and the reaction-controlled limit, where slow, intrinsic reaction kinetics rather than diffusion of the chain is rate limiting. In the diffusion-controlled limit, we find that the overall rate scales as N−b, where N is the number of monomers in the chain and b ≈ 2.2 for excluded volume chains. This value of the scaling exponent b is close to that derived from a simple approximate theory treating the dynamics of the chain end relative to the surface as one-dimensional diffusion in an effective potential. In the reaction-controlled limit, the value of the scaling exponent b is close to 1. We compare our findings with those for the related (and better studied) problem of end-to-end reactions within an unconstrained polymer chain and discuss their implications for electrochemical DNA sensors.

Understanding the rate at which various parts of a molecular chain come together to facilitate the folding of a biopolymer (e.g., a protein or RNA) into its functional form remains an elusive goal. Here we use experiments, simulations, and theory to study the kinetics of internal loop closure in disordered biopolymers such as single-stranded oligonucleotides and unfolded proteins. We present theoretical arguments and computer simulation data to show that the relationship between the timescale of internal loop formation and the positions of the monomers enclosing the loop can be recast in a form of a universal master dependence. We also perform experimental measurements of the loop closure times of single-stranded oligonucleotides and show that both these and previously reported internal loop closure kinetics of unfolded proteins are well described by this theoretically predicted dependence. Finally, we propose that experimental deviations from the master dependence can then be used as a sensitive probe of dynamical and structural order in unfolded proteins and other biopolymers.