Electrostatic steering enhances the rate of cAMP binding to phosphodiesterase: Brownian dynamics modeling

Removal of Metal from Carboxypeptidase A Proceeds via a Split Pathway: Implications for the General Mechanisms of Metalloenzyme Inactivation by Chelating Agents

The chelation of protein-bound metal ions is typically thought to follow either a dissociative (D) or an associative (A) path. While the former mechanism involves the spontaneous dissociation of the metal from the protein prior to chelation, the latter route is characterized by the formation of an intermediate protein-metal-chelator ternary complex. Using the prototypical zinc protease carboxypeptidase A (CPA) and a variety of charged and neutral chelating agents, we demonstrate that inactivation of the enzyme (and likely other metalloproteins) proceeds through a split pathway with contributions from both D- and A-type mechanisms. In the case of charged chelators such as ethylenediaminetetraacetic acid (EDTA), the proportions of both paths could be tuned over a wide range through variation of the chelator concentration and the ionic strength, I (from ∼95% A type at low I values to ∼5% at high I values). By measuring the EDTA concentration and time dependence of CPA inactivation and fitting the obtained kinetic data to a modified time-dependent inhibition model, we obtained the rate constants for the A and D paths (kinact and koff, respectively) and the inhibition constant (KI) for the formation of the CPA-Zn2+-EDTA ternary complex, indicating that the decreased contribution of the A-type mechanism at high ionic strengths originates from a large (40-fold; at I = 0.5 M) increase in KI. This observation might be related to a triarginine motif in CPA that electrostatically steers negatively charged substrates into the active site and may therefore also guide carboxylate-bearing chelators toward the Zn2+ ion.

On the possibility of the existence of orienting hydrodynamic steering effects in the kinetics of receptor-ligand association

In the vast majority of biologically relevant cases of receptor-ligand complex formation, the binding site of the receptor is a small part of its surface, and moreover, formation of a biologically active complex often requires a specific orientation of the ligand relative to the binding site. Before the formation of the initial form of the complex, only long-range, electrostatic and hydrodynamic interactions can act between the ligand approaching the binding site and the receptor. In this context, the question arises whether as a result of these interactions, there is a pre-orientation of the ligand towards the binding site, which to some extent would accelerate the formation of the complex. The role of electrostatic interactions in the orientation of the ligand relative to the binding site of the receptor is well documented. The analogous role of hydrodynamic interactions, although assessed as very significant by Brune and Kim (PNAS 91, 2930-2934, (1994)), is still debatable. In this article, I present the current state of knowledge on this subject and consider the possibilities of demonstrating the orienting effect of hydrodynamic interactions in the processes of receptor-ligand association, in an experimental way supported by computer simulations.

Searching for Hydrodynamic Orienting Effects in the Association of Tri-N-acetylglucosamine with Hen Egg-White Lysozyme

Using stopped-flow fluorometry, we determined rate constants for the formation of diffusional encounter complexes of tri-N-acetylglucosamine (NAG₃) with hen egg-white lysozyme (k_a^WT) and its double mutant Asp48Asn/Lys116Gln (k_a^MT). We defined binding anisotropy, κ ≡ (k_a^WT – k_a^MT)/(k_a^WT + k_a^MT), and determined its ionic strength dependence. Our goal was to check if this ionic strength dependence provides information about the orienting hydrodynamic effects in the ligand-binding process. We also computed ionic strength dependence of the binding anisotropy from Brownian dynamics simulations using simple models of the lysozyme–NAG₃ system. The results of our experiments indicate that in the case of lysozyme and NAG₃ such hydrodynamic orienting effects are rather negligible. On the other hand, the results of our Brownian dynamics simulations prove that there exist molecular systems for which such orienting effects are substantial. However, the ionic strength dependence of the rate constants for the wild-type and modified systems do not exhibit any qualitative features that would allow us to conclude the presence of hydrodynamic orienting effects from stopped-flow experiments alone. Nevertheless, the results of our simulations suggest the presence of hydrodynamic orienting effects in the receptor–ligand association when the anisotropy of binding depends on the solvent viscosity.

Mechanism of Dynamic Binding of Replication Protein A to ssDNA

Replication protein A (RPA) serves as a hub protein inside eukaryotic cells, where it coordinates crucial DNA metabolic processes and activates the DNA-damage response system. A characteristic feature of its action is to associate with single-stranded DNA (ssDNA) intermediates before handing them over to downstream proteins. The length of ssDNA intermediates differs for different pathways. This means that RPA must have mechanisms for selective processing of ssDNA intermediates based on their length, the knowledge of which is fundamental to elucidate when and how DNA repair and replication processes are symphonized. By employing extensive molecular dynamics simulations, we investigated the mechanism of binding of RPA to ssDNA of different lengths. We show that the binding involves dynamic equilibrium with a stable intermediate, the population of which increases with the length of ssDNA. The vital underlying factors are decoded through collective variable principal component analysis. It suggests a differently orchestrated set of interactions that define the action of RPA based on the length of ssDNA intermediates. We further estimated the association kinetics that matches excellently well with previous experimental studies and probed the diffusion mechanism of RPA to ssDNA. RPA diffuses on short ssDNA through progressive "bulge" formation. With long ssDNA, we observed a conformational change in ssDNA coupled with its binding to RPA in a cooperative fashion. This unanticipated binding mechanism successfully explains how the "short-lived", long ssDNA intermediates are processed quickly in vivo. This study thus reveals the molecular basis of several recent experimental observations related to RPA binding to ssDNA and provides novel insights into the RPA functioning in DNA repair and replication.

Protein-Protein Interactions of Phosphodiesterases

BACKGROUND

Phosphodiesterases (PDEs) are enzymes that play a key role in terminating cyclic nucleotides signalling by catalysing the hydrolysis of 3', 5'- cyclic adenosine monophosphate (cAMP) and/or 3', 5' cyclic guanosine monophosphate (cGMP), the second messengers within the cell that transport the signals produced by extracellular signalling molecules which are unable to get into the cells. However, PDEs are proteins which do not operate alone but in complexes that made up of a many proteins.

OBJECTIVE

This review highlights some of the general characteristics of PDEs and focuses mainly on the Protein-Protein Interactions (PPIs) of selected PDE enzymes. The objective is to review the role of PPIs in the specific mechanism for activation and thereby regulation of certain biological functions of PDEs.

METHODS

The article discusses some of the PPIs of selected PDEs as reported in recent scientific literature. These interactions are critical for understanding the biological role of the target PDE.

RESULTS

The PPIs have shown that each PDE has a specific mechanism for activation and thereby regulation a certain biological function.

CONCLUSION

Targeting of PDEs to specific regions of the cell is based on the interaction with other proteins where each PDE enzyme binds with specific protein(s) via PPIs.

Does Ionic Screening Lower Activation Barriers for Conformational Transitions in Proteins?

In this work, we investigated the kinetics of binding of hen egg-white lysozyme with tri- N-acetylglucosamine in aqueous solutions, at two values of pH, 3.2 and 11, as a function of ionic strength, by a stopped-flow method with tryptophyl fluorescence observation of the transients. We analyzed registered reaction progress curves by employing numerical integration of appropriate chemical master equations. We discriminated between several binding models and established that the process observed in experiments follows a two-step mechanism, composed of four elementary stages: diffusional formation of an encounter complex, dissociation of the encounter complex, conformational transition of the encounter complex to the final complex, and the reverse transformation, i.e., from the final complex to the encounter complex. We evaluated rate constants of these elementary stages and determined their dependencies on solution ionic strength. Regardless of solution pH, rate constants of both forward and reverse conformational transitions increase with an increasing ionic strength. This suggests that ionic screening of intramolecular electrostatic interactions may act to lower the activation barrier for conformational transition in proteins.

Brownian dynamic study of an enzyme metabolon in the TCA cycle: Substrate kinetics and channeling

Malate dehydrogenase (MDH) and citrate synthase (CS) are two pacemaking enzymes involved in the tricarboxylic acid (TCA) cycle. Oxaloacetate (OAA) molecules are the intermediate substrates that are transferred from the MDH to CS to carry out sequential catalysis. It is known that, to achieve a high flux of intermediate transport and reduce the probability of substrate leaking, a MDH-CS metabolon forms to enhance the OAA substrate channeling. In this study, we aim to understand the OAA channeling within possible MDH-CS metabolons that have different structural orientations in their complexes. Three MDH-CS metabolons from native bovine, wild-type porcine, and recombinant sources, published in recent work, were selected to calculate OAA transfer efficiency by Brownian dynamics (BD) simulations and to study, through electrostatic potential calculations, a possible role of charges that drive the substrate channeling. Our results show that an electrostatic channel is formed in the metabolons of native bovine and recombinant porcine enzymes, which guides the oppositely charged OAA molecules passing through the channel and enhances the transfer efficiency. However, the channeling probability in a suggested wild-type porcine metabolon conformation is reduced due to an extended diffusion length between the MDH and CS active sites, implying that the corresponding arrangements of MDH and CS result in the decrease of electrostatic steering between substrates and protein surface and then reduce the substrate transfer efficiency from one active site to another.

Hydrodynamic Steering in Protein Association Revisited: Surprisingly Minuscule Effects of Considerable Torques

We investigate the previously postulated hydrodynamic steering phenomenon, resulting from complication of molecular shapes, its magnitude and possible relevance for protein-ligand and protein-protein diffusional encounters, and the kinetics of diffusion-controlled association. We consider effects of hydrodynamic interactions in a prototypical model system consisting of a cleft enzyme and an elongated substrate, and real protein-protein complexes, that of barnase and barstar, and human growth hormone and its binding protein. The kinetics of diffusional encounters is evaluated on the basis of rigid-body Brownian dynamics simulations in which hydrodynamic interactions between molecules are modeled using the bead-shell method for detailed description of molecular surfaces, and the first-passage-time approach. We show that magnitudes of steering torques resulting from the hydrodynamic coupling of associating molecules, evaluated for the studied systems on the basis of the Stokes-Einstein type relations for arbitrarily shaped rigid bodies, are comparable with magnitudes of torques resulting from electrostatic interactions of binding partners. Surprisingly, however, unlike in the case of electrostatic torques that strongly affect the diffusional encounter, overall effects of hydrodynamic steering torques on the association kinetics, while clearly discernible in Brownian dynamics simulations, are rather minute. We explain this result as a consequence of the thermal agitation of the binding partners. Our finding is relevant for the general understanding of a wide spectrum of molecular processes in solution but there is also a more practical aspect to it if one considers the low level of shape detail of models that are usually employed to evaluate hydrodynamic interactions in particle-based Stokesian and Brownian dynamics simulations of multicomponent biomolecular systems. Results described in the current work justify, in part at least, such a low-resolution description.

Channeling of cAMP in PDE-PKA Complexes Promotes Signal Adaptation

Spatiotemporal control of the cAMP signaling pathway is governed by both hormonal stimulation of cAMP generation by adenylyl cyclases (activation phase) and cAMP hydrolysis by phosphodiesterases (PDEs) (termination phase). The termination phase is initiated by PDEs actively targeting the protein kinase A (PKA) R-subunit through formation of a PDE-PKAR-cyclic adenosine monophosphate (cAMP) complex (the termination complex). Our results using PDE8 as a model PDE, reveal that PDEs mediate active hydrolysis of cAMP bound to its receptor RIα by enhancing the enzymatic activity. This accelerated cAMP turnover occurs via formation of a stable PDE8-RIα complex, where the protein-protein interface forms peripheral contacts and the central ligand cements this ternary interaction. The basis for enhanced catalysis is active translocation of cAMP from its binding site on RIα to the hydrolysis site on PDE8 through direct "channeling." Our results reveal cAMP channeling in the PDE8-RIα complex and a molecular description of how this channel facilitates processive hydrolysis of unbound cAMP. Thus, unbound cAMP maintains the PDE8-RIα complex while being hydrolyzed, revealing an undiscovered mode for amplification of PKA activity by cAMP-mediated sequestration of the R-subunit by PDEs. This novel regulatory mode explains the paradox of cAMP signal amplification by accelerated PDE-mediated cAMP turnover. This highlights how target effector proteins of small-molecule ligands can promote enzyme-mediated ligand hydrolysis by scaffolding effects. Enhanced activity of the PDE8-RIα complex facilitates robust desensitization, allowing the cell to respond to dynamic levels of cAMP rather than steady-state levels. The PDE8-RIα complex represents a new class of PDE-based complexes for specific drug discovery targeting the cAMP signaling pathway.

Electric-Field-Assisted Anion-π Catalysis

This report focuses on the remote control of anion-π catalysis by electric fields. We have synthesized and immobilized anion-π catalysts to explore the addition reaction of malonic acid half thioesters to enolate acceptors on conductive indium tin oxide surfaces. Exposed to increasing electric fields, anion-π catalysts show an increase in activity and an inversion of selectivity. These changes originate from a more than 100-fold rate enhancement of the disfavored enolate addition reaction that coincides with an increase in selectivity of transition-state recognition by up to -14.8 kJ mol^-1. The addition of nitrate with strong π affinity nullified (IC₅₀ = 2.2 mM) the responsiveness of anion-π catalysts to electric fields. These results support that the polarization of the π-acidic naphthalenediimide surface in anion-π catalysts with electric fields increases the recognition of anionic intermediates and transition states on this polarized π surface, that is, the existence and relevance of electric-field-assisted anion-π catalysis.

Linking computation and experiments to study the role of charge-charge interactions in protein folding and stability

Over the past two decades there has been an increase in appreciation for the role of surface charge-charge interactions in protein folding and stability. The perception shifted from the belief that charge-charge interactions are not important for protein folding and stability to the near quantitative understanding of how these interactions shape the folding energy landscape. This led to the ability of computational approaches to rationally redesign surface charge-charge interactions to modulate thermodynamic properties of proteins. Here we summarize our progress in understanding the role of charge-charge interactions for protein stability using examples drawn from my own laboratory and touch upon unanswered questions.

Understanding ligand-receptor non-covalent binding kinetics using molecular modeling

Kinetic properties may serve as critical differentiators and predictors of drug efficacy and safety, in addition to the traditionally focused binding affinity. However the quantitative structure-kinetics relationship (QSKR) for modeling and ligand design is still poorly understood. This review provides an introduction to the kinetics of drug binding from a fundamental chemistry perspective. We focus on recent developments of computational tools and their applications to non-covalent binding kinetics.

Electrostatic Steering Accelerates C3d:CR2 Association

Electrostatic effects are ubiquitous in protein interactions and are found to be pervasive in the complement system as well. The interaction between complement fragment C3d and complement receptor 2 (CR2) has evolved to become a link between innate and adaptive immunity. Electrostatic interactions have been suggested to be the driving factor for the association of the C3d:CR2 complex. In this study, we investigate the effects of ionic strength and mutagenesis on the association of C3d:CR2 through Brownian dynamics simulations. We demonstrate that the formation of the C3d:CR2 complex is ionic strength-dependent, suggesting the presence of long-range electrostatic steering that accelerates the complex formation. Electrostatic steering occurs through the interaction of an acidic surface patch in C3d and the positively charged CR2 and is supported by the effects of mutations within the acidic patch of C3d that slow or diminish association. Our data are in agreement with previous experimental mutagenesis and binding studies and computational studies. Although the C3d acidic patch may be locally destabilizing because of unfavorable Coulombic interactions of like charges, it contributes to the acceleration of association. Therefore, acceleration of function through electrostatic steering takes precedence to stability. The site of interaction between C3d and CR2 has been the target for delivery of CR2-bound nanoparticle, antibody, and small molecule biomarkers, as well as potential therapeutics. A detailed knowledge of the physicochemical basis of C3d:CR2 association may be necessary to accelerate biomarker and drug discovery efforts.