Bartłomiej Zaborowski

A unified coarse-grained model of biological macromolecules based on mean-field multipole–multipole interactions

AbstractA unified coarse-grained model of three major classes of biological molecules—proteins, nucleic acids, and polysaccharides—has been developed. It is based on the observations that the repeated units of biopolymers (peptide groups, nucleic acid bases, sugar rings) are highly polar and their charge distributions can be represented crudely as point multipoles. The model is an extension of the united residue (UNRES) coarse-grained model of proteins developed previously in our laboratory. The respective force fields are defined as the potentials of mean force of biomacromolecules immersed in water, where all degrees of freedom not considered in the model have been averaged out. Reducing the representation to one center per polar interaction site leads to the representation of average site–site interactions as mean-field dipole–dipole interactions. Further expansion of the potentials of mean force of biopolymer chains into Kubo’s cluster-cumulant series leads to the appearance of mean-field dipole–dipole interactions, averaged in the context of local interactions within a biopolymer unit. These mean-field interactions account for the formation of regular structures encountered in biomacromolecules, e.g., α-helices and β-sheets in proteins, double helices in nucleic acids, and helicoidally packed structures in polysaccharides, which enables us to use a greatly reduced number of interacting sites without sacrificing the ability to reproduce the correct architecture. This reduction results in an extension of the simulation timescale by more than four orders of magnitude compared to the all-atom representation. Examples of the performance of the model are presented. FigureComponents of the Unified Coarse Grained Model (UCGM) of biological macromolecules

Analysis of Fluorescence Quenching of Coumarin Derivatives by 4-Hydroxy-TEMPO in Aqueous Solution

The fluorescence quenching of different coumarin derivatives (7-hydroxy-4-methylcoumarin, 5,7-dimethoxycoumarin, 7-amino-4-methyl-3-coumarinylacetic acid, 7-ethoxy-4-methylcoumarin, 7-methoxycoumarin, 7-hydroxycoumarin, 7-hydroxy-4-methyl-3-coumarinylacetic acid and 7-amino-4-methylcoumarin) by 4-hydroxy-TEMPO in aqueous solutions at the room temperature was studied with the use of UV–Vis absorption spectroscopy as well as a steady-state and time-resolved fluorescence spectroscopy. In order to understand the mechanism of quenching the absorption and fluorescence emission spectra of all coumarins along with fluorescence decays were recorded under the action of 4-hydroxy-TEMPO. The Stern-Volmer plots (both from time-averaged and time-resolved measurements) displayed no positive (upward) deviation from a linearity. The fluorescence quenching mechanism was found to be entirely dynamic, what was additionally confirmed by the registration of Stern-Volmer plots at different temperatures. The Stern-Volmer quenching constants and bimolecular quenching rate constants were obtained for all coumarins studied at the room temperature. The findings demonstrate the possibility of developing an analytical method for the quantitative determination of the free radicals’ scavenger, 4-hydroxy-TEMPO.

A Maximum-Likelihood Approach to Force-Field Calibration.

A new approach to the calibration of the force fields is proposed, in which the force-field parameters are obtained by maximum-likelihood fitting of the calculated conformational ensembles to the experimental ensembles of training system(s). The maximum-likelihood function is composed of logarithms of the Boltzmann probabilities of the experimental conformations, calculated with the current energy function. Because the theoretical distribution is given in the form of the simulated conformations only, the contributions from all of the simulated conformations, with Gaussian weights in the distances from a given experimental conformation, are added to give the contribution to the target function from this conformation. In contrast to earlier methods for force-field calibration, the approach does not suffer from the arbitrariness of dividing the decoy set into native-like and non-native structures; however, if such a division is made instead of using Gaussian weights, application of the maximum-likelihood method results in the well-known energy-gap maximization. The computational procedure consists of cycles of decoy generation and maximum-likelihood-function optimization, which are iterated until convergence is reached. The method was tested with Gaussian distributions and then applied to the physics-based coarse-grained UNRES force field for proteins. The NMR structures of the tryptophan cage, a small α-helical protein, determined at three temperatures (T = 280, 305, and 313 K) by Hałabis et al. ( J. Phys. Chem. B 2012 , 116 , 6898 - 6907 ), were used. Multiplexed replica-exchange molecular dynamics was used to generate the decoys. The iterative procedure exhibited steady convergence. Three variants of optimization were tried: optimization of the energy-term weights alone and use of the experimental ensemble of the folded protein only at T = 280 K (run 1); optimization of the energy-term weights and use of experimental ensembles at all three temperatures (run 2); and optimization of the energy-term weights and the coefficients of the torsional and multibody energy terms and use of experimental ensembles at all three temperatures (run 3). The force fields were subsequently tested with a set of 14 α-helical and two α + β proteins. Optimization run 1 resulted in better agreement with the experimental ensemble at T = 280 K compared with optimization run 2 and in comparable performance on the test set but poorer agreement of the calculated folding temperature with the experimental folding temperature. Optimization run 3 resulted in the best fit of the calculated ensembles to the experimental ones for the tryptophan cage but in much poorer performance on the training set, suggesting that use of a small α-helical protein for extensive force-field calibration resulted in overfitting of the data for this protein at the expense of transferability. The optimized force field resulting from run 2 was found to fold 13 of the 14 tested α-helical proteins and one small α + β protein with the correct topologies; the average structures of 10 of them were predicted with accuracies of about 5 Å C(α) root-mean-square deviation or better. Test simulations with an additional set of 12 α-helical proteins demonstrated that this force field performed better on α-helical proteins than the previous parametrizations of UNRES. The proposed approach is applicable to any problem of maximum-likelihood parameter estimation when the contributions to the maximum-likelihood function cannot be evaluated at the experimental points and the dimension of the configurational space is too high to construct histograms of the experimental distributions.

Fluorescence quenching of 7-amino-4-methylcoumarin by different TEMPO derivatives.

The fluorescence quenching of 7-amino-4-methylcoumarin by different TEMPO derivatives was studied in aqueous solutions with the use of steady-state, time-resolved fluorescence spectroscopy as well as UV-VIS absorption spectroscopy methods. In order to distinguish each TEMPO derivative from the others and to understand the mechanism of quenching, the absorption and fluorescence emission spectra as well as decays of the fluorescence of 7-amino-4-methylcoumarin were registered as a function of each TEMPO derivative concentration. There were no deviations from a linearity in the Stern-Volmer plots (determined from both, steady-state and time-resolved measurements). The fluorescence quenching mechanism was found to be entirely collisional, what was additionally confirmed by the registration of Stern-Volmer plots at 5 temperatures ranging from 15 to 55°C. Based on theoretical calculations of molecular radii and ionization potentials of all TEMPO derivatives the mechanism of electron transfer was rejected. The fluorescence quenching which was being studied seems to be diffusion-limited and caused by the increase of non-radiative processes, such as an internal conversion and an intersystem crossing. The Stern-Volmer quenching constants and bimolecular quenching constants were determined at the room temperature for all TEMPO derivatives studied. Among all TEMPO derivatives studied TEMPO-4-amino-4-carboxylic acid (TOAC) was found to be the most effective quencher of 7-amino-4-methylcoumarin fluorescence (kq for TOAC was approximately 1.5 higher than kq for other TEMPO compounds studied). The findings demonstrate the possibility of developing an analytical method for the quantitative determination of TOAC, which incorporation into membrane proteins may provide a direct detection of peptide backbone dynamics.

An analysis and evaluation of the WeFold collaborative for protein structure prediction and its pipelines in CASP11 and CASP12

Every two years groups worldwide participate in the Critical Assessment of Protein Structure Prediction (CASP) experiment to blindly test the strengths and weaknesses of their computational methods. CASP has significantly advanced the field but many hurdles still remain, which may require new ideas and collaborations. In 2012 a web-based effort called WeFold, was initiated to promote collaboration within the CASP community and attract researchers from other fields to contribute new ideas to CASP. Members of the WeFold coopetition (cooperation and competition) participated in CASP as individual teams, but also shared components of their methods to create hybrid pipelines and actively contributed to this effort. We assert that the scale and diversity of integrative prediction pipelines could not have been achieved by any individual lab or even by any collaboration among a few partners. The models contributed by the participating groups and generated by the pipelines are publicly available at the WeFold website providing a wealth of data that remains to be tapped. Here, we analyze the results of the 2014 and 2016 pipelines showing improvements according to the CASP assessment as well as areas that require further adjustments and research.

Use of the UNRES force field in template-assisted prediction of protein structures and the refinement of server models: Test with CASP12 targets

Knowledge-based methods are, at present, the most effective ones for the prediction of protein structures; however, their results heavily depend on the similarity of a target sequence to those of proteins with known structures. On the other hand, the physics-based methods, although still less accurate and more expensive to execute, are independent of databases and give reasonable results where the knowledge-based methods fail because of weak sequence similarity. Therefore, a plausible approach seems to be the use of knowledge-based methods to determine the sections of the structures that correspond to sufficient sequence similarity and physics-based methods to determine the remaining structure. By participating in the 12th Community Wide Experiment on the Critical Assessment of Techniques for Protein Structure Prediction (CASP12) as the KIAS-Gdansk group, we tested our recently developed hybrid approach, in which protein-structure prediction is carried out by using the physics-based UNRES coarse-grained energy function, with restraints derived from the server models. Best predictions among all groups were obtained for 2 targets and 80% of our models were in the upper 50% of the models submitted to CASP. Our method was also able to exclude, with about 70% confidence, the information from the servers that performed poorly on a given target. Moreover, the method resulted in the best models of 2 refinement targets and performed remarkably well on oligomeric targets.

The development of 1,3-diphenylisobenzofuran as a highly selective probe for the detection and quantitative determination of hydrogen peroxide

Abstract 1,3-Diphenylisobenzofuran (DPBF) has been developed as a selective probe for the detection and quantitative determination of hydrogen peroxide in samples containing different reactive nitrogen and oxygen species (RNOS). DPBF is a fluorescent probe which, for almost 20 years, was believed to react in a highly specific manner toward some reactive oxygen species (ROS) such as singlet oxygen and hydroxy, alkyloxy or alkylperoxy radicals. Under the action of these individuals DPBF has been rapidly transformed to 1,2-dibenzoylbenzene (DBB). In order to check if DPBF can act as a unique indicator of the total amount of different RNOS, as well as oxidative stress caused by an overproduction of these individuals, a series of experiments was carried out, in which DPBF reacted with peroxynitrite anion, superoxide anion, hydrogen peroxide, hypochlorite anion, and anions commonly present under biological conditions, namely nitrite and nitrate. In all cases, except for hydrogen peroxide, the product of the reaction is DBB. Only under the action of H2O2 9-hydroxyanthracen-10(9H)-one (oxanthrone) is formed. This product has been identified with the use of fluorescence spectroscopy, NMR spectroscopy, high performance liquid chromatography coupled with mass spectrometry, infrared spectroscopy, elemental analysis, and cyclic voltammetry (CV). A linear relationship was found between a decrease in the fluorescence intensity of DPBF and the concentration of hydrogen peroxide in the range of concentrations of 0.196–3.941 mM. DPBF responds to hydrogen peroxide in a very specific way with the limits of detection and quantitation of 88 and 122.8 μM, respectively. The kinetics of the reaction between DBBF and H2O2 was also studied.

A Novel Method for Force-Field Calibration Based on Maximum-Likelihood Approach and Thermal Unfolding Data

Calibration is the final and critical stage of the design of the force fields for proteins and other biological macromolecules. For proteins, the usual goal of this procedure is to optimize the force-field parameters to reproduce the native structures of selected training proteins. However, the resulting force fields are usually not sufficiently predictive, because only the structures of folded proteins are used. Thus, a force field is not sufficiently trained to distinguish folded structures from misfolded ones. In this work, we propose a novel approach, in which a force field is calibrated with the ensembles of structures determined by NMR at various temperatures that encompass the region of thermal unfolding. The method is based on applying the maximum-likelihood principle. Each conformation of the NMR-determined ensemble at a given temperature is an experimental point and the theoretical probability-density function is represented by a sum of Gaussians centred at the decoys from the corresponding ensembles generated by simulations; in this work the replica exchange molecular dynamics procedure was used. The maximum-likelihood function (-logL) is minimized using the current decoy set, then new decoys are generated with the optimized force-field parameters. The procedure is iterated until convergence. The method was applied to the physics-based coarse-grained UNRES force field developed in our laboratory. On the first attempt, NMR structures of a small alpha-helical protein, the tryptophan cage, were used. The resulting force field predicted correctly the structures of 13 out of 14 alpha-helical proteins with different helix-packing topology and size from 36 to 104 amino-acid residues. Results of the calibration of the UNRES force field with more proteins, including villin headpiece (alpha), the C-terminal fragment of the IGG protein (beta), and full-sequence design 1 (alpha+beta), will be presented.