Continuum solvation models, such as Poisson–Boltzmann and Generalized Born methods, have become increasingly popular tools for investigating the influence of electrostatics on biomolecular structure, energetics and dynamics. However, the use of such methods requires accurate and complete structural data as well as force field parameters such as atomic charges and radii. Unfortunately, the limiting step in continuum electrostatics calculations is often the addition of missing atomic coordinates to molecular structures from the Protein Data Bank and the assignment of parameters to biomolecular structures. To address this problem, we have developed the PDB2PQR web service (http://agave.wustl.edu/pdb2pqr/). This server automates many of the common tasks of preparing structures for continuum electrostatics calculations, including adding a limited number of missing heavy atoms to biomolecular structures, estimating titration states and protonating biomolecules in a manner consistent with favorable hydrogen bonding, assigning charge and radius parameters from a variety of force fields, and finally generating ‘PQR’ output compatible with several popular computational biology packages. This service is intended to facilitate the setup and execution of electrostatics calculations for both experts and non-experts and thereby broaden the accessibility to the biological community of continuum electrostatics analyses of biomolecular systems.

/pdf/pdb2pqr-an-automated-pipeline-for-the-setup-of-poisson-53luup701x.pdf

PDB2PQR: an automated pipeline for the setup of Poisson–Boltzmann electrostatics calculations

Molecular mechanics Poisson–Boltzmann surface area (MM-PBSA), a method to estimate interaction free energies, has been increasingly used in the study of biomolecular interactions. Recently, this method has also been applied as a scoring function in computational drug design. Here a new tool g_mmpbsa, which implements the MM-PBSA approach using subroutines written in-house or sourced from the GROMACS and APBS packages is described. g_mmpbsa was developed as part of the Open Source Drug Discovery (OSDD) consortium. Its aim is to integrate high-throughput molecular dynamics (MD) simulations with binding energy calculations. The tool provides options to select alternative atomic radii and different nonpolar solvation models including models based on the solvent accessible surface area (SASA), solvent accessible volume (SAV), and a model which contains both repulsive (SASA-SAV) and attractive components (described using a Weeks–Chandler–Andersen like integral method). We showcase the effectiveness of the tool ...

g_mmpbsa--a GROMACS tool for high-throughput MM-PBSA calculations.

The ClC chloride channels catalyse the selective flow of Cl- ions across cell membranes, thereby regulating electrical excitation in skeletal muscle and the flow of salt and water across epithelial barriers. Genetic defects in ClC Cl- channels underlie several familial muscle and kidney diseases. Here we present the X-ray structures of two prokaryotic ClC Cl- channels from Salmonella enterica serovar typhimurium and Escherichia coli at 3.0 and 3.5 A, respectively. Both structures reveal two identical pores, each pore being formed by a separate subunit contained within a homodimeric membrane protein. Individual subunits are composed of two roughly repeated halves that span the membrane with opposite orientations. This antiparallel architecture defines a selectivity filter in which a Cl- ion is stabilized by electrostatic interactions with alpha-helix dipoles and by chemical coordination with nitrogen atoms and hydroxyl groups. These findings provide a structural basis for further understanding the function of ClC Cl- channels, and establish the physical and chemical basis of their anion selectivity.

/pdf/x-ray-structure-of-a-clc-chloride-channel-at-3-0-a-reveals-1wk74mj1p6.pdf

X-ray structure of a ClC chloride channel at 3.0 A reveals the molecular basis of anion selectivity.

https://www.cell.com/trends/biochemical-sciences/pdf/S0968-0004(03)00273-1.pdf

TPR proteins: the versatile helix.

The recognition of specific DNA sequences by proteins is thought to depend on two types of mechanism: one that involves the formation of hydrogen bonds with specific bases, primarily in the major groove, and one involving sequence-dependent deformations of the DNA helix. By comprehensively analysing the three-dimensional structures of protein-DNA complexes, here we show that the binding of arginine residues to narrow minor grooves is a widely used mode for protein-DNA recognition. This readout mechanism exploits the phenomenon that narrow minor grooves strongly enhance the negative electrostatic potential of the DNA. The nucleosome core particle offers a prominent example of this effect. Minor-groove narrowing is often associated with the presence of A-tracts, AT-rich sequences that exclude the flexible TpA step. These findings indicate that the ability to detect local variations in DNA shape and electrostatic potential is a general mechanism that enables proteins to use information in the minor groove, which otherwise offers few opportunities for the formation of base-specific hydrogen bonds, to achieve DNA-binding specificity.

/pdf/the-role-of-dna-shape-in-protein-dna-recognition-3xfnu159dk.pdf

The role of DNA shape in protein–DNA recognition

This article describes a number of algorithms that are designed to improve both the efficiency and accuracy of finite difference solutions to the Poisson-Boltzmann equation (the FDPB method) and to extend its range of application. The algorithms are incorporated in the DelPhi program. The first algorithm involves an efficient and accurate semianalytical method to map the molecular surface of a molecule onto a three-dimensional lattice. This method constitutes a significant improvement over existing methods in terms of its combination of speed and accuracy. The DelPhi program has also been expanded to allow the definition of geometrical objects such as spheres, cylinders, cones, and parallelepipeds, which can be used to describe a system that may also include a standard atomic level depiction of molecules. Each object can have a different dielectric constant and a different surface or volume charge distribution. The improved definition of the surface leads to increased precision in the numerical solutions of the PB equation that are obtained. A further improvement in the precision of solvation energy calculations is obtained from a procedure that calculates induced surface charges from the FDPB solutions and then uses these charges in the calculation of reaction field energies. The program allows for finite difference grids of large dimension; currently a maximum of 571 3 can be used on molecules containing several thousand atoms and charges. As described elsewhere, DelPhi can also treat mixed salt systems containing mono- and divalent ions and provide electrostatic free energies as defined by the nonlinear PB equation.

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.

Alessandro Chiabrera

Papers

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.