Protein Electrostatic Potential Research Articles

Structural Dynamics of the Bacillus subtilis MntR Transcription Factor Is Locked by Mn2+ Binding.

Manganese (II) ions are essential for a variety of bacterial cellular processes. The transcription factor MntR is a metallosensor that regulates Mn2+ ion homeostasis in the bacterium Bacillus subtilis. Its DNA-binding affinity is increased by Mn2+ ion binding, allowing it to act as a transcriptional repressor of manganese import systems. Although experimentally well-researched, the molecular mechanism that regulates this process is still a puzzle. Computational simulations supported by circular dichroism (CD), differential scanning calorimetry (DSC) and native gel electrophoresis (native-PAGE) experiments were employed to study MntR structural and dynamical properties in the presence and absence of Mn2+ ions. The results of molecular dynamics (MD) simulations revealed that Mn2+ ion binding reduces the structural dynamics of the MntR protein and shifts the dynamic equilibrium towards the conformations adequate for DNA binding. Results of CD and DSC measurements support the computational results showing the change in helical content and stability of the MntR protein upon Mn2+ ion binding. Further, MD simulations show that Mn2+ binding induces polarization of the protein electrostatic potential, increasing the positive electrostatic potential of the DNA-binding helices in particular. In order to provide a deeper understanding of the changes in protein structure and dynamics due to Mn2+ binding, a mutant in which Mn2+ binding is mimicked by a cysteine bridge was constructed and also studied computationally and experimentally.

Spike protein of SARS-CoV-2 Omicron variant: An in-silico study evaluating spike interactions and immune evasion.

The fundamentals of the infectivity and immune evasion of the SARS-CoV-2 Omicron variant are not yet fully understood. Here, we carried out an in-silico study analyzing the spike protein, the protein electrostatic potential, and the potential immune evasion. The analysis was based on the structure of the spike protein from two SARS-CoV-2 variants, the original Wuhan and the Botswana (Omicron). The full-length genome sequences and protein sequences were obtained from databanks. The interaction of the spike proteins with the human Angiotensin Converting Enzyme 2 (ACE2) receptor was evaluated through the open-source software. The Immune Epitope Database was used to analyze the potential immune evasion of the viruses. Our data show that the Omicron spike protein resulted in 37 amino acid changes. The physicochemical properties of the spike had changed, and the electrostatic potentials differed between both variants. This resulted in a decrease in protein interactions, which does not establish a greater interaction with the ACE2 receptor. These changes compromise key receptor-binding motif residues in the SARS-CoV-2 spike protein that interact with neutralizing antibodies and ACE2. These mutations appear to confer enhanced properties of infectivity. The Omicron variant appears to be more effective at evading immune responses.

Effect of Protein–Protein Interactions on Translational Diffusion of Spheroidal Proteins

One of the commonly accepted approaches to estimate protein–protein interactions (PPI) in aqueous solutions is the analysis of their translational diffusion. The present review article observes a phenomenological approach to analyze PPI effects via concentration dependencies of self- and collective translational diffusion coefficient for several spheroidal proteins derived from the pulsed field gradient NMR (PFG NMR) and dynamic light scattering (DLS), respectively. These proteins are rigid globular α-chymotrypsin (ChTr) and human serum albumin (HSA), and partly disordered α-casein (α-CN) and β-lactoglobulin (β-Lg). The PPI analysis enabled us to reveal the dominance of intermolecular repulsion at low ionic strength of solution (0.003–0.01 M) for all studied proteins. The increase in the ionic strength to 0.1–1.0 M leads to the screening of protein charges, resulting in the decrease of the protein electrostatic potential. The increase of the van der Waals potential for ChTr and α-CN characterizes their propensity towards unstable weak attractive interactions. The decrease of van der Waals interactions for β-Lg is probably associated with the formation of stable oligomers by this protein. The PPI, estimated with the help of interaction potential and idealized spherical molecular geometry, are in good agreement with experimental data.

Assessment of the Components of the Electrostatic Potential of Proteins in Solution: Comparing Experiment and Theory.

In this work, the components of the protein electrostatic potentials in solution are analyzed with NMR paramagnetic relaxation enhancement experiments and compared with continuum solution theory, and multiscale simulations. To determine the contributions of the solution components, we analyze them at different ionic strengths from 0 to 745 mM. A theoretical approximation allows the determination of the electrostatic potential at a given proton without reference to the protein structure given the ratio of paramagnetic relaxation enhancements rates between a cationic and an anionic probe. The results derived from simulations show good agreement with experiment and simple continuum solvent theory for many of the residues. A discrepancy including a switch of sign of the electrostatic potential was observed for particular residues. By considering the components of the potential, we found the discrepancy is mainly caused by angular correlations of the probe molecules with these residues. The correction for the correlations allows a more accurate analysis of the experiments determining the electrostatic potential of proteins in solution.

How Can We Predict Accurate Electrochromic Shifts for Biochromophores? A Case Study on the Photosynthetic Reaction Center.

Protein-embedded chromophores are responsible for light harvesting, excitation energy transfer, and charge separation in photosynthesis. A critical part of the photosynthetic apparatus are reaction centers (RCs), which comprise groups of (bacterio)chlorophyll and (bacterio)pheophytin molecules that transform the excitation energy derived from light absorption into charge separation. The lowest excitation energies of individual pigments (site energies) are key for understanding photosynthetic systems, and form a prime target for quantum chemistry. A major theoretical challenge is to accurately describe the electrochromic (Stark) shifts in site energies produced by the inhomogeneous electric field of the protein matrix. Here, we present large-scale quantum mechanics/molecular mechanics calculations of electrochromic shifts for the RC chromophores of photosystem II (PSII) using various quantum chemical methods evaluated against the domain-based local pair natural orbital (DLPNO) implementation of the similarity-transformed equation of motion coupled cluster theory with single and double excitations (STEOM-CCSD). We show that certain range-separated density functionals (ωΒ97, ωΒ97X-V, ωΒ2PLYP, and LC-BLYP) correctly reproduce RC site energy shifts with time-dependent density functional theory (TD-DFT). The popular CAM-B3LYP functional underestimates the shifts and is not recommended. Global hybrid functionals are too insensitive to the environment and should be avoided, while nonhybrid functionals are strictly nonapplicable. Among the applicable approximate coupled cluster methods, the canonical versions of CC2 and ADC(2) were found to deviate significantly from the reference results both for the description of the lowest excited state and for the electrochromic shifts. By contrast, their spin-component-scaled (SCS) and particularly the scale-opposite-spin (SOS) variants compare well with the reference DLPNO-STEOM-CCSD and the best range-separated DFT methods. The emergence of RC excitation asymmetry is discussed in terms of intrinsic and protein electrostatic potentials. In addition, we evaluate a minimal structural scaffold of PSII, the D1–D2–CytB559 RC complex often employed in experimental studies, and show that it would have the same site energy distribution of RC chromophores as the full PSII supercomplex, but only under the unlikely conditions that the core protein organization and cofactor arrangement remain identical to those of the intact enzyme.

Brownian dynamics simulations of the interactions between lactate dehydrogenase (LDH) and G- or F-Actin. Part I: Muscle and heart homo-isoforms

Glycolytic enzymes may compartment in cells by binding to cytoskeletal structures. One potentially important step in compartmentation is the binding of the cytoskeletal protein F-actin by lactate dehydrogenase (LDH). Brownian dynamics (BD) simulations of LDH interacting with G- or F-actin provided first-encounter snapshots and relative binding free energies. Models of muscle (LDH-M4) and heart (LDH-H4) isoforms of lactate dehydrogenase for four species (rabbit, human, rat and pig) were examined. Strong electrostatic interactions were observed between both monomeric and polymerized actin for LDH-M4, but not LDH-H4. These electrostatic field potentials were examined in detail and explained the differences in complex formations between LDH-M4 and LDH-H4 with G- and F-actin. Complexes formed between LDH-M4 and G-actin involved residues found in the surface grooves between the A/D and B/C subunits. Several binding modes occurred between LDH-M4 and F-actin, the most frequent of which involved two subunits of actin interacting with two subunits of LDH-M4. LDH-M4 interactions with F-actin relied on 10 residues, nine of which were conserved between species. All of these residues except two were located in the positive grooves between A/D and B/C subunits, where the protein's electrostatic potential is enhanced by its quaternary structure. LDH-H4, however, does not show unique binding modes nor does it have a favorable relative binding free energy compared to LDH-M4. Thus, BD supports experimental findings [Ehmann, J. D.; Hultin, H. O. Arch. Biochem. Biophys.1973, 154, 471], that LDH-M4 binds actin, but LDH-H4 does not.

CpHMD-Then-QM/MM Identification of the Amino Acids Responsible for the Anabaena Sensory Rhodopsin pH-Dependent Electronic Absorption Spectrum.

Anabaena Sensory Rhodopsin (ASR), a microbial photoactive protein featuring the retinal chromophore in two different conformations, exhibits a pH-dependent electronic absorption spectrum. Using the recently developed CpHMD-then-QM/MM multiscale protocol applied to ASR embedded in a membrane model, the pH-induced changes in its maximum absorption wavelength have been reproduced and analyzed. While the acidic tiny red-shift is essentially correlated with the deprotonation of an aspartic acid located on the ASR extracellular side, the larger blue-shift experimentally reported at pH values larger than 5 involves a cluster of titrating residues sitting on the cytoplasmic side. The ASR pH-dependent spectrum is the consequence of the competitive stabilization of retinal ground and excited states by the protein electrostatic potential.

Electrostatic Complementarity as a Fast and Effective Tool to Optimize Binding and Selectivity of Protein-Ligand Complexes.

Electrostatic interactions between small molecules and their respective receptors are essential for molecular recognition and are also key contributors to the binding free energy. Assessing the electrostatic match of protein-ligand complexes therefore provides important insights into why ligands bind and what can be changed to improve binding. Ideally, the ligand and protein electrostatic potentials at the protein-ligand interaction interface should maximize their complementarity while minimizing desolvation penalties. In this work, we present a fast and efficient tool to calculate and visualize the electrostatic complementarity (EC) of protein-ligand complexes. We compiled benchmark sets demonstrating electrostatically driven structure-activity relationships (SAR) from literature data, including kinase, protein-protein interaction, and GPCR targets, and used these to demonstrate that the EC method can visualize, rationalize, and predict electrostatically driven ligand affinity changes and help to predict compound selectivity. The methodology presented here for the analysis of EC is a powerful and versatile tool for drug design.

On the control of the proton current in the voltage-gated proton channel Hv1

The nature of the action of voltage-activated proton transport proteins is a conundrum of great current interest. Here we approach this issue by exploring the action of Hv1, a voltage-gated proton channel found in different cells in humans and other organisms. Our study focuses on evaluating the free energy of transporting a proton through the channel, as well as the effect of the proton transfer through D112, in both the closed and open channel conformations. It is found that D112 allows a transported proton to bypass the electrostatic barrier of the open channel, while not being able to help in passing the barrier in the closed form. This reflects the change in position of the gating arginine residues relative to D112, upon voltage activation. Significantly, the effect of D112 accounts for the observed trend in selectivity by overcoming the electrostatic barrier at its highest point. Thus, the calculations provide a structure/function correlation for the Hv1 system. The present work also clarifies that the action of Hv1 is not controlled by a Grotthuss mechanism but, as is always the case, by the protein electrostatic potential at the rate-limiting barriers.

Effects of aligned α-helix peptide dipoles on experimental electrostatic potentials.

Aligned protein α-helix dipoles have been implicated in protein function and structure. The recent breakthroughs in high-resolution electron microscopy (EM) of macromolecules makes it possible to explore fundamental aspects of structural biology at the detailed molecular level. The electrostatic potential (ESP) generated by aligned protein α-helix dipole should be observable in high-resolution EM maps despite the fact that the effect may be partially screened by induced electric fields. Here, we show that aligned backbone dipoles in protein α-helices account for long-range features in the protein ESP functions. Our results are consistent with experimental EM maps and density functional theory calculations, including direct Fourier summation for proper calculation of the ESP due to the nonlocal nature of the ESP function from aligned dipoles and other partial atomic charges.

Biological effect of LOXL1 coding variants associated with pseudoexfoliation syndrome

Pseudoexfoliation (PEX) syndrome is a systemic disease involving the extracellular matrix. It increases the risk of glaucoma, an irreversible cause of blindness, and susceptibility to heart disease, stroke and hearing loss. Single nucleotide polymorphisms (SNPs) in the LOXL1 (Lysyl oxidase-like 1) gene are the major known genetic risk factor for PEX syndrome. Two coding SNPs, rs1048861 (G > T; Arg141Leu) and rs3825942 (G > A; Gly153Asp), in the LOXL1 gene are strongly associated with the disease risk in multiple populations worldwide. In the present study, we investigated functional effects of these SNPs on the LOXL1 protein. We show through molecular modelling that positions 141 and 153 are likely surface residues and hence possible recognition sites for protein–protein interactions; the Arg141Leu and Gly153Asp substitutions cause charge changes that would lead to local differences in protein electrostatic potential and in turn the potential to modify protein–protein interactions. In RFL-6 rat fetal lung fibroblast cells ectopically expressing the LOXL1 protein variants related to PEX (Arg141_Gly153, Arg141_Asp153 or Leu141_Gly153), immunoprecipitation of the secreted variants showed differences in their processing by endogenous proteins, possibly Bone morphogenetic protein-1 (BMP-1) that cleaves and leads to enzymatic activation of LOXL1. Immunofluorescence labelling of the ectopically expressed protein variants in RFL-6 cells showed no significant difference in their extracellular accumulation tendency. In conclusion, this is the first report of a biological effect of the coding SNPs in the LOXL1 gene associated with PEX syndrome, on the LOXL1 protein. The findings indicate that the disease associated coding variants themselves may be involved in the manifestation of PEX syndrome.

Abstracts of the 17th International Symposium on Bioluminescence and Chemiluminescence ‐ (ISBC 2012)

WARNING : The light-emitting molecular structures responsible for the chemiluminescence and fluorescence phenomena are not necessarily the same!

Theoretical study of the excited states of the photosynthetic reaction center in photosystem II: Electronic structure, interactions, and their origin

The excited states of the chlorophyll 6-mer in the photosystem II (PSII) reaction center (RC) were investigated theoretically using ab initio quantum chemical calculations, and the results are compared with those of the bacterial reaction center (bRC). A significant difference in the peak at the lowest energy in the absorption spectra arises from the structural asymmetry of the special pair (SP). The origin can be traced back to the structural difference in the CD helix. The low-lying excited states are characterized as a linear combination of the excited states of the chlorophyll monomers, which verifies the applicability of exciton theory. Analysis of the molecular interactions clearly explains the cause of the constructive/destructive interferences in the state transition moment. The protein electrostatic potential (ESP) decreases the energy of the charge-transfer (Chl D1 → Pheo D1) state. The ESP also localizes the HOMO distribution to the P D1 moiety and increases the ionization potential.

Excited States of Fluorescent Proteins, mKO and DsRed: Chromophore−Protein Electrostatic Interaction Behind the Color Variations

The emitting states of green fluorescent protein (GFP), monomeric Kusabira orange (mKO), and Discosoma red (DsRed) were studied using QM/MM and SAC-CI methods. By comparing the electronic structures among the green-, orange-, and red-emitting states as well as their electrostatic and quantum mechanical interactions within the protein cavity, the basic mechanisms for determining emission colors have been clarified. We found that the orange and red emissions of mKO and DsRed, respectively, result from cancellation between two effects, the pi skeleton extension (red shift) and protein electrostatic potential (blue shift). The extension of the pi skeleton enhances the intramolecular charge-transfer character of the transition, which makes the fluorescence energy more sensitive to the protein's electrostatic potential. On the basis of this mechanism, we predicted amino acid mutations that could red shift the emission energy of DsRed. A novel single amino acid mutation, which was examined computationally, reduced the DsRed emission energy from 2.14 (579 nm) to 1.95 eV (636 nm), which is approaching near-infrared fluorescence.

Cross-species analysis of the glycolytic pathway by comparison of molecular interaction fields

The electrostatic potential of an enzyme is a key determinant of its substrate interactions and catalytic turnover. Here we invoke comparative analysis of protein electrostatic potentials, along with sequence and structural analysis, to classify and characterize all the enzymes in an entire pathway across a set of different organisms. The electrostatic potentials of the enzymes from the glycolytic pathway of 11 eukaryotes were analyzed by qPIPSA (quantitative protein interaction property similarity analysis). The comparison allows the functional assignment of neuron-specific isoforms of triosephosphate isomerase from zebrafish, the identification of unusual protein surface interaction properties of the mosquito glucose-6-phosphate isomerase and the functional annotation of ATP-dependent phosphofructokinases and cofactor-dependent phosphoglycerate mutases from plants. We here show that plants possess two parallel pathways to convert glucose. One is similar to glycolysis in humans, the other is specialized to let plants adapt to their environmental conditions. We use differences in electrostatic potentials to estimate kinetic parameters for the triosephosphate isomerases from nine species for which published parameters are not available. Along the core glycolytic pathway, phosphoglycerate mutase displays the most conserved electrostatic potential. The largest cross-species variations are found for glucose-6-phosphate isomerase, enolase and fructose-1,6-bisphosphate aldolase. The extent of conservation of electrostatic potentials along the pathway is consistent with the absence of a single rate-limiting step in glycolysis.

WebPIPSA: a web server for the comparison of protein interaction properties

Protein molecular interaction fields are key determinants of protein functionality. PIPSA (Protein Interaction Property Similarity Analysis) is a procedure to compare and analyze protein molecular interaction fields, such as the electrostatic potential. PIPSA may assist in protein functional assignment, classification of proteins, the comparison of binding properties and the estimation of enzyme kinetic parameters. webPIPSA is a web server that enables the use of PIPSA to compare and analyze protein electrostatic potentials. While PIPSA can be run with downloadable software (see http://projects.eml.org/mcm/software/pipsa), webPIPSA extends and simplifies a PIPSA run. This allows non-expert users to perform PIPSA for their protein datasets. With input protein coordinates, the superposition of protein structures, as well as the computation and analysis of electrostatic potentials, is automated. The results are provided as electrostatic similarity matrices from an all-pairwise comparison of the proteins which can be subjected to clustering and visualized as epograms (tree-like diagrams showing electrostatic potential differences) or heat maps. webPIPSA is freely available at: http://pipsa.eml.org.

Optimization of cutting schemes for the evaluation of molecular electrostatic potentials in proteins via Moving-Domain QM/MM

This work presents new developments of the moving-domain QM/MM (MoD-QM/MM) method for modeling protein electrostatic potentials. The underlying goal of the method is to map the electronic density of a specific protein configuration into a point-charge distribution. Important modifications of the general strategy of the MoD-QM/MM method involve new partitioning and fitting schemes and the incorporation of dynamic effects via a single-step free energy perturbation approach (FEP). Selection of moderately sized QM domains partitioned between C (alpha) and C (from C=O), with incorporation of delocalization of electrons over neighboring domains, results in a marked improvement of the calculated molecular electrostatic potential (MEP). More importantly, we show that the evaluation of the electrostatic potential can be carried out on a dynamic framework by evaluating the free energy difference between a non-polarized MEP and a polarized MEP. A simplified form of the potassium ion channel protein Gramicidin-A from Bacillus brevis is used as the model system for the calculation of MEP.

A Self-Consistent Space-Domain Decomposition Method for QM/MM Computations of Protein Electrostatic Potentials.

This paper introduces a self-consistent computational protocol for modeling protein electrostatic potentials according to static point-charge model distributions. The protocol involves a simple space-domain decomposition scheme where individual molecular domains are modeled as Quantum-Mechanical (QM) layers embedded in the otherwise classical Molecular-Mechanics (MM) protein environment. ElectroStatic-Potential (ESP) atomic charges of the constituent molecular domains are computed, to account for mutual polarization effects, and iterated until obtaining a self-consistent point-charge model of the protein electrostatic potential. The novel protocol achieves quantitative agreement with full QM calculations in the description of electrostatic potentials of small polypeptides where polarization effects are significant, showing a remarkable improvement relative to the corresponding electrostatic potentials obtained with popular MM force fields. The capabilities of the method are demonstrated in several applications, including calculations of the electrostatic potential in the potassium channel protein and the description of protein-protein electrostatic interactions.

Electrostatic Potential of Nucleotide-free Protein is Sufficient for Discrimination Between Adenine and Guanine-specific Binding Sites

Despite sharing many common features, adenine-binding and guanine-binding sites in proteins often show a clear preference for the cognate over the non-cognate ligand. We have analyzed electrostatic potential (ESP) patterns at adenine and guanine-binding sites of a large number of non-redundant proteins where each binding site was first annotated as adenine/guanine-specific or non-specific from a survey of primary literature. We show that more than 90% of ESP variance at the binding sites is accounted for by only two principal component ESP vectors, each aligned to molecular dipoles of adenine and guanine. Projected on these principal component vectors, the adenine/guanine-specific and non-specific binding sites, including adenine-containing dinucleotides, show non-overlapping distributions. Adenine or guanine specificities of the binding sites also show high correlation with the corresponding electrostatic replacement (cognate by non-cognate ligand) energies. High correlation coefficients (0.94 for 35 adenine-binding sites and 1.0 for 20 guanine-binding sites) were obtained when adenine/guanine specificities were predicted using the replacement energies. Our results demonstrate that ligand-free protein ESP is an excellent indicator for discrimination between adenine and guanine-specific binding sites and that ESP of ligand-free protein can be used as a tool to annotate known and putative purine-binding sites in proteins as adenine or guanine-specific.

Species dependence of enzyme-substrate encounter rates for triose phosphate isomerases

Triose phosphate isomerase (TIM) is a diffusion-controlled enzyme whose rate is limited by the diffusional encounter of the negatively charged substrate glyceraldehyde 3-phosphate (GAP) with the homodimeric enzyme's active sites. Translational and orientational steering of GAP toward the active sites by the electrostatic field of chicken muscle TIM has been observed in previous Brownian dynamics (BD) simulations. Here we report simulations of the association of GAP with TIMs from four species with net charges at pH 7 varying from -12e to +12e. Computed second-order rate constants are in good agreement with experimental data. The BD simulations and computation of average Boltzmann factors of substrate-protein interaction energies show that the protein electrostatic potential enhances the rates for all the enzymes. There is much less variation in the computed rates than might be expected on the basis of the net charges. Comparison of the electrostatic potentials by means of similarity indices shows that this is due to conservation of the local electrostatic potentials around the active sites which are the primary determinants of electrostatic steering of the substrate.