DelPhi

Last updated
DelPhi
The electrostatic field lines of TRAP, trp RNA binding attenuation protein (PDB ID, 2EXS).jpg
The electrostatic field lines of TRAP---trp RNA binding attenuation protein (PDB ID: 2EXS). The blue and red colors represent the positive and negative polarities, respectively. The field lines indicate the directions of the electrostatic forces surrounding the protein.
Original author(s) Barry Honig
Developer(s) DelPhi Development Team
Operating system Linux, Mac OS X, Microsoft Windows
Website compbio.clemson.edu/delphi

DelPhi is a scientific application which calculates electrostatic potentials in and around macromolecules and the corresponding electrostatic energies. It incorporates the effects of ionic strength mediated screening by evaluating the Poisson-Boltzmann equation at a finite number of points within a three-dimensional grid box. DelPhi is commonly used in protein science to visualize variations in electrostatics along a protein or other macromolecular surface and to calculate the electrostatic components of various energies. [1]

Contents

Development

One of the main problems in modeling the electrostatic potential of biological macromolecules is that they exist in water at a given ionic strength and that they have an irregular shape. Analytical solutions of the corresponding Poisson-Boltzmann Equation (PBE) are not available for such cases and the distribution of the potential can be found only numerically. DelPhi, developed in Professor Barry Honig's lab in 1986, was the first PBE solver used by many researchers. The widespread popularity of DelPhi is due to its speed, accuracy (calculation of the electrostatic free energy is only slightly dependent on the resolution of the grid) and the ability to handle extremely high grid dimensions.

Features

Additional features such as assigning different dielectric constants to different regions of space, smooth Gaussian-based dielectric distribution function, [2] modeling geometric objects and charge distributions, and treating systems containing mixed salt solutions also attracted many researchers. In addition to the typical potential map, DelPhi can generate and output the calculated distribution of either the dielectric constant or ion concentration, providing the biomedical community with extra tools for their research. [3] [4] [5]

Pdb files are typically used as input for DelPhi calculations. Other required inputs are an atomic radii file and a charge file . [6] Binary Potential files as output from DelPhi can be viewed in most molecular viewers such as UCSF Chimera, Jmol, and VMD, and can either be mapped onto surfaces or visualized at a fixed cutoff. [7]

Versions

Delphi distribution comes as a sequential as well as parallelized codes, runs on Linux, Mac OS X and Microsoft Windows systems and the source code is available in Fortran 95 and C++ programming languages. DELPHI is also implemented into an accessible web-server. [8] DELPHI has also been utilized to build a server that predicts pKa's of biological macromolecules such as proteins, RNAs and DNAs which can be accessed via web. [9]

DelPhi v.7 is distributed in four versions:

  1. IRIX version, compiled under IRIX 6.5 Operating System, 32bits, using f77 and cc compilers.
  2. IRIX version, compiled under IRIX 6.5 Operating System, 64bits, using f77 and cc compilers.
  3. LINUX version, compiled under Red Hat 7.1, kernel 2.4.2 Operating System, using GNU gfortran compilers,
  4. PC version, compiled under Windows Operating System, using Microsoft Developer Studio C++ and Fortran compilers.

Their way of working is very similar; however, unexpected differences may appear due to different numerical precision or to the porting of the software to different architectures. For example, the elapsed time in the PC version is not calculated at present.

Each distribution contains one executable (named delphi or delphi.exe), the source codes (with corresponding makefile when needed), and some worked examples.

See also

Related Research Articles

Poissons equation Expression frequently encountered in mathematical physics, generalization of Laplaces equation.

Poisson's equation is an elliptic partial differential equation of broad utility in theoretical physics. For example, the solution to Poisson's equation is the potential field caused by a given electric charge or mass density distribution; with the potential field known, one can then calculate electrostatic or gravitational (force) field. It is a generalization of Laplace's equation, which is also frequently seen in physics. The equation is named after French mathematician and physicist Siméon Denis Poisson.

Chemistry at Harvard Macromolecular Mechanics (CHARMM) is the name of a widely used set of force fields for molecular dynamics, and the name for the molecular dynamics simulation and analysis computer software package associated with them. The CHARMM Development Project involves a worldwide network of developers working with Martin Karplus and his group at Harvard to develop and maintain the CHARMM program. Licenses for this software are available, for a fee, to people and groups working in academia.

Electrostatics Study of stationary electric charge

Electrostatics is a branch of physics that studies electric charges at rest.

Molecular mechanics Use of classical mechanics to model molecular systems

Molecular mechanics uses classical mechanics to model molecular systems. The Born–Oppenheimer approximation is assumed valid and the potential energy of all systems is calculated as a function of the nuclear coordinates using force fields. Molecular mechanics can be used to study molecule systems ranging in size and complexity from small to large biological systems or material assemblies with many thousands to millions of atoms.

Visual Molecular Dynamics

Visual Molecular Dynamics (VMD) is a molecular modelling and visualization computer program. VMD is developed as mainly a tool to view and analyze the results of molecular dynamics simulations. It also includes tools for working with volumetric data, sequence data, and arbitrary graphics objects. Molecular scenes can be exported to external rendering tools such as POV-Ray, RenderMan, Tachyon, Virtual Reality Modeling Language (VRML), and many others. Users can run their own Tcl and Python scripts within VMD as it includes embedded Tcl and Python interpreters. VMD runs on Unix, Apple Mac macOS, and Microsoft Windows. VMD is available to non-commercial users under a distribution-specific license which permits both use of the program and modification of its source code, at no charge.

Electric potential energy

Electric potential energy, is a potential energy that results from conservative Coulomb forces and is associated with the configuration of a particular set of point charges within a defined system. An object may have electric potential energy by virtue of two key elements: its own electric charge and its relative position to other electrically charged objects.

Force field (chemistry) Concept on molecular modeling

In the context of chemistry and molecular modelling, a force field is a computational method that is used to estimate the forces between atoms within molecules and also between molecules. More precisely, the force field refers to the functional form and parameter sets used to calculate the potential energy of a system of atoms or coarse-grained particles in molecular mechanics, molecular dynamics, or Monte Carlo simulations. The parameters for a chosen energy function may be derived from experiments in physics and chemistry, calculations in quantum mechanics, or both. Force fields are interatomic potentials and utilize the same concept as force fields in classical physics, with the difference that the force field parameters in chemistry describe the energy landscape, from which the acting forces on every particle are derived as a gradient of the potential energy with respect to the particle coordinates.

In computational biology, protein pKa calculations are used to estimate the pKa values of amino acids as they exist within proteins. These calculations complement the pKa values reported for amino acids in their free state, and are used frequently within the fields of molecular modeling, structural bioinformatics, and computational biology.

The Poisson–Boltzmann equation is a useful equation in many settings, whether it be to understand physiological interfaces, polymer science, electron interactions in a semiconductor, or more. It aims to describe the distribution of the electric potential in solution in the direction normal to a charged surface. This distribution is important to determine how the electrostatic interactions will affect the molecules in solution. The Poisson–Boltzmann equation is derived via mean-field assumptions. From the Poisson–Boltzmann equation many other equations have been derived with a number of different assumptions.

Implicit solvation is a method to represent solvent as a continuous medium instead of individual “explicit” solvent molecules, most often used in molecular dynamics simulations and in other applications of molecular mechanics. The method is often applied to estimate free energy of solute-solvent interactions in structural and chemical processes, such as folding or conformational transitions of proteins, DNA, RNA, and polysaccharides, association of biological macromolecules with ligands, or transport of drugs across biological membranes.

Periodic boundary conditions

Periodic boundary conditions (PBCs) are a set of boundary conditions which are often chosen for approximating a large (infinite) system by using a small part called a unit cell. PBCs are often used in computer simulations and mathematical models. The topology of two-dimensional PBC is equal to that of a world map of some video games; the geometry of the unit cell satisfies perfect two-dimensional tiling, and when an object passes through one side of the unit cell, it re-appears on the opposite side with the same velocity. In topological terms, the space made by two-dimensional PBCs can be thought of as being mapped onto a torus (compactification). The large systems approximated by PBCs consist of an infinite number of unit cells. In computer simulations, one of these is the original simulation box, and others are copies called images. During the simulation, only the properties of the original simulation box need to be recorded and propagated. The minimum-image convention is a common form of PBC particle bookkeeping in which each individual particle in the simulation interacts with the closest image of the remaining particles in the system.

The hybrid QM/MM approach is a molecular simulation method that combines the strengths of ab initio QM calculations (accuracy) and MM (speed) approaches, thus allowing for the study of chemical processes in solution and in proteins. The QM/MM approach was introduced in the 1976 paper of Warshel and Levitt. They, along with Martin Karplus, won the 2013 Nobel Prize in Chemistry for "the development of multiscale models for complex chemical systems".

Anthony Nicholls (physicist)

Anthony (Ant) Nicholls is a physicist and entrepreneur from Plympton, Plymouth, England.

The Alexander Hollaender Award in Biophysics is awarded by the U.S. National Academy of Sciences "for outstanding contributions in biophysics". Named in honor of Alexander Hollaender, it has been awarded every three years since 1998.

Biology Monte Carlo methods (BioMOCA) have been developed at the University of Illinois at Urbana-Champaign to simulate ion transport in an electrolyte environment through ion channels or nano-pores embedded in membranes. It is a 3-D particle-based Monte Carlo simulator for analyzing and studying the ion transport problem in ion channel systems or similar nanopores in wet/biological environments. The system simulated consists of a protein forming an ion channel (or an artificial nanopores like a Carbon Nano Tube, CNT), with a membrane (i.e. lipid bilayer) that separates two ion baths on either side. BioMOCA is based on two methodologies, namely the Boltzmann transport Monte Carlo (BTMC) and particle-particle-particle-mesh (P3M). The first one uses Monte Carlo method to solve the Boltzmann equation, while the later splits the electrostatic forces into short-range and long-range components.

Adsorption is the accumulation and adhesion of molecules, atoms, ions, or larger particles to a surface, but without surface penetration occurring. The adsorption of larger biomolecules such as proteins is of high physiological relevance, and as such they adsorb with different mechanisms than their molecular or atomic analogs. Some of the major driving forces behind protein adsorption include: surface energy, intermolecular forces, hydrophobicity, and ionic or electrostatic interaction. By knowing how these factors affect protein adsorption, they can then be manipulated by machining, alloying, and other engineering techniques to select for the most optimal performance in biomedical or physiological applications.

Discovery Studio is a suite of software for simulating small molecule and macromolecule systems. It is developed and distributed by Dassault Systemes BIOVIA.

Polyelectrolytes are charged polymers capable of stabilizing colloidal emulsions through electrostatic interactions. Their effectiveness can be dependent on molecular weight, pH, solvent polarity, ionic strength, and the hydrophilic-lipophilic balance (HLB). Stabilized emulsions are useful in many industrial processes, including deflocculation, drug delivery, petroleum waste treatment, and food technology.

A depletion force is an effective attractive force that arises between large colloidal particles that are suspended in a dilute solution of depletants, which are smaller solutes that are preferentially excluded from the vicinity of the large particles. One of the earliest reports of depletion forces that lead to particle coagulation is that of Bondy, who observed the separation or 'creaming' of rubber latex upon addition of polymer depletant molecules to solution. More generally, depletants can include polymers, micelles, osmolytes, ink, mud, or paint dispersed in a continuous phase.

Barry H. Honig is an American biochemist, molecular biophysicist, and computational biophysicist, who develops theoretical methods and computer software for "analyzing the structure and function of biological macromolecules."

References

  1. Rohs R, West SM, Sosinsky A, Liu P, Mann RS, Honig B (October 2009). "The role of DNA shape in protein-DNA recognition". Nature. 461 (7268): 1248–53. doi:10.1038/nature08473. PMC   2793086 . PMID   19865164.
  2. Li L, Li C, Zhang Z, Alexov E (2013). "On the Dielectric "Constant" of Proteins: Smooth Dielectric Function for Macromolecular Modeling and Its Implementation in DelPhi". Journal of Chemical Theory and Computation. 9 (4): 2126–2136. doi:10.1021/ct400065j. PMC   3622359 . PMID   23585741.
  3. "Software:DelPhi". Honig Lab. Columbia University.
  4. Rocchia W, Sridharan S, Nicholls A, Alexov E, Chiabrera A, Honig B (2002). "Rapid grid-based construction of the molecular surface and the use of induced surface charge to calculate reaction field energies: applications to the molecular systems and geometric objects". Journal of Computational Chemistry. 23 (1): 128–37. doi:10.1002/jcc.1161. PMID   11913378.
  5. Li L, Li C, Sarkar S, Zhang J, Witham S, Zhang Z, Wang L, Smith N, Petukh M, Alexov E (2012). "DelPhi: a comprehensive suite for DelPhi software and associated resources". BMC Biophysics. 5: 9. doi:10.1186/2046-1682-5-9. PMC   3463482 . PMID   22583952.
  6. "DelPhi Parameter Files".
  7. DelPhi manual
  8. "DelPhi Webserver".
  9. "DelPhiPka Webserver".