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.

Ligand-Induced structural changes in TEM‑1 probed by molecular dynamics and relative binding free energy calculations

The TEM family of enzymes has had a crucial impact on the pharmaceutical industry due to their important role in antibiotic resistance. Even with the latest technologies in structural biology and genomics, no 3D structure of a TEM- 1/antibiotic complex is known previous to acylation. Therefore, the comprehension of their capability in acylate antibiotics is based on the protein macromolecular structure uncomplexed. In this work, molecular docking, molecular dynamic simulations, and relative free energy calculations were applied in order to get a comprehensive and thorough analysis of TEM-1/ampicillin and TEM-1/amoxicillin complexes. We described the complexes and analyzed the effect of ligand binding on the overall structure. We clearly demonstrate that the key residues involved in the stability of the ligand (hot-spots) vary with the nature of the ligand. Structural effects such as (i) the distances between interfacial residues (Ser70−Oγ and Lys73−Nζ, Lys73−Nζ and Ser130−Oγ, and Ser70−Oγ−Ser130−Oγ), (ii) side chain rotamer variation (Tyr105 and Glu240), and (iii) the presence of conserved waters can be also influenced by ligand binding. This study supports the hypothesis that TEM-1 suffers structural modifications upon ligand binding.

Comparison between computational alanine scanning and per-residue binding free energy decomposition for protein-protein association using MM-GBSA: application to the TCR-p-MHC complex.

Recognition by the T-cell receptor (TCR) of immunogenic peptides (p) presented by Class I major histocompatibility complexes (MHC) is the key event in the immune response against virus-infected cells or tumor cells. A study of the 2C TCR/SIYR/H-2K(b) system using a computational alanine scanning and a much faster binding free energy decomposition based on the Molecular Mechanics-Generalized Born Surface Area (MM-GBSA) method is presented. The results show that the TCR-p-MHC binding free energy decomposition using this approach and including entropic terms provides a detailed and reliable description of the interactions between the molecules at an atomistic level. Comparison of the decomposition results with experimentally determined activity differences for alanine mutants yields a correlation of 0.67 when the entropy is neglected and 0.72 when the entropy is taken into account. Similarly, comparison of experimental activities with variations in binding free energies determined by computational alanine scanning yields correlations of 0.72 and 0.74 when the entropy is neglected or taken into account, respectively. Some key interactions for the TCR-p-MHC binding are analyzed and some possible side chains replacements are proposed in the context of TCR protein engineering. In addition, a comparison of the two theoretical approaches for estimating the role of each side chain in the complexation is given, and a new ad hoc approach to decompose the vibrational entropy term into atomic contributions, the linear decomposition of the vibrational entropy (LDVE), is introduced. The latter allows the rapid calculation of the entropic contribution of interesting side chains to the binding. This new method is based on the idea that the most important contributions to the vibrational entropy of a molecule originate from residues that contribute most to the vibrational amplitude of the normal modes. The LDVE approach is shown to provide results very similar to those of the exact but highly computationally demanding method.

MM-GBSA binding free energy decomposition and T cell receptor engineering.

Recognition by the T-cell receptor (TCR) of immunogenic peptides (p) presented by class I major histocompatibility complexes (MHC) is the key event in the immune response against virus infected cells or tumor cells. The major determinant of T cell activation is the affinity of the TCR for the peptide-MHC complex, though kinetic parameters are also important. A study of the 2C TCR/SIYR/H-2Kb system using a binding free energy decomposition (BFED) based on the MM-GBSA approach had been performed to assess the performance of the approach on this system. The results showed that the TCR-p-MHC BFED including entropic terms provides a detailed and reliable description of the energetics of the interaction (Zoete and Michielin, 2007). Based on these results, we have developed a new approach to design sequence modifications for a TCR recognizing the human leukocyte antigen (HLA)-A2 restricted tumor epitope NY-ESO-1. NY-ESO-1 is a cancer testis antigen expressed not only in melanoma, but also on several other types of cancers. It has been observed at high frequencies in melanoma patients with unusually positive clinical outcome and, therefore, represents an interesting target for adoptive transfer with modified TCR. Sequence modifications of TCR potentially increasing the affinity for this epitope have been proposed and tested in vitro. T cells expressing some of the proposed TCR mutants showed better T cell functionality, with improved killing of peptide-loaded T2 cells and better proliferative capacity compared to the wild type TCR expressing cells. These results open the door of rational TCR design for adoptive transfer cancer therapy.

Binding free energy differences in a TCR-peptide-MHC complex induced by a peptide mutation: a simulation analysis.

Recognition by the T-cell receptor (TCR) of immunogenic peptides presented by class I major histocompatibility complexes (MHCs) is the determining event in the specific cellular immune response against virus-infected cells or tumor cells. It is of great interest, therefore, to elucidate the molecular principles upon which the selectivity of a TCR is based. These principles can in turn be used to design therapeutic approaches, such as peptide-based immunotherapies of cancer. In this study, free energy simulation methods are used to analyze the binding free energy difference of a particular TCR (A6) for a wild-type peptide (Tax) and a mutant peptide (Tax P6A), both presented in HLA A2. The computed free energy difference is 2.9 kcal/mol, in good agreement with the experimental value. This makes possible the use of the simulation results for obtaining an understanding of the origin of the free energy difference which was not available from the experimental results. A free energy component analysis makes possible the decomposition of the free energy difference between the binding of the wild-type and mutant peptide into its components. Of particular interest is the fact that better solvation of the mutant peptide when bound to the MHC molecule is an important contribution to the greater affinity of the TCR for the latter. The results make possible identification of the residues of the TCR which are important for the selectivity. This provides an understanding of the molecular principles that govern the recognition. The possibility of using free energy simulations in designing peptide derivatives for cancer immunotherapy is briefly discussed.

Resummation of infrared divergences in the free energy of spin-two fields

We derive a closed form expression for the sum of all the infrared divergent contributions to the free energy of a gas of gravitons. An important ingredient of our calculation is the use of a gauge fixing procedure such that the graviton propagator becomes both traceless and transverse. This has been shown to be possible, in a previous work, using a general gauge fixing procedure, in the context of the lowest order expansion of the Einstein-Hilbert action, describing noninteracting spin-two fields. In order to encompass the problems involving thermal loops, such as the resummation of the free energy, in the present work, we have extended this procedure to the situations when the interactions are taken into account.

Free energy of Lorentz-violating QED at high temperature

In this paper we study the one-and two-loop contribution to the free energy in QED with Lorentz symmetry breaking introduced via constant CPT-even Lorentz-breaking parameters at the high temperature limit. We find the impact of the Lorentz-violating term for the free energy and carry out a numerical estimation for the Lorentz-breaking parameter.

Scaling of the elastic contribution to the surface free energy of a nematic liquid crystal on a sawtoothed substrate

We characterize the elastic contribution to the surface free energy of a nematic liquid crystal in the presence of a sawtooth substrate. Our findings are based on numerical minimization of the Landau-de Gennes model and analytical calculations on the Frank-Oseen theory. The nucleation of disclination lines (characterized by non-half-integer winding numbers) in the wedges and apexes of the substrate induces a leading order proportional to q ln q to the elastic contribution to the surface free-energy density, with q being the wave number associated with the substrate periodicity.

Density dependence of the symmetry free energy of hot nuclei

The density and excitation energy dependence of symmetry energy and symmetry free energy for finite nuclei are calculated microscopically in a microcanonical framework, taking into account thermal and expansion effects. A finite-range momentum and density-dependent two-body effective interaction is employed for this purpose. The role of mass, isospin, and equation of state (EOS) on these quantities is also investigated; our calculated results are in consonance with the available experimental data.

The Gibbs free energy of formation of a glass alloy

A simple expression for the Gibbs free energy of formation of a pure component or a eutectic alloy glass, relative to the stable crystalline phase (or phases) at the same temperature is deduced by use of thermodynamic arguments. The expression obtained is supposed to apply to both monocomponent and multicomponent liquid alloys that might become glasses from the supercooled liquid state, irrespective of the critical cooling rate needed to avoid crystallization.

On the Ewald method and the quartic Helmholtz free energy of an anharmonic metallic crystal

Thesis (M. Sc.) - Brock University, 1978.

The Landau free energy of hard ellipses obtained from microscopic simulations

Systems of two-dimensional hard ellipses of varying aspect ratios and packing fractions are studied by Monte Carlo simulations in the generalised canonical ensemble. From this microscopic model, we extract a coarse-grained macroscopic Landau-de Gennes free energy as a function of packing fraction and orientational order parameter. We separate the free energy into the ideal orientational entropy of non-interacting two-dimensional spins and an excess free energy associated with excluded volume interactions. We further explore the isotropic-nematic phase transition using our empirical expression for the free energy and find that the nature of the phase transition is continuous for the aspect ratios we studied.

Effect of collagen matrix saturation on the surface free energy of dentin using different agents

The surface free energy of conditioned-dentin is one of the factors that interfere with monomeric infiltration of the interfibrillar spaces. Saturation of the tooth matrix with different substances may modulate this energy and, consequently, the wettability of the dentin. To evaluate the influence of different substances used to saturate conditioned-dentin on surface free energy (SFE) of this substrate. Dentin blocks (4 × 7 × 1 mm, n = 6/ group), obtained from the roots of bovine incisors, were etched using phosphoric acid for 15 seconds, rinsed and gently dried. The surfaces were treated for 60 seconds with: ultra-purified water (H20-control); ethanol (EtOH), acetone (ACT), chlorhexidine (CHX), ethylenediaminetetraacetic acid (EDTA); or sodium hypochlorite (NaOCl). The tooth surfaces were once again dried with absorbent paper and prepared for SFE evaluation using three standards: water, formamide and bromonaphthalene. Analysis of variance (ANOVA) and Dunnet's tests (a = 0.05) were applied to the data. Ethylenediaminetetraacetic acid was the only substance that caused a change to the contact angle for the standards water and formamide, while only EtOH influenced the angles formed between formamide and the dentin surface. None of the substances exerted a significant effect for bromonaphtha-lene. In comparison to the control, only EDTA and NaOCl altered both polar components of the SFE. Total SFE was increased by saturation of the collagen matrix by EDTA and reduced when NaOCl was used. Saturation of the collagen matrix by EDTA and EtOH changed the surface free energy of the dentin. In addition, the use of NaOCl negatively interfered with the properties evaluated. The increase of surface free energy and wettability of the dentin surface would allow higher penetration of the the adhesive system, which would be of importance to the clinical success of resin-dentin union.

Evaluation of the surface free energy of plant surfaces: toward standardizing the procedure

Plant surfaces have been found to have a major chemical and physical heterogeneity and play a key protecting role against multiple stress factors. During the last decade, there is a raising interest in examining plant surface properties for the development of biomimetic materials. Contact angle measurement of different liquids is a common tool for characterizing synthetic materials, which is just beginning to be applied to plant surfaces. However, some studies performed with polymers and other materials showed that for the same surface, different surface free energy values may be obtained depending on the number and nature of the test liquids analyzed, materials' properties, and surface free energy calculation methods employed. For 3 rough and 3 rather smooth plant materials, we calculated their surface free energy using 2 or 3 test liquids and 3 different calculation methods. Regardless of the degree of surface roughness, the methods based on 2 test liquids often led to the under- or over-estimation of surface free energies as compared to the results derived from the 3-Liquids method. Given the major chemical and structural diversity of plant surfaces, it is concluded that 3 different liquids must be considered for characterizing materials of unknown physico-chemical properties, which may significantly differ in terms of polar and dispersive interactions. Since there are just few surface free energy data of plant surfaces with the aim of standardizing the calculation procedure and interpretation of the results among for instance, different species, organs, or phenological states, we suggest the use of 3 liquids and the mean surface tension values provided in this study.

Free energy of burying hydrophobic residues in the interface between protein subunits

We have obtained an experimental estimate of the free energy change associated with variations at the interface between protein subunits, a subject that has raised considerable interest since the concept of accessible surface area was introduced by Lee and Richards [Lee, B. & Richards, F. M. (1971) J. Mol. Biol. 55, 379–400]. We determined by analytical ultracentrifugation the dimer–tetramer equilibrium constant of five single and three double mutants of human Hb. One mutation is at the stationary α1β1 interface, and all of the others are at the sliding α1β2 interface where cleavage of the tetramer into dimers and ligand-linked allosteric changes are known to occur. A surprisingly good linear correlation between the change in the free energy of association of the mutants and the change in buried hydrophobic surface area was obtained, after corrections for the energetic cost of losing steric complementarity at the αβ dimer interface. The slope yields an interface stabilization free energy of −15 ± 1.2 cal/mol upon burial of 1 Å2 of hydrophobic surface, in very good agreement with the theoretical estimate given by Eisenberg and McLachlan [Eisenberg, D. & McLachlan, A. D. (1986) Nature (London) 319, 199–203].