Charge profile analysis reveals that activation of pro-apoptotic regulators Bax and Bak relies on charge transfer mediated allosteric regulation

NAChRDB: A Web Resource of Structure-Function Annotations to Unravel the Allostery of Nicotinic Acetylcholine Receptors

Nicotinic acetylcholine receptors (nAChRs) comprise a large and ancient family of allosteric ion channels mediating synaptic transmission. The vast knowledge about nAChRs has become difficult to navigate. NAChRDB is a web-accessible resource of curated residue-level functional annotations of neuromuscular nAChRs. Interactive three-dimensional (3D) visualization and sequence alignment give further context to this rich and growing collection of experimental observations and computational predictions. NAChRDB is freely available at https://crocodile.ncbr.muni.cz/Apps/NAChRDB/, with interactive tutorials and regular updates to the content and web interface. No installation or user registration is required. NAChRDB is accessible through any modern internet browser on desktops and mobile devices. By providing immediate and systematic access to practical knowledge gained through decades of research, NAChRDB represents a powerful educational tool and helps guide discovery by revealing gaps in current knowledge and aiding the interpretation of results of molecular and structural biology experiments or computational studies.

Measuring how two proteins affect each other's net charge in a crowded environment

Theory predicts that the net charge (Z) of a protein can be altered by the net charge of a neighboring protein as the two approach one another below the Debye length. This type of charge regulation suggests that a protein's charge and perhaps function might be affected by neighboring proteins without direct binding. Charge regulation during protein crowding has never been directly measured due to analytical challenges. Here, we show that lysine specific protein crosslinkers (NHS ester-Staudinger pairs) can be used to mimic crowding by linking two non-interacting proteins at a maximal distance of ~7.9 Å. The net charge of the regioisomeric dimers and preceding monomers can then be determined with lysine-acyl "protein charge ladders" and capillary electrophoresis. As a proof of concept, we covalently linked myoglobin (Z_monomer  = -0.43 ± 0.01) and α-lactalbumin (Z_monomer  = -4.63 ± 0.05). Amide hydrogen/deuterium exchange and circular dichroism spectroscopy demonstrated that crosslinking did not significantly alter the structure of either protein or result in direct binding (thus mimicking crowding). Ultimately, capillary electrophoretic analysis of the dimeric charge ladder detected a change in charge of ΔZ = -0.04 ± 0.09 upon crowding by this pair (Z_dimer  = -5.10 ± 0.07). These small values of ΔZ are not necessarily general to protein crowding (qualitatively or quantitatively) but will vary per protein size, charge, and solvent conditions.

Atomic Charge Calculator II: web-based tool for the calculation of partial atomic charges

Partial atomic charges serve as a simple model for the electrostatic distribution of a molecule that drives its interactions with its surroundings. Since partial atomic charges are frequently used in computational chemistry, chemoinformatics and bioinformatics, many computational approaches for calculating them have been introduced. The most applicable are fast and reasonably accurate empirical charge calculation approaches. Here, we introduce Atomic Charge Calculator II (ACC II), a web application that enables the calculation of partial atomic charges via all the main empirical approaches and for all types of molecules. ACC II implements 17 empirical charge calculation methods, including the highly cited (QEq, EEM), the recently published (EQeq, EQeq+C), and the old but still often used (PEOE). ACC II enables the fast calculation of charges even for large macromolecular structures. The web server also offers charge visualization, courtesy of the powerful LiteMol viewer. The calculation setup of ACC II is very straightforward and enables the quick calculation of high-quality partial charges. The application is available at https://acc2.ncbr.muni.cz.

Transient Unfolding and Long-Range Interactions in Viral BCL2 M11 Enable Binding to the BECN1 BH3 Domain

Viral BCL2 proteins (vBCL2s) help to sustain chronic infection of host proteins to inhibit apoptosis and autophagy. However, details of conformational changes in vBCL2s that enable binding to BH3Ds remain unknown. Using all-atom, multiple microsecond-long molecular dynamic simulations (totaling 17 μs) of the murine γ-herpesvirus 68 vBCL2 (M11), and statistical inference techniques, we show that regions of M11 transiently unfold and refold upon binding of the BH3D. Further, we show that this partial unfolding/refolding within M11 is mediated by a network of hydrophobic interactions, which includes residues that are 10 Å away from the BH3D binding cleft. We experimentally validate the role of these hydrophobic interactions by quantifying the impact of mutating these residues on binding to the Beclin1/BECN1 BH3D, demonstrating that these mutations adversely affect both protein stability and binding. To our knowledge, this is the first study detailing the binding-associated conformational changes and presence of long-range interactions within vBCL2s.

Iterative Solvers for Empirical Partial Atomic Charges: Breaking the Curse of Cubic Numerical Complexity

Rational drug design involves a vast amount of computations to get thermodynamically reliable results and often relies on atomic charges as a means to model electrostatic interactions within the system. Computational inefficiency often hampers the development of new and wider dissemination of the known methods; thus, any source to speed up the calculations without a sacrifice in quality is warranted. At the heart of many empirical methods of calculating atomic charges is the solution of a system of linear algebraic equations (SLAE). The classical method of solving SLAE-the Gauss method-has in general case a cubic computational complexity. It is shown that the use of iterative methods for solving SLAE, characteristic to typical empirical atomic charge calculation methods, makes it possible to significantly reduce the amount of calculations and to obtain a computational complexity approaching a quadratic one. Despite the fact that this phenomenon is well-known in numerical methods, iterative solvers surprisingly do not seem to have been systematically applied to calculation of atomic charges via empirical schemes. Another finding is the relative values of the matrix elements, determined by the physical grounds of the interactions within the empirical system, generally lead to SLAE's with well-defined matrices, suited to use with iterative solvers to fasten computation compared to using the noniterative solvers. This finding broadens the applicability range of atomic charges obtained with empirical methods for such cases as, e.g., account of polarizability via "on-the-fly" recalculation of charges in changing surroundings within the force fields in molecular dynamics settings.

Mitochondrial VDAC2 and cell homeostasis: highlighting hidden structural features and unique functionalities

Voltage-dependent anion channels (VDACs) are the gateway to mitochondrial processes, interlinking the cytosolic and mitochondrial compartments. The mitochondrion acts as a storehouse for cytochrome c, the effector of apoptosis, and hence VDACs become intricately involved in the apoptotic pathway. Isoform 1 of VDAC is abundant in the outer mitochondrial membrane of many cell types, while isoform 2 is the preferred channel in specialized cells including brain and some cancer cells. The primary role of VDACs is metabolite flux. The pro- and anti-apoptotic role of VDAC1 and VDAC2, respectively, are secondary, and are influenced by external factors and interacting proteins. Herein, we focus on the less-studied VDAC2, and shed light on its unique functions and features. VDAC2, along with sharing many of its functions with VDAC1, such as metabolite and Ca²⁺ transport, also has many delineating functions. VDAC2 is closely engaged in the gametogenesis and steroidogenesis pathways and in protection from oxidative stress as well as in neurodegenerative diseases like Alzheimer's and epilepsy. A closer examination of the functional pathways of VDACs indicates that the unique functions of VDAC2 are a result of the different interactome of this isoform. We couple functional differences to the structural and biophysical evidence obtained for the VDACs, and present a testament of why the two VDAC isoforms with >90% sequence similarity, are functionally diverse. Based on these differences, we suggest that the VDAC isoforms now be considered as paralogs. An in-depth understanding of VDAC2 will help us to design better biomolecule targets for cancer and neurodegenerative diseases.

AtomicChargeCalculator: interactive web-based calculation of atomic charges in large biomolecular complexes and drug-like molecules

Background

Partial atomic charges are a well-established concept, useful in understanding and modeling the chemical behavior of molecules, from simple compounds, to large biomolecular complexes with many reactive sites.

Results

This paper introduces AtomicChargeCalculator (ACC), a web-based application for the calculation and analysis of atomic charges which respond to changes in molecular conformation and chemical environment. ACC relies on an empirical method to rapidly compute atomic charges with accuracy comparable to quantum mechanical approaches. Due to its efficient implementation, ACC can handle any type of molecular system, regardless of size and chemical complexity, from drug-like molecules to biomacromolecular complexes with hundreds of thousands of atoms. ACC writes out atomic charges into common molecular structure files, and offers interactive facilities for statistical analysis and comparison of the results, in both tabular and graphical form.

Conclusions

Due to high customizability and speed, easy streamlining and the unified platform for calculation and analysis, ACC caters to all fields of life sciences, from drug design to nanocarriers. ACC is freely available via the Internet at http://ncbr.muni.cz/ACC.

Electronic supplementary material

The online version of this article (doi:10.1186/s13321-015-0099-x) contains supplementary material, which is available to authorized users.

The indirect activation model of mitochondrial outer membrane permeabilisation (MOMP) initiation requires a trade-off between robustness in the absence of and sensitivity in the presence of stress

Mitochondrial outer membrane permeabilisation (MOMP) is a hallmark of apoptosis and is controlled by antagonising members of the BCL 2 protein family. However, whether the effector proteins for MOMP, BAX and BAK, require a separate activation step (direct activation model) or if the sole inhibition of otherwise active effectors by anti-apoptotic BCL 2 proteins suffices (indirect activation model) is controversial. To address this question, we recently identified that, in most cells, effector proteins were more abundant than their inhibitors. We now employ systems modelling to elucidate fundamental consequences of this imbalance for the indirect activation model by integrating its key features into the recently validated implementation of the direct activation model. First, by modelling constitutively active BAK, several cell lines were predicted to be unstable against MOMP even in the absence of stress. This instability could be rectified by assuming BAK pores to be subject to degradation, which however demanded BAK to underlie a high and biologically unreasonable protein turnover. Stability against MOMP in the absence of stress was similarly reconstituted by assuming BAK oligomerisation to be weak, but predicted some cells to be MOMP resistant even under high stress. Assuming also BAX to be constitutively active required highly effective BAX re-translocation to the cytosol along with weak BAX oligomerisation, which again rendered some cells insensitive to induce MOMP in the presence of stress. Our results suggest that an activation step of the two effectors is required to allow stability against MOMP in the absence of stress and execution of MOMP under stress.

Accelerated Molecular Dynamics Simulations with the AMOEBA Polarizable Force Field on Graphics Processing Units

The accelerated molecular dynamics (aMD) method has recently been shown to enhance the sampling of biomolecules in molecular dynamics (MD) simulations, often by several orders of magnitude. Here, we describe an implementation of the aMD method for the OpenMM application layer that takes full advantage of graphics processing units (GPUs) computing. The aMD method is shown to work in combination with the AMOEBA polarizable force field (AMOEBA-aMD), allowing the simulation of long time-scale events with a polarizable force field. Benchmarks are provided to show that the AMOEBA-aMD method is efficiently implemented and produces accurate results in its standard parametrization. For the BPTI protein, we demonstrate that the protein structure described with AMOEBA remains stable even on the extended time scales accessed at high levels of accelerations. For the DNA repair metalloenzyme endonuclease IV, we show that the use of the AMOEBA force field is a significant improvement over fixed charged models for describing the enzyme active-site. The new AMOEBA-aMD method is publicly available (http://wiki.simtk.org/openmm/VirtualRepository) and promises to be interesting for studying complex systems that can benefit from both the use of a polarizable force field and enhanced sampling.