Protocol for MM/PBSA molecular dynamics simulations of proteins

Continuum solvent models have been employed in past years for understanding processes such as protein folding or biomolecular association. In the last decade, several attempts have been made to merge atomic detail molecular dynamics simulations with solvent continuum models. Among continuum models, the Poisson-Boltzmann solvent accessible surface area model is one of the oldest and most fundamental. Notwithstanding its wide usage for simulation of biomolecular electrostatic potential, the Poisson-Boltzmann equation has been very seldom used to obtain solvation forces for molecular dynamics simulation. We propose here a fast and reliable methodology to implement continuum forces in standard molecular mechanics and dynamics algorithms. Results for a totally unrestrained 1 ns molecular dynamics simulation of a small protein are quantitatively similar to results obtained by explicit solvent molecular dynamics simulations.     

Electrostatic contribution to the bending of DNA

A model is derived that accounts for the short-range electrostatic contribution to the bending of DNA molecule in solution and in complexes with proteins in terms of the non-linear Poisson-Boltzmann equation. We defined that the short-range electrostatic interactions depend on the changes of the polyion surface charge density under deformation, while the long-range interactions depend on the bending-induced changes in distances between each two points along the polyion axis. After an appropriate simplification of the Poisson-Boltzmann equation, the short-range term is calculated separately giving the lower limit for the electrostatic contribution to the DNA persistence length. The result is compared with the theoretical approaches developed earlier [M. Fixman, J. Chem. Phys. 76 (1982) 6346; M. Le Bret, J. Chem. Phys. 76 (1982) 6243] and with the experimental data. The conclusion is made that the results of Fixman-Le Bret, which took into account both types of the electrostatic interactions for a uniformly bent polyion, give the upper limit for the electrostatic persistence length at low ionic strength, and the actual behavior of the DNA persistence length lies between two theoretical limits. Only the short-range term is significant at moderate-to-high ionic strength where our results coincide with the predictions of Fixman-Le Bret. The bending of DNA on the protein surface that is accompanied by an asymmetric neutralization of the DNA charge is also analyzed. In this case, the electrostatic bending energy gives a significant favorite contribution to the total bending energy of DNA. Important implications to the mechanisms of DNA-protein interactions, particularly in the nucleosome particle, are discussed.     

APBSmem: A Graphical Interface for Electrostatic Calculations at the Membrane

Electrostatic forces are one of the primary determinants of molecular interactions. They help guide the folding of proteins, increase the binding of one protein to another and facilitate protein-DNA and protein-ligand binding. A popular method for computing the electrostatic properties of biological systems is to numerically solve the Poisson-Boltzmann (PB) equation, and there are several easy-to-use software packages available that solve the PB equation for soluble proteins. Here we present a freely available program, called APBSmem, for carrying out these calculations in the presence of a membrane. The Adaptive Poisson-Boltzmann Solver (APBS) is used as a back-end for solving the PB equation, and a Java-based graphical user interface (GUI) coordinates a set of routines that introduce the influence of the membrane, determine its placement relative to the protein, and set the membrane potential. The software Jmol is embedded in the GUI to visualize the protein inserted in the membrane before the calculation and the electrostatic potential after completing the computation. We expect that the ease with which the GUI allows one to carry out these calculations will make this software a useful resource for experimenters and computational researchers alike. Three examples of membrane protein electrostatic calculations are carried out to illustrate how to use APBSmem and to highlight the different quantities of interest that can be calculated.     

Brownian dynamics simulation of protein solutions: structural and dynamical properties

The study of solutions of biomacromolecules provides an important basis for understanding the behavior of many fundamental cellular processes, such as protein folding, self-assembly, biochemical reactions, and signal transduction. Here, we describe a Brownian dynamics simulation procedure and its validation for the study of the dynamic and structural properties of protein solutions. In the model used, the proteins are treated as atomically detailed rigid bodies moving in a continuum solvent. The protein-protein interaction forces are described by the sum of electrostatic interaction, electrostatic desolvation, nonpolar desolvation, and soft-core repulsion terms. The linearized Poisson-Boltzmann equation is solved to compute electrostatic terms. Simulations of homogeneous solutions of three different proteins with varying concentrations, pH, and ionic strength were performed. The results were compared to experimental data and theoretical values in terms of long-time self-diffusion coefficients, second virial coefficients, and structure factors. The results agree with the experimental trends and, in many cases, experimental values are reproduced quantitatively. There are no parameters specific to certain protein types in the interaction model, and hence the model should be applicable to the simulation of the behavior of mixtures of macromolecules in cell-like crowded environments.     

Energetics of charge-charge interactions in proteins

Electrostatic interactions between pairs of atoms in proteins are calculated with a model based on the linearized Poisson-Boltzmann equation. The equation is solved accurately by a method that takes into account the detailed shape of the protein. This paper presents applications to several systems. Experimental data for the interaction of ionized residues with an active site histidine in subtilisin BPN' allow the model to be tested, using various assumptions for the electrical properties of the protein and solvent. The electrostatic stabilization of the active site thiolate of rhodanese is analyzed, with attention to the influence of alpha-helices. Finally, relationships between electrostatic potential and charge-charge distance are reported for large and small globular proteins. The above results are compared with those of simpler electrostatic models, including Coulomb's law with both a distance-dependent dielectric constant (epsilon = R) and a fixed dielectric constant (epsilon = 2), and Tanford-Kirkwood theory. The primary conclusions are as follows: 1) The Poisson-Boltzmann model agrees with the subtilisin data over a range of ionic strengths; 2) two alpha-helices generate a large potential in the active site of rhodanese; 3) epsilon = R overestimates weak electrostatic interactions but yields relatively good results for strong ones; 4) Tanford-Kirkwood theory is a useful approximation to detailed solutions of the linearized Poisson-Boltzmann equation in globular proteins; and 5) the modified Tanford-Kirkwood theory over-screens the measured electrostatic interactions in subtilisin.     

Contact pairs of RNA with magnesium ions-electrostatics beyond the Poisson-Boltzmann equation

The electrostatic interaction of RNA with its aqueous environment is most relevant for defining macromolecular structure and biological function. The attractive interaction of phosphate groups in the RNA backbone with ions in the water environment leads to the accumulation of positively charged ions in the first few hydration layers around RNA. Electrostatics of this ion atmosphere and the resulting ion concentration profiles have been described by solutions of the nonlinear Poisson-Boltzmann equation and atomistic molecular dynamics (MD) simulations. Much less is known on contact pairs of RNA phosphate groups with ions at the RNA surface, regarding their abundance, molecular geometry, and role in defining RNA structure. Here, we present a combined theoretical and experimental study of interactions of a short RNA duplex with magnesium (Mg²⁺) ions. MD simulations covering a microsecond time range give detailed hydration geometries as well as electrostatics and spatial arrangements of phosphate-Mg²⁺ pairs, including both pairs in direct contact and separated by a single water layer. The theoretical predictions are benchmarked by linear infrared absorption and nonlinear two-dimensional infrared spectra of the asymmetric phosphate stretch vibration which probes both local interaction geometries and electric fields. Contact pairs of phosphate groups and Mg²⁺ ions are identified via their impact on the vibrational frequency position and line shape. A quantitative analysis of infrared spectra for a range of Mg²⁺-excess concentrations and comparison with fluorescence titration measurements shows that on average 20–30% of the Mg²⁺ ions interacting with the RNA duplex form contact pairs. The experimental and MD results are in good agreement. In contrast, calculations based on the nonlinear Poisson-Boltzmann equation fail in describing the ion arrangement, molecular electrostatic potential, and local electric field strengths correctly. Our results underline the importance of local electric field mapping and molecular-level simulations to correctly account for the electrostatics at the RNA-water interface.     

Protein Electrostatics: Rapid Multigrid-Based Newton Algorithm for Solution of the Full Nonlinear Poisson-Boltzmann Equation

Abstract

A new method for solving the full nonlinear Poisson-Boltzmann equation is outlined. This method is robust and efficient, and uses a combination of the multigrid and inexact Newton algorithms. The novelty of this approach lies in the appropriate combination of the two methods, neither of which by themselves are capable of solving the nonlinear problem accurately. Features of the Poisson-Boltzmann equation are fully exploited by each component of the hybrid algorithm to provide robustness and speed. The advantages inherent in this method increase with the size of the problem. The efficacy of the method is illustrated by calculations of the electrostatic potential around the enzyme Superoxide Dismutase. The CPU time required to solve the full nonlinear equation is less than half that needed for a conjugate gradient solution of the corresponding linearized Poisson-Boltzmann equation. The solutions reveal that the field around the active sites is significantly reduced as compared to that obtained by solving the corresponding linearized Poisson-Boltzmann equation. This new method for the nonlinear Poisson-Boltzmann equation will enable fast and accurate solutions of large protein electrostatics problems.     

蓖麻毒素A链突变体(MRTA)的分子设计

利用同源模建的方法,借助分子力学优化,分子力学模拟退火设计构建了删除部分氨基酸序列的蓖麻毒素Ａ链突变体。采用泊松－玻尔兹曼方程对比分析了蓖麻毒素Ａ链与ＭＲＴＡ表现静电势分布,研究了ＲＴＡ与ＭＲＴＡ蛋白表面静电性质;     

Surface electrostatic effects in oligonucleotide microarrays: control and optimization of binding thermodynamics

We present a theoretical thermodynamic framework for the design of more efficient oligonucleotide microarrays. A general thermodynamic relation is derived to describe the electrostatic surface effects on the binding of the assayed biomolecule to a surface-tethered molecular probe. The relation is applied to analyze how the nucleic acid target, the oligonuleotide probe, and their DNA duplex electrostatic interactions with the surface affect the hybridization on DNA arrays. Taking advantage of a closed form exact solution of the linear Poisson-Boltzmann equation for a charged ion-penetrable sphere in electrolyte solution interacting with a plane wall, we study the effects of the surface and solution conditions. Binding free energy is found as a function of the surface material, dielectric or metal, the surface charge density, linker molecule length, temperature, and added salt content. The charge or electric potential of the dielectric or metal surface, respectively, is shown to dominate the hybridization, especially at low added salt or short linker length. We predict that substantial enhancement of sensitivity, selectivity, and reliability of microarrays can be achieved by control of the surface conditions. As examples, we discuss how to overcome two limitations of current technologies: nonequal sensitivity of the probes with different GC and AT bases content, and poor match/mismatch discrimination. In addition, we suggest the design of microarray conditions where the tested nucleic acid is unfolded, thus making possible the screening of a larger sequence with single nucleotide resolution. These promising findings are discussed and further experimental tests suggested.     

蓖麻毒素A链突变体(MRTA)的分子设计

利用同源模建的方法,借助分子力学优化、分子动力学模拟退火设计构建了删除部分氨基酸序列的蓖麻毒素Ａ链突变体（ＭＲＴＡ）。采用泊松—玻尔兹曼方程对比分析了蓖麻毒素Ａ链（ＲＴＡ）与ＭＲＴＡ表观静电势分布,研究ＲＴＡ与ＭＲＴＡ蛋白表面静电性质;通过半经验量子化学ＡＭ１与分子力学结合方法探讨ＲＴＡ与ＭＲＴＡ功能域氨基酸前线分子轨道性质、能级分布,从理论上预测ＭＲＴＡ功能活性     

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.     

Poisson-Boltzmann analysis of the lambda repressor-operator interaction.

A theoretical study of the ion atmosphere contribution to the binding free energy of the lambda repressor-operator complex is presented. The finite-difference form of the Poisson-Boltzmann equation was solved to calculate the electrostatic interaction energy of the amino-terminal domain of the lambda repressor with a 9 or 45 base pair oligonucleotide. Calculations were performed at various distances between repressor and operator as well as at different salt concentrations to determine ion atmosphere contributions to the total electrostatic interaction. Details in the distribution of charges on DNA and protein atoms had a strong influence on the calculated total interaction energies. In contrast, the calculated salt contributions are relatively insensitive to changes in the details of the charge distribution. The results indicate that the ion atmosphere contribution favors association at all protein-DNA distances studied. The theoretical number of ions released upon repressor-operator binding appears to be in reasonable agreement with experimental data.     

Calculating the electrostatic properties of RNA provides new insights into molecular interactions and function.

Solutions to the nonlinear Poisson-Boltzmann equation were used to obtain the electrostatic potentials of RNA molecules that have known three-dimensional structures. The results are described in terms of isopotential contours and surface electrostatic potential maps. Both representations have unexpected features: 'cavities' within isopotential contours and areas of enhanced negative potential on molecular surfaces. Intriguingly, the sites of unusual electrostatic features correspond to functionally important regions, suggesting that electrostatic properties play a key role in RNA recognition and stabilization. These calculations reveal that the electrostatic potentials generated by RNA molecules have a variety of functionally important characteristics that cannot be discerned by simple visual inspection of the molecular structure.     

Concepts and protocols for electrostatic free energies

Electrostatic free energies play an essential role in numerous biomolecular processes occurring in solution. Difficulties arise when the long-range Coulomb interaction is computed for idealised infinite simulation models with periodic boundary conditions. To maintain a neutral simulation box and a finite per-box energy, a neutralising charge density or ‘gellium’ is commonly used, leading to a mean box potential that is constrained to be rigorously equal to zero at all times. Thus, in considering quantities such as ion solvation free energy, the potential drop to move from solvent into the usual, gas phase reference state is missing. In fact, for an infinite molecular system, the electrostatic potential itself is not uniquely defined, but takes the form of an infinite series that is only conditionally convergent. This leads to several possible computational conventions that give different values for the potential and field, all mathematically valid. For experimentally measurable quantities, however, unique results are obtained when sufficiently large simulation boxes are utilised. These concepts are detailed, as well as a fundamental, linear response theoretical framework that provides qualitative understanding of the physical processes involved, especially dielectric relaxation of the environment in response to a new solute charge. Illustrative applications to ligand binding and biomolecular electron transfer are described.     

Modeling salt-mediated electrostatics of macromolecules: the discrete surface charge optimization algorithm and its application to the nucleosome

Much progress has been achieved on quantitative assessment of electrostatic interactions on the all-atom level by molecular mechanics and dynamics, as well as on the macroscopic level by models of continuum solvation. Bridging of the two representations-an area of active research-is necessary for studying integrated functions of large systems of biological importance. Following perspectives of both discrete (N-body) interaction and continuum solvation, we present a new algorithm, DiSCO (Discrete Surface Charge Optimization), for economically describing the electrostatic field predicted by Poisson-Boltzmann theory using a discrete set of Debye-Hückel charges distributed on a virtual surface enclosing the macromolecule. The procedure in DiSCO relies on the linear behavior of the Poisson-Boltzmann equation in the far zone; thus contributions from a number of molecules may be superimposed, and the electrostatic potential, or equivalently the electrostatic field, may be quickly and efficiently approximated by the summation of contributions from the set of charges. The desired accuracy of this approximation is achieved by minimizing the difference between the Poisson-Boltzmann electrostatic field and that produced by the linearized Debye-Hückel approximation using our truncated Newton optimization package. DiSCO is applied here to describe the salt-dependent electrostatic environment of the nucleosome core particle in terms of several hundred surface charges. This representation forms the basis for modeling-by dynamic simulations (or Monte Carlo)-the folding of chromatin. DiSCO can be applied more generally to many macromolecular systems whose size and complexity warrant a model resolution between the all-atom and macroscopic levels.     

A free energy analysis of nucleic acid base stacking in aqueous solution.

This paper reports a theoretical study of the free energy contributions to nucleic acid base stacking in aqueous solution. Electrostatic interactions are treated by using the finite difference Poisson-Boltzmann method and nonpolar effects are treated with explicit calculation of van der Waals interactions and/or free energy-surface area relationships. Although for some pairs of bases there is a favorable Coulombic interaction in the stacked conformation, generally the net effect of electrostatic interactions is to oppose stacking. This result is caused by the loss of favorable base-solvent electrostatic interactions, that accompany the partial removal of polar atoms from water in the stacked conformation. Nonpolar interactions, involving the hydrophobic effect and enhancement of van der Waals interactions caused by close-packing, drive stacking. The calculations qualitatively reproduce the experimental dependence of stacking free energy on purine-pyrimidine composition.     

Nonlinear poisson equation for heterogeneous media

The Poisson equation is a widely accepted model for electrostatic analysis. However, the Poisson equation is derived based on electric polarizations in a linear, isotropic, and homogeneous dielectric medium. This article introduces a nonlinear Poisson equation to take into consideration of hyperpolarization effects due to intensive charges and possible nonlinear, anisotropic, and heterogeneous media. Variational principle is utilized to derive the nonlinear Poisson model from an electrostatic energy functional. To apply the proposed nonlinear Poisson equation for the solvation analysis, we also construct a nonpolar solvation energy functional based on the nonlinear Poisson equation by using the geometric measure theory. At a fixed temperature, the proposed nonlinear Poisson theory is extensively validated by the electrostatic analysis of the Kirkwood model and a set of 20 proteins, and the solvation analysis of a set of 17 small molecules whose experimental measurements are also available for a comparison. Moreover, the nonlinear Poisson equation is further applied to the solvation analysis of 21 compounds at different temperatures. Numerical results are compared to theoretical prediction, experimental measurements, and those obtained from other theoretical methods in the literature. A good agreement between our results and experimental data as well as theoretical results suggests that the proposed nonlinear Poisson model is a potentially useful model for electrostatic analysis involving hyperpolarization effects.