Understanding and Manipulating Electrostatic Fields at the Protein-Protein Interface Using Vibrational Spectroscopy and Continuum Electrostatics Calculations

Determining the Electrostatic Contributions of GTPase-GEF Complexes on Interfacial Drug Binding Specificity: A Case Study of a Protein-Drug-Protein Complex

Understanding the factors that contribute to specificity of protein-protein interactions allows for design of orthosteric small molecules. Within this environment, a small molecule requires both structural and electrostatic complementarity. While the structural contribution to protein-drug-protein specificity is well characterized, electrostatic contributions require more study. To this end, we used a series of protein complexes involving Arf1 bound to guanine nucleotide exchange factors (GEFs) that are sensitive or resistant to the small molecule brefeldin A (BFA). By comparing BFA-sensitive Arf1-Gea1p and Arf1-ARNO with different combinations of four BFA sensitizing ARNO mutations (ARNOwt, ARNO1M, ARNO3M, and ARNO4M), we describe how electrostatic environments at each interface guide BFA binding specificity. We labeled Arf1 with cyanocysteine at several interfacial sites and measured by nitrile adsorption frequencies to map changes in electric field at each interface using the linear Stark equation. Temperature dependence of nitrile vibrational spectra was used to investigate differences in hydrogen bonding environments. These comparisons showed that interfacial electric field at the surface of Arf1 varied substantially depending on the GEF. The greatest differences were seen between Arf1-ARNOwt and Arf1-ARNO4M, suggesting a greater change in electric field is required for BFA binding to Arf1-ARNO. Additionally, rigidity of the interface of the Arf1-ARNO complex correlated strongly with BFA sensitivity, indicating that flexible interfaces are sensitive to disruption upon orthosteric small molecule binding. These findings demonstrate a qualitatively consistent electrostatic environment for Arf1 binding and more subtle differences preventing BFA specificity. We discuss how these results will guide improved design of other small molecules that can target protein-protein interfaces.

Electrostatic Impact of Brefeldin A on Thiocyanate Probes Surrounding the Interface of Arf1-BFA-ARNO4M, a Protein-Drug-Protein Complex

Protein-protein interactions regulate many cellular processes, making them ideal drug candidates. Design of such drugs, however, is hindered by a lack of understanding of the factors that contribute to the interaction specificity. Specific protein-protein complexes possess both structural and electrostatic complementarity, and while structural complementarity of protein complexes has been extensively investigated, fundamental understanding of the complicated networks of electrostatic interactions at these interfaces is lacking, thus hindering the rational design of orthosterically binding small molecules. To better understand the electrostatic interactions at protein interfaces and how a small molecule could contribute to and fit within that environment, we used a model protein-drug-protein system, Arf1-BFA-ARNO4M, to investigate how small molecule brefeldin A (BFA) perturbs the Arf1-ARNO4M interface. By using nitrile probe labeled Arf1 sites and measuring vibrational Stark effects as well as temperature dependent infrared shifts, we measured changes in the electric field and hydrogen bonding at this interface upon BFA binding. At all five probe locations of Arf1, we found that the vibrational shifts resulting from BFA binding corroborate trends found in Poisson-Boltzmann calculations of surface potentials of Arf1-ARNO4M and Arf1-BFA-ARNO4M, where BFA contributes negative electrostatic potential to the protein interface. The data also corroborate previous hypotheses about the mechanism of interfacial binding and confirm that alternating patches of hydrophobic and polar interactions lead to BFA binding specificity. These findings demonstrate the impact of BFA on this protein-protein interface and have implications for the design of other interfacial drug candidates.

Quantification of the electric field inside protein active sites and fullerenes

While electrostatic interactions are exceedingly accountable for biological functions, no simple method exists to directly estimate or measure the electrostatic field in protein active sites. The electrostatic field inside the protein is generally inferred from the shift in the vibrational stretching frequencies of nitrile and thionitrile probes at the active sites through several painstaking and time-consuming experiments like vibrational Stark effect spectroscopy (VSS). Here we present a simple, fast, and reliable methodology, which can efficiently predict the vibrational Stark tuning rates (VSRs) of a large variety of probes within 10% error of the reported experimental data. Our methodology is based on geometry optimization and frequency calculations in the presence of an external electric field to predict the accurate VSR of newly designed nitrile/thionitrile probes. A priori information of VSRs is useful for difficult experiments such as catalytic/enzymatic study and in structural biology. We also applied our methodology successfully to estimate the electric field inside fullerenes and nano-onions, which is encouraging for researchers to adopt it for further applications in materials science and supramolecular chemistry.

Protein electrostatics: From computational and structural analysis to discovery of functional fingerprints and biotechnological design

Computationally driven engineering of proteins aims to allow them to withstand an extended range of conditions and to mediate modified or novel functions. Therefore, it is crucial to the biotechnological industry, to biomedicine and to afford new challenges in environmental sciences, such as biocatalysis for green chemistry and bioremediation. In order to achieve these goals, it is important to clarify molecular mechanisms underlying proteins stability and modulating their interactions. So far, much attention has been given to hydrophobic and polar packing interactions and stability of the protein core. In contrast, the role of electrostatics and, in particular, of surface interactions has received less attention. However, electrostatics plays a pivotal role along the whole life cycle of a protein, since early folding steps to maturation, and it is involved in the regulation of protein localization and interactions with other cellular or artificial molecules. Short- and long-range electrostatic interactions, together with other forces, provide essential guidance cues in molecular and macromolecular assembly. We report here on methods for computing protein electrostatics and for individual or comparative analysis able to sort proteins by electrostatic similarity. Then, we provide examples of electrostatic analysis and fingerprints in natural protein evolution and in biotechnological design, in fields as diverse as biocatalysis, antibody and nanobody engineering, drug design and delivery, molecular virology, nanotechnology and regenerative medicine.

A Critical Evaluation of Vibrational Stark Effect (VSE) Probes with the Local Vibrational Mode Theory

Over the past two decades, the vibrational Stark effect has become an important tool to measure and analyze the in situ electric field strength in various chemical environments with infrared spectroscopy. The underlying assumption of this effect is that the normal stretching mode of a target bond such as CO or CN of a reporter molecule (termed vibrational Stark effect probe) is localized and free from mass-coupling from other internal coordinates, so that its frequency shift directly reflects the influence of the vicinal electric field. However, the validity of this essential assumption has never been assessed. Given the fact that normal modes are generally delocalized because of mass-coupling, this analysis was overdue. Therefore, we carried out a comprehensive evaluation of 68 vibrational Stark effect probes and candidates to quantify the degree to which their target normal vibration of probe bond stretching is decoupled from local vibrations driven by other internal coordinates. The unique tool we used is the local mode analysis originally introduced by Konkoli and Cremer, in particular the decomposition of normal modes into local mode contributions. Based on our results, we recommend 31 polyatomic molecules with localized target bonds as ideal vibrational Stark effect probe candidates.

Vibrational Lifetime of the SCN Protein Label in H2O and D2O Reports Site-Specific Solvation and Structure Changes During PYP's Photocycle

The application of vibrational labels such as thiocyanate  (-S-C≡N) for studying protein structure and dynamics is thriving. Absorption spectroscopy is usually employed to obtain wavenumber and line shape of the label. An observable of great significance might be the vibrational lifetime, which can be obtained by pump probe or 2D-IR spectroscopy. Due to the insulating effect of the heavy sulfur atom in the case of the SCN label, the lifetime of the C≡N oscillator is expected to be particularly sensitive to its surrounding as it is not dominated by through-bond relaxation. We therefore investigate the vibrational lifetime of the SCN label at various positions in the blue light sensor protein Photoactive Yellow Protein (PYP) in the ground state and signaling state of the photoreceptor. We find that the vibrational lifetime of the C≡N stretching mode is strongly affected both by its protein environment and by the degree of exposure to the solvent. Even for label positions where the line shape and wavenumber observed by FTIR are barely changing upon activation of the photoreceptor, we find that the lifetime can change considerably. To obtain an unambiguous measure for the solvent exposure of the labeled site, we show that it is imperative to compare the lifetimes in H₂O and D₂O. Importantly, the lifetimes shorten in H₂O as compared to D₂O for water exposed labels, while they stay largely the same for buried labels. We quantify this effect by defining a solvent exclusion coefficient (SEC). The response of the label's vibrational lifetime to its solvent exposure renders it a suitable universal probe for protein investigations. This applies even to systems that are otherwise hard to address, such as transient or short-lived states, which could be created during a protein's working cycle (as here in PYP) or during protein folding. It is also applicable to flexible systems (intrinsically disordered proteins), protein-protein and protein-membrane interactions.

Following local light-induced structure changes and dynamics of the photoreceptor PYP with the thiocyanate IR label

Photoactive Yellow Protein (PYP) is a bacterial blue light receptor that enters a photocycle after excitation. The intermediate states are formed on time scales ranging from femtoseconds up to hundreds of milliseconds, after which the signaling state with a lifetime of about 1 s is reached. To investigate structural changes and dynamics, we incorporated the SCN IR label at distinct positions of the photoreceptor via cysteine mutation and cyanylation. FT-IR measurements of the SCN label at different sites of the well-established dark state structure of PYP characterized the spectral response of the label to differences in the environment. Under constant blue light irradiation, we observed the formation of the signaling state with significant changes of wavenumber and lineshape of the SCN bands. Thereby we deduced light-induced structural changes in the local environment of the labels. These results were supported by molecular dynamics simulations on PYP providing the solvent accessible surface area (SASA) at the different positions. To follow protein dynamics via the SCN label during the photocycle, we performed step-scan FT-IR measurements with a time resolution of 10 μs. Global analysis yielded similar time constants of τ1 = 70 μs, τ2 = 640 μs, and τ3 > 20 ms for the wild type and τ1 = 36 μs, τ2 = 530 μs, and τ3 > 20 ms for the SCN-labeled mutant PYP-A44C*, a mutant which provided a sufficiently large SCN difference signal to measure step-scan FT-IR spectra. In comparison to the protein (amide, E46) and chromophore bands the dynamics of the SCN label show a different behavior. This result indicates that the local kinetics sensed by the label are different from the global protein kinetics.

New Insights into Quinine-DNA Binding Using Raman Spectroscopy and Molecular Dynamics Simulations

Quinine's ability to bind DNA and potentially inhibit transcription and translation has been examined as a mode of action for its antimalarial activity. UV absorption and fluorescence-based studies have lacked the chemical specificity to develop an unambiguous molecular-level picture of the binding interaction. To address this, we use Raman spectroscopy and molecular dynamics (MD) to investigate quinine-DNA interactions. We demonstrate that quinine's strongest Raman band in the fingerprint region, which derives from a symmetric stretching mode of the quinoline ring, is highly sensitive to the local chemical environment and pH. The frequency shifts observed for this mode in solvents of varying polarity can be explained in terms of the Stark effect using a simple Onsager solvation model, indicating that the vibration reports on the local electrostatic environment. However, specific chemical interactions between the quinoline ring and its environment, such as hydrogen bonding and π-stacking, perturb the frequency of this mode in a more complicated but predictable manner. We use this vibration as a spectroscopic probe to investigate the binding interaction between quinine and DNA. We find that, when the quinoline ring is protonated, quinine weakly intercalates into DNA by forming π-stacking interactions with the base pairs. The Raman spectra indicate that quinine can intercalate into DNA with a ratio reaching up to roughly one molecule per 25 base pairs. Our results are confirmed by MD simulations, which also show that the quinoline ring adopts a t-shaped π-stacking geometry with the DNA base pairs, whereas the quinuclidine head group weakly interacts with the phosphate backbone in the minor groove. We expect that the spectral correlations determined here will enable future studies to probe quinine's antimalarial activities, such as disrupting hemozoin biocrystallization, which is hypothesized to be, among other things, one of its primary modes of action against Plasmodium parasites.

Long-Range Modulations of Electric Fields in Proteins

Electrostatic interactions are essential for controlling the protein structure and function. Whereas so far experimental and theoretical efforts focused on the effect of local electrostatics, this work aims at elucidating the long-range modulation of electric fields in proteins upon binding to charged surfaces. The study is based on cytochrome c (Cytc) variants carrying nitrile reporters for the vibrational Stark effect that are incorporated into the protein via genetic engineering and chemical modification. The Cytc variants were thoroughly characterized with respect to possible structural perturbations due to labeling. For the proteins in solution, the relative hydrogen bond occupancy and the calculated electric fields, both obtained from molecular dynamics (MD) simulations, and the experimental nitrile stretching frequencies were used to develop a relationship for separating hydrogen-bonding and non-hydrogen-bonding electric field effects. This relationship provides an excellent description for the stable Cytc variants in solution. For the proteins bound to Au electrodes coated with charged self-assembled monolayers (SAMs), the underlying MD simulations can only account for the electric field changes Δ E_ads due to the formation of the electrostatic SAM-Cytc complexes but not for the additional contribution, Δ E_int, representing the consequences of the potential drops over the electrode/SAM/protein interfaces. Both Δ E_ads and Δ E_int, determined at distances between 20 and 30 Å with respect to the SAM surface, are comparable in magnitude to the non-hydrogen-bonding electric field in the unbound protein. This long-range modulation of the internal electric field may be of functional relevance for proteins in complexes with partner proteins (Δ E_ads) and attached to membranes (Δ E_ads + Δ E_int).

Pandemic Avian Influenza and Intra/Interhaemagglutinin Subtype Electrostatic Variation among Viruses Isolated from Avian, Mammalian, and Human Hosts

Host jump can result in deadly pandemic events when avian influenza A viruses broaden their host specificity and become able to infect mammals, including humans. Haemagglutinin—the major capsid protein in influenza A viruses—is subjected to high rate mutations, of which several occur at its “head”: the receptor-binding domain that mediates specific binding to host cell receptors. Such surface-changing mutations may lead to antigenically novel influenza A viruses hence in pandemics by host jump and in vaccine escape by antigenic drift. Changes in haemagglutinin surface electrostatics have been recently associated with antigenic drift and with clades evolution and spreading in H5N1 and H9N2 viruses. We performed a comparative analysis of haemagglutinin surface electrostatics to investigate clustering and eventual fingerprints among representative pandemic (H5 and H7) and nonpandemic (H4 and H6) avian influenza viral subtypes. We observed preferential sorting of viruses isolated from mammalian/human hosts among these electrostatic clusters of a subtype; however, sorting was not “100% specific” to the different clusters. Therefore, electrostatic fingerprints can help in understanding, but they cannot explain alone the host jumping mechanism.

Exploring local solvation environments of a heme protein using the spectroscopic reporter 4-cyano-l-phenylalanine

The vibrational reporter unnatural amino acid (UAA) 4-cyano-l-phenylalanine (pCNF) was genetically incorporated individually at three sites (5, 36, and 78) in the heme protein Caldanaerobacter subterraneus H-NOX to probe local hydration environments. The UAA pCNF was incorporated site-specifically using an engineered, orthogonal tRNA synthetase in E. coli. The ability of all of the pCNF-containing H-NOX proteins to form the ferrous CO, NO, or O₂ ligated and unligated states was confirmed with UV-Vis spectroscopy. The solvation state at each site of the three sites of pCNF incorporation was assessed using temperature-dependent infrared spectroscopy. Specifically, the frequency-temperature line slope (FTLS) method was utilized to show that the nitrile group at site 36 was fully solvated and the nitrile group at site 78 was de-solvated (buried) in the heme pocket. The nitrile group at site 5 was found to be partially solvated suggesting that the nitrile group was involved in moderate strength hydrogen bonds. These results were confirmed by the determination of the X-ray crystal structure of the H-NOX protein construct containing pCNF at site 5.

Conformational Ensembles of Calmodulin Revealed by Nonperturbing Site-Specific Vibrational Probe Groups

Seven native residues on the regulatory protein calmodulin, including three key methionine residues, were replaced (one by one) by the vibrational probe amino acid cyanylated cysteine, which has a unique CN stretching vibration that reports on its local environment. Almost no perturbation was caused by this probe at any of the seven sites, as reported by CD spectra of calcium-bound and apo calmodulin and binding thermodynamics for the formation of a complex between calmodulin and a canonical target peptide from skeletal muscle myosin light chain kinase measured by isothermal titration. The surprising lack of perturbation suggests that this probe group could be applied directly in many protein–protein binding interfaces. The infrared absorption bands for the probe groups reported many dramatic changes in the probes’ local environments as CaM went from apo- to calcium-saturated to target peptide-bound conditions, including large frequency shifts and a variety of line shapes from narrow (interpreted as a rigid and invariant local environment) to symmetric to broad and asymmetric (likely from multiple coexisting and dynamically exchanging structures). The fast intrinsic time scale of infrared spectroscopy means that the line shapes report directly on site-specific details of calmodulin’s variable structural distribution. Though quantitative interpretation of the probe line shapes depends on a direct connection between simulated ensembles and experimental data that does not yet exist, formation of such a connection to data such as that reported here would provide a new way to evaluate conformational ensembles from data that directly contains the structural distribution. The calmodulin probe sites developed here will also be useful in evaluating the binding mode of calmodulin with many uncharacterized regulatory targets.

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.

Dangerous liaisons: anion-induced protonation in phosphate–polyamine interactions and their implications for the charge states of biologically relevant surfaces

The specific binding of phosphates on surface amino groups has complex consequences on the protonation equilibriums and effective surface charges with serious implications for further biochemically relevant interactions.

Identification of Ion-Pair Structures in Solution by Vibrational Stark Effects

Ion pairing is a fundamental consideration in many areas of chemistry and has implications in a wide range of sciences and technologies that include batteries and organic photovoltaics. Ions in solution are known to inhabit multiple possible states, including free ions (FI), contact ion pairs (CIP), and solvent-separated ion pairs (SSIP). However, in solutions of organic radicals and nonmetal electrolytes, it is often difficult to distinguish between these states. In the first part of this work, we report evidence for the formation of SSIPs in low-polarity solvents and distinct measurements of CIP, SSIP, and FI, by using the ν(C≡N) infrared (IR) band of a nitrile-substituted fluorene radical anion. Use of time-resolved IR detection following pulse radiolysis allowed us to unambiguously assign the peak of the FI. In the presence of nonmetal electrolytes, two distinct red-shifted peaks were observed and assigned to the CIP and SSIP. The assignments are interpreted in the framework of the vibrational Stark effect (VSE) and are supported by (1) the solvent dependence of ion-pair populations, (2) the observation of a cryptand-separated sodium ion pair that mimics the formation of SSIPs, and (3) electronic structure calculations. In the second part of this work, we show that a blue-shift of the ν(C≡N) IR band due to the VSE can be induced in a nitrile-substituted fluorene radical anion by covalently tethering it to a metal-chelating ligand that forms an intramolecular ion pair upon reduction and complexation with sodium ion. This adds support to the conclusion that the shift in IR absorptions by ion pairing originates from the VSE. These results combined show that we can identify ion-pair structures by using the VSE, including the existence of SSIPs in a low-polarity solvent.

Vibrational solvatochromism of nitrile infrared probes: beyond the vibrational Stark dipole approach

Systematic probing of local environments around biopolymers is important for understanding their functions.