Importance of explicit salt ions for protein stability in molecular dynamics simulation.

The accurate and efficient treatment of electrostatic interactions is one of the challenging problems of molecular dynamics simulation. Truncation procedures such as switching or shifting energies or forces lead to artifacts and significantly reduced accuracy. The particle mesh Ewald (PME) method is one approach to overcome these problems by providing a computationally efficient means of calculating all long-range electrostatic interactions in a periodic simulation box by use of fast Fourier transformation techniques. For the application of the PME method to the simulation of a protein with a net charge in aqueous solution, counterions are added to neutralize the system. The usual procedure is to add charge-balancing counterions close to charged residues to neutralize the protein surface. In the present article, we show that for MD simulation of a small protein of marginal stability, the YAP-WW domain, explicit modeling of 0.2 M ionic strength (in addition to the charge-balancing counterions) is necessary to maintain a stable protein structure. Without explicit ions throughout the periodic simulation box, the charge-balancing counterions on the protein surface diffuse away from the protein, resulting in destruction of the beta-sheet secondary structure of the WW domain.     

Peptide Binding to a PDZ Domain by Electrostatic Steering via Nonnative Salt Bridges

We have captured the binding of a peptide to a PDZ domain by unbiased molecular dynamics simulations. Analysis of the trajectories reveals on-pathway encounter complex formation, which is driven by electrostatic interactions between negatively charged carboxylate groups in the peptide and positively charged side chains surrounding the binding site. In contrast, the final stereospecific complex, which matches the crystal structure, features completely different interactions, namely the burial of the hydrophobic side chain of the peptide C-terminal residue and backbone hydrogen bonds. The simulations show that nonnative salt bridges stabilize kinetically the encounter complex during binding. Unbinding follows the inverse sequence of events with the same nonnative salt bridges in the encounter complex. Thus, in contrast to protein folding, which is driven by native interactions, the binding of charged peptides can be steered by nonnative interactions, which might be a general mechanism, e.g., in the recognition of histone tails by bromodomains.     

Interactions of the pleckstrin homology domain with phosphatidylinositol phosphate and membranes: characterization via molecular dynamics simulations

The mechanism of interaction of pleckstrin homology (PH) domains with phosphatidylinositol 4,5-bisphosphate (PIP 2)-containing lipid bilayers remains uncertain. While crystallographic studies have emphasized PH-inositol 1,4,5-trisphosphate (IP 3) interactions, biophysical studies indicate a degree of less specific protein-bilayer interactions. We have used molecular dynamics simulations to characterize the interactions of the PH domain from phospholipase C-delta1 with IP 3 and with PIP 2, the latter in lipid bilayers and in detergent micelles. Simulations of the PH domain in water reveal a reduction in protein flexibility when IP 3 is bound. Simulations of the PH domain bound to PIP 2 in lipid bilayers indicate a tightening of ligand-protein interactions relative to the PH-IP 3 complex, alongside formation of H-bonds between PH side chains and lipid (PC) headgroups, and a degree of penetration of hydrophobic side chains into the core of the bilayer. Comparison with simulations of the PH-bound domain to a PC bilayer in the absence of PIP 2 suggests that the presence of PIP 2 increases the extent of PH-membrane interactions. Thus, comparative molecular dynamics simulations reveal how a PI-binding domain undergoes changes in conformational dynamics on binding to a PIP 2-containing membrane and how interactions additional to those with the PI headgroup are formed.     

Structure of a PH domain from the C. elegans muscle protein UNC-89 suggests a novel function

BACKGROUND: Pleckstrin homology (PH) domains constitute a structurally conserved family present in many signaling and regulatory proteins. PH domains have been shown to bind to phospholipids, and many function in membrane targeting. They generally have a strong electrostatic polarization and interact with negatively charged phospholipids via the positive pole. On the basis of electrostatic modeling, however, we have previously identified a class of PH domains with a predominantly negative charge and predicted that these domains recognize other targets. Here, we report the first experimental structure of such a PH domain. RESULTS: The structure of the PH domain from Caenorhabditis elegans muscle protein UNC-89 has been determined by heteronuclear NMR. The domain adopts the classic PH fold, but has an unusual closed conformation of the "inositol binding loops. This creates a small opening to a deep hydrophobic pocket lined with negative charges on one side, and provides a molecular explanation for the lack of association with inositol-1,4,5-triphosphate. As predicted, the PH domain of UNC-89 has a strongly negative overall electrostatic potential. Modeling the Dbl homology (DH)-linked PH domains from the C. elegans genome shows that a large proportion of these modules are negatively charged. CONCLUSIONS: We present the first structure of a PH domain with a strong negative overall electrostatic potential. The presence of a deep pocket lined with negative charges suggests that the domain binds to ligands other than acidic phospholipids. The abundance of this class of PH domain in the C. elegans genome suggests a prominent role in mediating protein-protein interactions.     

Salt effects on hydrophobic‐core formation in folding of a helical miniprotein studied by molecular dynamics simulations

We have investigated effects of salt ions on folding events of a helical miniprotein chicken villin headpiece subdomain HP36. Low concentrations of ions alter electrostatic interactions between charged groups of a protein and can change the populations of conformers. Here, we compare two data sets of folding simulations of HP36 in explicit water solvent with or without ions. For efficient sampling of the conformational space of HP36, the multicanonical replica‐exchange molecular dynamics method was employed. Our analyses suggest that salt alters salt‐bridging nature of the protein at later stages of folding at room temperature. Especially, more nonnative, nonlocal salt bridges are formed at near‐native conformations in pure water. Our analyses also show that such salt‐bridge formation hinders the fully native hydrophobic‐core packing at the final stages of folding. Proteins 2014; 82:933–943. © 2013 Wiley Periodicals, Inc.     

Open-Boundary Molecular Dynamics of a DNA Molecule in a Hybrid Explicit/Implicit Salt Solution

The composition and electrolyte concentration of the aqueous bathing environment have important consequences for many biological processes and can profoundly affect the behavior of biomolecules. Nevertheless, because of computational limitations, many molecular simulations of biophysical systems can be performed only at specific ionic conditions: either at nominally zero salt concentration, i.e., including only counterions enforcing the system’s electroneutrality, or at excessive salt concentrations. Here, we introduce an efficient molecular dynamics simulation approach for an atomistic DNA molecule at realistic physiological ionic conditions. The simulations are performed by employing the open-boundary molecular dynamics method that allows for simulation of open systems that can exchange mass and linear momentum with the environment. In our open-boundary molecular dynamics approach, the computational burden is drastically alleviated by embedding the DNA molecule in a mixed explicit/implicit salt-bathing solution. In the explicit domain, the water molecules and ions are both overtly present in the system, whereas in the implicit water domain, only the ions are explicitly present and the water is described as a continuous dielectric medium. Water molecules are inserted and deleted into/from the system in the intermediate buffer domain that acts as a water reservoir to the explicit domain, with both water molecules and ions free to enter or leave the explicit domain. Our approach is general and allows for efficient molecular simulations of biomolecules solvated in bathing salt solutions at any ionic strength condition.     

Do salt bridges stabilize proteins? A continuum electrostatic analysis

The electrostatic contribution to the free energy of folding was calculated for 21 salt bridges in 9 protein X-ray crystal structures using a continuum electrostatic approach with the DELPHI computer-program package. The majority (17) were found to be electrostatically destabilizing; the average free energy change, which is analogous to mutation of salt bridging side chains to hydrophobic isosteres, was calculated to be 3.5 kcal/mol. This is fundamentally different from stability measurements using pKa shifts, which effectively measure the strength of a salt bridge relative to 1 or more charged hydrogen bonds. The calculated effect was due to a large, unfavorable desolvation contribution that was not fully compensated by favorable interactions within the salt bridge and between salt-bridge partners and other polar and charged groups in the folded protein. Some of the salt bridges were studied in further detail to determine the effect of the choice of values for atomic radii, internal protein dielectric constant, and ionic strength used in the calculations. Increased ionic strength resulted in little or no change in calculated stability for 3 of 4 salt bridges over a range of 0.1-0.9 M. The results suggest that mutation of salt bridges, particularly those that are buried, to "hydrophobic bridges" (that pack at least as well as wild type) can result in proteins with increased stability. Due to the large penalty for burying uncompensated ionizable groups, salt bridges could help to limit the number of low free energy conformations of a molecule or complex and thus play a role in determining specificity (i.e., the uniqueness of a protein fold or protein-ligand binding geometry).     

Reaction within the coulomb-cage; science in retrospect

The Coulomb-cage is defined as the space where the electrostatic interaction between two bodies is more intensive than the thermal energy (k_BT). For small molecule, the Coulomb-cage is a small sphere, extending only few water molecules towards the bulk and its radius is sensitive to the ionic strength of the solution. For charged proteins or membranal structures, the Coulomb-cage can engulf large fraction of the surface and provides a preferred pathway for ion propagation along the surface. Similarly, electrostatic potential at the inner space of a channel can form preferential trajectories passage for ions. The dynamics of ions inside the Coulomb-cage of ions was formulated by the studies of proton-anion recombination of excited photoacids. In the present article, we recount the study of intra- Coulomb-cage reaction taking place on the surface of macro-molecular bodies like micelles, membranes, proteins and intra-protein cavities. The study progressed stepwise, tracing the dynamics of a proton ejected from a photo-acid molecule located at defined sites (on membrane, inter-membrane space, active site of enzyme, inside Large Pore Channels etc.). Accumulation of experimental observations encouraged us to study of the reaction mechanism by molecular dynamics simulations of ions within the Coulomb-cage of proteins surface or inside large pores. The intra-Coulomb-cage proton transfer events follows closely the fine structure of the electrostatic field inside the cage and reflects the shape of nearby dielectric boundaries, the temporal ordering of the solvent molecules and the structural fluctuations of the charged side chains. The article sums some 40 years of research, which in retrospect clarifies the intra-Coulomb-cage reaction mechanism.     

Contribution of arginine-glutamate salt bridges to helix stability

Peptide side chain interactions were studied by molecular dynamics simulation using explicit solvent on a peptide with the sequence AAARAAAAEAAEAAAARA. Three different protonation states of the glutamic acid side chains were simulated for four 20 ns runs each, a total simulation time of 240 ns. Two different salt bridge geometries were observed and the preferred geometry was found to depend on Glu — Arg residue spacing. Stable charge clusters were also observed, particularly in the fully charged peptide. Salt bridges were selectively interrupted upon protonation, with concomitant changes in secondary structure. The fully charged peptide was highly helical between residues 9 and 13, although protonation increased helicity near the N-terminus. The contribution of salt bridges to helix stability therefore depends on both position and relative position of charged residues within a sequence.     

The role of positively charged amino acids and electrostatic interactions in the complex of U1A protein and U1 hairpin II RNA

Previous kinetic investigations of the N-terminal RNA recognition motif (RRM) domain of spliceosomal protein U1A, interacting with its RNA target U1 hairpin II, provided experimental evidence for a ‘lure and lock’ model of binding in which electrostatic interactions first guide the RNA to the protein, and close range interactions then lock the two molecules together. To further investigate the ‘lure’ step, here we examined the electrostatic roles of two sets of positively charged amino acids in U1A that do not make hydrogen bonds to the RNA: Lys20, Lys22 and Lys23 close to the RNA-binding site, and Arg7, Lys60 and Arg70, located on ‘top’ of the RRM domain, away from the RNA. Surface plasmon resonance-based kinetic studies, supplemented with salt dependence experiments and molecular dynamics simulation, indicate that Lys20 predominantly plays a role in association, while nearby residues Lys22 and Lys23 appear to be at least as important for complex stability. In contrast, kinetic analyses of residues away from the RNA indicate that they have a minimal effect on association and stability. Thus, well-positioned positively charged residues can be important for both initial complex formation and complex maintenance, illustrating the multiple roles of electrostatic interactions in protein–RNA complexes.     

Effect of ion binding on protein transport through ultrafiltration membranes

Electrostatic interactions can have a significant impact on protein transmission through semipermeable membranes. Experimental data for the transport of bovine serum albumin (BSA) through a polyethersulfone ultrafiltration membrane were obtained in different salt solutions over a range of pH and salt concentrations. Net BSA charge under the same conditions was evaluated from mobility data measured by capillary electrophoresis. The results show that specific ionic composition, in addition to solution pH and ionic strength, can strongly affect the rate of protein transport through semipermeable ultrafiltration membranes. The effects of different ions on BSA sieving are due primarily to differences in ion binding to the protein, which leads to significant differences in the net protein charge at a given pH and ionic strength. This effect could be described in terms of an effective protein radius, which accounts for the electrostatic exclusion of the charged protein from the membrane pores. These results provide important insights into the nature of the electrostatic interactions in membrane systems.     

Deciphering the factors responsible for the stability of a GFP variant resistant to alkaline pH using molecular dynamics simulations

Charged amino acids having ionizable side chains play crucial roles in maintaining the solubility and stability of a protein. These charged amino acids are mostly exposed on protein surface and participate in electrostatic interactions with neighboring charged amino acids as well as with solvent. Therefore, the change in the solvent pH affects the protein stability in most cases. Previously, we reported a GFP variant, GFP14R having 14 surface lysines replaced with arginines, that showed enhanced stability under alkaline pH. Here, we analyzed the factors that contribute to the stability of the GFP14R under alkaline pH quantitatively using molecular dynamics simulations. Protonation state of the charged amino acids of GFP14R and control GFP under neutral pH and alkaline pH were modeled, and molecular dynamics simulations were performed. This comparative analysis revealed that the GFP14R with more arginine frequency on the surface maintained the stability under both pH conditions without much change in their salt-bridge interactions as well as the hydrogen bond interactions with solvent. On the other hand, these interactions were significantly reduced for the control GFP under alkaline pH due to the deprotonated lysine side chains. These results suggest that the advantageous property of arginine over lysine can be considered one of the parameter for the protein stability engineering under alkaline pH conditions.     

An electrostatic/hydrogen bond switch as the basis for the specific interaction of phosphatidic acid with proteins

Phosphatidic acid (PA) is a minor but important phospholipid that, through specific interactions with proteins, plays a central role in several key cellular processes. The simple yet unique structure of PA, carrying just a phosphomonoester head group, suggests an important role for interactions with the positively charged essential residues in these proteins. We analyzed by solid-state magic angle spinning 31P NMR and molecular dynamics simulations the interaction of low concentrations of PA in model membranes with positively charged side chains of membrane-interacting peptides. Surprisingly, lysine and arginine residues increase the charge of PA, predominantly by forming hydrogen bonds with the phosphate of PA, thereby stabilizing the protein-lipid interaction. Our results demonstrate that this electrostatic/hydrogen bond switch turns the phosphate of PA into an effective and preferred docking site for lysine and arginine residues. In combination with the special packing properties of PA, PA may well be nature's preferred membrane lipid for interfacial insertion of positively charged membrane protein domains.     

Classification of protein sequences by homology modeling and quantitative analysis of electrostatic similarity

Protein electrostatics plays a key role in ligand binding and protein-protein interactions. Therefore, similarities or dissimilarities in electrostatic potentials can be used as indicators of similarities or dissimilarities in protein function. We here describe a method to compare the electrostatic properties within protein families objectively and quantitatively. Three-dimensional structures are built from database sequences by comparative modeling. Molecular potentials are then computed for these with a continuum solvation model by finite difference solution of the Poisson-Boltzmann equation or analytically as a multipole expansion that permits rapid comparison of very large datasets. This approach is applied to 104 members of the Pleckstrin homology (PH) domain family. The deviation of the potentials of the homology models from those of the corresponding experimental structures is comparable to the variation of the potential in an ensemble of structures from nuclear magnetic resonance data or between snapshots from a molecular dynamics simulation. For this dataset, the results for analysis of the full electrostatic potential and the analysis using only monopole and dipole terms are very similar. The electrostatic properties of the PH domains are generally conserved despite the extreme sequence divergence in this family. Notable exceptions from this conservation are seen for PH domains linked to a Db1 homology (DH) domain and in proteins with internal PH domain repeats.     

How does protein phosphorylation regulate photosynthesis?

Phosphorylation of light-harvesting antenna proteins redirects absorbed light energy between reaction centres of photosynthetic membranes. A generally accepted explanation for this is that electrostatic forces drive the more negatively charged, phosphorylated antenna proteins between membrane domains that differ in surface charge. However, structural studies on soluble phosphoproteins indicate that phosphorylated amino acid side chains have specific effects on molecular recognition, by ligand blocking or by intramolecular interactions which alter protein structure. These studies suggest alternative mechanisms for phosphorylation in control of pairwise protein-protein interactions in biological membranes. Thus, in photosynthesis, the surface charge model is only one possible interpretation.     

Molecular dynamics study of the interaction between fatty acid binding proteins with palmitate mini-micelles

The fatty acid binding proteins (FAPBs) function as intracellular carriers of fatty acid (FA) and related compounds. During the digestion of lipids, the local concentration of FA exceeds their critical micellar concentration; the excess ratio of FA/FABP can be as high as ~1,000/1, consequently building micelles. Considering that the micelle formation is a rapid process, the FABP must be able to remove the mini-micelle. In this study, we describe the results of molecular dynamics simulations of liver basic FABP (Lb-FABP), carried out in the presence of ~20 mM palmitate ions, all in the presence of explicit water and at ionic strength of ~100 mM, approximating physiological conditions. The Lb-FABP appears to react, along with a free FA, with mini-micelle creating a stable complex (on the time scale of the simulations), which is attached to the anti-portal domain of the protein. The complex may be formed by the stepwise addition of free FA or through the interaction of a pre-formed mini-micelle with the free protein. The driving force of the mini-micelle-FABP complex is a combination of electrostatic attraction between the negative carboxylates of the mini-micelle with the positive charge of the N terminal amine residues and Lennard-Jones FA–protein interactions. The preferred tendency of the mini-micelle to react with the anti-portal domain retains the α-helixes of the portal region free for its electrostatic interaction with the membrane, ensuring a rapid unloading of the cargo on the membrane.     

Exposure of hydrophobic core in human prion protein pathogenic mutant H187R

Pathogenesis studies have revealed that H187R mutation of human prion protein (huPrP) is related to GSS type of TSE diseases. Its pathogenic mechanism is still unclear. We here studied the globular domain of this mutant protein by molecular dynamics simulations. Compared to the wide-type protein, the mutant has similar dynamics and stability profiles in our simulation. Conformational rearrangements are concentrated around the mutation site, due to the introduction the positively charged side chain of Arg187. The strong electrostatic repulsion between Arg156 and Arg187 drives both side chains away from their original positions, leaving its hydrophobic core to be solvent accessible. Such a unfavorable conformational change may destabilize the mutant protein and make it more susceptible to unfolding.     

Steric selectivity in Na channels arising from protein polarization and mobile side chains

Monte Carlo simulations of equilibrium selectivity of Na channels with a DEKA locus are performed over a range of radius R and protein dielectric coefficient epsilon(p). Selectivity arises from the balance of electrostatic forces and steric repulsion by excluded volume of ions and side chains of the channel protein in the highly concentrated and charged (approximately 30 M) selectivity filter resembling an ionic liquid. Ions and structural side chains are described as mobile charged hard spheres that assume positions of minimal free energy. Water is a dielectric continuum. Size selectivity (ratio of Na+ occupancy to K+ occupancy) and charge selectivity (Na+ to Ca2+) are computed in concentrations as low as 10(-5) M Ca2+. In general, small R reduces ion occupancy and favors Na+ over K+ because of steric repulsion. Small epsilon(p) increases occupancy and favors Na+ over Ca2+ because protein polarization amplifies the pore's net charge. Size selectivity depends on R and is independent of epsilon(p); charge selectivity depends on both R and epsilon(p). Thus, small R and epsilon(p) make an efficient Na channel that excludes K+ and Ca2+ while maximizing Na+ occupancy. Selectivity properties depend on interactions that cannot be described by qualitative or verbal models or by quantitative models with a fixed free energy landscape.     

How Does a Voltage Sensor Interact with a Lipid Bilayer? Simulations of a Potassium Channel Domain

The nature of voltage sensing by voltage-activated ion channels is a key problem in membrane protein structural biology. The way in which the voltage-sensor (VS) domain interacts with its membrane environment remains unclear. In particular, the known structures of Kv channels do not readily explain how a positively charged S4 helix is able to stably span a lipid bilayer. Extended (2 x 50 ns) molecular dynamics simulations of the high-resolution structure of the isolated VS domain from the archaebacterial potassium channel KvAP, embedded in zwitterionic and in anionic lipid bilayers, have been used to explore VS/lipid interactions at atomic resolution. The simulations reveal penetration of water into the center of the VS and bilayer. Furthermore, there is significant local deformation of the lipid bilayer by interactions between lipid phosphate groups and arginine side chains of S4. As a consequence of this, the electrostatic field is "focused" across the center of the bilayer.