Electrostatic enhancement of diffusion-controlled protein-protein association: comparison of theory and experiment on barnase and barstar

Computational electrostatic engineering of nanobodies for enhanced SARS-CoV-2 receptor binding domain recognition

This study presents a novel computational approach for engineering nanobodies (Nbs) for improved interaction with receptor-binding domain (RBD) of the SARS-CoV-2 spike protein. Using Protein Structure Reliability reports, RBD (7VYR_R) was selected and refined for subsequent Nb-RBD interactions. By leveraging electrostatic complementarity (EC) analysis, we engineered and characterized five Electrostatically Complementary Nbs (ECSb1-ECSb5) based on the CeVICA library's SR6c3 Nb. Through targeted modifications in the complementarity-determining regions (CDR) and framework regions (FR), we optimized electrostatic interactions to improve binding affinity and specificity. The engineered Nbs (ECSb3, ECSb4, and ECSb5) demonstrated high binding specificity for AS3, CA1, and CA2 epitopes. Interestingly, ECSb1 and ECSb2 selectively engaged with AS3 and CA1 instead of AS1 and AS2, respectively, due to a preference for residues that conferred superior binding complementarities. Furthermore, ECSbs significantly outperformed SR6c3 Nb in MM/GBSA results, notably, ECSb4 and ECSb3 exhibited superior binding free energies of -182.58 kcal.mol-1 and -119.07 kcal.mol-1, respectively, compared to SR6c3 (-105.50 kcal.mol-1). ECSbs exhibited significantly higher thermostability (100.4-148.3 kcal·mol⁻1) compared to SR6c3 (62.6 kcal·mol⁻1). Similarly, enhanced electrostatic complementarity was also observed for ECSb4-RBD and ECSb3-RBD (0.305 and 0.390, respectively) relative to SR6c3-RBD (0.233). Surface analyses confirmed optimized electrostatic patches and reduced aggregation propensity in the engineered Nb. This integrated EC and structural engineering approach successfully developed engineered Nbs with enhanced binding specificity, increased thermostability, and reduced aggregation, laying the groundwork for novel therapeutic applications targeting the SARS-CoV-2 spike protein.

Optimal Dielectric Boundary for Binding Free Energy Estimates in the Implicit Solvent

Accuracy of binding free energy calculations utilizing implicit solvent models is critically affected by parameters of the underlying dielectric boundary, specifically, the atomic and water probe radii. Here, a multidimensional optimization pipeline is used to find optimal atomic radii, specifically for binding calculations in the implicit solvent. To reduce overfitting, the optimization target includes separate, weighted contributions from both binding and hydration free energies. The resulting five-parameter radii set, OPT_BIND5D, is evaluated against experiment for binding free energies of 20 host-guest (H-G) systems, unrelated to the types of structures used in the training. The resulting accuracy for this H-G test set (root mean square error of 2.03 kcal/mol, mean signed error of -0.13 kcal/mol, mean absolute error of 1.68 kcal/mol, and Pearson's correlation of r = 0.79 with the experimental values) is on par with what can be expected from the fixed charge explicit solvent models. Best agreement with the experiment is achieved when the implicit salt concentration is set equal or close to the experimental conditions.

Exploring the significance of potassium homeostasis in copper ion binding to human αB-Crystallin

αB-Crystallin (αB-Cry) is a small heat shock protein known for its protective role, with an adaptable structure that responds to environmental changes through oligomeric dynamics. Cu(II) ions are crucial for cellular processes but excessive amounts are linked to diseases like cataracts and neurodegeneration. This study investigated how optimal and detrimental Cu(II) concentrations affect αB-Cry oligomers and their chaperone activity, within the potassium-regulated ionic-strength environment. Techniques including isothermal titration calorimetry, differential scanning calorimetry, fluorescence spectroscopy, inductively coupled plasma atomic emission spectroscopy, cyclic voltammetry, dynamic light scattering, circular dichroism, and MTT assay were employed and complemented by computational methods. Results showed that potassium ions affected αB-Cry's structure, promoting Cu(II) binding at multiple sites and scavenging ability, and inhibiting ion redox reactions. Low concentrations of Cu(II), through modifications of oligomeric interfaces, induce regulation of surface charge and hydrophobicity, resulting in an increase in chaperone activity. Subunit dynamics were regulated, maintaining stable interfaces, thereby inhibiting further aggregation and allowing the functional reversion to oligomers after stress. High Cu(II) disrupted charge/hydrophobicity balance, sewing sizable oligomers together through subunit-subunit interactions, suppressing oligomer dissociation, and reducing chaperone efficiency. This study offers insights into how Cu(II) and potassium ions influence αB-Cry, advancing our understanding of Cu(II)-related diseases.

Influence of molecular rebinding on the reaction rate of complex formation

We simulated Brownian diffusion and reaction-diffusion processes to study the influence of molecular rebinding on the reaction rates of bimolecular reactions. We found that the number of rebinding events, N_reb, is proportional to the target's size and inversely proportional to the diffusion coefficient D and simulation time-step Δt. We found the proportionality constant close to π^-1/2. We confirmed that N_reb is defined as a ratio of the activation-limited rate constant k_a to the diffusion-limited rate constant, k_D. We provide the formula describing the reactivity coefficient κ, modelling the transient-native complex transition for the activation-controlled reaction rates. We show that κ is proportional to (D/Δt)^1/2. Finally, we apply our rebinding-including reaction rate model to the real reactions of photoacid dissociation and protein association. Based on literature data for both types of reactions, we found the Δt time-scale. We show that for the photodissociation of a proton, the Δt is equal to 171 ± 18 fs and the average number of rebinding events is approximately equal to 40. For proteins, Δt is of the order of 100 ps with around 20 rebinding events. In both cases the timescale is similar to the timescale of fluctuation of the solvent molecules surrounding the reactants; vibrations and bending in the case of photoacid dissociation and diffusional motion for proteins.

Shark Antibody Variable Domains Rigidify Upon Affinity Maturation-Understanding the Potential of Shark Immunoglobulins as Therapeutics

Sharks and other cartilaginous fish are the phylogenetically oldest living organisms that have antibodies as part of their adaptive immune system. As part of their humoral adaptive immune response, they produce an immunoglobulin, the so-called immunoglobulin new antigen receptor (IgNAR), a heavy-chain only antibody. The variable domain of an IgNAR, also known as V NAR , binds the antigen as an independent soluble domain. In this study, we structurally and dynamically characterized the affinity maturation mechanism of the germline and somatically matured (PBLA8) V NAR to better understand their function and their applicability as therapeutics. We observed a substantial rigidification upon affinity maturation, which is accompanied by a higher number of contacts, thereby contributing to the decrease in flexibility. Considering the static x-ray structures, the observed rigidification is not obvious, as especially the mutated residues undergo conformational changes during the simulation, resulting in an even stronger network of stabilizing interactions. Additionally, the simulations of the V NAR in complex with the hen egg-white lysozyme show that the V NAR antibodies evidently follow the concept of conformational selection, as the binding-competent state already preexisted even without the presence of the antigen. To have a more detailed description of antibody-antigen recognition, we also present here the binding/unbinding mechanism between the hen egg-white lysozyme and both the germline and matured V NAR s. Upon maturation, we observed a substantial increase in the resulting dissociation-free energy barrier. Furthermore, we were able to kinetically and thermodynamically describe the binding process and did not only identify a two-step binding mechanism, but we also found a strong population shift upon affinity maturation toward the native binding pose.

Discrete-state stochastic kinetic models for target DNA search by proteins: Theory and experimental applications

To perform their functions, transcription factors and DNA-repair/modifying enzymes randomly search DNA in order to locate their specific targets on DNA. Discrete-state stochastic kinetic models have been developed to explain how the efficiency of the search process is influenced by the molecular properties of proteins and DNA as well as by other factors such as molecular crowding. These theoretical models not only offer explanations on the relation of microscopic processes to macroscopic behavior of proteins, but also facilitate the analysis and interpretation of experimental data. In this review article, we provide an overview on discrete-state stochastic kinetic models and explain how these models can be applied to experimental investigations using stopped-flow, single-molecule, nuclear magnetic resonance (NMR), and other biophysical and biochemical methods.

PAINT using proteins: A new brush for super-resolution artists

PAINT (points accumulation for imaging in nanoscale topography) refers to methods that achieve the sparse temporal labeling required for super-resolution imaging by using transient interactions between a biomolecule of interest and a fluorophore. There have been a variety of different implementations of this method since it was first described in 2006. Recent papers illustrate how transient peptide-protein interactions, rather than small molecule binding or DNA oligonucleotide duplex formation, can be employed to perform PAINT-based single molecule localization microscopy (SMLM). We discuss the different approaches to PAINT using peptide and protein interactions, and their applications in vitro and in vivo. We highlight the important parameters to consider when selecting suitable peptide-protein interaction pairs for such studies. We also note the opportunities for protein scientists to apply their expertise in guiding the choice of peptide and protein pairs that are used. Finally, we discuss the potential for expanding super-resolution imaging methods based on transient peptide-protein interactions, including the development of simultaneous multicolor imaging of multiple proteins and the study of very high and very low abundance proteins in live cells.

Influence of electrostatic forces on the association kinetics and conformational ensemble of an intrinsically disordered protein

Recent work has revealed that the association of a disordered region of a protein with a folded binding partner can occur as rapidly as association between two folded proteins. This is the case for the phosphatase calcineurin (CaN) and its association with its activator calmodulin. Calmodulin binds to the intrinsically disordered regulatory domain of CaN. Previous studies have shown that electrostatic steering can accelerate the binding of folded proteins with disordered ligands. Given that electrostatic forces are strong determinants of disordered protein ensembles, the relationship between electrostatics, conformational ensembles, and quaternary interactions is unclear. Here, we employ experimental approaches to explore the impact of electrostatic interactions on the association of calmodulin with the disordered regulatory region of CaN. We find that estimated association rate constants of calmodulin with our chosen calmodulin-substrates are within the diffusion-limited regime. The association rates are dependent on the ionic strength, indicating that favorable electrostatic forces increase the rate of association. Further, we show that charged amino acids outside the calmodulin-binding site modulate the binding rate. Conformational ensembles obtained from computer simulations suggest that electrostatic interactions within the regulatory domain might bias the conformational ensemble such that the calmodulin binding region is readily accessible. Given the prevalence of charged residues in disordered protein chains, our findings are likely relevant to many protein-protein interactions.

Co-localization and confinement of ecto-nucleotidases modulate extracellular adenosine nucleotide distributions

Nucleotides comprise small molecules that perform critical signaling roles in biological systems. Adenosine-based nucleotides, including adenosine tri-, di-, and mono-phosphate, are controlled through their rapid degradation by diphosphohydrolases and ecto-nucleotidases (NDAs). The interplay between nucleotide signaling and degradation is especially important in synapses formed between cells, which create signaling 'nanodomains'. Within these 'nanodomains', charged nucleotides interact with densely-packed membranes and biomolecules. While the contributions of electrostatic and steric interactions within such nanodomains are known to shape diffusion-limited reaction rates, less is understood about how these factors control the kinetics of nucleotidase activity. To quantify these factors, we utilized reaction-diffusion numerical simulations of 1) adenosine triphosphate (ATP) hydrolysis into adenosine monophosphate (AMP) and 2) AMP into adenosine (Ado) via two representative nucleotidases, CD39 and CD73. We evaluate these sequentially-coupled reactions in nanodomain geometries representative of extracellular synapses, within which we localize the nucleotidases. With this model, we find that 1) nucleotidase confinement reduces reaction rates relative to an open (bulk) system, 2) the rates of AMP and ADO formation are accelerated by restricting the diffusion of substrates away from the enzymes, and 3) nucleotidase co-localization and the presence of complementary (positive) charges to ATP enhance reaction rates, though the impact of these contributions on nucleotide pools depends on the degree to which the membrane competes for substrates. As a result, these contributions integratively control the relative concentrations and distributions of ATP and its metabolites within the junctional space. Altogether, our studies suggest that CD39 and CD73 nucleotidase activity within junctional spaces can exploit their confinement and favorable electrostatic interactions to finely control nucleotide signaling.

Molecular Dynamics Ensemble Refinement of Intrinsically Disordered Peptides According to Deconvoluted Spectra from Circular Dichroism

We have developed a computational method of atomistically refining the structural ensemble of intrinsically disordered peptides (IDPs) facilitated by experimental measurements using circular dichroism spectroscopy (CD). A major challenge surrounding this approach stems from the deconvolution of experimental CD spectra into secondary structure features of the IDP ensemble. Currently available algorithms for CD deconvolution were designed to analyze the spectra of proteins with stable secondary structures. Herein, our work aims to minimize any bias from the peptide deconvolution analysis by implementing a non-negative linear least-squares fitting algorithm in conjunction with a CD reference data set that contains soluble and denatured proteins (SDP48). The non-negative linear least-squares method yields the best results for deconvolution of proteins with higher disordered content than currently available methods, according to a validation analysis of a set of protein spectra with Protein Data Bank entries. We subsequently used this analysis to deconvolute our experimental CD data to refine our computational model of the peptide secondary structure ensemble produced by all-atom molecular dynamics simulations with implicit solvent. We applied this approach to determine the ensemble structures of a set of short IDPs, that mimic the calmodulin binding domain of calcium/calmodulin-dependent protein kinase II and its 1-amino-acid and 3-amino-acid mutants. Our study offers a, to our knowledge, novel way to solve the ensemble secondary structures of IDPs in solution, which is important to advance the understanding of their roles in regulating signaling pathways through the formation of complexes with multiple partners.

Transition path dynamics in the binding of intrinsically disordered proteins: A simulation study

Association of proteins and other biopolymers is a ubiquitous process in living systems. Recent single-molecule measurements probe the dynamics of association in unprecedented detail by measuring the properties of association transition paths, i.e., short segments of molecular trajectories between the time the proteins are close enough to interact and the formation of the final complex. Interpretation of such measurements requires adequate models for describing the dynamics of experimental observables. In an effort to develop such models, here we report a simulation study of the association dynamics of two oppositely charged, disordered polymers. We mimic experimental measurements by monitoring intermonomer distances, which we treat as "experimental reaction coordinates." While the dynamics of the distance between the centers of mass of the molecules is found to be memoryless and diffusive, the dynamics of the experimental reaction coordinates displays significant memory and can be described by a generalized Langevin equation with a memory kernel. We compute the most commonly measured property of transition paths, the distribution of the transition path time, and show that, despite the non-Markovianity of the underlying dynamics, it is well approximated as one-dimensional diffusion in the potential of mean force provided that an apparent value of the diffusion coefficient is used. This apparent value is intermediate between the slow (low frequency) and fast (high frequency) limits of the memory kernel. We have further studied how the mean transition path time depends on the ionic strength and found only weak dependence despite strong electrostatic attraction between the polymers.

Population balance modeling of assembly formation; accounting for the H-NS protein oligomerization and DNA binding mechanisms

The results obtained in this study demonstrate H-NS oligomerization in solution is not sufficient in forming large structures of H-NS-DNA assemblies indicating oligomerization along bacterial DNA is crucial. We present a comprehensive and novel multi-scale modeling scheme to study protein-protein and protein-DNA interactions, with the focus on population balance modeling of assemblies formed by the H-NS protein and bacterial DNA. The model helps understanding how the concentration of H-NS influences the formation of H-NS-DNA assemblies. The presented model covers the two underlying molecular mechanisms involved in assembly formation, H-NS oligomerization in solution and H-NS-DNA binding. Within protein oligomerization interactions, protein molecules are entitled to dimerize, propagate, recombine and deform (break). In addition, DNA binding and unbinding interactions are included to account for formation of filaments along DNA. All mechanisms have been modeled with their associated forward (formation) and backward (deformation) interactions. The results obtained agree well with former experimental studies.

Comparative analysis of plastocyanin-cytochrome f complex formation in higher plants, green algae and cyanobacteria

Mechanisms of the complex formation between plastocyanin and cytochrome f in higher plants (Spinacia oleracea and Brassica rapa), green microalgae Chlamydomonas reinhardtii and two species of cyanobacteria (Phormidium laminosum and Nostoc sp.) were investigated using combined Brownian and molecular dynamics simulations and hierarchical cluster analysis. In higher plants and green algae, electrostatic interactions force plastocyanin molecule close to the heme of cytochrome f. In the subsequent rotation of plastocyanin molecule around the point of electrostatic contact in the vicinity of cytochrome f, copper (Cu) atom approaches cytochrome heme forming a stable configuration where cytochrome f molecule behaves as a rather rigid body without conformational changes. In Nostoc plastocyanin molecule approaches cytochrome f in a different orientation (head-on) where the stabilization of the plastocyanin-cytochrome f complex is accompanied by the conformational changes of the G188E189D190 loop that stabilizes the whole complex. In cyanobacterium P. laminosum, electrostatic preorientation of the approaching molecules was not detected, thus indicating that random motions rather than long-range electrostatic interactions are responsible for the proper mutual orientation. We demonstrated that despite the structural similarity of the investigated electron transport proteins in different photosynthetic organisms, the complexity of molecular mechanisms of the complex formation increases in the following sequence: non-heterocystous cyanobacteria - heterocystous cyanobacteria - green algae - flowering plants.

Electrostatic Steering Enables Flow-Activated Von Willebrand Factor to Bind Platelet Glycoprotein, Revealed by Single-Molecule Stretching and Imaging

Von Willebrand factor (VWF), a large multimeric blood protein, senses changes in shear stress during bleeding and responds by binding platelets to plug ruptures in the vessel wall. Molecular mechanisms underlying this dynamic process are difficult to uncover using standard approaches due to the challenge of applying mechanical forces while monitoring structure and activity. By combining single-molecule fluorescence imaging with high-pressure, rapidly switching microfluidics, we reveal the key role of electrostatic steering in accelerating the binding between flow-activated VWF and GPIbα, and in rapidly immobilizing platelets under flow. We measure the elongation and tension-dependent activation of individual VWF multimers under a range of ionic strengths and pH levels, and find that the association rate is enhanced by 4 orders of magnitude by electrostatic steering. Under supraphysiologic salt concentrations, strong electrostatic screening dramatically decreases platelet binding to VWF in flow, revealing the critical role of electrostatic attraction in VWF-platelet binding during bleeding.

Electrostatic control of calcineurin's intrinsically-disordered regulatory domain binding to calmodulin

Calcineurin (CaN) is a serine/threonine phosphatase that regulates a variety of physiological and pathophysiological processes in mammalian tissue. The calcineurin (CaN) regulatory domain (RD) is responsible for regulating the enzyme's phosphatase activity, and is believed to be highly-disordered when inhibiting CaN, but undergoes a disorder-to-order transition upon diffusion-limited binding with the regulatory protein calmodulin (CaM). The prevalence of polar and charged amino acids in the regulatory domain (RD) suggests electrostatic interactions are involved in mediating calmodulin (CaM) binding, yet the lack of atomistic-resolution data for the bound complex has stymied efforts to probe how the RD sequence controls its conformational ensemble and long-range attractions contribute to target protein binding. In the present study, we investigated via computational modeling the extent to which electrostatics and structural disorder facilitate CaM/CaN association kinetics. Specifically, we examined several RD constructs that contain the CaM binding region (CAMBR) to characterize the roles of electrostatics versus conformational diversity in controlling diffusion-limited association rates, via microsecond-scale molecular dynamics (MD) and Brownian dynamic (BD) simulations. Our results indicate that the RD amino acid composition and sequence length influence both the dynamic availability of conformations amenable to CaM binding, as well as long-range electrostatic interactions to steer association. These findings provide intriguing insight into the interplay between conformational diversity and electrostatically-driven protein-protein association involving CaN, which are likely to extend to wide-ranging diffusion-limited processes regulated by intrinsically-disordered proteins.

Diffusion-limited association of disordered protein by non-native electrostatic interactions

Intrinsically disordered proteins (IDPs) usually fold during binding to target proteins. In contrast to interactions between folded proteins, this additional folding step makes the binding process more complex. Understanding the mechanism of coupled binding and folding of IDPs requires analysis of binding pathways that involve formation of the transient complex (TC). However, experimental characterization of TC is challenging because it only appears for a very brief period during binding. Here, we use single-molecule fluorescence spectroscopy to investigate the mechanism of diffusion-limited association of an IDP. A large enhancement of the association rate is observed due to the stabilization of TC by non-native electrostatic interactions. Moreover, photon-by-photon analysis reveals that the lifetime of TC for IDP binding is at least two orders of magnitude longer than that for binding of two folded proteins. This result suggests the long lifetime of TC is generally required for folding of IDPs during binding processes.

Intrinsically disordered proteins (IDPs) usually fold during binding to target proteins which involves the formation of a transient complex (TC). Here authors use single-molecule FRET to show that the lifetime of TC for IDP binding is very long due to the stabilization by non-native electrostatic interactions, which makes fast association possible.

Stopped-Flow Kinetic Techniques for Studying Binding Reactions of Intrinsically Disordered Proteins

Intrinsically disordered proteins are abundant in signaling processes such as transcription. Suitable binding and unbinding rates of proteins with their partners are critical for allowing them to perform their biological roles. Understanding how these are achieved, and indeed designing strategies for intervening or modulating related biological processes, therefore requires kinetic studies. In this chapter, we describe stopped-flow-based methods for determining association and dissociation rate constants for pairs of macromolecular binding partners. We describe how to select the simplest appropriate model to describe the interaction, and highlight cases where it is possible to distinguish between induced fit and conformational selection binding mechanisms. Finally, we go on to describe methods for examining the role of electrostatic forces in binding processes, and for describing the transition state for binding processes that have folding associated with them.

Ribosome surface properties may impose limits on the nature of the cytoplasmic proteome

Much of the molecular motion in the cytoplasm is diffusive, which possibly limits the tempo of processes. We studied the dependence of protein mobility on protein surface properties and ionic strength. We used surface-modified fluorescent proteins (FPs) and determined their translational diffusion coefficients (D) in the cytoplasm of Escherichia coli, Lactococcus lactis and Haloferax volcanii. We find that in E. coli D depends on the net charge and its distribution over the protein, with positive proteins diffusing up to 100-fold slower than negative ones. This effect is weaker in L. lactis and Hfx. volcanii due to electrostatic screening. The decrease in mobility is probably caused by interaction of positive FPs with ribosomes as shown in in vivo diffusion measurements and confirmed in vitro with purified ribosomes. Ribosome surface properties may thus limit the composition of the cytoplasmic proteome. This finding lays bare a paradox in the functioning of prokaryotic (endo)symbionts.

Kinetic mechanism of controlled Fab-arm exchange for the formation of bispecific immunoglobulin G1 antibodies

Bispecific antibodies (bsAbs) combine the antigen specificities of two distinct Abs and demonstrate therapeutic promise based on novel mechanisms of action. Among the many platforms for creating bsAbs, controlled Fab-arm exchange (cFAE) has proven useful based on minimal changes to native Ab structure and the simplicity with which bsAbs can be formed from two parental Abs. Despite a published protocol for cFAE and its widespread use in the pharmaceutical industry, the reaction mechanism has not been determined. Knowledge of the mechanism could lead to improved yields of bsAb at faster rates as well as foster adoption of process control. In this work, a combination of Förster resonance energy transfer (FRET), nonreducing SDS-PAGE, and strategic mutation of the Ab hinge region was employed to identify and characterize the individual steps of cFAE. Fluorescence correlation spectroscopy (FCS) was used to determine the affinity of parental (homodimer) and bispecific (heterodimer) interactions within the C_H3 domain, further clarifying the thermodynamic basis for bsAb formation. The result is a clear sequence of events with rate constants that vary with experimental conditions, where dissociation of the K409R parental Ab into half-Ab controls the rate of the reaction.

Affinity of IDPs to their targets is modulated by ion-specific changes in kinetics and residual structure

Intrinsically disordered proteins (IDPs) are characterized by a lack of defined structure. Instead, they populate ensembles of rapidly interconverting conformations with marginal structural stabilities. Changes in solution conditions such as temperature and crowding agents consequently affect IDPs more than their folded counterparts. Here we reveal that the residual structure content of IDPs is modulated both by ionic strength and by the type of ions present in solution. We show that these ion-specific structural changes result in binding affinity shifts of up to sixfold, which happen through alteration of both association and dissociation rates. These effects follow the Hofmeister series, but unlike the well-established effects on the stability of folded proteins, they already occur at low, hypotonic concentrations of salt. We attribute this sensitivity to the marginal stability of IDPs, which could have physiological implications given the role of IDPs in signaling, the asymmetric ion profiles of different cellular compartments, and the role of ions in biology.

Electrostatic forces govern the binding mechanism of intrinsically disordered histone chaperones

A unified picture to understand the protein recognition and function must include the native binding complex structure ensembles and the underlying binding mechanisms involved in specific biological processes. However, quantifications of both binding complex structures and dynamical mechanisms are still challenging for IDP. In this study, we have investigated the underlying molecular mechanism of the chaperone Chz1 and histone H2A.Z-H2B association by equilibrium and kinetic stopped-flow fluorescence spectroscopy. The dependence of free energy and kinetic rate constant on electrolyte mean activity coefficient and urea concentration are uncovered. Our results indicate a previous unseen binding kinetic intermediate. An initial conformation selection step of Chz1 is also revealed before the formation of this intermediate state. Based on these observations, a mixed mechanism of three steps including both conformation selection and induced fit is proposed. By combination of the ion- and denaturant-induced experiments, we demonstrate that electrostatic forces play a dominant role in the recognition of bipolar charged intrinsically disordered protein Chz1 to its preferred partner H2A.Z-H2B. Both the intra-chain and inter-chain electrostatic interactions have direct impacts on the native collapsed structure and binding mechanism.

Kinetic characterization of apoptotic Ras signaling through Nore1-MST1 complex formation

Ras-mediated apoptotic signaling is expected to be mediated via Rassf-MST complexes, but the system has been poorly characterized in vitro until now. Here we demonstrate that active H-Ras, Nore1A and MST1 form a stable ternary complex in vitro without other external factors, Nore1A interacting simultaneously with H-Ras and MST1 via its RBD and SARAH domain, respectively. Moreover, our data show for the first time that the SARAH domain of Nore1A plays a role in the Nore1A binding to H-Ras. Finally, we analyze the relation between the electrostatic and hydrophobic forces and kinetic constants of the Nore1A - H-Ras complex.

Role of non-native electrostatic interactions in the coupled folding and binding of PUMA with Mcl-1

PUMA, which belongs to the BH3-only protein family, is an intrinsically disordered protein (IDP). It binds to its cellular partner Mcl-1 through its BH3 motif, which folds upon binding into an α helix. We have applied a structure-based coarse-grained model, with an explicit Debye-Hückel charge model, to probe the importance of electrostatic interactions both in the early and the later stages of this model coupled folding and binding process. This model was carefully calibrated with the experimental data on helical content and affinity, and shown to be consistent with previously published experimental data on binding rate changes with respect to ionic strength. We find that intramolecular electrostatic interactions influence the unbound states of PUMA only marginally. Our results further suggest that intermolecular electrostatic interactions, and in particular non-native electrostatic interactions, are involved in formation of the initial encounter complex. We are able to reveal the binding mechanism in more detail than is possible using experimental data alone however, and in particular we uncover the role of non-native electrostatic interactions. We highlight the potential importance of such electrostatic interactions for describing the binding reactions of IDPs. Such approaches could be used to provide predictions for the results of mutational studies.

The conformational feasibility for the formation of reaching dimer in ASV and HIV integrase: a molecular dynamics study

Retroviral integrases are reported to form alternate dimer assemblies like the core-core dimer and reaching dimer. The core-core dimer is stabilized predominantly by an extensive interface between two catalytic core domains. The reaching dimer is stabilized by N-terminal domains that reach to form intermolecular interfaces with the other subunit's core and C-terminal domains (CTD), as well as CTD-CTD interactions. In this study, molecular dynamics (MD), Brownian dynamics (BD) simulations, and free energy analyses, were performed to elucidate determinants for the stability of the reaching dimer forms of full-length Avian Sarcoma Virus (ASV) and Human Immunodeficiency Virus (HIV) IN, and to examine the role of the C-tails (the last ~16-18 residues at the C-termini) in their structural dynamics. The dynamics of an HIV reaching dimer derived from small angle X-ray scattering and protein crosslinking data, was compared with the dynamics of a core-core dimer model derived from combining the crystal structures of two-domain fragments. The results showed that the core domains in the ASV reaching dimer express free dynamics, whereas those in the HIV reaching dimer are highly stable. BD simulations suggest a higher rate of association for the HIV core-core dimer than the reaching dimer. The predicted stability of these dimers was therefore ranked in the following order: ASV reaching dimer < HIV reaching dimer < composite core-core dimer. Analyses of MD trajectories have suggested residues that are critical for intermolecular contacts in each reaching dimer. Tests of these predictions and insights gained from these analyses could reveal a potential pathway for the association and dissociation of full-length IN multimers.

Forces Driving Chaperone Action

It is still unclear what molecular forces drive chaperone-mediated protein folding. Here, we obtain a detailed mechanistic understanding of the forces that dictate the four key steps of chaperone-client interaction: initial binding, complex stabilization, folding, and release. Contrary to the common belief that chaperones recognize unfolding intermediates by their hydrophobic nature, we discover that the model chaperone Spy uses long-range electrostatic interactions to rapidly bind to its unfolded client protein Im7. Short-range hydrophobic interactions follow, which serve to stabilize the complex. Hydrophobic collapse of the client protein then drives its folding. By burying hydrophobic residues in its core, the client's affinity to Spy decreases, which causes client release. By allowing the client to fold itself, Spy circumvents the need for client-specific folding instructions. This mechanism might help explain how chaperones can facilitate the folding of various unrelated proteins.

Botulinum Neurotoxin Serotype A Recognizes Its Protein Receptor SV2 by a Different Mechanism than Botulinum Neurotoxin B Synaptotagmin

Botulinum neurotoxins (BoNTs) exhibit extraordinary potency due to their exquisite neurospecificity, which is achieved by dual binding to complex polysialo-gangliosides and synaptic vesicle proteins. The luminal domain 4 (LD4) of the three synaptic vesicle glycoprotein 2 isoforms, SV2A‐C, identified as protein receptors for the most relevant serotype BoNT/A, binds within the 50 kDa cell binding domain H_C of BoNT/A. Here, we deciphered the BoNT/A‐SV2 interactions in more detail. In pull down assays, the binding of H_CA to SV2-LD4 isoforms decreases from SV2C >> SV2A > SV2B. A binding constant of 200 nM was determined for BoNT/A to rat SV2C-LD4 in GST pull down assay. A similar binding constant was determined by surface plasmon resonance for H_CA to rat SV2C and to human SV2C, the latter being slightly lower due to the substitution L563F in LD4. At pH 5, as measured in acidic synaptic vesicles, the binding constant of H_CA to hSV2C is increased more than 10-fold. Circular dichroism spectroscopy reveals that the quadrilateral helix of SV2C-LD4 already exists in solution prior to BoNT/A binding. Hence, the BoNT/A‐SV2C interaction is of different nature compared to BoNT/B‐Syt-II. In particular, the preexistence of the quadrilateral β-sheet helix of SV2 and its pH-dependent binding to BoNT/A via backbone–backbone interactions constitute major differences. Knowledge of the molecular details of BoNT/A‐SV2 interactions drives the development of high affinity peptides to counteract BoNT/A intoxications or to capture functional BoNT/A variants in innovative detection systems for botulism diagnostic.

Energetics of Hydrophilic Protein-Protein Association and the Role of Water

Hydrophilic protein-protein interfaces constitute a major part of all protein-protein interfaces and are thus of great importance. However, the quantitative characterization of their association is still an ongoing challenge and the driving force behind their association remains poorly characterized. Here, we have addressed the association of hydrophilic proteins and the role of water by means of extensive molecular dynamics simulations in explicit water using three well studied protein complexes; Barnase-Barstar, cytochrome c-cytochrome c peroxidase, and the N-terminal domain of enzyme I-histidine-containing phosphocarrier. The one-dimensional free energy profiles obtained from umbrella sampling simulations are downhill or, in other words, barrierless. Using these one-dimensional free energy profiles, the computed standard free energies of binding are -12.7 ± 1.1 kcal/mol, -9.4 ± 0.7 kcal/mol, and -8.4 ± 1.9 kcal/mol that are in reasonable to very good agreement with the experimental values of -19.6 kcal/mol, -8.8 kcal/mol, and -7.8 kcal/mol. As expected, analysis of the confined water between the hydrophilic complex partners shows that the density and the orientational order parameter deviate noticeably from the bulk values, especially at close separations of the confining proteins.

Free Energy Diagram for the Heterogeneous Enzymatic Hydrolysis of Glycosidic Bonds in Cellulose

Kinetic and thermodynamic data have been analyzed according to transition state theory and a simplified reaction scheme for the enzymatic hydrolysis of insoluble cellulose. For the cellobiohydrolase Cel7A from Hypocrea jecorina (Trichoderma reesei), we were able to measure or collect relevant values for all stable and activated complexes defined by the reaction scheme and hence propose a free energy diagram for the full heterogeneous process. For other Cel7A enzymes, including variants with and without carbohydrate binding module (CBM), we obtained activation parameters for the association and dissociation of the enzyme-substrate complex. The results showed that the kinetics of enzyme-substrate association (i.e. formation of the Michaelis complex) was almost entirely entropy-controlled and that the activation entropy corresponded approximately to the loss of translational and rotational degrees of freedom of the dissolved enzyme. This implied that the transition state occurred early in the path where the enzyme has lost these degrees of freedom but not yet established extensive contact interactions in the binding tunnel. For dissociation, a similar analysis suggested that the transition state was late in the path where most enzyme-substrate contacts were broken. Activation enthalpies revealed that the rate of dissociation was far more temperature-sensitive than the rates of both association and the inner catalytic cycle. Comparisons of one- and two-domain variants showed that the CBM had no influence on the transition state for association but increased the free energy barrier for dissociation. Hence, the CBM appeared to promote the stability of the complex by delaying dissociation rather than accelerating association.

Structural and Biochemical Basis for the Inhibitory Effect of Liprin-α3 on Mouse Diaphanous 1 (mDia1) Function

Diaphanous-related formins are eukaryotic actin nucleation factors regulated by an autoinhibitory interaction between the N-terminal RhoGTPase-binding domain (mDiaN) and the C-terminal Diaphanous-autoregulatory domain (DAD). Although the activation of formins by Rho proteins is well characterized, its inactivation is only marginally understood. Recently, liprin-α3 was shown to interact with mDia1. Overexpression of liprin-α3 resulted in a reduction of the cellular actin filament content. The molecular mechanisms of how liprin-α3 exerts this effect and counteracts mDia1 activation by RhoA are unknown. Here, we functionally and structurally define a minimal liprin-α3 core region, sufficient to recapitulate the liprin-α3 determined mDia1-respective cellular functions. We show that liprin-α3 alters the interaction kinetics and thermodynamics of mDiaN with RhoA·GTP and DAD. RhoA displaces liprin-α3 allosterically, whereas DAD competes with liprin-α3 for a highly overlapping binding site on mDiaN. Liprin-α3 regulates actin polymerization by lowering the regulatory potency of RhoA and DAD on mDiaN. We present a model of a mechanistically unexplored and new aspect of mDiaN regulation by liprin-α3.

Electrostatic interactions in the binding pathway of a transient protein complex studied by NMR and isothermal titration calorimetry

Much of our knowledge of protein binding pathways is derived from extremely stable complexes that interact very tightly, with lifetimes of hours to days. Much less is known about weaker interactions and transient complexes because these are challenging to characterize experimentally. Nevertheless, these types of interactions are ubiquitous in living systems. The combination of NMR relaxation dispersion Carr-Purcell-Meiboom-Gill (CPMG) experiments and isothermal titration calorimetry allows the quantification of rapid binding kinetics for complexes with submillisecond lifetimes that are difficult to study using conventional techniques. We have used this approach to investigate the binding pathway of the Src homology 3 (SH3) domain from the Fyn tyrosine kinase, which forms complexes with peptide targets whose lifetimes are on the order of about a millisecond. Long range electrostatic interactions have been shown to play a critical role in the binding pathways of tightly binding complexes. The role of electrostatics in the binding pathways of transient complexes is less well understood. Similarly to previously studied tight complexes, we find that SH3 domain association rates are enhanced by long range electrostatics, whereas short range interactions are formed late in the docking process. However, the extent of electrostatic association rate enhancement is several orders of magnitudes less, whereas the electrostatic-free basal association rate is significantly greater. Thus, the SH3 domain is far less reliant on electrostatic enhancement to achieve rapid association kinetics than are previously studied systems. This suggests that there may be overall differences in the role played by electrostatics in the binding pathways of extremely stable versus transient complexes.

Allostery within a transcription coactivator is predominantly mediated through dissociation rate constants

The kinase-inducible domain interacting (KIX) domain of CREB binding protein binds to multiple intrinsically disordered transcription factors in vivo at two distinct sites on its surface. Several reports have been made of allosteric communication between these two sites in this well-characterized model system. In this work, we have performed fluorescence stopped-flow measurements to investigate the kinetics of binding of five KIX binding proteins. We find that they all have similar association and dissociation rate constants for complex formation, despite their wide range of intrinsic helical propensities. Furthermore, by careful arrangement of pseudofirst-order conditions, we have been able to show that both association and dissociation rate constants are decreased when a partner is bound at the alternative site. These decreases suggest that positive allosteric effects are not mediated by structural changes in binding sites but rather, through a more general mechanism, largely mediated through dissociation, which we propose is largely related to changes in the flexibility of the KIX domain itself.

Positive and negative impacts of nonspecific sites during target location by a sequence-specific DNA-binding protein: origin of the optimal search at physiological ionic strength

The inducible transcription factor Egr-1, which recognizes a 9-bp target DNA sequence via three zinc-finger domains, rapidly activates particular genes upon cellular stimuli such as neuronal signals and vascular stresses. Here, using the stopped-flow fluorescence method, we measured the target search kinetics of the Egr-1 zinc-finger protein at various ionic strengths between 40 and 400 mM KCl and found the most efficient search at 150 mM KCl. We further investigated the kinetics of intersegment transfer, dissociation, and sliding of this protein on DNA at distinct concentrations of KCl. Our data suggest that Egr-1's kinetic properties are well suited for efficient scanning of chromosomal DNA in vivo. Based on a newly developed theory, we analyzed the origin of the optimal search efficiency at physiological ionic strength. Target association is accelerated by nonspecific binding to nearby sites and subsequent sliding to the target as well as by intersegment transfer. Although these effects are stronger at lower ionic strengths, such conditions also favor trapping of the protein at distant nonspecific sites, decelerating the target association. Our data demonstrate that Egr-1 achieves the optimal search at physiological ionic strength through a compromise between the positive and negative impacts of nonspecific interactions with DNA.

Inhibition of antithrombin and bovine serum albumin native state aggregation by heparin

Protein native state aggregation, a major problem in pharmaceutical and biological processes, has been addressed pharmacologically by the addition of protein-binding excipients. Heparin (Hp), a highly sulfated polysaccharide, interacts with numerous proteins with moderate to high affinity, but reports about its effect on protein aggregation are contradictory. We studied the pH dependence of the aggregation of antithrombin (AT) and bovine serum albumin (BSA) in the presence and absence of heparin. High-precision turbidimetry showed strong aggregation for both AT and BSA in I = 10 mM NaCl, conditions at which electrostatically driven Hp binding and aggregation both occur, with more obvious aggregation of heparin-free AT appearing as larger aggregate size. Aggregation of AT was dramatically inhibited at Hp: protein 6:1 (mole ratio); however, the effect at 0.5:1 Hp:protein was greater for BSA. Frontal analysis capillary electrophoresis showed a much larger equilibrium association constant Kobs between Hp and AT, in accord with the onset of Hp binding at a higher pH; both effects are explained by the higher charge density of the positive domain for AT as revealed by modeling with DelPhi. The corresponding modeling images showed that these domains persist at high salt only for AT, consistent with the 160-fold drop in Kobs at 100 mM salt for BSA-Hp binding. The smaller inhibition effect for AT arises from the tendency of its uncomplexed monomer to form larger aggregates more rapidly, but the stronger binding of Hp to AT does not facilitate Hp-induced aggregate dissolution which occurs more readily for BSA. This can be attributed to the higher density of AT aggregates evidenced by higher fractal dimensions. Differences between inhibition and reversal by Hp arise because the former may depend on the stage at which Hp enters the aggregation process and the latter on aggregate size and morphology.

Fixed charges in the gel matrix of sensor chips and dissociation in diffusion gradients influence the detection of fast protein-protein interactions

In molecular interaction studies based on surface plasmon resonance (SPR) measurements, the ligand is often immobilized in a thin carboxydextran gel matrix. Here we investigated the influence of the charged gel on the results of such SPR measurements. At physiological ionic strength, analytes with a net charge of more than about 5 are considerably enriched or depleted due to the Donnan potential under commonly applied experimental conditions. Below physiological ionic strength, enrichment was found to be even stronger than predicted by Donnan theory. The influence of the gel matrix on the apparent binding is prevented in competition experiments, in which SPR measurements are only used to discriminate between free and complexed analyte while the interaction between analyte and ligand is studied in solution. However, if the analyte-ligand interaction is very fast, thermodynamic equilibrium is disturbed near the interface where free analyte binds to the immobilized ligand due to mass transport limitation. Consequently, the soluble analyte-ligand complex dissociates, which results in an overestimation of free analyte. In experiments of calmodulin binding to fragments of the KCNH1 ion channel protein this mass-transport-induced dissociation led to a systematic underestimation of the affinity. We conclude that the insufficient discrimination between the true analyte-ligand binding and the complex interactions of the analyte with the gel phase may result in systematic errors. The theoretical framework for recognizing and avoiding such errors is provided.

Electrostatically accelerated encounter and folding for facile recognition of intrinsically disordered proteins

Achieving facile specific recognition is essential for intrinsically disordered proteins (IDPs) that are involved in cellular signaling and regulation. Consideration of the physical time scales of protein folding and diffusion-limited protein-protein encounter has suggested that the frequent requirement of protein folding for specific IDP recognition could lead to kinetic bottlenecks. How IDPs overcome such potential kinetic bottlenecks to viably function in signaling and regulation in general is poorly understood. Our recent computational and experimental study of cell-cycle regulator p27 (Ganguly et al., J. Mol. Biol. (2012)) demonstrated that long-range electrostatic forces exerted on enriched charges of IDPs could accelerate protein-protein encounter via "electrostatic steering" and at the same time promote "folding-competent" encounter topologies to enhance the efficiency of IDP folding upon encounter. Here, we further investigated the coupled binding and folding mechanisms and the roles of electrostatic forces in the formation of three IDP complexes with more complex folded topologies. The surface electrostatic potentials of these complexes lack prominent features like those observed for the p27/Cdk2/cyclin A complex to directly suggest the ability of electrostatic forces to facilitate folding upon encounter. Nonetheless, similar electrostatically accelerated encounter and folding mechanisms were consistently predicted for all three complexes using topology-based coarse-grained simulations. Together with our previous analysis of charge distributions in known IDP complexes, our results support a prevalent role of electrostatic interactions in promoting efficient coupled binding and folding for facile specific recognition. These results also suggest that there is likely a co-evolution of IDP folded topology, charge characteristics, and coupled binding and folding mechanisms, driven at least partially by the need to achieve fast association kinetics for cellular signaling and regulation.

Heparin decamer bridges a growth factor and an oligolysine by different charge-driven interactions

Full-length heparin is widely used in tissue engineering applications due its multiple protein-binding sites that allow it to retain growth factor affinity while associating with oligopeptide components of the tissue scaffold. However, the extent to which oligopeptide coupling interferes with cognate protein binding is difficult to predict. To investigate such simultaneous interactions, we examined a well-defined ternary system comprised of acidic fibroblast growth factor (FGF), tetralysine (K4), with a heparin decamer (dp10) acting as a noncovalent coupler. Electrospray ionization mass spectrometry was used to assess binding affinities and complex stoichiometries as a function of ionic strength for dp10·K4 and FGF·dp10. The ionic strength dependence of K4·dp10 formation is qualitatively consistent with binding driven by the release of condensed counterions previously suggested for native heparin with divalent oligopeptides (Mascotti, D. P.; Lohman, T. M. Biochemistry 1995, 34, 2908-2915). On the other hand, FGF binding displays more complex ionic strength dependence, with higher salt resistance. Remarkably, dp10 that can bind two FGF molecules can only bind one tetralysine. The limited binding of K4 to dp10 suggests that the tetralysine might not block growth factor binding, and the 1:1:1 ternary complex is indeed observed. The analysis of mass distribution of the bound dp10 chains in FGF·dp10, FGF2·dp10, and FGF·dp10·K4 complexes indicated that higher degrees of dp10 sulfation promote the formation of FGF2·dp10 and FGF·dp10·K4. Thus, the selectivity of appropriately chosen short heparin chains could be used to modulate growth factor sequestration and release in a way not feasible with heterogeneous native heparin. In support of this, human hepatocellular carcinoma cells (HEP3Bs) treated with FGF·dp10·K4 were found to exhibit biological activity similar to cells treated with FGF.

Remarkably fast coupled folding and binding of the intrinsically disordered transactivation domain of cMyb to CBP KIX

Association rates for interactions between folded proteins have been investigated extensively, allowing the development of computational and theoretical prediction methods. Less is known about association rates for complexes where one or more partner is initially disordered, despite much speculation about how they may compare to those for folded proteins. We have attached a fluorophore to the N-terminus of the 25 amino acid cMyb peptide used previously in NMR and equilibrium studies (termed FITC-cMyb), and used this to monitor the kinetics of its interaction with the KIX protein. We have investigated the ionic strength and temperature dependence of the kinetics, and conclude that the association process is extremely fast, apparently exceeding the rates predicted by formulations applicable to interactions between pairs of folded proteins. This is despite the fact that not all collisions result in complex formation (there is an observable activation energy for the association process). We propose that this is partially a result of the disordered nature of the FITC-cMyb peptide itself.

Hydrodynamic effects on the relative rotational velocity of associating proteins

Hydrodynamic steering effects on the barnase-barstar association were studied through the analysis of the relative rotational velocity of the proteins. We considered the two proteins approaching each other in response to their electrostatic attraction and employed a method that accounts for the long-range and many-body character of the hydrodynamic interactions, as well as the complicated shapes of the proteins. Hydrodynamic steering effects were clearly seen when attractive forces were applied to the geometric centers of the proteins (resulting in zero torques) and the attraction acted along the line that connects centers of geometry of proteins in their crystallographic complex. When we rotated barstar relative to barnase around this line by an angle in the range from -90° to 60°, the rotational velocity arising solely from hydrodynamic interactions restored the orientation of the proteins in the crystal structure. However, because, in reality, both electrostatic forces and torques act on the proteins and these forces and torques depend on the protein-protein distance and the relative orientation of the binding partners, we also investigated more realistic situations employing continuum electrostatics calculations based on atomistic protein models. Overall, we conclude that hydrodynamic interactions aid barnase and barstar in assuming a proper relative orientation upon complex formation.

Tyrosine replacement of PSGL-1 reduces association kinetics with P- and L-selectin on the cell membrane

Binding of selectins to P-selectin glycoprotein ligand-1 (PSGL-1) mediates tethering and rolling of leukocytes on the endothelium during inflammation. Previous measurements obtained with a flow-chamber assay have shown that mutations of three tyrosines at the PSGL-1 N-terminus (Y46, Y48, and Y51) increase the reverse rates and their sensitivity to the force of bonds with P- and L-selectin. However, the effects of these mutations on the binding affinities and forward rates have not been studied. We quantified these effects by using an adhesion frequency assay to measure two-dimensional affinity and kinetic rates at zero force. Wild-type PSGL-1 has 2.2- to 8.5-fold higher binding affinities for P- and L-selectin than PSGL-1 mutants with two of three tyrosines substituted by phenylalanines, and 9.6- to 49-fold higher affinities than the PSGL-1 mutant with all three tyrosines replaced. In descending order, the affinity decreased from wild-type to Y48/51F, Y46/51F, Y46/48F, and Y46/48/51F. The affinity differences were attributed to major changes in the forward rate and minor changes in the reverse rate, suggesting that these tyrosines regulate the accessibility of PSGL-1 to P- and L-selectin via electrostatic interactions, which is supported by molecular-dynamics simulations. Our results provide insights into the structure-function relationship of receptor-ligand binding at a single-residue level.