Accelerating Steered Molecular Dynamics: Toward Smaller Velocities in Forced Unfolding Simulations

Accumulation of Hydrogen Bonds and van der Waals Interactions Determines Force Response between Two Parallel Cellulose Chains: Steered Molecular Dynamics Simulations

We investigated the response forces between two parallel cellulose chains during the shearing and tearing processes by using steered molecular dynamics simulations. It was found that there are two logarithmic dependencies between response force and pulling speed in shearing processes but only one in tearing, according to Bell's equation by fitting the f-ln v curve. The mechanism is that there are 2-fold interactions determining the force response between two parallel cellulose chains resisting chain separation during a shearing process. Our results indicate that hydrogen bonds dominate the interchain interactions in the fast pull mode (FPM) for shearing, while van der Waals interactions dominate in the slow pull mode (SPM). For tearing, the one-by-one breaking of hydrogen bonds and van der Waals interactions plays a main role.

Plasticity of the proteasome-targeting signal Fat10 enhances substrate degradation

The proteasome controls levels of most cellular proteins, and its activity is regulated under stress, quiescence, and inflammation. However, factors determining the proteasomal degradation rate remain poorly understood. Proteasome substrates are conjugated with small proteins (tags) like ubiquitin and Fat10 to target them to the proteasome. It is unclear if the structural plasticity of proteasome-targeting tags can influence substrate degradation. Fat10 is upregulated during inflammation, and its substrates undergo rapid proteasomal degradation. We report that the degradation rate of Fat10 substrates critically depends on the structural plasticity of Fat10. While the ubiquitin tag is recycled at the proteasome, Fat10 is degraded with the substrate. Our results suggest significantly lower thermodynamic stability and faster mechanical unfolding in Fat10 compared to ubiquitin. Long-range salt bridges are absent in the Fat10 structure, creating a plastic protein with partially unstructured regions suitable for proteasome engagement. Fat10 plasticity destabilizes substrates significantly and creates partially unstructured regions in the substrate to enhance degradation. NMR-relaxation-derived order parameters and temperature dependence of chemical shifts identify the Fat10-induced partially unstructured regions in the substrate, which correlated excellently to Fat10-substrate contacts, suggesting that the tag-substrate collision destabilizes the substrate. These results highlight a strong dependence of proteasomal degradation on the structural plasticity and thermodynamic properties of the proteasome-targeting tags.

The Thermodynamic Properties of Fat10ylated Proteins Are Regulated by the Fat10ylation Site

Degradation of proteins by the proteasome is crucial in regulating their levels in the cell. Post-translational modifications, such as ubiquitylation and Fat10ylation, trigger proteasomal degradation of the substrate proteins. While ubiquitylation regulates multiple cellular pathways, Fat10ylation functions explicitly in the inflammatory response pathway. At the proteasome, ubiquitin is recycled after being cleaved from the substrate, while Fat10 is degraded simultaneously with its substrate. Although the thermodynamic properties of the substrate are critical for effective proteasomal degradation, they remain poorly understood for the Fat10-proteasome pathway. We studied the thermodynamic properties of the Fat10∼substrate conjugate to uncover mechanistic details of the pathway. First, the mechanical unfolding of Fat10∼substrate was studied by molecular dynamics simulations, which suggested that the unfolding pathway and unfolding energy of the substrate depend on the site of Fat10 modification. We also investigated different pathways for the entry of the Fat10∼substrate into the proteasome core. Our analysis supports a model where the entry of Fat10, followed by the substrate, is the energetically preferred pathway. Further, we studied Fat10's effect on the thermodynamic properties of distinct substrates, considering their size, flexibility, and surface properties. The results uncovered significant entropic destabilization of substrates due to Fat10ylation, particularly in smaller substrates. For larger substrates, multi-monoFat10ylation is necessary to induce destabilization. Our study further reveals that Fat10 modification at negative patches on substrate surfaces is essential for optimal destabilization and subsequent degradation. These findings provide atomistic insights into the degradation mechanisms in the Fat10 proteasome pathway with potential implications for therapeutic interventions.

The Mechanism of Channel Opening of Anion Channelrhodopsin GtACR1: A Molecular Dynamics Simulation

Guillardia theta anion channelrhodopsin 1 (GtACR1) is a widely used inhibitor of optogenetics with unique conductance mechanisms and photochemistry. However, the molecular mechanism of light-gated anion conduction is poorly understood without a crystal structure for the intermediate state. In this study, we built the dark-state model based on the crystal structure of retinal and isomerized the model by twisting the C12-C13=C14-C15 dihedral step by step using molecular dynamics simulation. The conformational changes revealed the all-trans to 13-cis photoisomerization of the retinal chromophore cannot open the channel. There is no water influx, and a pre-opened K-like intermediate after photoisomerization of retinal is formed. During the opening of the ion channel, proton transfer occurs between E68 and D234. Steered molecular dynamics (SMD) and umbrella sampling indicated that the E68 and D234 were the key residues for chloride-ion conducting. We propose a revised channel opening pathway model of GtACR1 after analyzing (de)protonation of E68 and D234. Reprotonation of D234 will result in two different early L intermediates, named L1-like and L1‘-like, which correspond to the L1 and L1‘ intermediates reported in a recent study. Simulation results showed that L1-like may convert by parallel paths into L1‘-like and L2-like states. This model provides conformational details for the intermediate as well.

Molecular interpretation of single-molecule force spectroscopy experiments with computational approaches

Single molecule force-spectroscopy techniques have granted access to unprecedented molecular-scale details about biochemical and biological mechanisms. However, the interpretation of the experimental data is often challenging. Computational and simulation approaches (all-atom steered MD simulations in particular) are key to provide molecular details about the associated mechanisms, to help test different hypotheses and to predict experimental results. In this review, particular recent efforts directed towards the molecular interpretation of single-molecule force spectroscopy experiments on proteins and protein-related systems (often in close collaboration with experimental groups) will be presented. These results will be discussed in the broader context of the field, highlighting the recent achievements and the ongoing challenges for computational biophysicists and biochemists. In particular, I will focus on the input gained from molecular simulations approaches to rationalize the origin of the unfolded protein elasticity and the protein conformational behavior under force, to understand how force denaturation differs from chemical, thermal or shear unfolding, and to unravel the molecular details of unfolding events for a variety of systems. I will also discuss the use of models based on Langevin dynamics on a 1-D free-energy surface to understand the effect of protein segmentation on the work exerted by a force, or, at the other end of the spectrum of computational techniques, how quantum calculations can help to understand the reactivity of disulfide bridges exposed to force.

Recent Advances and Emerging Challenges in the Molecular Modeling of Mechanobiological Processes

Many biological processes result from the effect of mechanical forces on macromolecular structures and on their interactions. In particular, the cell shape, motion, and differentiation directly depend on mechanical stimuli from the extracellular matrix or from neighboring cells. The development of experimental techniques that can measure and characterize the tiny forces acting at the cellular scale and down to the single-molecule, biomolecular level has enabled access to unprecedented details about the involved mechanisms. However, because the experimental observables often do not provide a direct atomistic picture of the corresponding phenomena, particle-based simulations performed at various scales are instrumental in complementing these experiments and in providing a molecular interpretation. Here, we will review the recent key achievements in the field, and we will highlight and discuss the many technical challenges these simulations are facing, as well as suggest future directions for improvement.

Effect of Loading on the Adhesion and Frictional Characteristics of Top Layer Articular Cartilage Nanoscale Contact: A Molecular Dynamics Study

Articular cartilage is a water-lubricated naturally occurring biological interface imparting unique mechanical and ultralow frictional properties in bone joints. Although the material of cartilage, synovial fluid composition, and their lubricating modes and properties have been extensively investigated at various scales experimentally, there is still a lack of understanding of load bearing, adhesion, and friction mechanisms of the cartilage-cartilage interface from an atomistic perspective under heavy loads. In this study, the effect of loading on adhesion and frictional behavior in articular cartilage is investigated with a proposed atomistic model for top layer cartilage-cartilage contact in unhydrated conditions using molecular dynamics (MD) simulations. Pull-off tests reveal that cohesive interactions occur at the interface due to formation of heavily interpenetrated atomistic sites leading to stretching and localized pulling of fragments during sliding. Sliding tests show that friction is load- and direction-dependent with the coefficient of friction (COF) obtained in the range of 0.20-0.75 at the interface for sliding in parallel and perpendicular directions to the collagen axis. These values are in good agreement with earlier nanoscale experimental results reported for the top layer cartilage-cartilage interface. The COF reduces with an increase in load and tends to be higher for the parallel sliding case than for the perpendicular case owing to the presence of the constant number of H-bonds. Overall, this work contributes toward understanding sliding in unhydrated biointerfaces, which is the precursor of wear, and provides insights into implant research.

Mechanical Unfolding of Spectrin Repeats Induces Water-Molecule Ordering

Mechanical processes are involved at many stages of the development of living cells, and often external forces applied to a biomolecule result in its unfolding. Although our knowledge of the unfolding mechanisms and the magnitude of the forces involved has evolved, the role that water molecules play in the mechanical unfolding of biomolecules has not yet been fully elucidated. To this end, we investigated with steered molecular dynamics simulations the mechanical unfolding of dystrophin's spectrin repeat 1 and related the changes in the protein's structure to the ordering of the surrounding water molecules. Our results indicate that upon mechanically induced unfolding of the protein, the solvent molecules become more ordered and increase their average number of hydrogen bonds. In addition, the unfolded structures originating from mechanical pulling expose an increasing amount of the hydrophobic residues to the solvent molecules, and the uncoiled regions adapt a convex surface with a small radius of curvature. As a result, the solvent molecules reorganize around the protein's small protrusions in structurally ordered waters that are characteristic of the so-called "small-molecule regime," which allows water to maintain a high hydrogen bond count at the expense of an increased structural order. We also determined that the response of water to structural changes in the protein is localized to the specific regions of the protein that undergo unfolding. These results indicate that water plays an important role in the mechanically induced unfolding of biomolecules. Our findings may prove relevant to the ever-growing interest in understanding macromolecular crowding in living cells and their effects on protein folding, and suggest that the hydration layer may be exploited as a means for short-range allosteric communication.

Molecular Dynamics Simulation of Transmembrane Transport of Chloride Ions in Mutants of Channelrhodopsin

Channelrhodopsins (ChRs) are light-gated transmembrane cation channels which are widely used for optogenetic technology. Replacing glutamate located at the central gate of the ion channel with positively charged amino acid residues will reverse ion selectivity and allow anion conduction. The structures and properties of the ion channel, the transport of chloride, and potential of mean force (PMF) of the chimera protein (C1C2) and its mutants, EK-TC, ER-TC and iChloC, were investigated by molecular dynamics simulation. The results show that the five-fold mutation in E122Q-E129R-E140S-D195N-T198C (iChloC) increases the flexibility of the transmembrane channel protein better than the double mutations in EK-TC and ER-TC, and results in an expanded ion channel pore size and decreased steric resistance. The iChloC mutant was also found to have a higher affinity for chloride ions and, based on surface electrostatic potential analysis, provides a favorable electrostatic environment for anion conduction. The PMF free energy curves revealed that high affinity Cl^- binding sites are generated near the central gate of the three mutant proteins. The energy barriers for the EK-TC and ER-TC were found to be much higher than that of iChloC. The results suggest that the transmembrane ion channel of iChloC protein is better at facilitating the capture and transport of chloride ions.

Formation Mechanism of Ion Channel in Channelrhodopsin-2: Molecular Dynamics Simulation and Steering Molecular Dynamics Simulations

Channelrhodopsin-2 (ChR2) is a light-activated and non-selective cationic channel protein that can be easily expressed in specific neurons to control neuronal activity by light. Although ChR2 has been extensively used as an optogenetic tool in neuroscience research, the molecular mechanism of cation channel formation following retinal photoisomerization in ChR2 is not well understood. In this paper, studies of the closed and opened state ChR2 structures are presented. The formation of the cationic channel is elucidated in atomic detail using molecular dynamics simulations on the all-trans-retinal (ChR2-trans) configuration of ChR2 and its isomerization products, 13-cis-retinal (ChR2-cis) configuration, respectively. Photoisomerization of the retinal-chromophore causes the destruction of interactions among the crucial residues (e.g., E90, E82, N258, and R268) around the channel and the extended H-bond network mediated by numerous water molecules, which opens the pore. Steering molecular dynamics (SMD) simulations show that the electrostatic interactions at the binding sites in intracellular gate (ICG) and central gate (CG) can influence the transmembrane transport of Na⁺ in ChR2-cis obviously. Potential of mean force (PMF) constructed by SMD and umbrella sampling also found the existing energy wells at these two binding sites during the transportation of Na⁺. These wells partly hinder the penetration of Na⁺ into cytoplasm through the ion channel. This investigation provides a theoretical insight on the formation mechanism of ion channels and the mechanism of ion permeation.

Energy Landscape of Negatively Charged BSA Adsorbed on a Negatively Charged Silica Surface

We study the energy landscape of the negatively charged protein bovine serum albumin adsorbed on a negatively charged silica surface at pH 7. We use fully atomistic molecular dynamics (MD) and steered MD (SMD) to probe the energy of adsorption and the pathway for the surface diffusion of the protein and its associated activation energy. We find an adsorption energy ∼1.2 eV, which implies that adsorption is irreversible even on experimental time scales of hours. In contrast, the activation energy for surface diffusion is ∼0.4 eV so that it is observable on the MD simulation time scale of 100 ns. This analysis paves the way for a more detailed understanding of how a protein layer forms on biomaterial surfaces, even when the protein and surface share the same electrical polarity.

α-Helix Unwinding as Force Buffer in Spectrins

Spectrins are cytoskeletal proteins located at the inner face of the plasma membrane, making connections between membrane anchors and the actin cortex, and between actin filaments. Spectrins share a common structure forming a bundle of 3 α-helices and play a major role during cell deformation. Here, we used high-speed force spectroscopy and steered molecular dynamics simulations to understand the mechanical stability of spectrin, revealing a molecular force buffering function. We find that spectrin acts as a soft spring at short extensions (70-100 Å). Under continuous external stretching, its α-helices unwind, leading to a viscous mechanical response over larger extensions (100-300 Å), represented by a constant-force plateau in force/extension curves. This viscous force buffering emerges from a quasi-equilibrium competition between disruption and re-formation of α-helical hydrogen bonds. Our results suggest that, in contrast to β-sheet proteins, which unfold in a catastrophic event, α-helical spectrins dominantly unwind, providing a viscous force buffer over extensions about 5 times their folded length.

Soft Vibrational Modes Predict Breaking Events during Force-Induced Protein Unfolding

We investigate the correlation between soft vibrational modes and unfolding events in simulated force spectroscopy of proteins. Unfolding trajectories are obtained from molecular dynamics simulations of a Gō model of a monomer of a mutant of superoxide dismutase 1 protein containing all heavy atoms in the protein, and a normal mode analysis is performed based on the anisotropic network model. We show that a softness map constructed from the superposition of the amplitudes of localized soft modes correlates with unfolding events at different stages of the unfolding process. Soft residues are up to eight times more likely to undergo disruption of native structure than the average amino acid. The memory of the softness map is retained for extensions of up to several nanometers, but decorrelates more rapidly during force drops.

Accelerated Molecular Dynamics Study of the Effects of Surface Hydrophilicity on Protein Adsorption

The adsorption of streptavidin is studied on two surfaces, graphite and titanium dioxide, using accelerated molecular dynamics. Adsorption on graphite leads to strong conformational changes while the protein spreads out over the surface. Interestingly, also adsorption on the highly hydrophilic rutile surface induces considerable spreading of the protein. We pin down the cause for this unfolding to the interaction of the protein with the ordered water layers above the rutile surface. For special orientations, the protein penetrates the ordered water layers and comes into direct contact with the surface where the positively charged amino acids settle in places adjacent to the negatively charged top surface atom layer of rutile. We conclude that for both surface materials studied, streptavidin changes its conformation so strongly that it loses its potential for binding biotin. Our results are in good qualitative agreement with available experimental studies.

Forced Desorption of Bovine Serum Albumin and Lysozyme from Graphite: Insights from Molecular Dynamics Simulation

We use molecular dynamics (MD) simulation to study the adsorption and desorption of two widely different proteins, bovine serum albumin (BSA) and lysozyme, on a graphite surface. The adsorption is modeled using accelerated MD to allow the proteins to find optimum conformations on the surface. Our results demonstrate that the "hard protein" lysozyme retains much of its secondary structure during adsorption, whereas BSA loses it almost completely. BSA has a considerably larger adsorption energy compared to that of lysozyme, which does not scale with chain length. Desorption simulations are carried out using classical steered MD. The BSA molecule becomes fully unzipped during pull-off, whereas several helices survive this process in lysozyme. The unzipping process shows up in the force-distance curve of BSA as a series of peaks, whereas only a single or few, depending on protein orientation, force peaks occur for lysozyme. The maximum desorption force is larger for BSA than for lysozyme, but only by a factor of about 2.3.