Force-based Analysis of Multidimensional Energy Landscapes: Application of Dynamic Force Spectroscopy and Steered Molecular Dynamics Simulations to an Antibody Fragment–Peptide Complex

Structural basis of the tarocystatin inhibitory mechanism against papain

Plant pathogens pose a severe threat to global food security by compromising the availability, quality, and safety of crops for human and animal consumption. Given the urgent need for alternatives to chemical pesticides, natural inhibitors of phytopathogenic proteases represent promising biopesticides. Tarocystatin has been characterized in Colocasia esculenta as a defense protein against phytopathogenic nematodes and fungi. Despite its biotechnological potential, few studies describe its mechanical, structural, and energetic properties. In this study, we characterized the inhibitory mechanism of tarocystatin against papain using a computational biophysics approach. Through extensive molecular dynamics (MD) and steered molecular dynamics (SMD) simulations, we explored the dynamic, energetic, structural, and mechanical basis of tarocystatin and its specific binding to papain. Our results suggest that the stability of the complex is characterized by a lack of conformational rearrangements, showing invariability in its secondary structure across all MD replicas. Additionally, electrostatic analysis revealed a high complementarity of the tarocystatin-papain complex, which was later corroborated by the hydrogen-bond network established at the complex interface, explaining its strong inhibitory capacity. Moreover, we determined that the substrate-competitive inhibitory mechanism is due to the binding ability of conserved motifs in tarocystatin, which efficiently interact with the catalytic active site of papain. This was also confirmed through SMD, where we observed that the N-terminal region acts as a spring to prevent the dissociation of the complex under external pulling forces. Overall, our study is the first to provide a comprehensive exploration of the biophysical properties of the tarocystatin-papain complex, offering insights into the tarocystatin's inhibition mechanism. These results lay the foundation for future development of tarocystatin-based antifungal alternatives, as well as for exploring its inhibitory activity in other pathogens or enhancing its efficacy through molecular engineering.

Mechanical control of antigen detection and discrimination by T and B cell receptors

The adaptive immune response is orchestrated by just two cell types, T cells and B cells. Both cells possess the remarkable ability to recognize virtually any antigen through their respective antigen receptors-the T cell receptor (TCR) and B cell receptor (BCR). Despite extensive investigations into the biochemical signaling events triggered by antigen recognition in these cells, our ability to predict or control the outcome of T and B cell activation remains elusive. This challenge is compounded by the sensitivity of T and B cells to the biophysical properties of antigens and the cells presenting them-a phenomenon we are just beginning to understand. Recent insights underscore the central role of mechanical forces in this process, governing the conformation, signaling activity, and spatial organization of TCRs and BCRs within the cell membrane, ultimately eliciting distinct cellular responses. Traditionally, T cells and B cells have been studied independently, with researchers working in parallel to decipher the mechanisms of activation. While these investigations have unveiled many overlaps in how these cell types sense and respond to antigens, notable differences exist. To fully grasp their biology and harness it for therapeutic purposes, these distinctions must be considered. This review compares and contrasts the TCR and BCR, placing emphasis on the role of mechanical force in regulating the activity of both receptors to shape cellular and humoral adaptive immune responses.

Connections between efficient control and spontaneous transitions in an Ising model

A system can be driven between metastable configurations by a time-dependent driving protocol, which uses external control parameters to change the potential energy of the system. Here we investigate the correspondence between driving protocols that are designed to minimize work and the spontaneous transition paths of the system in the absence of driving. We study the spin-inversion reaction in a 2D Ising model, quantifying the timing of each spin flip and heat flow to the system during both a minimum-work protocol and a spontaneous transition. The general order of spin flips during the transition mechanism is preserved between the processes, despite the coarseness of control parameters that are unable to reproduce more detailed features of the spontaneous mechanism. Additionally, external control parameters provide energy to each system component to compensate changes in internal energy, showing how control parameters are tuned during a minimum-work protocol to counteract underlying energetic features. This paper supports a correspondence between minimum-work protocols and spontaneous transition mechanisms.

Multidimensional minimum-work control of a 2D Ising model

A system's configurational state can be manipulated using dynamic variation of control parameters, such as temperature, pressure, or magnetic field; for finite-duration driving, excess work is required above the equilibrium free-energy change. Minimum-work protocols in multidimensional control-parameter space have potential to significantly reduce work relative to one-dimensional control. By numerically minimizing a linear-response approximation to the excess work, we design protocols in control-parameter spaces of a 2D Ising model that efficiently drive the system from the all-down to all-up configuration. We find that such designed multidimensional protocols take advantage of more flexible control to avoid control-parameter regions of high system resistance, heterogeneously input and extract work to make use of system relaxation, and flatten the energy landscape, making accessible many configurations that would otherwise have prohibitively high energy and thus decreasing spin correlations. Relative to one-dimensional protocols, this speeds up the rate-limiting spin-inversion reaction, thereby keeping the system significantly closer to equilibrium for a wide range of protocol durations, and significantly reducing resistance and hence work.

Next Generation Methods for Single-Molecule Force Spectroscopy on Polyproteins and Receptor-Ligand Complexes

Single-molecule force spectroscopy with the atomic force microscope provides molecular level insights into protein function, allowing researchers to reconstruct energy landscapes and understand functional mechanisms in biology. With steadily advancing methods, this technique has greatly accelerated our understanding of force transduction, mechanical deformation, and mechanostability within single- and multi-domain polyproteins, and receptor-ligand complexes. In this focused review, we summarize the state of the art in terms of methodology and highlight recent methodological improvements for AFM-SMFS experiments, including developments in surface chemistry, considerations for protein engineering, as well as theory and algorithms for data analysis. We hope that by condensing and disseminating these methods, they can assist the community in improving data yield, reliability, and throughput and thereby enhance the information that researchers can extract from such experiments. These leading edge methods for AFM-SMFS will serve as a groundwork for researchers cognizant of its current limitations who seek to improve the technique in the future for in-depth studies of molecular biomechanics.

Mechanobiology: protein refolding under force

The application of direct force to a protein enables to probe wide regions of its energy surface through conformational transitions as unfolding, extending, recoiling, collapsing, and refolding. While unfolding under force typically displayed a two-state behavior, refolding under force, from highly extended unfolded states, displayed a more complex behavior. The first recording of protein refolding at a force quench step displayed an initial rapid elastic recoil, followed by a plateau phase at some extension, concluding with a collapse to a final state, at which refolding occurred. These findings stirred a lively discussion, which led to further experimental and theoretical investigation of this behavior. It was demonstrated that the polymeric chain of the unfolded protein is required to fully collapse to a globular conformation for the maturation of native structure. This behavior was modeled using one-dimensional free energy landscape over the end-to-end length reaction coordinate, the collective measured variable. However, at low forces, conformational space is not well captured by such models, and using two-dimensional energy surfaces provides further insight into the dynamics of this process. This work reviews the main concepts of protein refolding under constant force, which is essential for understanding how mechanotransducing proteins operate in vivo.

Designed inhibitors with hetero linkers for gastric proton pump H+,K+-ATPase: Steered molecular dynamics and metadynamics studies

Acid suppressant SCH28080 and its derivatives reversibly reduce acid secretion activity of the H⁺,K⁺-ATPase in a K⁺ competitive manner. The results on homologation of the SCH28080 by varying the linker chain length suggested the improvement in efficacy. However, the pharmacokinetic studies reveal that the hydrophobic nature of the CH₂ linker units may not help it to function as a better acid suppressant. We have exploited the role of linker unit to enhance the efficacy of such reversible acid suppressant drug molecules using hetero linker, i.e., disulfide and peroxy linkers. The logarithm of partition coefficient defined for a drug molecule relates to the partition coefficient, which allows the optimum solubility characteristics to reach the active site. The logarithm of partition coefficient calculated for the designed inhibitors suggests that inhibitors would possibly reach the active site in sufficient concentration like in the case of SCH28080. The steered molecular dynamics studies have revealed that the Inhibitor-1 with disulfide linker unit is more stable at the active site due to greater noncovalent interactions compared to the SCH28080. Centre of mass distance analysis suggests that the Cysteine-813 amino acid residue selectively plays an important role in the inhibition of H⁺,K⁺-ATPase for Inhibitor-1. Furthermore, the quantum chemical calculations with M11L/6-31+G(d,p) level of theory have been performed to account the noncovalent interactions responsible for the stabilization of inhibitor molecules in the active site gorge of the gastric proton pump at different time scale. The hydrogen bonding and hydrophobic interaction studies corroborate the center of mass distance analysis as well. Well-tempered metadynamics free energy surface and center of mass separation analysis for the Inhibitor-1 is in good agreement with the steered molecular dynamics results. The torsional angle of the linker units seems to be crucial for better efficacy of drug molecules. The torsional angle of linker units of SCH28080 (COCH₂C) and of Inhibitor 1 (CSSC) prefers to lie within ∼60°-90° for a longer time during the simulations, whereas, the peroxy linker (COOC) of Inhibitor 2 prefers to adopt ∼120-160°. Therefore, it appears that the smaller torsion angle of linker units can achieve better interactions with the active site residues of H⁺,K⁺-ATPase to inhibit the acid secretion activity. The reversible drug molecules with disulfide linker unit would be a promising candidate as proton pump antagonist to H⁺,K⁺-ATPase.

Lagrangian descriptors of driven chemical reaction manifolds

The persistence of a transition state structure in systems driven by time-dependent environments allows the application of modern reaction rate theories to solution-phase and nonequilibrium chemical reactions. However, identifying this structure is problematic in driven systems and has been limited by theories built on series expansion about a saddle point. Recently, it has been shown that to obtain formally exact rates for reactions in thermal environments, a transition state trajectory must be constructed. Here, using optimized Lagrangian descriptors [G. T. Craven and R. Hernandez, Phys. Rev. Lett. 115, 148301 (2015)PRLTAO0031-900710.1103/PhysRevLett.115.148301], we obtain this so-called distinguished trajectory and the associated moving reaction manifolds on model energy surfaces subject to various driving and dissipative conditions. In particular, we demonstrate that this is exact for harmonic barriers in one dimension and this verification gives impetus to the application of Lagrangian descriptor-based methods in diverse classes of chemical reactions. The development of these objects is paramount in the theory of reaction dynamics as the transition state structure and its underlying network of manifolds directly dictate reactivity and selectivity.

Post-Translational Sortase-Mediated Attachment of High-Strength Force Spectroscopy Handles

Single-molecule force spectroscopy greatly benefits from site-specific surface immobilization and specific probing with a functionalized cantilever. Here, we describe a streamlined approach to such experiments by covalently attaching mechanically stable receptors onto proteins of interest (POI) to improve pickup efficiency and specificity. This platform provides improved throughput, allows precise control over the pulling geometry, and allows for multiple constructs to be probed with the same ligand-modified cantilever. We employ two orthogonal enzymatic ligation reactions [sortase and phosphopantetheinyl transferase (Sfp)] to covalently immobilize POI to a pegylated surface and to subsequently ligate the POI to a mechanically stable dockerin domain at the protein's C-terminus for use as a high-strength pulling handle. Our configuration permits expression and folding of the POI to proceed independently from the mechanically stable receptor used for specific probing and requires only two short terminal peptide sequences (i.e., ybbR-tag and sortase C-tag). We applied this system successfully to proteins expressed using in vitro transcription and translation reactions without a protein purification step and to purified proteins expressed in Escherichia coli.

Revealing the importance of linkers in K-series oxime reactivators for tabun-inhibited AChE using quantum chemical, docking and SMD studies

Inhibition of acetylcholinesterase (AChE) with organophosphorus compounds has a detrimental effect on human life. Oxime K203 seems to be one of the promising reactivators for tabun-inhibited AChE than (K027, K127, and K628). These reactivators differ only in the linker units between the two pyridinium rings. The conformational analyses performed with quantum chemical RHF/6-31G* level for K027, K127, K203 and K628 showed that the minimum energy conformers have different orientations of the active and peripheral pyridinium rings for these reactivator molecules. K203 with (-CH₂-CH=CH-CH₂-) linker unit possesses more open conformation compared to the other reactivators. Such orientation of K203 experiences favorable interaction with the surrounding residues of catalytic anionic site (CAS) and peripheral anionic site (PAS) of tabun-inhibited AChE. From the steered molecular dynamics simulations, it has been observed that the oxygen atom of the oxime group of K203 reactivator approaches nearest to the P-atom of the SUN203 (3.75 Å) at lower time scales (less than ~1000 ps) as compared to the other reactivators. K203 experiences less number of hydrophobic interaction with the PAS residues which is suggested to be an important factor for the efficient reactivation process. In addition, K203 crates large number of H-bonding with CAS residues SUN203, Phe295, Tyr337, Phe338 and His447. K203 barely changes its conformation during the SMD simulation process and hence the energy penalty to adopt any other conformation is minimal in this case as compared to the other reactivators. The molecular mechanics and Poisson-Boltzmann surface area binding energies obtained for the interaction of K203 inside the gorge of tabun inhibited AChE is substantially higher (-290.2 kcal/mol) than the corresponding K628 reactivator (-260.4 kcal/mol), which also possess unsaturated aromatic linker unit.

Simulated Force Quench Dynamics Shows GB1 Protein Is Not a Two State Folder

Single molecule force spectroscopy is a useful technique for investigating mechanically induced protein unfolding and refolding under reduced forces by monitoring the end-to-end distance of the protein. The data is often interpreted via a "two-state" model based on the assumption that the end-to-end distance alone is a good reaction coordinate and the thermodynamic behavior is then ascribed to the free energy as a function of this one reaction coordinate. In this paper, we determined the free energy surface (PMF) of GB1 protein from atomistic simulations in explicit solvent under different applied forces as a function of two collective variables (the end-to-end-distance, and the fraction of native contacts ρ). The calculated 2-d free energy surfaces exhibited several distinct states, or basins, mostly visible along the ρ coordinate. Brownian dynamics (BD) simulations on the smoothed free energy surface show that the protein visits a metastable molten globule state and is thus a three state folder, not the two state folder inferred using the end-to-end distance as the sole reaction coordinate. This study lends support to recent experiments that suggest that GB1 is not a two-state folder.

Single versus dual-binding conformations in cellulosomal cohesin-dockerin complexes

Cohesins and dockerins are complementary interacting protein modules that form stable and highly specific receptor-ligand complexes. They play a crucial role in the assembly of cellulose-degrading multi-enzyme complexes called cellulosomes and have potential applicability in several technology areas, including biomass conversion processes. Here, we describe several exceptional properties of cohesin-dockerin complexes, including their tenacious biochemical affinity, remarkably high mechanostability and a dual-binding mode of recognition that is contrary to the conventional lock-and-key model of receptor-ligand interactions. We focus on structural aspects of the dual mode of cohesin-dockerin binding, highlighting recent single-molecule analysis techniques for its explicit characterization.

The reactivation of tabun-inhibited mutant AChE with Ortho-7: steered molecular dynamics and quantum chemical studies

A highly toxic nerve agent, tabun, can inhibit acetylcholinesterase (AChE) at cholinergic sites, which leads to serious cardiovascular complications, respiratory compromise and death. We have examined the structural features of the tabun-conjugated AChE complex with an oxime reactivator, Ortho-7, to provide a strategy for designing new and efficient reactivators. Mutation of mAChE within the choline binding site by Y337A and F338A and its interaction with Ortho-7 has been investigated using steered molecular dynamics (SMD) and quantum chemical methods. The overall study shows that after mutagenesis (Y337A), the reactivator can approach more freely towards the phosphorylated active site of serine without any significant steric hindrance in the presence of tabun compared to the wild type and double mutant. Furthermore, the poor binding of Ortho-7 with the peripheral residues of mAChE in the case of the single mutant compared to that of the wild-type and double mutant (Y337A/F338A) can contribute to better efficacy in the former case. Ortho-7 has formed a greater number of hydrogen bonds with the active site surrounding residues His447 and Phe295 in the case of the single mutant (Y337A), and that stabilizes the drug molecule for an effective reactivation process. The DFT M05-2X/6-31+G(d) level of theory shows that the binding energy of Ortho-7 with the single mutant (Y337A) is energetically more preferred (-19.8 kcal mol(-1)) than the wild-type (-8.1 kcal mol(-1)) and double mutant (Y337A/F338A) (-16.0 kcal mol(-1)). The study reveals that both the orientation of the oxime reactivator for nucleophilic attack and the stabilization of the reactivator at the active site would be crucial for the design of an efficient reactivator.

Deconstructing field-induced ketene isomerization through Lagrangian descriptors

Phase space contours (shown in color) constructed using the method of Lagrangian descriptors resolve the separatrices governing state transitions on the reaction-path potential energy surface (shown in white) for field-induced ketene isomerization.

Can hydroxylamine be a more potent nucleophile for the reactivation of tabun-inhibited AChE than prototype oxime drugs? An answer derived from quantum chemical and steered molecular dynamics studies

Organophosphorus nerve agents are highly toxic compounds which strongly inhibit acetylcholinesterase (AChE) in the blood and in the central nervous system (CNS). Tabun is one of the highly toxic organophosphorus (OP) compounds and is resistant to many oxime drugs formulated for the reactivation of AChE. The reactivation mechanism of tabun-conjugated AChE with various drugs has been examined with density functional theory and ab initio quantum chemical calculations. The presence of a lone-pair located on the amidic group resists the nucleophilic attack at the phosphorus center of the tabun-conjugated AChE. We have shown that the newly designed drug candidate N-(pyridin-2-yl)hydroxylamine, at the MP2/6-31+G*//M05-2X/6-31G* level in the aqueous phase with the polarizable continuum solvation model (PCM), is more effective in reactivating the tabun-conjugated AChE than typical oxime drugs. The rate determining activation barrier with N-(pyridin-2-yl)hydroxylamine was found to be ∼1.7 kcal mol(-1), which is 7.2 kcal mol(-1) lower than the charged oxime trimedoxime (one of the most efficient reactivators in tabun poisonings). The greater nucleophilicity index (ω(-)) and higher CHelpG charge of pyridinylhydroxylamine compared to TMB4 support this observation. Furthermore, we have also examined the reactivation process of tabun-inhibited AChE with some other bis-quaternary oxime drug candidates such as methoxime (MMB4) and obidoxime. The docking analysis suggests that charged bis-quaternary pyridinium oximes have greater binding affinity inside the active-site gorge of AChE compared to the neutral pyridinylhydroxylamine. The peripheral ligand attached to the neutral pyridinylhydroxylamine enhanced the binding with the aromatic residues in the active-site gorge of AChE through effective π-π interactions. Steered molecular dynamics (SMD) simulations have also been performed with the charged oxime (TMB4) and the neutral hydroxylamine. From protein-drug interaction parameters (rupture force profiles, hydrogen bonds, hydrophobic interactions), geometry and the orientation of the drug candidates, the hydroxylamine is suggested to orchestrate the reactivation process better than TMB4. Furthermore, the calculated log P values show the effective penetration of the neutral drug candidate through the blood-brain barrier. The toxicity measurements and the IC50 values (a measure of the intrinsic affinity toward AChE) suggest that the pyridinylhydroxylamine compound could have similar toxic behavior compared to the prototype oxime antidotes used for reactivation purposes. The newly designed pyridinylhydroxylamine drug candidate can be an effective antidote both kinetically and structurally to reactivate the tabun-inhibited enzyme.

Lagrangian Descriptors of Thermalized Transition States on Time-Varying Energy Surfaces

Thermalized chemical reactions driven under dynamical load are characteristic of activated dynamics for arbitrary nonautonomous systems. Recent generalizations of transition state theory to obtain formally exact rates have required the construction of a time-dependent transition state trajectory. Here, we show that Lagrangian descriptors can be used to obtain this structure directly. By developing a phase space separatrix that is void of recrossings, these constructs allow for the principal criterion in the implementation of modern rate theories to be satisfied. Thus, the reactive flux over a time-varying barrier can be determined without ambiguity in chemical reactions. The generality of the formalism suggests that this approach is applicable to any activated system subjected to arbitrary driving and thermal fluctuations.

Nanomechanical force transducers for biomolecular and intracellular measurements: is there room to shrink and why do it?

Over the past couple of decades there has been a tremendous amount of progress on the development of ultrasensitive nanomechanical instruments, which has enabled scientists to peer for the first time into the mechanical world of biomolecular systems. Currently, work-horse instruments such as the atomic force microscope and optical/magnetic tweezers have provided the resolution necessary to extract quantitative force data from various molecular systems down to the femtonewton range, but it remains difficult to access the intracellular environment with these analytical tools as they have fairly large sizes and complicated feedback systems. This review is focused on highlighting some of the major milestones and discoveries in the field of biomolecular mechanics that have been made possible by the development of advanced atomic force microscope and tweezer techniques as well as on introducing emerging state-of-the-art nanomechanical force transducers that are addressing the size limitations presented by these standard tools. We will first briefly cover the basic setup and operation of these instruments, and then focus heavily on summarizing advances in in vitro force studies at both the molecular and cellular level. The last part of this review will include strategies for shrinking down the size of force transducers and provide insight into why this may be important for gaining a more complete understanding of cellular activity and function.

Quantum chemical and steered molecular dynamics studies for one pot solution to reactivate aged acetylcholinesterase with alkylator oxime

Dimethyl(pyridin-2-yl)sulfonium based oxime has been designed to reverse the aging process of organophosphorus inhibited AChE and to reactivate the aged-AChE adduct. We have employed DFT M05-2X/6-31G(∗) level of theory in aqueous phase with polarizable continuum solvation model (PCM) for the methylation of phosphonate ester monoanion of the soman-aged adduct. The calculated free energy of activation for the methyl transfer process from designed dimethyl(phenyl)sulfonium (1) to aged AChE-OP adduct occurs with a barrier of 31.4kcal/mol at M05-2X/6-31G(∗) level of theory, which is 6.4kcal/mol lower compared to the aging process signifies the preferential reversal process to recover the aged AChE-OP adduct. The pyridine ring containing alkylated sulfonium species, dimethyl(pyridin-2-yl)sulfonium (2), reduced the free energy of activation by 4.4kcal/mol compared to the previously reported alkylating agent N-methyl-2-methoxypyridinium species (A) for the alkylation of aged AChE-OP adduct. The free enzyme can be liberated from the inhibited acetylcholinesterase with the sulfonium compound decorated with an oxime group to avoid the administration of oxime drugs separately. The calculated potential energy surfaces show that the oxime based sulfonium compound (3) can effectively methylate the aged phosphonate ester monoanion of soman aged-adduct. The calculated global reactivity descriptors of the oxime 3 also shed light on this observation. To gain better understanding for protein drug interaction as well as the unbinding and conformational changes of the drug candidate in the active site of cholinesterase, steered molecular dynamics (SMD) simulation with 3 has been performed. Through a protein-drug interaction parameters (rupture force profiles, hydrogen bonds, hydrophobic interactions), geometrical and the orientation of drug-like candidate, the oxime 3 suggests to orchestrate the better reactivation process. The docking studies have been performed with 3 in the aged AChE and BChE to obtain the initial geometry of the SMD studies. The docking methods adopted in this study have been verified with the available crystal geometry of 1-methyl-2-(pentafluorobenzyloxyimino)pyridinium compound in aged soman inhibited human BChE (PDB code: 4B0P). The computational study suggests that the newly designed oxime is a potential candidate to reactivate the aged-AChE adduct.

Investigating receptor-ligand systems of the cellulosome with AFM-based single-molecule force spectroscopy

Cellulosomes are discrete multienzyme complexes used by a subset of anaerobic bacteria and fungi to digest lignocellulosic substrates. Assembly of the enzymes onto the noncatalytic scaffold protein is directed by interactions among a family of related receptor-ligand pairs comprising interacting cohesin and dockerin modules. The extremely strong binding between cohesin and dockerin modules results in dissociation constants in the low picomolar to nanomolar range, which may hamper accurate off-rate measurements with conventional bulk methods. Single-molecule force spectroscopy (SMFS) with the atomic force microscope measures the response of individual biomolecules to force, and in contrast to other single-molecule manipulation methods (i.e. optical tweezers), is optimal for studying high-affinity receptor-ligand interactions because of its ability to probe the high-force regime (>120 pN). Here we present our complete protocol for studying cellulosomal protein assemblies at the single-molecule level. Using a protein topology derived from the native cellulosome, we worked with enzyme-dockerin and carbohydrate binding module-cohesin (CBM-cohesin) fusion proteins, each with an accessible free thiol group at an engineered cysteine residue. We present our site-specific surface immobilization protocol, along with our measurement and data analysis procedure for obtaining detailed binding parameters for the high-affinity complex. We demonstrate how to quantify single subdomain unfolding forces, complex rupture forces, kinetic off-rates, and potential widths of the binding well. The successful application of these methods in characterizing the cohesin-dockerin interaction responsible for assembly of multidomain cellulolytic complexes is further described.

PROBING αIIbβ3: LIGAND INTERACTIONS BY DYNAMIC FORCE SPECTROSCOPY AND SURFACE PLASMON RESONANCE

The interaction between platelet integrin αIIbβ3 and fibrin(ogen) plays a key role in blood clot formation and stability. Integrin antagonists, a class of pharmaceuticals used to prevent and treat cardiovascular disease, are designed to competitively interfere with this process. However, the energetics of the integrin-drug binding are not fully understood, potentially hampering further development of this class of pharmaceuticals. We integrated dynamic force spectroscopy (DFS) and surface plasmon resonance (SPR) to probe the energetics of complex formation between αIIbβ3 and cHarGD, a cyclic peptide integrin antagonist. Analysis of αIIbβ3:cHarGD DFS rupture force data at pulling rates of 14 000 pN/s, 42 000 pN/s and 70 000 pN/s yielded koff = 0.02-0.09 s-1, a dissociation energy barrier [Formula: see text] = 22-29 kJ/mol, and a potential well width x-1 = 0.5-0.8 nm. SPR kinetic data yielded an association rate constant kon = 7 × 103 L/mol-s and a dissociation rate constant koff = 10-2 s-1, followed by a slower stabilization step (τ ~ 400 s). Both DFS and SPR detected minimal interactions between αIIbβ3 and cHarGA demonstrating a key role for electrostatic interactions between the ligand aspartate and the integrin metal ion-dependent adhesion site (MIDAS). Our work provides new insights into the energy landscape of αIIbβ3's interactions with pharmacological and physiological ligands.

Atomic force microscopy-based molecular studies on the recognition of immunogenic chlorinated ovalbumin by macrophage receptors

This report presents simple and reliable approach developed to study the specific recognition events between chlorinated ovalbumin (OVA) and macrophages using atomic force microscopy (AFM). Thanks to the elimination of nonspecific adhesion, the interactions of the native and chlorinated OVA with a membrane of macrophages could be quantified using exclusively the so-called adhesion frequency (AF). The proposed system not only enabled the application of AFM-based force measurements for such poorly defined ligand-receptor pairs but also significantly improved both the acquisition and the processing of the data. The proteins were immobilized on the gold-coated AFM tips from the aqueous solutions containing charged thiol adsorbates. Such surface dilution of the proteins ensured the presence of single or just a few macromolecules at the tip-surface contact. The formation of negatively charged monolayer on the tip dramatically limited its nonspecific interactions with the macrophage surface. In such systems, AF was used as a measure of the recognition events even if the interaction forces varied significantly for sets of measurements. The system with the native OVA, a weak immunogen, showed only negligible AF compared with 85% measured for the immunogenic chlorinated OVA. The AF values varied with the tip-macrophage contact time and loading velocity. Blocking of the receptors by the chlorinated OVA was also confirmed. The developed approach can be also used to study other ligand-receptor interactions in poorly defined biological systems with intrinsically broad distribution of the rupture forces, thus opening new fields for AFM-based recognition on molecular level.

Delineation of the unbinding pathway of α-conotoxin ImI from the α7 nicotinic acetylcholine receptor

α-Conotoxins potently and specifically inhibit isoforms of nicotinic acetylcholine receptors (nAChRs) and are used as molecular probes and as drugs or drug leads. Interactions occurring during binding and unbinding events are linked to binding kinetics, and knowledge of these interactions could help in the development of α-conotoxins as drugs. Here, the unbinding process for the prototypical α-conotoxin ImI/α7-nAChR system was investigated theoretically, and three exit routes were identified using random accelerated molecular dynamics simulations. The route involving the smallest conformation perturbation was further divided into three subpathways, which were studied using steered molecular dynamics simulations. Of the three subpathways, two had better experimental support and lower potential of mean force, indicating that they might be sampled more frequently. Additionally, these subpathways were supported by previous experimental studies. Several pairwise interactions, including a cation-π interaction and charge and hydrogen bond interactions, were identified as potentially playing important roles in the unbinding event.

Strength of hydrogen bond network takes crucial roles in the dissociation process of inhibitors from the HIV-1 protease binding pocket

To understand the underlying mechanisms of significant differences in dissociation rate constant among different inhibitors for HIV-1 protease, we performed steered molecular dynamics (SMD) simulations to analyze the entire dissociation processes of inhibitors from the binding pocket of protease at atomistic details. We found that the strength of hydrogen bond network between inhibitor and the protease takes crucial roles in the dissociation process. We showed that the hydrogen bond network in the cyclic urea inhibitors AHA001/XK263 is less stable than that of the approved inhibitor ABT538 because of their large differences in the structures of the networks. In the cyclic urea inhibitor bound complex, the hydrogen bonds often distribute at the flap tips and the active site. In contrast, there are additional accessorial hydrogen bonds formed at the lateral sides of the flaps and the active site in the ABT538 bound complex, which take crucial roles in stabilizing the hydrogen bond network. In addition, the water molecule W301 also plays important roles in stabilizing the hydrogen bond network through its flexible movement by acting as a collision buffer and helping the rebinding of hydrogen bonds at the flap tips. Because of its high stability, the hydrogen bond network of ABT538 complex can work together with the hydrophobic clusters to resist the dissociation, resulting in much lower dissociation rate constant than those of cyclic urea inhibitor complexes. This study may provide useful guidelines for design of novel potent inhibitors with optimized interactions.

Thiol-based, site-specific and covalent immobilization of biomolecules for single-molecule experiments

The success of single-molecule (SM) experiments critically depends on the functional immobilization of the biomolecule(s) to be studied. With the continuing trend of combining SM fluorescence with SM force experiments, methods are required that are suitable for both types of measurements. We describe a general protocol for the site-specific and covalent coupling of any type of biomolecule that can be prepared with a free thiol group. The protocol uses a poly(ethylene glycol) (PEG) spacer, which carries an N-hydroxy succinimide (NHS) group on one end and a maleimide group on the other. After reacting the NHS group with an amino-functionalized surface, the relatively stable but highly reactive maleimide group allows the coupling of the biomolecule. This protocol provides surfaces with low fluorescence background, low nonspecific binding and a large number of reactive sites. Surfaces containing immobilized biomolecules can be obtained within 6 h.

Single-molecule force spectroscopy: a method for quantitative analysis of ligand–receptor interactions

The quantitative analysis of molecular interactions is of high interest in medical research. Most methods for the investigation of ligand–receptor complexes deal with huge ensembles of biomolecules, but often neglect interactions with low affinity or small subpopulations with different binding properties. Single-molecule force spectroscopy offers fascinating possibilities for the quantitative analysis of ligand–receptor interactions in a wide affinity range and the sensitivity to detect point mutations. Furthermore, this technique allows one to address questions about the related binding energy landscape. In this article, we introduce single-molecule force spectroscopy with a focus on novel developments in both data analysis and theoretical models for the technique. We also demonstrate two examples of the capabilities of this method.

Structural insights into substrate recognition and product expulsion in CTX-M enzymes

beta-Lactamase-mediated resistance to beta-lactam antibiotics poses a major threat to our antibiotic armamentarium. Among beta-lactamases, a significant threat comes from enzymes that hydrolyze extended-spectrum cephalosporins such as cefotaxime. Among the enzymes that exhibit this phenotype, the CTX-M family is found worldwide. These enzymes have a small active site, which makes it difficult to explain how they hydrolyze the bulky extended-spectrum cephalosporins into the binding site. We investigated noncovalent substrate recognition and product release in CTX-M enzymes using steered molecular dynamics simulation and X-ray diffraction. An arginine residue located far from the binding site favors the capture and tracking of substrates during entrance into the catalytic pocket. We show that the accommodation of extended-spectrum cephalosporins by CTX-M enzymes induced subtle changes in the active site and established a high density of electrostatic interactions. Interestingly, the product of the catalytic reaction initiates its own release because of steric hindrances and electrostatic repulsions. This suggests that there exists a general mechanism for product release for all members of the beta-lactamase family and probably for most carboxypeptidases.

Grading the commercial optical biosensor literature-Class of 2008: 'The Mighty Binders'

Optical biosensor technology continues to be the method of choice for label-free, real-time interaction analysis. But when it comes to improving the quality of the biosensor literature, education should be fundamental. Of the 1413 articles published in 2008, less than 30% would pass the requirements for high-school chemistry. To teach by example, we spotlight 10 papers that illustrate how to implement the technology properly. Then we grade every paper published in 2008 on a scale from A to F and outline what features make a biosensor article fabulous, middling or abysmal. To help improve the quality of published data, we focus on a few experimental, analysis and presentation mistakes that are alarmingly common. With the literature as a guide, we want to ensure that no user is left behind.

Differential binding of bispyridinium oxime drugs with acetylcholinesterase

AIM

To performe a time-dependent topographical delineation of protein-drug interactions to gain molecular insight into the supremacy of Ortho-7 over HI-6 in reactivating tabun-conjugated mouse acetylcholinesterase (mAChE).

METHODS

We conducted all-atom steered molecular dynamics simulations of the two protein-drug complexes. Through a host of protein-drug interaction parameters (rupture force profiles, hydrogen bonds, water bridges, hydrophobic interactions), geometrical, and orientation ordering of the drugs, we monitored the enzyme's response during the release of the drugs from its active-site.

RESULTS

The results show the preferential binding of the drugs with the enzyme. The pyridinium ring of HI-6 shows excellent complementary binding with the peripheral anionic site, whereas one of two identical pyridinium rings of Ortho-7 has excellent binding compatibility in the enzyme active-site where it can orchestrate the reactivation process. We found that the active pyridinium ring of HI-6 undergoes a complete turn along the active site axis, directed away from the active-site region during the course of the simulation.

CONCLUSION

Due to excellent cooperative binding of Ortho-7, as rendered by several cation-pi interactions with the active-site gorge of the enzyme, Ortho-7 may be a more efficient reactivator than HI-6. Our work supports the growing body of evidence that the efficacy of the drugs is due to the differential bindings of the oximes with AChE and can aid to the rational design of oxime drugs.

The effect of different force applications on the protein-protein complex Barnase-Barstar

Steered molecular dynamics simulations are a tool to examine the energy landscape of protein-protein complexes by applying external forces. Here, we analyze the influence of the velocity and geometry of the probing forces on a protein complex using this tool. With steered molecular dynamics, we probe the stability of the protein-protein complex Barnase-Barstar. The individual proteins are mechanically labile. The Barnase-Barstar binding site is more stable than the folds of the individual proteins. By using different force protocols, we observe a variety of responses of the system to the applied tension.

Controlled manipulation of bacteriophages using single-virus force spectroscopy

A method is described for the site-directed manipulation of single filamentous bacteriophages, by using phage display technology and atomic force microscopy. f1 filamentous bacteriophages were genetically engineered to display His-tags on their pIX tail. Following adsorption on nitrilotriacetate-terminated surfaces, force spectroscopy with tips bearing monoclonal anti-pIII antibodies was used to pull on individual phages via their pIII head. Analysis of the force-extension profiles revealed that upon pulling, the phages are progressively straightened into an extended orientation until rupture of the anti-pIII/pIII complex. The single-virus manipulation technique presented here provides new opportunities for understanding the forces driving cell-virus and material-virus interactions, and for characterizing the binding properties of polypeptide sequences or proteins selected by the phage display technology.

Triggering enzymatic activity with force

Integrating single molecule force spectroscopy with fluorescence-based techniques allows the manipulation of an enzyme with a periodic stretching and relaxation protocol while simultaneously monitoring its catalytic activity. After releasing the stretching force we observe a higher probability for enzymatic activity at a time of 1.7 s. A detailed theoretical analysis reveals that the relaxation from the force-induced enzyme conformation to the observed active conformation follows a cascade reaction with several steps and a free energy difference of at least 8 k(B)T. Our study clearly points out the direct influence of force on enzymatic activity and opens up a new way to study and manipulate (bio)catalytic reactions at the single molecule level.