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

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.

Molecular Mechanism of Interaction between DNA Aptamer and Receptor-Binding Domain of Severe Acute Respiratory Syndrome Coronavirus 2 Variants Revealed by Steered Molecular Dynamics Simulations

The ongoing SARS-CoV-2 pandemic has underscored the urgent need for versatile and rapidly deployable antiviral strategies. While vaccines have been pivotal in controlling the spread of the virus, the emergence of new variants continues to pose significant challenges to global health. Here, our study focuses on a novel approach to antiviral therapy using DNA aptamers, short oligonucleotides with high specificity and affinity for their targets, as potential inhibitors against the spike protein of SARS-CoV-2 variants Omicron and JN.1. Our research utilizes steered molecular dynamics (SMD) simulations to elucidate the binding mechanisms of a specifically designed DNA aptamer, AM032-4, to the receptor-binding domain (RBD) of the aforementioned variants. The simulations reveal detailed molecular insights into the aptamer-RBD interaction, demonstrating the aptamer's potential to maintain effective binding in the face of rapid viral evolution. Our work not only demonstrates the dynamic interaction between aptamer-RBD for possible antiviral therapy but also introduces a computational method to study aptamer-protein interactions.

Xenobiotic Bi3+ Coordination by Cysteine-Rich Metallothionein-3 Reveals a Cooperatively Formed Thiolate-Sharing Bi2S5 Cluster

The field of designing artificial metalloproteins has yet to effectively tackle the incorporation of multimetal clusters, which is a key component of natural metalloproteins, such as metallothioneins (MTs) and calmodulin. MT is a physiological, essential, cysteine-rich metalloprotein that binds to a variety of metals but is only known to form metal-thiolate clusters with Cd2+, Zn2+, and Cu+. Bismuth is a xenobiotic metal and a component of metallodrugs used to treat gastric ulcers and cancer, as well as an emerging metal used in industrial practices. Electrospray ionization mass spectrometry, UV-visible spectroscopy, and extended X-ray absorption fine structure spectroscopy were used to probe the Bi3+ binding site structures in apo-MT3 (brain-located MT) at pH 7.4 and 2 and provide the complete set of binding affinities. We discovered the highly cooperative formation of a novel Bi3+ species, Bi2MT3, under physiological conditions, where each Bi3+ ion is coordinated by three cysteinyl thiolates, with one of the thiolates bridging between the two Bi3+ ions. This cluster structure was associated with a strong visible region absorption band, which was disrupted by the addition of Zn2+ and reversibly disrupted by acidification and increased temperature. This is the first reported presence of bridging cysteines for a xenobiotic metal in MT3 and the Bi2MT structure is the first Bi cluster found in a metalloprotein.

Quantification of carboxylate-bridged di-zinc site stability in protein due ferri by single-molecule force spectroscopy

Carboxylate-bridged diiron proteins belong to a protein family involved in different physiological processes. These proteins share the conservative EXXH motif, which provides the carboxylate bridge and is critical for metal binding. Here, we choose de novo-designed single-chain due ferri protein (DFsc), a four-helical protein with two EXXH motifs as a model protein, to study the stability of the carboxylate-bridged di-metal binding site. The mechanical and kinetic properties of the di-Zn site in DFsc were obtained by atomic force microscopy-based single-molecule force spectroscopy. Zn-DFsc showed a considerable rupture force of ~200 pN, while the apo-protein is mechanically labile. In addition, multiple rupture pathways were observed with different probabilities, indicating the importance of the EXXH-based carboxylate-bridged metal site. These results demonstrate carboxylate-bridged di-metal site is mechanically stable and improve our understanding of this important type of metalloprotein.

Enzymatic Protein Immobilization on Amino-Functionalized Nanoparticles

The immobilization of proteins on nanoparticles has received much attention in recent years. Among different approaches, enzymatic protein immobilization shows unique advantages because of its site-specific connection. OaAEP1 is a recently engineered peptide ligase which can specifically recognize an N-terminal GL residue (NH₂-Gly-Leu) and a C-terminal NGL amino acid residue (Asn-Gly-Leu-COOH) and ligates them efficiently. Herein, we report OaAEP1-mediated protein immobilization on synthetic magnetic nanoparticles. Our work showed that OaAEP1 could mediate C-terminal site-specific protein immobilization on the amino-functionalized Fe₃O₄ nanoparticles. Our work demonstrates a new method for site-specific protein immobilization on nanoparticles.

One-step asparaginyl endopeptidase (OaAEP1)-based protein immobilization for single-molecule force spectroscopy

Enzymatic protein ligation has become the most powerful and widely used method for high-precision atomic force microscopy single-molecule force spectroscopy (AFM-SMFS) study of protein mechanics. However, this methodology typically requires the functionalization of the glass surface with a corresponding peptide sequence/tag for enzymatic recognition and multiple steps are needed. Thus, it is time-consuming and a high level of experience is needed for reliable results. To solve this problem, we simplified the procedure using two strategies both based on asparaginyl endopeptidase (AEP). First, we designed a heterobifunctional peptide-based crosslinker, GL-peptide-propargylglycine, which links to an N 3-functionalized surface via the click reaction. Then, the target protein with a C-terminal NGL sequence can be immobilized via the AEP-mediated ligation. Furthermore, we took advantage of the direct ligation between primary amino in a small molecule and protein with C-terminal NGL by AEP. Thus, the target protein can be immobilized on an amino-functionalized surface via AEP in one step. Both approaches were successfully applied to the AFM-SMFS study of eGFP, showing consistent single-molecule results.

Single-Molecule Force Spectroscopy Reveals Stability of mitoNEET and its [2Fe2Se] Cluster in Weakly Acidic and Basic Solutions

The outer mitochondrial membrane protein mitoNEET (mNT) is a recently identified iron-sulfur protein containing a unique Fe₂ S₂ (His)₁ (Cys)₃ metal cluster with a single Fe-N(His87) coordinating bond. This labile Fe-N bond led to multiple unfolding/rupture pathways of mNT and its cluster by atomic force microscopy-based single-molecule force spectroscopy (AFM-SMFS), one of most common tools for characterizing the molecular mechanics. Although previous ensemble studies showed that this labile Fe-N(His) bond is essential for protein function, they also indicated that the protein and its [2Fe2S] cluster are stable under acidic conditions. Thus, we applied AFM-SMFS to measure the stability of mNT and its cluster at pH values of 6, 7, and 8. Indeed, all previous multiple unfolding pathways of mNT were still observed. Moreover, single-molecule measurements revealed that the stabilities of the protein and the [2Fe2S] cluster are consistent at these pH values with only ≈20 pN force differences. Thus, we found that the behavior of the protein is consistent in both weakly acidic and basic solutions despite a labile Fe-N bond.

Enzymatic Protein-Protein Conjugation through Internal Site Verified at the Single-Molecule Level

Enzymes are widely used for protein ligation because of their efficient and site-specific connections under mild conditions. However, many enzymatic ligations are restricted to connections between protein termini while protein-protein conjugation at a specific internal site is limited. Previous work has found that Sortase A (SrtA) conjugates small molecules/peptides to a pilin protein at an internal lysine site via an isopeptide bond. Herein, we apply this noncanonical ligation property of SrtA for protein-protein conjugation at a designed YPKH site. Both a small protein domain, I27, and a large protein, GFP, were ligated at the designed internal site. Moreover, besides characterization by classic methods at the ensemble level, the specific ligation site at the interior YPKH motif is unambiguously verified by atomic force microscopy-based single-molecule force spectroscopy, showing the characteristic unfolding signature at the single-molecule level. Finally, steered molecular dynamics simulations also agreed with the results.

Facile Synthesis of Peptide-Conjugated Gold Nanoclusters with Different Lengths

The synthesis of ultra-small gold nanoclusters (Au NCs) with sizes down to 2 nm has received increasing interest due to their unique optical and electronic properties. Like many peptide-coated gold nanospheres synthesized before, modified gold nanoclusters with peptide conjugation are potentially significant in biomedical and catalytic fields. Here, we explore whether such small-sized gold nanoclusters can be conjugated with peptides also and characterize them using atomic force microscopy. Using a long and flexible elastin-like polypeptide (ELP)₂₀ as the conjugated peptide, (ELP)₂₀-Au NCs was successfully synthesized via a one-pot synthesis method. The unique optical and electronic properties of gold nanoclusters are still preserved, while a much larger size was obtained as expected due to the peptide conjugation. In addition, a short and rigid peptide (EAAAK)₃ was conjugated to the gold nanoclusters. Their Yong’s modulus was characterized using atomic force microscopy (AFM). Moreover, the coated peptide on the nanoclusters was pulled using AFM-based single molecule-force spectroscopy (SMFS), showing expected properties as one of the first force spectroscopy experiments on peptide-coated nanoclusters. Our results pave the way for further modification of nanoclusters based on the conjugated peptides and show a new method to characterize these materials using AFM-SMFS.

An Integrated Mass Spectrometry and Molecular Dynamics Simulations Approach Reveals the Spatial Organization Impact of Metal-Binding Sites on the Stability of Metal-Depleted Metallothionein-2 Species

Mammalian metallothioneins (MTs) are a group of cysteine-rich proteins that bind metal ions in two α- and β-domains and represent a major cellular Zn(II)/Cu(I) buffering system in the cell. At cellular free Zn(II) concentrations (10-11-10-9 M), MTs do not exist in fully loaded forms with seven Zn(II)-bound ions (Zn7MTs). Instead, MTs exist as partially metal-depleted species (Zn4-6MT) because their Zn(II) binding affinities are on the nano- to picomolar range comparable to the concentrations of cellular Zn(II). The mode of action of MTs remains poorly understood, and thus, the aim of this study is to characterize the mechanism of Zn(II) (un)binding to MTs, the thermodynamic properties of the Zn1-6MT2 species, and their mechanostability properties. To this end, native mass spectrometry (MS) and label-free quantitative bottom-up and top-down MS in combination with steered molecular dynamics simulations, well-tempered metadynamics (WT-MetaD), and parallel-bias WT-MetaD (amounting to 3.5 μs) were integrated to unravel the chemical coordination of Zn(II) in all Zn1-6MT2 species and to explain the differences in binding affinities of Zn(II) ions to MTs. Differences are found to be the result of the degree of water participation in MT (un)folding and the hyper-reactive character of Cys21 and Cys29 residues. The thermodynamics properties of Zn(II) (un)binding to MT2 are found to differ from those of Cd(II), justifying their distinctive roles. The potential of this integrated strategy in the investigation of numerous unexplored metalloproteins is attested by the results highlighted in the present study.

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.

From folding to function: complex macromolecular reactions unraveled one-by-one with optical tweezers

Single-molecule manipulation with optical tweezers has uncovered macromolecular behaviour hidden to other experimental techniques. Recent instrumental improvements have made it possible to expand the range of systems accessible to optical tweezers. Beyond focusing on the folding and structural changes of isolated single molecules, optical tweezers studies have evolved into unraveling the basic principles of complex molecular processes such as co-translational folding on the ribosome, kinase activation dynamics, ligand-receptor binding, chaperone-assisted protein folding, and even dynamics of intrinsically disordered proteins (IDPs). In this mini-review, we illustrate the methodological principles of optical tweezers before highlighting recent advances in studying complex protein conformational dynamics - from protein synthesis to physiological function - as well as emerging future issues that are beginning to be addressed with novel approaches.

Effects of Ligand Binding on the Energy Landscape of Acyl-CoA-Binding Protein

Binding of ligands is often crucial for function yet the effects of ligand binding on the mechanical stability and energy landscape of proteins are incompletely understood. Here, we use a combination of single-molecule optical tweezers and MD simulations to investigate the effect of ligand binding on the energy landscape of acyl-coenzyme A (CoA)-binding protein (ACBP). ACBP is a topologically simple and highly conserved four-α-helix bundle protein that acts as an intracellular transporter and buffer for fatty-acyl-CoA and is active in membrane assembly. We have previously described the behavior of ACBP under tension, revealing a highly extended transition state (TS) located almost halfway between the unfolded and native states. Here, we performed force-ramp and force-jump experiments, in combination with advanced statistical analysis, to show that octanoyl-CoA binding increases the activation free energy for the unfolding reaction of ACBP without affecting the position of the transition state along the reaction coordinate. It follows that ligand binding enhances the mechanical resistance and thermodynamic stability of the protein, without changing its mechanical compliance. Steered molecular dynamics simulations allowed us to rationalize the results in terms of key interactions that octanoyl-CoA establishes with the four α-helices of ACBP and showed that the unfolding pathway is marginally affected by the ligand. The results show that ligand-induced mechanical stabilization effects can be complex and may prove useful for the rational design of stabilizing ligands.

Blue light-induced low mechanical stability of ruthenium-based coordination bonds: an AFM-based single-molecule force spectroscopy study

A HA–Ru^II complex was conjugated to a hyaluronan polymer through amide bonds. In AFM experiments using the “multi-fishhook” approach, the cantilever tip made contact with the polymeric molecule, resulting in stretching, indicated by sawtooth-like force-extension curves.