Mechanical Unfolding of Macromolecules Coupled to Bond Dissociation

Multiscale modeling of unfolding and bond dissociation of rubredoxin metalloprotein

Mechanical properties of proteins that have a crucial effect on their operation. This study used a molecular dynamics simulation package to investigate rubredoxin unfolding on the atomic scale. Different simulation techniques were applied, and due to the dissociation of covalent/hydrogen bonds, this protein demonstrates several intermediate states in force-extension behavior. A conceptual model based on the cohesive finite element method was developed to consider the intermediate damages that occur during unfolding. This model is based on force-displacement curves derived from molecular dynamics results. The proposed conceptual model is designed to accurately identify bond rupture points and determine the associated forces. This is achieved by conducting a thorough comparison between molecular dynamics and cohesive finite element results. The utilization of a viscoelastic cohesive zone model allows for the consideration of loading rate effects. This rate-dependent model can be further developed and integrated into the multiscale modeling of large assemblies of metalloproteins, providing a comprehensive understanding of mechanical behavior while maintaining a reduced computational cost.

Interdomain Linker Effect on the Mechanical Stability of Ig Domains in Titin

Titin is the largest protein in humans, composed of more than one hundred immunoglobulin (Ig) domains, and plays a critical role in muscle’s passive elasticity. Thus, the molecular design of this giant polyprotein is responsible for its mechanical function. Interestingly, most of these Ig domains are connected directly with very few interdomain residues/linker, which suggests such a design is necessary for its mechanical stability. To understand this design, we chose six representative Ig domains in titin and added nine glycine residues (9G) as an artificial interdomain linker between these Ig domains. We measured their mechanical stabilities using atomic force microscopy-based single-molecule force spectroscopy (AFM-SMFS) and compared them to the natural sequence. The AFM results showed that the linker affected the mechanical stability of Ig domains. The linker mostly reduces its mechanical stability to a moderate extent, but the opposite situation can happen. Thus, this effect is very complex and may depend on each particular domain’s property.

Direct Measurements of the Cobalt-thiolate Bonds Strength in Rubredoxin by Single-Molecule Force Spectroscopy

Cobalt is a trace transition metal. Although it is not abundant on earth, tens of cobalt-containing proteins exist in life. Moreover, the characteristic spectrum of Co(II) ion makes it a powerful probe for the characterization of metal-binding proteins through the formation of cobalt-ligand bonds. Since most of these natural and artificial cobalt-containing proteins are stable, we believe that these cobalt-ligand bonds in the protein system are also mechanically stable. To prove this, we used atomic force microscopy-based single-molecule force spectroscopy (AFM-SMFS) to directly measure the rupture force of Co(II)-thiolate bond in Co-substituted rubredoxin (CoRD). By combining the chemical denature/renature method for building metalloprotein and cysteine coupling-based polyprotein construction strategy, we successfully prepared the polyprotein sample (CoRD) n suitable for single-molecule study. Thus, we quantified the strength of Co(II)-thiolate bonds in  rubredoxin with a rupture force of ~140 pN, revealing that the bond is a stable chemical bond. In addition, the Co-S bond is more labile than the Zn-S bond in proteins, similar to the result from the metal-competing titration experiment.

Highly Dynamic Polynuclear Metal Cluster Revealed in a Single Metallothionein Molecule

Human metallothionein (MT) is a small-size yet efficient metal-binding protein, playing an essential role in metal homeostasis and heavy metal detoxification. MT contains two domains, each forming a polynuclear metal cluster with an exquisite hexatomic ring structure. The apoprotein is intrinsically disordered, which may strongly influence the clusters and the metal-thiolate (M-S) bonds, leading to a highly dynamic structure. However, these features are challenging to identify due to the transient nature of these species. The individual signal from dynamic conformations with different states of the cluster and M-S bond will be averaged and blurred in classic ensemble measurement. To circumvent these problems, we combined a single-molecule approach and multiscale molecular simulations to investigate the rupture mechanism and chemical stability of the metal cluster by a single MT molecule, focusing on the Zn₄S₁₁ cluster in the α domain upon unfolding. Unusual multiple unfolding pathways and intermediates are observed for both domains, corresponding to different combinations of M-S bond rupture. None of the pathways is clearly preferred suggesting that unfolding proceeds from the distribution of protein conformational substates with similar M-S bond strengths. Simulations indicate that the metal cluster may rearrange, forming and breaking metal-thiolate bonds even when MT is folded independently of large protein backbone reconfiguration. Thus, a highly dynamic polynuclear metal cluster with multiple conformational states is revealed in MT, responsible for the binding promiscuity and diverse cellular functions of this metal-carrier protein.

Role of the copper ion in pseudoazurin during the mechanical unfolding process

Metalloproteins require the corresponding metal cofactors to exert their proper function. The presence of metal cofactors in the metalloprotein makes it more difficult to investigate its folding and unfolding process. In this study, we employed atomic-force-microscopy-based single-molecule force spectroscopy to reveal the unfolding process of pseudoazurin (PAZ) that belongs to blue copper proteins. Our study shows that holo-PAZ requires a higher rupture force for mechanical unfolding comparing with the apo-PAZ. This result demonstrates that the copper atom not only enables PAZ access to transfer electron, but should also have an influence on its stability. The results also suggest that the electronic configuration of the metal cofactors has a striking effect on the strength of the organometallic bonds. Moreover, the results also reveal that there is an intermediate state during the unfolding process of PAZ. This study provides insight into the characteristics of metalloproteins and leads to a better knowledge of their interaction at the individual molecule level.

Single-Molecule Desorption Studies of Poly(acrylic acid) at Electrolyte/Oxide/TiAlN Interfaces

The presented studies correlate the surface chemistry of electrochemically oxidized TiAlN hard coatings with the desorption forces of poly(acrylic acid) (PAA) at the electrolyte/oxide/TiAlN interface. Atomic force microscopy (AFM)-based single-molecule force spectroscopy (SMFS) was performed at different pH values to investigate surface chemistry-induced changes in desorption force. The chemical state was characterized by X-ray photoemission spectroscopy and electrochemical analysis. The results show that the desorption forces continuously decrease with increasing pH in the range from pH 5 to 9. The comparison of the desorption forces on rf-sputtered titanium dioxide and aluminum oxide films shows that the electrochemically oxidized surface of TiAlN, in agreement with the revealed surface composition, shows interfacial adhesive properties in contact with PAA and water that resemble a pure titanium oxide layer. Load rate-dependent measurements were performed to analyze both the free energy barrier and the transition state distance.

The rupture mechanism of rubredoxin is more complex than previously thought

The surprisingly low rupture force and remarkable mechanical anisotropy of rubredoxin have been known for several years. Exploiting the first combination of steered molecular dynamics and the quantum chemical Judgement of Energy DIstribution (JEDI) analysis, the common belief that hydrogen bonds between neighboring amino acid backbones and the sulfur atoms of the central FeS₄ unit in rubredoxin determine the low mechanical resistance of the protein is invalidated. The distribution of strain energy in the central part of rubredoxin is elucidated in real-time with unprecedented detail, giving important insights into the mechanical unfolding pathway of rubredoxin. While structural anisotropy as well as the contribution of angle bendings in the FeS₄ unit have a significant influence on the mechanical properties of rubredoxin, these factors are insufficient to explain the experimentally observed low rupture force. Instead, the rupture mechanism of rubredoxin is far more complex than previously thought and requires more than just a hydrogen bond network.

From Brazil to Germany: Challenges and Advantages

Moving to a new country, with a different culture and a new environment, is not an easy decision. In this perspective, I present some reasons that made me, a Brazilian computational biochemist, move abroad to do postdoctoral research and some of the challenges I faced before and after moving.

How strong are hydrogen bonds in the peptide model?

We provide a methodology based on the steered molecular dynamics simulations and dynamic force spectroscopy calculations to determine the kinetic and energetic characteristics of hydrogen bonds.

Modeling the Hydrolysis of Iron-Sulfur Clusters

Iron-sulfur (FeS) clusters are essential metal cofactors involved in a wide variety of biological functions. Their catalytic efficiency, biosynthesis, and regulation depend on FeS stability in aqueous solution. Here, molecular modeling is used to investigate the hydrolysis of an oxidized (ferric) mononuclear FeS cluster by bare dissociation and water substitution mechanisms in neutral and acidic solution. First, approximate electronic structure descriptions of FeS reactions by density functional theory are validated against high-level wave function CCSD(T) calculations. Solvation contributions are included by an all-atom model with hybrid quantum chemical/molecular mechanical (QM/MM) potentials and enhanced sampling molecular dynamics simulations. The free energy profile obtained for FeS cluster hydrolysis indicates that the hybrid functional M06 together with an implicit solvent correction capture the most important aspects of FeS cluster reactivity in aqueous solution. Then, 20 reaction channels leading to two consecutive Fe-S bond ruptures were explored with this calibrated model. For all protonation states, nucleophilic substitution with concerted bond breaking and forming to iron is the preferred mechanism, both kinetic and thermodynamically. In neutral solution, proton transfer from water to the sulfur leaving group is also concerted. Dissociative reactions show higher barriers and will not be relevant for FeS reactivity when exposed to solvent. These hydrolysis mechanisms may help to explain the stability and catalytic mechanisms of FeS clusters of multiple sizes and proteins.

Single molecule force spectroscopy of a streptomycin-binding RNA aptamer: An out-of-equilibrium molecular dynamics study

Here, we investigate the unfolding behavior of a streptomycin-binding ribonucleic acid (RNA) aptamer under application of force in shear geometry. Using Langevin out-of-equilibrium simulations to emulate the single-molecule force spectroscopy (SMFS) experiment, we were able to understand the hierarchical unfolding process that occurs in the RNA molecule under application of stretching force and the influence of streptomycin modifying this unfolding. Subsequently, the application of the Jarzynski equality to the force profiles obtained in the pulling simulations shows that the free energies for individual systems and the difference of unfolding free energy upon streptomycin binding to the RNA free aptamer are in fair agreement with the experimental values, obtained through SMFS by Nick et al. [J. Phys. Chem. B 120, 6479 (2016)].

Multistep Protein Unfolding Scenarios from the Rupture of a Complex Metal Cluster Cd3S9

Protein (un)folding is a complex and essential process. With the rapid development of single-molecule techniques, we can detect multiple and transient proteins (un)folding pathways/intermediates. However, the observation of multiple multistep (>2) unfolding scenarios for a single protein domain remains limited. Here, we chose metalloprotein with relatively stable and multiple metal-ligand coordination bonds as a system for such a purpose. Using AFM-based single-molecule force spectroscopy (SMFS), we successfully demonstrated the complex and multistep protein unfolding scenarios of the β-domain of a human protein metallothionein-3 (MT). MT is a protein of ~60 amino acids (aa) in length with 20 cysteines for various metal binding, and the β-domain (βMT) is of ~30 aa with an M3S9 metal cluster. We detected four different types of three-step protein unfolding scenarios from the Cd-βMT, which can be possibly explained by the rupture of Cd-S bonds in the complex Cd3S9 metal cluster. In addition, complex unfolding scenarios with four rupture peaks were observed. The Cd-S bonds ruptured in both single bond and multiple bonds modes. Our results provide not only evidence for multistep protein unfolding phenomena but also reveal unique properties of metalloprotein system using single-molecule AFM.

Mechanical Unfolding Pathway of the High-Potential Iron-Sulfur Protein Revealed by Single-Molecule Atomic Force Microscopy: Toward a General Unfolding Mechanism for Iron-sulfur Proteins

High-potential iron-sulfur proteins (HiPIPs) are an important class of metalloproteins with a [4Fe-4S] cluster coordinated by four cysteine residues. Distinct from other iron-sulfur proteins, the cluster in HiPIP has a high reduction potential, making it an essential electron carrier in bacterial photosynthesis. Here, we combined single-molecule atomic force microscopy and protein engineering techniques to investigate the mechanical unfolding mechanism of HiPIP from Chromatium tepidum (cHiPIP). We found that cHiPIP unfolds in a two-step fashion with the protein sequence sequestered by the iron-sulfur center as a stable unfolding intermediate state. The rupture of the iron-sulfur center of cHiPIP proceeds in two distinct parallel pathways; one pathway involves the concurrent rupture of multiple iron-thiolate bonds, and the other one involves the sequential rupture of the iron-thiolate bonds. This mechanistic information was further confirmed by mutational studies. We found that the rupture of the iron-thiolate bonds in reduced and oxidized cHiPIP occurred in the range of 150-180 pN at a pulling speed of 400 nm/s, similar to that measured for iron-thiolate bonds in rubredoxin and ferredoxin. Our results may have important implications for understanding the general unfolding mechanism governing iron-sulfur proteins, as well as the mechanism governing the mechanical rupture of the iron-sulfur center.

Cell-cell adhesion interface: orthogonal and parallel forces from contraction, protrusion, and retraction

The epithelial lateral membrane plays a central role in the integration of intercellular signals and, by doing so, is a principal determinant in the emerging properties of epithelial tissues. Mechanical force, when applied to the lateral cell-cell interface, can modulate the strength of adhesion and influence intercellular dynamics. Yet the relationship between mechanical force and epithelial cell behavior is complex and not completely understood. This commentary aims to provide an investigative look at the usage of cellular forces at the epithelial cell-cell adhesion interface.

Simulated and experimental force spectroscopy of lysozyme on silica

The force spectra of proteins detaching from oxide surfaces measured by atomic force microscopy (AFM) often present complex patterns of peaks, which are difficult to correlate with individual bond-breaking events at the atomic scale. In this work we rationalize experimental AFM force spectra of hen-egg-white lysozyme detaching from silica by means of all-atom steered molecular dynamics (SMD) simulations. In particular, we demonstrate that the native tertiary structure of lysozyme is preserved if, and only if, its four intramolecular disulfide bridges are intact. Otherwise, the protein pulled off the surface undergoes severe unfolding, which is well captured by SMD simulations in explicit solvent. Implicit solvent simulations, on the contrary, wrongly predict protein unfolding even in the presence of S-S bridges, due to the lack of additional structural stabilization provided by the water's hydrogen-bond network within and surrounding the protein. On the basis of our combined experimental and theoretical findings, we infer that the rugged force spectra characteristic of lysozyme/silica interfaces are not due to the successive breaking of internal disulfide bonds leading to partial unfolding events. Rather, they reflect the detachment of several molecules bound to the same AFM tip, each anchored to the surface via multiple hydrogen and ionic bonds.