Halogen bonding in ligand-receptor systems in the framework of classical force fields

Halogen bond is an important non-covalent interaction which is receiving a growing attention in the study of protein-ligand complexes. Many drugs are halogenated molecules and it has been recently shown that many halogenated ligands establish halogen bonds with biomolecules. As the halogen bond nature is due to an anisotropy of the electrostatic potential around halogen atoms, it is not possible to use traditional force fields based on a set of atom-centred charges to study halogen bonds in biomolecules. We show that the introduction of pseudo-atoms on halogens permits us to correctly describe the anisotropy of the electrostatic potential and to perform molecular dynamics simulations on complexes of proteins with halogenated ligands that reproduce experimental values. The results are compared with crystallographic data and with hybrid quantum mechanics/molecular mechanics calculations.     

Temperature-responsive chromatography for the separation of biomolecules

Temperature-responsive chromatography for the separation of biomolecules utilizing poly(N-isopropylacrylamide) (PNIPAAm) and its copolymer-modified stationary phase is performed with an aqueous mobile phase without using organic solvent. The surface properties and function of the stationary phase are controlled by external temperature changes without changing the mobile-phase composition. This analytical system is based on nonspecific adsorption by the reversible transition of a hydrophilic-hydrophobic PNIPAAm-grafted surface. The driving force for retention is hydrophobic interaction between the solute molecules and the hydrophobized polymer chains on the stationary phase surface. The separation of the biomolecules, such as nucleotides and proteins was achieved by a dual temperature- and pH-responsive chromatography system. The electrostatic and hydrophobic interactions could be modulated simultaneously with the temperature in an aqueous mobile phase, thus the separation system would have potential applications in the separation of biomolecules. Additionally, chromatographic matrices prepared by a surface-initiated atom transfer radical polymerization (ATRP) exhibit a strong interaction with analytes, because the polymerization procedure forms a densely packed polymer, called a polymer brush, on the surfaces. The copolymer brush grafted surfaces prepared by ATRP was an effective tool for separating basic biomolecules by modulating the electrostatic and hydrophobic interactions. Applications of thermally responsive columns for the separations of biomolecules are reviewed here.     

Highly accurate biomolecular electrostatics in continuum dielectric environments

Implicit solvent models based on the Poisson-Boltzmann (PB) equation are frequently used to describe the interactions of a biomolecule with its dielectric continuum environment. A novel, highly accurate Poisson-Boltzmann solver is developed based on the matched interface and boundary (MIB) method, which rigorously enforces the continuity conditions of both the electrostatic potential and its flux at the molecular surface. The MIB based PB solver attains much better convergence rates as a function of mesh size compared to conventional finite difference and finite element based PB solvers. Consequently, highly accurate electrostatic potentials and solvation energies are obtained at coarse mesh sizes. In the context of biomolecular electrostatic calculations it is demonstrated that the MIB method generates substantially more accurate solutions of the PB equation than other established methods, thus providing a new level of reference values for such models. Initial results also indicate that the MIB method can significantly improve the quality of electrostatic surface potentials of biomolecules that are frequently used in the study of biomolecular interactions based on experimental structures.     

Monte Carlo-based linear Poisson-Boltzmann approach makes accurate salt-dependent solvation free energy predictions possible

The prediction of salt-mediated electrostatic effects with high accuracy is highly desirable since many biological processes where biomolecules such as peptides and proteins are key players can be modulated by adjusting the salt concentration of the cellular milieu. With this goal in mind, we present a novel implicit-solvent based linear Poisson-Boltzmann (PB) solver that provides very accurate nonspecific salt-dependent electrostatic properties of biomolecular systems. To solve the linear PB equation by the Monte Carlo method, we use information from the simulation of random walks in the physical space. Due to inherent properties of the statistical simulation method, we are able to account for subtle geometric features in the biomolecular model, treat continuity and outer boundary conditions and interior point charges exactly, and compute electrostatic properties at different salt concentrations in a single PB calculation. These features of the Monte Carlo-based linear PB formulation make it possible to predict the salt-dependent electrostatic properties of biomolecules with very high accuracy. To illustrate the efficiency of our approach, we compute the salt-dependent electrostatic solvation free energies of arginine-rich RNA-binding peptides and compare these Monte Carlo-based PB predictions with computational results obtained using the more mature deterministic numerical methods.     

Application of the thin-shell formulation to the numerical modeling of Stern layer in biomolecular electrostatics

In this article, the thin-shell formulation is applied to efficiently modeling the Stern layer within computational algorithms oriented toward the boundary element solution of the linearized Poisson-Boltzmann equation. The attention is focused on the calculation of the electrostatic potential in proximity to a biomolecule immersed in an electrolyte medium. Following the proposed approach, the Stern layer is made to collapse to a zero-thickness region (two-dimensional surface) with interface conditions linking the electrostatic potential over the molecular and the bulk ion accessible surfaces. Advantages lie in the limitation of divergent integral problems and in the halving of the unknown number, with a significant impact on computational time and memory requirements when modeling large biomolecules.     

An empirical potential for interactions of large molecules: Application to binding of dipeptides to DNA

An empirical potential is introduced that is suitable for calculation of the interaction enegies of biomolecules with thousands of atoms. The potential consists of electrostatic, repulsion, and dispersion energy terms. The approach used for parametrization of the potential is entirely different from that used with other existing potentials. Namely, all the terms were parametrized independently to retain their physical significance. The sum of the electrostatic and repulsion terms mimic the SCF interaction energy calculated using Huzinaga's minimal basis set MINI -1. The dispersion energy is very important and is usually predominant for the interactions of large (polar) molecules in the gas phase as well as in the liquid phase.     

Light‐Triggered Capture and Release of DNA and Proteins by Host–Guest Binding and Electrostatic Interaction

The development of an effective and general delivery method that can be applied to a large variety of structurally diverse biomolecules remains a bottleneck in modern drug therapy. Herein, we present a supramolecular system for the dynamic trapping and light‐stimulated release of both DNA and proteins. Self‐assembled ternary complexes act as nanoscale carriers, comprising vesicles of amphiphilic cyclodextrin, the target biomolecules and linker molecules with an azobenzene unit and a charged functionality. The non‐covalent linker binds to the cyclodextrin by host–guest complexation with the azobenzene. Proteins or DNA are then bound to the functionalized vesicles through multivalent electrostatic attraction. The photoresponse of the host–guest complex allows a light‐induced switch from the multivalent state that can bind the biomolecules to the low‐affinity state of the free linker, thereby providing external control over the cargo release. The major advantage of this delivery approach is the wide variety of targets that can be addressed by multivalent electrostatic interaction, which we demonstrate on four types of DNA and six different proteins.     

Reliable treatment of electrostatics in combined QM/MM simulation of macromolecules

A robust approach for dealing with electrostatic interactions for spherical boundary conditions has been implemented in the QM/MM framework. The development was based on the generalized solvent boundary potential (GSBP) method proposed by Im et al. [J. Chem. Phys. 114, 2924 (2001)], and the specific implementation was applied to the self-consistent-charge density-functional tight-binding approach as the quantum mechanics (QM) level, although extension to other QM methods is straightforward. Compared to the popular stochastic boundary-condition scheme, the new protocol offers a balanced treatment between quantum mechanics/molecular mechanics (QM/MM) and MM/MM interactions; it also includes the effect of the bulk solvent and macromolecule atoms outside of the microscopic region at the Poisson-Boltzmann level. The new method was illustrated with application to the enzyme human carbonic anhydrase II and compared to stochastic boundary-condition simulations using different electrostatic treatments. The GSBP-based QM/MM simulations were most consistent with available experimental data, while conventional stochastic boundary simulations yielded various artifacts depending on different electrostatic models. The results highlight the importance of carefully treating electrostatics in QM/MM simulations of biomolecules and suggest that the commonly used truncation schemes should be avoided in QM/MM simulations, especially in simulations that involve extensive conformational samplings. The development of the GSBP-based QM/MM protocol has opened up the exciting possibility of studying chemical events in very complex biomolecular systems in a multiscale framework.     

A semi-implicit solvent model for the simulation of peptides and proteins

We present a new model of biomolecules hydration based on macroscopic electrostatic theory, that can both describe the microscopic details of solvent-solute interactions and allow for an efficient evaluation of the electrostatic hydration free energy. This semi-implicit model considers the solvent as an ensemble of polarizable pseudoparticles whose induced dipole describe both the electronic and orientational solvent polarization. In the presented version of the model, there is no mutual dipolar interaction between the particles, and they only interact through short-ranged Lennard-Jones interactions. The model has been integrated into a molecular dynamics code, and offers the possibility to simulate efficiently the conformational evolution of biomolecules. It is able to provide estimations of the electrostatic solvation free energy within short time windows during the simulation. It has been applied to the study of two small peptides, the octaalanine and the N-terminal helix of ribonuclease A, and two proteins, the bovine pancreatic trypsin inhibitor and the B1 immunoglobin-binding domain of streptococcal protein G. Molecular dynamics simulations of these biomolecules, using a slightly modified Amber force field, provide stable and meaningful trajectories in overall agreement with experiments and all-atom simulations. Correlations with respect to Poisson-Boltzmann electrostatic solvation free energies are also presented to discuss the parameterization of the model and its consequences.     

A multipole-based water potential with implicit polarization for biomolecular simulations

A new water potential, DMIP (distributed multipoles, implicit polarization), is constructed using distributed multipoles to describe the electrostatic interactions, while accounting for polarization implicitly. In this procedure, small clusters are randomly sampled from atomistic simulations of bulk water using the AMOEBA (Ren and Ponder, J Comput Chem 2002, 23, 1497) potential. The multipole moments of the central water in each cluster are obtained from ab initio densities for each cluster, and the moments are then averaged over all clusters. Properties of bulk water calculated using DMIP compare favorably with existing data from AMOEBA simulations and experiment, with a conservative estimate of reduction in compute time of roughly 40%. The implicit force-field is also shown to work compatibly with existing polarizable multipole-based force-fields for biomolecules.     

Scanning Kelvin nanoprobe detection in materials science and biochemical analysis

The Kelvin nanoprobe is an extremely sensitive instrument capable of discerning subtle molecular interactions using vibrating electromagnetic and acoustic fields. It is based on the measurement of a fundamental material property, the work function. Modulation of this substrate parameter is caused by the adsorption or desorption of molecules, oxidation, corrosion, contamination, mechanical stress, illumination, temperature changes, electrostatic charging, surface treatment, attached dipolar structures and/or the immobilization of biomolecules. The present article explains the general principles of the method and offers an indication of the wide range of possible applications, with an emphasis on potential use in the biotechnological arena.     

Ultrafast UV spectroscopy: from a local to a global view of dynamical processes in macromolecules

The aim of this Perspective article is to cover recent developments in the application of femtosecond UV spectroscopy to understand molecular dynamics, and outlining potential future directions in this area. With several examples from recent literature the added-value of these techniques and their capability to study in real time changes in structure, dynamics and electrostatic fields of macromolecules in a site-specific fashion, as well as to uncover concerted dynamics in biomolecules, will be shown and discussed. The emerging fields of UV pulse-shaping techniques and UV optical nonlinear spectroscopies will be discussed to outline their potential to generate a novel family of coherent nonlinear spectroscopies for spectroscopic and microscopic applications.     

A nano-confined charged layer defies the principle of electrostatic interaction

The reactivity between two charged molecules and the activity of charged biomolecules are mainly governed by the principle of electrostatic interaction, i.e., like charges repel and opposite charges attract. In the present study it is shown that the principle of electrostatic interaction is violated in the nano-confined biomimetic environment. Thus a positively charged molecule shows more preference to a positively charged surface compared to a negatively charged surface.     

Accuracy assessment of semiempirical molecular electrostatic potential of proteins

 Accurate electrostatic maps of proteins are of great importance in research of protein interaction with ligands, solvent media, drugs, and other biomolecules. The large size of real-life proteins imposes severe limitations on computational methods one can use for obtaining the electrostatic map. Well-known accurate second-order M?ller–Plesset and density functional theory methods are not routinely applicable to systems larger than several hundred atoms. Conventional semiempirical tools, as less resource demanding ones, could be an attractive solution but they do not yield sufficiently accurate calculation results with reference to protein systems, as our analysis demonstrates. The present work performs a thorough analysis of the accuracy issues of the modified neglect of differential overlap type semiempirical Hamiltonians AM1 and PM3 on example of the calculation of the molecular electrostatic potential and the dipole moment of natural amino acids. Real capabilities and limitations of these methods with application to protein modeling are discussed. Received: 26 April 2002 / Accepted: 19 September 2002 / Published online: 14 February 2003     

Aqueous chromatographic system for separation of biomolecules using thermoresponsive polymer modified stationary phase

We have investigated a new method for HPLC using packing materials modified with a functional polymer, such as thermoresponsive poly(N-isopropylacrylamide) (PNIPAAm). PNIPAAm-modified silica exhibits temperature-controlled hydrophilic-hydrophobic surface property changes in aqueous systems. Temperature-responsive chromatography is performed with an aqueous mobile phase without using an organic solvent. We designed ternary copolymers of NIPAAm introduced 2-(dimethyl-amino) ethyl methacrylate (DMAEMA) as a cationic monomer and butyl methacrylate (BMA) as a hydrophobic monomer. A cationic thermoresponsive hydrogel grafted surface would produce an alterable stationary phase with both thermally regulated hydrophobicity and charge density for separation of bioactive compounds. In this study, we achieved successful separation of lysozyme without the loss of bioactivity by temperature-responsive chromatography. The electrostatic and hydrophobic interactions could be modulated simultaneously with the temperature in an aqueous mobile phase, thus the separation system would have potential applications in the separation of biomolecules.     

三种金属硫蛋白聚合物静电效应的研究

考察了三种金属硫蛋白(大鼠金属硫蛋白亚型II,兔肝金属硫蛋白亚型I和兔肝金属硫蛋白亚型II)的单体、二聚体和三聚体在pH为5.6-8.5和10.6两种缓冲条件下的静电势分布。其中大鼠金属硫蛋白亚型II的结构直接来自于晶体数据,兔肝金属硫蛋白亚型I和II的结构则通过同源蛋白模型搭建。三种金属硫蛋白的静电势通过有限差分方法求解Poisson-Boltzmann方程得到。对于三种金属硫蛋白的二聚体,pH为5.6-8.5时,单体和单体之间的静电势分布具有明显的互补性;但pH≥10.6时,这种互补性会大大削弱。对于三种金属硫蛋白的三聚体,单体和二聚体之间主要表现为静电排斥,而且pH在10.6下的静电排斥力明显强于pH为5.6-8.5时的静电排斥。     

An efficient second-order poisson–boltzmann method

Immersed interface method (IIM) is a promising high-accuracy numerical scheme for the Poisson–Boltzmann model that has been widely used to study electrostatic interactions in biomolecules. However, the IIM suffers from instability and slow convergence for typical applications. In this study, we introduced both analytical interface and surface regulation into IIM to address these issues. The analytical interface setup leads to better accuracy and its convergence closely follows a quadratic manner as predicted by theory. The surface regulation further speeds up the convergence for nontrivial biomolecules. In addition, uncertainties of the numerical energies for tested systems are also reduced by about half. More interestingly, the analytical setup significantly improves the linear solver efficiency and stability by generating more precise and better-conditioned linear systems. Finally, we implemented the bottleneck linear system solver on GPUs to further improve the efficiency of the method, so it can be widely used for practical biomolecular applications. © 2019 Wiley Periodicals, Inc.     

Biosensors based on immobilization of biomolecules by electrogenerated polymer films

The concept and potentialities of electrochemical procedures of biomolecule immobilization are described. The entrapment of biomolecules within electropolymerized films consists of the application of an appropriate potential to an electrode soaked in an aqueous solution containing monomer and biomolecules. This method of biosensor construction is compared with a two-step procedure based on the adsorption of an aqueous amphiphilic pyrrole monomer-biomolecule mixture on an electrode followed by the electropolymerization of the adsorbed monomers. Another approach is based on the electrogeneration of polymer films functionalized by specific groups allowing subsequently the attachment of biomolecules. The immobilization of biomolecules on these films by covalent binding or noncovalent interactions is described.     

A hexapeptide motif that electrostatically binds to the surface of titanium

A hexapeptide motif, RKLPDA, that recognizes the surface of titanium was identified from a peptide phage selection. Mutational analyses showed electrostatic interaction was a major factor in the binding, paving a new way to modify the surface of Ti with active biomolecules.