Analysis of thermodynamic non-ideality in terms of protein solvation

The effects of thermodynamic non-ideality on the forms of sedimentation equilibrium distributions for several isoelectric proteins have been analysed on the statistical-mechanical basis of excluded volume to obtain an estimate of the extent of protein solvation. Values of the effective solvation. parameter delta are reported for ellipsoidal as well as spherical models of the proteins, taken to be rigid, impenetrable macromolecular structures. The dependence of the effective solvated radius upon protein molecular mass exhibits reasonable agreement with the relationship calculated for a model in which the unsolvated protein molecule is surrounded by a 0.52-nm solvation shell. Although the observation that this shell thickness corresponds to a double layer of water molecules may be of questionable relevance to mechanistic interpretation of protein hydration, it augurs well for the assignment of magnitudes to the second virial coefficients of putative complexes in the quantitative characterization of protein-protein interactions under conditions where effects of thermodynamic non-ideality cannot justifiably be neglected. (C) 2001 Elsevier Science B.V. All rights reserved.

Molecular dynamics simulations of proteins: can the explicit water model be varied?

In molecular mechanics simulations of biological systems, the solvation water is typically represented by a default water model which is an integral part of the force field. Indeed, protein nonbonding parameters are chosen in order to obtain a balance between water-water and protein-water interactions and hence a reliable description of protein solvation. However, less attention has been paid to the question of whether the water model provides a reliable description of the water properties under the chosen simulation conditions, for which more accurate water models often exist. Here we consider the case of the CHARMM protein force field, which was parametrized for use with a modified TIP3P model. Using quantum mechanical and molecular mechanical calculations, we investigate whether the CHARMM force field can be used with other water models: TIP4P and TIP5P. Solvation properties of N-methylacetamide (NMA), other small solute molecules, and a small protein are examined. The results indicate differences in binding energies and minimum energy geometries, especially for TIP5P, but the overall description of solvation is found to be similar for all models tested. The results provide an indication that molecular mechanics simulations with the CHARMM force field can be performed with water models other than TIP3P, thus enabling an improved description of the solvent water properties.

Influence of Electromagnetic Fields On Biological Signalling: An Experimental and Theoretical Approach

The primary goals of this study were to develop a cell-free in vitro assay for the assessment of nonthermal electromagnetic (EMF) bioeffects and to develop theoretical models in accord with current experimental observations. Based upon the hypothesis that EMF effects operate by modulating Ca2+/CaM binding, an in vitro nitric oxide (NO) synthesis assay was developed to assess the effects of a pulsed radiofrequency (PRF) signal used for treatment of postoperative pain and edema. No effects of PRF on NO synthesis were observed. Effects of PRF on Ca2+/CaM binding were also assessed using a Ca2+-selective electrode, also yielding no EMF Ca2+/CaM binding. However, a PRF effect was observed on the interaction of hemoglobin (Hb) with tetrahydrobiopterin, leading to the development of an in vitro Hb deoxygenation assay, showing a reduction in the rate of Hb deoxygenation for exposures to both PRF and a static magnetic field (SMF). Structural studies using pyranine fluorescence, Gd3+ vibronic sideband luminescence and attenuated total reflectance Fourier transform infrared (ATR-FTIR) spectroscopy were conducted in order to ascertain the mechanism of this EMF effect on Hb. Also, the effect of SMF on Hb oxygen saturation (SO2) was assessed under gas-controlled conditions. These studies showed no definitive changes in protein/solvation structure or SO2 under equilibrium conditions, suggesting the need for real-time instrumentation or other means of observing out-of-equilibrium Hb dynamics. Theoretical models were developed for EMF transduction, effects on ion binding, neuronal spike timing, and dynamics of Hb deoxygenation. The EMF sensitivity and simplicity of the Hb deoxygenation assay suggest a new tool to further establish basic biophysical EMF transduction mechanisms. If an EMF-induced increase in the rate of deoxygenation can be demonstrated in vivo, then enhancement of oxygen delivery may be a new therapeutic method by which clinically relevant EMF-mediated enhancement of growth and repair processes can occur.

Deuterium isotope effect on the intramolecular electron transfer in Pseudomonas aeruginosa azurin

Intramolecular electron transfer in azurin in water and deuterium oxide has been studied over a broad temperature range. The kinetic deuterium isotope effect, kH/kD, is smaller than unity (0.7 at 298 K), primarily caused by the different activation entropies in water (−56.5 J K−1 mol−1) and in deuterium oxide (−35.7 J K−1 mol−1). This difference suggests a role for distinct protein solvation in the two media, which is supported by the results of voltammetric measurements: the reduction potential (E0′) of Cu2+/+ at 298 K is 10 mV more positive in D2O than in H2O. The temperature dependence of E0′ is also different, yielding entropy changes of −57 J K−1 mol−1 in water and −84 J K−1 mol−1 in deuterium oxide. The driving force difference of 10 mV is in keeping with the kinetic isotope effect, but the contribution to ΔS‡ from the temperature dependence of E0′ is positive rather than negative. Isotope effects are, however, also inherent in the nuclear reorganization Gibbs free energy and in the tunneling factor for the electron transfer process. A slightly larger thermal protein expansion in H2O than in D2O (0.001 nm K−1) is sufficient both to account for the activation entropy difference and to compensate for the different temperature dependencies of E0′. Thus, differences in driving force and thermal expansion appear as the most straightforward rationale for the observed isotope effect.

Use of the FACTS solvation model for protein-ligand docking calculations. Application to EADock.

Protein-ligand docking has made important progress during the last decade and has become a powerful tool for drug development, opening the way to virtual high throughput screening and in silico structure-based ligand design. Despite the flattering picture that has been drawn, recent publications have shown that the docking problem is far from being solved, and that more developments are still needed to achieve high successful prediction rates and accuracy. Introducing an accurate description of the solvation effect upon binding is thought to be essential to achieve this goal. In particular, EADock uses the Generalized Born Molecular Volume 2 (GBMV2) solvent model, which has been shown to reproduce accurately the desolvation energies calculated by solving the Poisson equation. Here, the implementation of the Fast Analytical Continuum Treatment of Solvation (FACTS) as an implicit solvation model in small molecules docking calculations has been assessed using the EADock docking program. Our results strongly support the use of FACTS for docking. The success rates of EADock/FACTS and EADock/GBMV2 are similar, i.e. around 75% for local docking and 65% for blind docking. However, these results come at a much lower computational cost: FACTS is 10 times faster than GBMV2 in calculating the total electrostatic energy, and allows a speed up of EADock by a factor of 4. This study also supports the EADock development strategy relying on the CHARMM package for energy calculations, which enables straightforward implementation and testing of the latest developments in the field of Molecular Modeling.

Direct measurement of salt-bridge solvation energies using a peptide model system: implications for protein stability.

The solvation energies of salt bridges formed between the terminal carboxyl of the host pentapeptide AcWL- X-LL and the side chains of Arg or Lys in the guest (X) position have been measured. The energies were derived from octanol-to-buffer transfer free energies determined between pH 1 and pH 9. 13C NMR measurements show that the salt bridges form in the octanol phase, but not in the buffer phase, when the side chains and the terminal carboxyl group are charged. The free energy of salt-bridge formation in octanol is approximately -4 kcal/mol (1 cal = 4.184 J), which is equal to or slightly larger than the sum of the solvation energies of noninteracting pairs of charged side chains. This is about one-half the free energy that would result from replacing a charge pair in octanol with a pair of hydrophobic residues of moderate size. Therefore, salt bridging in octanol can change the favorable aqueous solvation energy of a pair of oppositely charged residues to neutral or slightly unfavorable but cannot provide the same free energy decrease as hydrophobic residues. This is consistent with recent computational and experimental studies of protein stability.

Design of bio-inspired ionic liquids for protein stabilisation

Ionic Liquids (ILs) consist in organic salts that are liquid at/or near room temperature. Since ILs are entirely composed of ions, the formation of ion pairs is expected to be one essential feature for describing solvation in ILs. In recent years, protein - ionic liquid (P-IL) interactions have been the subject of intensive studies mainly because of their capability to promote folding/unfolding of proteins. However, the ion pairs and their lifetimes in ILs in P-IL thematic is dismissed, since the action of ILs is therefore the result of a subtle equilibrium between anion-cation interaction, ion-solvent and ion-protein interaction. The work developed in this thesis innovates in this thematic, once the design of ILs for protein stabilisation was bio-inspired in the high concentration of organic charged metabolites found in cell milieu. Although this perception is overlooked, those combined concentrations have been estimated to be ~300 mM among the macromolecules at concentrations exceeding 300 g/L (macromolecular crowding) and transient ion-pair can naturally occur with a potential specific biological role. Hence the main objective of this work is to develop new bio-ILs with a detectable ion-pair and understand its effects on protein structure and stability, under crowding environment, using advanced NMR techniques and calorimetric techniques. The choline-glutamate ([Ch][Glu]) IL was synthesized and characterized. The ion-pair was detected in water solutions using mainly the selective NOE NMR technique. Through the same technique, it was possible to detect a similar ion-pair promotion under synthetic and natural crowding environments. Using NMR spectroscopy (protein diffusion, HSQC experiments, and hydrogen-deuterium exchange) and differential scanning calorimetry (DSC), the model protein GB1 (production and purification in isotopic enrichment media) it was studied in the presence of [Ch][Glu] under macromolecular crowding conditions (PEG, BSA, lysozyme). Under dilute condition, it is possible to assert that the [Ch][Glu] induces a preferential hydration by weak and non-specific interactions, which leads to a significant stabilisation. On the other hand, under crowding environment, the [Ch][Glu] ion pair is promoted, destabilising the protein by favourable weak hydrophobic interactions , which disrupt the hydration layer of the protein. However, this capability can mitigates the effect of protein crowders. Overall, this work explored the ion-pair existence and its consequences on proteins in conditions similar to cell milieu. In this way, the charged metabolites found in cell can be understood as key for protein stabilisation.

Protein-ligand binding free energy estimation using molecular mechanics and continuum electrostatics. Application to HIV-1 protease inhibitors.

A method is proposed for the estimation of absolute binding free energy of interaction between proteins and ligands. Conformational sampling of the protein-ligand complex is performed by molecular dynamics (MD) in vacuo and the solvent effect is calculated a posteriori by solving the Poisson or the Poisson-Boltzmann equation for selected frames of the trajectory. The binding free energy is written as a linear combination of the buried surface upon complexation, SASbur, the electrostatic interaction energy between the ligand and the protein, Eelec, and the difference of the solvation free energies of the complex and the isolated ligand and protein, deltaGsolv. The method uses the buried surface upon complexation to account for the non-polar contribution to the binding free energy because it is less sensitive to the details of the structure than the van der Waals interaction energy. The parameters of the method are developed for a training set of 16 HIV-1 protease-inhibitor complexes of known 3D structure. A correlation coefficient of 0.91 was obtained with an unsigned mean error of 0.8 kcal/mol. When applied to a set of 25 HIV-1 protease-inhibitor complexes of unknown 3D structures, the method provides a satisfactory correlation between the calculated binding free energy and the experimental pIC5o without reparametrization.

Attracting cavities for docking. Replacing the rough energy landscape of the protein by a smooth attracting landscape.

Molecular docking is a computational approach for predicting the most probable position of ligands in the binding sites of macromolecules and constitutes the cornerstone of structure-based computer-aided drug design. Here, we present a new algorithm called Attracting Cavities that allows molecular docking to be performed by simple energy minimizations only. The approach consists in transiently replacing the rough potential energy hypersurface of the protein by a smooth attracting potential driving the ligands into protein cavities. The actual protein energy landscape is reintroduced in a second step to refine the ligand position. The scoring function of Attracting Cavities is based on the CHARMM force field and the FACTS solvation model. The approach was tested on the 85 experimental ligand-protein structures included in the Astex diverse set and achieved a success rate of 80% in reproducing the experimental binding mode starting from a completely randomized ligand conformer. The algorithm thus compares favorably with current state-of-the-art docking programs. © 2015 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

Function of the hydration layer around an antifreeze protein revealed by atomistic molecular dynamics simulations

Atomistic molecular dynamics simulations are used to investigate the mechanism by which the antifreeze protein from the spruce budworm, Choristoneura fumiferana, binds to ice. Comparison of structural and dynamic properties of the water around the three faces of the triangular prism-shaped protein in aqueous solution reveals that at low temperature the water structure is ordered and the dynamics slowed down around the ice-binding face of the protein, with a disordering effect observed around the other two faces. These results suggest a dual role for the solvation water around the protein. The preconfigured solvation shell around the ice-binding face is involved in the initial recognition and binding of the antifreeze protein to ice by lowering the barrier for binding and consolidation of the protein:ice interaction surface. Thus, the antifreeze protein can bind to the molecularly rough ice surface by becoming actively involved in the formation of its own binding site. Also, the disruption of water structure around the rest of the protein helps prevent the adsorbed protein becoming covered by further ice growth.

Computational study of ligand binding in lipid transfer proteins: Structures, interfaces, and free energies of protein-lipid complexes

Plant nonspecific lipid transfer proteins (nsLTPs) bind a wide variety of lipids, which allows them to perform disparate functions. Recent reports on their multifunctionality in plant growth processes have posed new questions on the versatile binding abilities of these proteins. The lack of binding specificity has been customarily explained in qualitative terms on the basis of a supposed structural flexibility and nonspecificity of hydrophobic protein-ligand interactions. We present here a computational study of protein-ligand complexes formed between five nsLTPs and seven lipids bound in two different ways in every receptor protein. After optimizing geometries inmolecular dynamics calculations, we computed Poisson- Boltzmann electrostatic potentials, solvation energies, properties of the protein-ligand interfaces, and estimates of binding free energies of the resulting complexes. Our results provide the first quantitative information on the ligand abilities of nsLTPs, shed new light into protein-lipid interactions, and reveal new features which supplement commonly held assumptions on their lack of binding specificity.

Temperature and length scale dependence of hydrophobic effects and their possible implications for protein folding

The Lum–Chandler–Weeks theory of hydrophobicity [Lum, K., Chandler, D. & Weeks, J. D. (1999) J. Phys. Chem. 103, 4570–4577] is applied to treat the temperature dependence of hydrophobic solvation in water. The application illustrates how the temperature dependence for hydrophobic surfaces extending less than 1 nm differs significantly from that for surfaces extending more than 1 nm. The latter is the result of water depletion, a collective effect, that appears at length scales of 1 nm and larger. Because of the contrasting behaviors at small and large length scales, hydrophobicity by itself can explain the variable behavior of entropies of protein folding.

Native protein sequences are close to optimal for their structures

How large is the volume of sequence space that is compatible with a given protein structure? Starting from random sequences, low free energy sequences were generated for 108 protein backbone structures by using a Monte Carlo optimization procedure and a free energy function based primarily on Lennard–Jones packing interactions and the Lazaridis–Karplus implicit solvation model. Remarkably, in the designed sequences 51% of the core residues and 27% of all residues were identical to the amino acids in the corresponding positions in the native sequences. The lowest free energy sequences obtained for ensembles of native-like backbone structures were also similar to the native sequence. Furthermore, both the individual residue frequencies and the covariances between pairs of positions observed in the very large SH3 domain family were recapitulated in core sequences designed for SH3 domain structures. Taken together, these results suggest that the volume of sequence space optimal for a protein structure is surprisingly restricted to a region around the native sequence.

Taurine Supplementation Increases K(atp) Channel Protein Content, Improving Ca2+ Handling And Insulin Secretion In Islets From Malnourished Mice Fed On A High-fat Diet.

Pancreatic β-cells are highly sensitive to suboptimal or excess nutrients, as occurs in protein-malnutrition and obesity. Taurine (Tau) improves insulin secretion in response to nutrients and depolarizing agents. Here, we assessed the expression and function of Cav and KATP channels in islets from malnourished mice fed on a high-fat diet (HFD) and supplemented with Tau. Weaned mice received a normal (C) or a low-protein diet (R) for 6 weeks. Half of each group were fed a HFD for 8 weeks without (CH, RH) or with 5% Tau since weaning (CHT, RHT). Isolated islets from R mice showed lower insulin release with glucose and depolarizing stimuli. In CH islets, insulin secretion was increased and this was associated with enhanced KATP inhibition and Cav activity. RH islets secreted less insulin at high K(+) concentration and showed enhanced KATP activity. Tau supplementation normalized K(+)-induced secretion and enhanced glucose-induced Ca(2+) influx in RHT islets. R islets presented lower Ca(2+) influx in response to tolbutamide, and higher protein content and activity of the Kir6.2 subunit of the KATP. Tau increased the protein content of the α1.2 subunit of the Cav channels and the SNARE proteins SNAP-25 and Synt-1 in CHT islets, whereas in RHT, Kir6.2 and Synt-1 proteins were increased. In conclusion, impaired islet function in R islets is related to higher content and activity of the KATP channels. Tau treatment enhanced RHT islet secretory capacity by improving the protein expression and inhibition of the KATP channels and enhancing Synt-1 islet content.

A Putative Otu Domain-containing Protein 1 Deubiquitinating Enzyme Is Differentially Expressed In Thyroid Cancer And Identifies Less-aggressive Tumours.

This study aimed to identify novel biomarkers for thyroid carcinoma diagnosis and prognosis. We have constructed a human single-chain variable fragment (scFv) antibody library that was selected against tumour thyroid cells using the BRASIL method (biopanning and rapid analysis of selective interactive ligands) and phage display technology. One highly reactive clone, scFv-C1, with specific binding to papillary thyroid tumour proteins was confirmed by ELISA, which was further tested against a tissue microarray that comprised of 229 thyroid tissues, including: 110 carcinomas (38 papillary thyroid carcinomas (PTCs), 42 follicular carcinomas, 30 follicular variants of PTC), 18 normal thyroid tissues, 49 nodular goitres (NG) and 52 follicular adenomas. The scFv-C1 was able to distinguish carcinomas from benign lesions (P=0.0001) and reacted preferentially against T1 and T2 tumour stages (P=0.0108). We have further identified an OTU domain-containing protein 1, DUBA-7 deubiquitinating enzyme as the scFv-binding antigen using two-dimensional polyacrylamide gel electrophoresis and mass spectrometry. The strategy of screening and identifying a cell-surface-binding antibody against thyroid tissues was highly effective and resulted in a useful biomarker that recognises malignancy among thyroid nodules and may help identify lower-risk cases that can benefit from less-aggressive management.