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 Adaptive Poisson-Boltzmann Solver (APBS) software was developed to solve the equations of continuum electrostatics for large biomolecular assemblages that has provided impact in the study of a broad range of chemical, biological, and biomedical applications. APBS addresses three key technology challenges for understanding solvation and electrostatics in biomedical applications: accurate and efficient models for biomolecular solvation and electrostatics, robust and scalable software for applying those theories to biomolecular systems, and mechanisms for sharing and analyzing biomolecular electrostatics data in the scientific community. To address new research applications and advancing computational capabilities, we have continually updated APBS and its suite of accompanying software since its release in 2001. In this manuscript, we discuss the models and capabilities that have recently been implemented within the APBS software package including: a Poisson-Boltzmann analytical and a semi-analytical solver, an optimized boundary element solver, a geometry-based geometric flow solvation model, a graph theory based algorithm for determining p$K_a$ values, and an improved web-based visualization tool for viewing electrostatics.

Improvements to the APBS biomolecular solvation software suite

Molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) and molecular mechanics generalized Born surface area (MM/GBSA) are arguably very popular methods for binding free energy prediction since they are more accurate than most scoring functions of molecular docking and less computationally demanding than alchemical free energy methods. MM/PBSA and MM/GBSA have been widely used in biomolecular studies such as protein folding, protein-ligand binding, protein-protein interaction, etc. In this review, methods to adjust the polar solvation energy and to improve the performance of MM/PBSA and MM/GBSA calculations are reviewed and discussed. The latest applications of MM/GBSA and MM/PBSA in drug design are also presented. This review intends to provide readers with guidance for practically applying MM/PBSA and MM/GBSA in drug design and related research fields.

End-Point Binding Free Energy Calculation with MM/PBSA and MM/GBSA: Strategies and Applications in Drug Design

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the current coronavirus disease 2019 (COVID-19) pandemic. A major virulence factor of SARS-CoVs is the nonstructural protein 1 (Nsp1), which suppresses host gene expression by ribosome association. Here, we show that Nsp1 from SARS-CoV-2 binds to the 40S ribosomal subunit, resulting in shutdown of messenger RNA (mRNA) translation both in vitro and in cells. Structural analysis by cryo-electron microscopy of in vitro-reconstituted Nsp1-40S and various native Nsp1-40S and -80S complexes revealed that the Nsp1 C terminus binds to and obstructs the mRNA entry tunnel. Thereby, Nsp1 effectively blocks retinoic acid-inducible gene I-dependent innate immune responses that would otherwise facilitate clearance of the infection. Thus, the structural characterization of the inhibitory mechanism of Nsp1 may aid structure-based drug design against SARS-CoV-2.

https://science.sciencemag.org/content/sci/early/2020/07/16/science.abc8665.full.pdf

Structural basis for translational shutdown and immune evasion by the Nsp1 protein of SARS-CoV-2.

The Adaptive Poisson-Boltzmann Solver (APBS) software was developed to solve the equations of continuum electrostatics for large biomolecular assemblages that have provided impact in the study of a broad range of chemical, biological, and biomedical applications. APBS addresses the three key technology challenges for understanding solvation and electrostatics in biomedical applications: accurate and efficient models for biomolecular solvation and electrostatics, robust and scalable software for applying those theories to biomolecular systems, and mechanisms for sharing and analyzing biomolecular electrostatics data in the scientific community. To address new research applications and advancing computational capabilities, we have continually updated APBS and its suite of accompanying software since its release in 2001. In this article, we discuss the models and capabilities that have recently been implemented within the APBS software package including a Poisson-Boltzmann analytical and a semi-analytical solver, an optimized boundary element solver, a geometry-based geometric flow solvation model, a graph theory-based algorithm for determining pKa values, and an improved web-based visualization tool for viewing electrostatics.

https://onlinelibrary.wiley.com/doi/pdf/10.1002/pro.3280

Improvements to the APBS biomolecular solvation software suite.

Background: Accurate modeling of electrostatic potential and corresponding energies becomes increasingly important for understanding properties of biological macromolecules and their complexes. However, this is not an easy task due to the irregular shape of biological entities and the presence of water and mobile ions. Results: Here we report a comprehensive suite for the well-known Poisson-Boltzmann solver, DelPhi, enriched with additional features to facilitate DelPhi usage. The suite allows for easy download of both DelPhi executable files and source code along with a makefile for local installations. The users can obtain the DelPhi manual and parameter files required for the corresponding investigation. Non-experienced researchers can download examples containing all necessary data to carry out DelPhi runs on a set of selected examples illustrating various DelPhi features and demonstrating DelPhi’s accuracy against analytical solutions. Conclusions: DelPhi suite offers not only the DelPhi executable and sources files, examples and parameter files, but also provides links to third party developed resources either utilizing DelPhi or providing plugins for DelPhi. In addition, the users and developers are offered a forum to share ideas, resolve issues, report bugs and seek help with respect to the DelPhi package. The resource is available free of charge for academic users from URL: http:// compbio.clemson.edu/DelPhi.php.

/pdf/delphi-a-comprehensive-suite-for-delphi-software-and-2facqdjp72.pdf

DelPhi: a comprehensive suite for DelPhi software and associated resources

The role of electrostatics in protein–protein interactions and binding is reviewed in this paper. A brief outline of the computational modeling, in the framework of continuum electrostatics, is presented and the basic electrostatic effects occurring upon the formation of the complex are discussed. The effect of the salt concentration and pH of the water phase on protein–protein binding free energy is demonstrated which indicates that the increase of the salt concentration tends to weaken the binding, an observation that is attributed to the optimization of the charge–charge interactions across the interface. It is pointed out that the pH-optimum (pH of optimal binding affinity) varies among the protein–protein complexes, and perhaps is a result of their adaptation to particular subcellular compartments. The similarities and differences between hetero- and homo-complexes are outlined and discussed with respect to the binding mode and charge complementarity.

/pdf/on-the-role-of-electrostatics-in-protein-protein-3j4m2z7un3.pdf

On the role of electrostatics in protein–protein interactions

Chloride intracellular channel 2 (CLIC2) protein is a member of the glutathione transferase class of proteins. Its' only known function is the regulation of ryanodine receptor (RyR) intracellular Ca(2+) release channels. These RyR proteins play a major role in the regulation of Ca(2+) signaling in many cells. Utilizing exome capture and deep sequencing of genes on the X-chromosome, we have identified a mutation in CLIC2 (c.303C>G, p.H101Q) which is associated with X-linked intellectual disability (ID), atrial fibrillation, cardiomegaly, congestive heart failure (CHF), some somatic features and seizures. Functional studies of the H101Q variant indicated that it stimulated rather than inhibited the action of RyR channels, with channels remaining open for longer times and potentially amplifying Ca(2+) signals dependent on RyR channel activity. The overly active RyRs in cardiac and skeletal muscle cells and neuronal cells would result in abnormal cardiac function and trigger post-synaptic pathways and neurotransmitter release. The presence of both cardiomegaly and CHF in the two affected males and atrial fibrillation in one are consistent with abnormal RyR2 channel function. Since the dysfunction of RyR2 channels in the brain via 'leaky mutations' can result in mild developmental delay and seizures, our data also suggest a vital role for the CLIC2 protein in maintaining normal cognitive function via its interaction with RyRs in the brain. Therefore, our patients appear to suffer from a new channelopathy comprised of ID, seizures and cardiac problems because of enhanced Ca(2+) release through RyRs in neuronal cells and cardiac muscle cells.

/pdf/an-x-linked-channelopathy-with-cardiomegaly-due-to-a-clic2-4gtr4qz0fr.pdf

An X-linked channelopathy with cardiomegaly due to a CLIC2 mutation enhancing ryanodine receptor channel activity

Large-scale next generation resequencing of X chromosome genes identified a missense mutation in the CLIC2 gene on Xq28 in a male with X-linked intellectual disability (XLID) and not found in healthy individuals. At the same time, numerous nsSNPs (nonsynonomous SNP) have been reported in the CLIC2 gene in healthy individuals indicating that the CLIC2 protein can tolerate amino acid substitutions and be fully functional. To test the possibility that p.H101Q is a disease-causing mutation, we performed in silico simulations to calculate the effects of the p.H101Q mutation on CLIC2 stability, dynamics, and ionization states while comparing the effects obtained for presumably harmless nsSNPs. It was found that p.H101Q, in contrast with other nsSNPs, (a) lessens the flexibility of the joint loop which is important for the normal function of CLIC2, (b) makes the overall 3D structure of CLIC2 more stable and thus reduces the possibility of the large conformational change expected to occur when CLIC2 moves from a soluble to membrane form, and (c) removes the positively charged residue, H101, which may be important for the membrane association of CLIC2. The results of in silico modeling, in conjunction with the polymorphism analysis, suggest that p.H101Q may be a disease-causing mutation, the first one suggested in the CLIC family.

/pdf/a-missense-mutation-in-clic2-associated-with-intellectual-4k92a1nd8r.pdf

A missense mutation in CLIC2 associated with intellectual disability is predicted by in silico modeling to affect protein stability and dynamics

Here we report a web server, the DelPhi web server, which utilizes DelPhi program to calculate electrostatic energies and the corresponding electrostatic potential and ionic distributions, and dielectric map. The server provides extra services to fix structural defects, as missing atoms in the structural file and allows for generation of missing hydrogen atoms. The hydrogen placement and the corresponding DelPhi calculations can be done with user selected force field parameters being either Charmm22, Amber98 or OPLS. Upon completion of the calculations, the user is given option to download fixed and protonated structural file, together with the parameter and Delphi output files for further analysis. Utilizing Jmol viewer, the user can see the corresponding structural file, to manipulate it and to change the presentation. In addition, if the potential map is requested to be calculated, the potential can be mapped onto the molecule surface. The DelPhi web server is available from http://compbio.clemson.edu/delphi_webserver.

/pdf/delphi-web-server-a-comprehensive-online-suite-for-220lvib4n3.pdf

Shawn Witham

Papers

DelPhi: a comprehensive suite for DelPhi software and associated resources

On the role of electrostatics in protein–protein interactions

An X-linked channelopathy with cardiomegaly due to a CLIC2 mutation enhancing ryanodine receptor channel activity

A missense mutation in CLIC2 associated with intellectual disability is predicted by in silico modeling to affect protein stability and dynamics

DelPhi Web Server: A comprehensive online suite for electrostatic calculations of biological macromolecules and their complexes