Thermodynamic aspects of cell spreading on solid substrata

To verify the validity of thermodynamic approaches to the prediction of cellular behavior, cell spreading of three different cell types on solid substrata was determined in vitro. Solid substrata as well as cell types were selected on the basis of their surface free energies, calculated from contact angle measurements. The surface free energies of the solid substrata ranged from 18–116 erg cm⁻². To measure contact angles on cells, a technique was developed in which a multilayer of cells was deposited on a filter and air dried. Cell surface free energies ranged from 60 erg cm⁻² for fibroblasts, and 57 for smooth muscle cells, to 91 for HeLa epithelial cells. After adsorption of serum proteins, cell surface free energies of all three cell types converged to approx 74 erg cm⁻². The spreading of these cell types from RPMI 1640 medium on the various solid substrata showed that both in the presence and in the absence of serum proteins in the medium, cells spread poorly on low energy substrata (Y _s <50 erg cm⁻²), whereas good cell spreading was observed on the higher energy substrata. Calculations of the interfacial free energy of adhesion (ΔF _adh) show that ΔF _adh decreases with increasingY _s , and equals zero around 45 erg cm⁻² for all three cell types in the presence of serum proteins and for HeLa epithelium cells in the absence of serum proteins. This explains the spreading of these cells on the various substrata upon a thermodynamic basis. The results clearly show that substratum surface free energy has a predictive value with respect to cell spreading in vitro, both in the presence and absence of serum proteins. It is noted, however, that interfacial thermodynamics fail to explain the behavior of fibroblasts and smooth muscle cells in the absence of serum proteins, most likely because of the relatively high surface charges of these two cell types.     

平衡型氨基酸透析液对维持性血液透析患者血浆游离氨基酸含量的影响

目的:研究一种平衡型氨基酸透析液,探讨其对维持性血透患者血游离氨基酸含量的影响。方法:实验采用自身对照设计,应用蛋白水解法测定、分析应用氨基酸透析液血透前后患者血浆氨基酸含量,并同时测定应用碳酸盐透析液患者及健康人血浆游离氨基酸作为比较。结果:经碳酸盐透析后,患者血浆大多数游离氨基酸和总氨基酸含量显著降低。给予平衡氨基酸透析液能不同程度改善患者血浆游离氨基酸的含量。结论:平衡型氨基酸透析液能减少血中部分氨基酸的丢失。     

QM/MM free energy simulations: recent progress and challenges

Due to the higher computational cost relative to pure molecular mechanical (MM) simulations, hybrid quantum mechanical/molecular mechanical (QM/MM) free energy simulations particularly require a careful consideration of balancing computational cost and accuracy. Here, we review several recent developments in free energy methods most relevant to QM/MM simulations and discuss several topics motivated by these developments using simple but informative examples that involve processes in water. For chemical reactions, we highlight the value of invoking enhanced sampling technique (e.g. replica-exchange) in umbrella sampling calculations and the value of including collective environmental variables (e.g. hydration level) in metadynamics simulations; we also illustrate the sensitivity of string calculations, especially free energy along the path, to various parameters in the computation. Alchemical free energy simulations with a specific thermodynamic cycle are used to probe the effect of including the first solvation shell into the QM region when computing solvation free energies. For cases where high-level QM/MM potential functions are needed, we analyse two different approaches: the QM/MM–MFEP method of Yang and co-workers and perturbative correction to low-level QM/MM free energy results. For the examples analysed here, both approaches seem productive although care needs to be exercised when analysing the perturbative corrections.     

Surface free energy and interaction of Staphylococcus epidermidis with biomaterials

The adhesion of twenty nine Staphylococcus epidermidis strains to teflon, polyethylene, polycarbonate and bovine pericardium was studied in vitro and examined in relation to the surface free energies of both bacteria and biomaterials. All S. epidermidis strains had similar surface free energies, close to that of water, and adhered better to the materials with analogous surface free energies. There was a significant correlation (Kendall's Tau B = 1000) of biomaterial's surface free energy with the number of adhering bacteria. This correlation is inverse (Kendall's Tau B = -1000) when surface hydrophobicity is considered instead of surface free energy. This indicates that in Staphylococcus epidermidis adherence to biomaterials is inversely correlated to the surface hydrophobicity of the last, being so just the opposite of that occurring with other bacteria.     

Physicochemical characterization ofEscherichia coli

EightEscherichia coli strains were characterized by determining their adhesion to xylene, surface free energy, zeta potential, relative surface charge, and their chemical composition. The latter was done by applying X-ray photoelectron spectroscopy (XPS) and infrared spectroscopy (IR). No relationship between the adhesion to xylene and the water contact angles of these strains was found. Three strains had significantly lower surface free energies than the other strains. Surface free energies were either obtained from polar and dispersion parts or from Lifshitz-van der Waals and acid/base parts of the surface free energy. A correlation (r=0.97) between the polar parts and the electron-donor contributions to the acid/base part of the surface free energy was found. The zeta potentials of all strains, measured as a function of pH (2–11), were negative. Depending on the zeta potential as a function of pH, three groups were recognized among the strains tested. A relationship (r=0.84) was found between the acid/base component of the surface free energy and the zeta potential measured at pH=7.4. There was no correlation between results of XPS and IR studies. Data from the literature of XPS and IR studies of the gram-positive staphylococci and streptococci were compared with data from the gram-negativeE. coli used in this study. It appeared that in these three groups of bacteria, the polysaccharide content detected by IR corresponded well with the oxygen-to-carbon ratio detected by XPS.     

Insights into protein-protein binding by binding free energy calculation and free energy decomposition for the Ras-Raf and Ras-RalGDS complexes

Absolute binding free energy calculations and free energy decompositions are presented for the protein-protein complexes H-Ras/C-Raf1 and H-Ras/RalGDS. Ras is a central switch in the regulation of cell proliferation and differentiation. In our study, we investigate the capability of the molecular mechanics (MM)-generalized Born surface area (GBSA) approach to estimate absolute binding free energies for the protein-protein complexes. Averaging gas-phase energies, solvation free energies, and entropic contributions over snapshots extracted from trajectories of the unbound proteins and the complexes, calculated binding free energies (Ras-Raf: -15.0(+/-6.3)kcal mol(-1); Ras-RalGDS: -19.5(+/-5.9)kcal mol(-1)) are in fair agreement with experimentally determined values (-9.6 kcal mol(-1); -8.4 kcal mol(-1)), if appropriate ionic strength is taken into account. Structural determinants of the binding affinity of Ras-Raf and Ras-RalGDS are identified by means of free energy decomposition. For the first time, computationally inexpensive generalized Born (GB) calculations are applied in this context to partition solvation free energies along with gas-phase energies between residues of both binding partners. For selected residues, in addition, entropic contributions are estimated by classical statistical mechanics. Comparison of the decomposition results with experimentally determined binding free energy differences for alanine mutants of interface residues yielded correlations with r(2)=0.55 and 0.46 for Ras-Raf and Ras-RalGDS, respectively. Extension of the decomposition reveals residues as far apart as 25A from the binding epitope that can contribute significantly to binding free energy. These "hotspots" are found to show large atomic fluctuations in the unbound proteins, indicating that they reside in structurally less stable regions. Furthermore, hotspot residues experience a significantly larger-than-average decrease in local fluctuations upon complex formation. Finally, by calculating a pair-wise decomposition of interactions, interaction pathways originating in the binding epitope of Raf are found that protrude through the protein structure towards the loop L1. This explains the finding of a conformational change in this region upon complex formation with Ras, and it may trigger a larger structural change in Raf, which is considered to be necessary for activation of the effector by Ras.     

The basis of the hydrophobic effect

The property of a molecule that most reliably determines the magnitude of the hydrophobic effect that it will experience is the number of hydrogen–carbon bonds it contains not the accessible surface area of its nonpolar portions. This conclusion follows from an examination of the standard free energies of transfer of alkanes, alkenes, alkadienes, and arenes from water to hexadecane. When the standard free energies of transfer for hydrocarbons in these different classes are plotted as a function of the number of hydrogen–carbon bonds they contain, all of the data fall upon the same line. These standard free energies of transfer are also directly proportional to the number of hydrogen–carbon bonds the hydrocarbons contain. When the same standard free energies of transfer are plotted as a function of the accessible surface areas of the hydrocarbons, the data do not fall upon the same line nor are the standard free energies of transfer directly proportional to the accessible surface areas. An examination of the standard free energies of transfer for the different classes of hydrocarbons from the gas phase to water and from the gas phase to hexadecane reinforces the conclusion that the number of hydrogen–carbon bonds in a molecule rather than its accessible surface area is the basis of the hydrophobic effect. Consequently, estimates of the contribution of different functional groups to the hydrophobic effect providing the free energy of folding of a molecule of protein or providing the free energy of dissociation for the association of two proteins or the association of a ligand with a protein should be made by counting the number of hydrogen–carbon bonds excluded from water rather than computing the accessible surface areas excluded from water.     

Thermodynamically feasible kinetic models of reaction networks

The dynamics of biological reaction networks are strongly constrained by thermodynamics. An holistic understanding of their behavior and regulation requires mathematical models that observe these constraints. However, kinetic models may easily violate the constraints imposed by the principle of detailed balance, if no special care is taken. Detailed balance demands that in thermodynamic equilibrium all fluxes vanish. We introduce a thermodynamic-kinetic modeling (TKM) formalism that adapts the concepts of potentials and forces from irreversible thermodynamics to kinetic modeling. In the proposed formalism, the thermokinetic potential of a compound is proportional to its concentration. The proportionality factor is a compound-specific parameter called capacity. The thermokinetic force of a reaction is a function of the potentials. Every reaction has a resistance that is the ratio of thermokinetic force and reaction rate. For mass-action type kinetics, the resistances are constant. Since it relies on the thermodynamic concept of potentials and forces, the TKM formalism structurally observes detailed balance for all values of capacities and resistances. Thus, it provides an easy way to formulate physically feasible, kinetic models of biological reaction networks. The TKM formalism is useful for modeling large biological networks that are subject to many detailed balance relations.     

Estimation of the free energy of hard-sphere crystals via a free-volume approach

The free energies of the face-centred (FCC) and base-centred cubic (BCC) hard-sphere (HS) crystals have been estimated via the free-volume approach. We present two free-volume equations for the FCC and BCC HS crystals, which are different from those predicted by Velasco et al. [Langmuir 14(19) (1998), 5652–5655], and the equations exhibit more accuracy than Velasco et al.'s equations. The limitation of using the free-volume approach was assessed by comparing with true free energies obtained by other well-known methods, and possible improvement and application are discussed.     

Investigation on the binding mechanism of loratinib with the c-ros oncogene 1 (ROS1) receptor tyrosine kinase via molecular dynamics simulation and binding free energy calculations

The c-ros oncogene 1 (ROS1) has proven to be an important cancer target for the treatment of various human cancers. The anaplastic lymphoma kinase inhibitor crizotinib has been granted approval for the treatment of patients with ROS1 positive metastatic non-small-cell lung cancer by the Food and Drug Administration on 2016. However, serious resistance due to the secondary mutation of glycine 2032 to arginine (G2032R) was developed in clinical studies. Loratinib (PF-06463922), a macrocyclic analog of crizotinib, showed significantly improved inhibitory activity against wild–type (WT) ROS1 and ROS1^G2032R mutant. To provide insights into the inhibition mechanism, molecular dynamics simulations and free energy calculations were carried out for the complexes of loratinib with WT and G2032R mutated ROS1. The apo-ROS1^WT and apo-ROS1^G2032R systems showed similar RMSF distributions, while ROS1^G2032R-loratinib showed significantly higher than that of WT ROS1-loratinib, which revealed that the binding of loratinib to ROS1^G2032R significantly interfered the ?uctuation of protein. Calculations of binding free energies indicate that G2032R mutation significantly reduces the binding affinity of loratinib for ROS1, which arose mostly from the increase of conformation entropy and the decrease of solvation energy. Furthermore, detailed per-residue binding free energies highlighted the increased and decreased contributions of some residues in the G2032R mutated systems. The present study revealed the detailed inhibitory mechanism of loratinib as potent WT and G2032R mutated ROS1 inhibitor, which was expected to provide a basis for rational drug design.     

A computational analysis of pyrazole-based inhibitors binding to Hsp90 using molecular dynamics simulation and the MM-GBSA method

Owing to the key role of heat-shock protein 90 (Hsp90) in the evolution, development and disease pathogenesis of cancer, it has been an important target for anti-cancer chemotherapy over the years. A five-nanosecond molecular dynamics simulation combined with the calculation of the binding free energy was carried out to investigate the binding mechanisms of three Hsp90 inhibitors 4BH, 2E1 and 2D9 to Hsp90. The binding free energy of each complex was computed using the molecular mechanics–generalised Born surface area method. Detailed binding free energies between each inhibitor and residues of Hsp90 were calculated using a per-residue basis decomposition method. The detailed inhibitor–residue interaction provides insights into binding mechanisms and in-depth understanding of the structure–affinity relationship. This study suggests that van der Waals energy is primarily responsible for driving the binding of the inhibitors to Hsp90, and the three inhibitors bind to Hsp90 in a similar binding mode. However, a substituent in 2D9 leads to higher binding free energy than the other two inhibitors. These data may assist in designing new potent drugs to combat cancer.     

Insights into the functional role of protonation states in the HIV-1 protease-BEA369 complex: molecular dynamics simulations and free energy calculations

The molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method combined with molecular dynamics (MD) simulations were used to investigate the functional role of protonation in human immunodeficiency virus type 1 (HIV-1) protease complexed with the inhibitor BEA369. Our results demonstrate that protonation of two aspartic acids (Asp25/Asp25′) has a strong influence on the dynamics behavior of the complex, the binding free energy of BEA369, and inhibitor–residue interactions. Relative binding free energies calculated using the MM-PBSA method show that protonation of Asp25 results in the strongest binding of BEA369 to HIV-1 protease. Inhibitor–residue interactions computed by the theory of free energy decomposition also indicate that protonation of Asp25 has the most favorable effect on binding of BEA369. In addition, hydrogen-bond analysis based on the trajectories of the MD simulations shows that protonation of Asp25 strongly influences the water-mediated link of a conserved water molecule, Wat301. We expect that the results of this study will contribute significantly to binding calculations for BEA369, and to the design of high affinity inhibitors.     

Computational analysis of binding of P1 variants to trypsin

The binding of P1 variants of bovine pancreatic trypsin inhibitor (BPTI) to trypsin has been investigated by means of molecular dynamics simulations. The specific interaction formed between the amino acid at the primary binding (P1) position of the binding loop of BPTI and the specificity pocket of trypsin was estimated by use of the linear interaction energy (LIE) method. Calculations for 13 of the naturally occurring amino acids at the P1 position were carried out, and the results obtained were found to correlate well with the experimental binding free energies. The LIE calculations rank the majority of the 13 variants correctly according to the experimental association energies and the mean error between calculated and experimental binding free energies is only 0.38 kcal/mole, excluding the Glu and Asp variants, which are associated with some uncertainties regarding protonation and the possible presence of counter-ions. The three-dimensional structures of the complex with three of the P1 variants (Asn, Tyr, and Ser) included in this study have not at present been solved by any experimental techniques and, therefore, were modeled on the basis of experimental data from P1 variants of similar size. Average structures were calculated from the MD simulations, from which specific interactions explaining the broad variation in association energies were identified. The present study also shows that explicit treatment of the complex water-mediated hydrogen bonding network at the protein-protein interface is of crucial importance for obtaining reliable binding free energies. The successful reproduction of relative binding energies shows that this type of methodology can be very useful as an aid in rational design and redesign of biologically active macromolecules.     

Calculation of binding affinities of HIV-1 RT and β-secretase inhibitors using the linear interaction energy method with explicit and continuum solvation approaches

The linear interaction energy (LIE) approach has been applied to estimate the binding free energies of representative sets of HIV-1 RT and β-Secretase inhibitors, using both molecular dynamics (MD) and tethered energy minimization sampling protocols with the OPLS-AA potential, using a range of solvation methodologies. Generalized Born (GB), ‘shell’ and periodic boundary condition (PBC) solvation were used, the latter with reaction field (RF) electrostatics. Poisson-Boltzmann (PB) and GB continuum electrostatics schemes were applied to the simulation trajectories for each solvation type to estimate the electrostatic ligand-water interaction energy in both the free and bound states. Reasonable agreement of the LIE predictions was obtained with respect to experimental binding free energy estimates for both systems: for instance, ‘PB’ fits on MD trajectories carried out with PBC solvation and RF electrostatics led to models with standard errors of 1.11 and 1.03 kcal mol⁻¹ and coefficients of determination, r ² of 0.76 and 0.75 for the HIV-1 RT and β-Secretase sets. However, it was also found that results from MD sampling using PBC solvation provided only slightly better fits than from simulations using shell or Born solvation or tethered energy minimization sampling. Figure Evolution of the running averages for compound H11 (binding to HIV-1RT) of the bound state ligand-water and ligand-protein interaction energies. The ligand-water electrostatic terms are twice the corresponding GB and PB electrostatic solvation free energies. The ligand-receptor van der Waals and Coulombic interaction energies are also shown, in addition to the ligand-water van der Waals interaction term. The terms were calculated (without application of a cut-off) from a trajectory sampled under PBC solvation with reaction field electrostatics Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Calculation of conformational free energies by confinement simulations in explicit water with implicit desolvation

Abstract

The confinement method is a robust and conceptually simple free energy simulation method that allows the calculation of conformational free energy differences between highly dissimilar states. Application of the method to explicitly solvated systems requires a multi-stage simulation protocol for the calculation of desolvation free energies. Here we show that these desolvation free energies can be readily obtained from an implicit treatment, which is simpler and less costly. The accuracy and robustness of this protocol was shown by the calculation of conformational free energy differences of a series of explicitly solvated test systems. Given the accuracy and ease by which these free energy differences were obtained, the confinement method is promising for the treatment of conformational changes in large and complex systems.     

A structural and free energy analysis of Ag complexes to five small peptides

The complexes of Ag⁺ with the peptides MetGly, ProGly, GlyPro, GlyHis and GlyProAla were investigated using hybrid density functional theory at the B3LYP/DZVP level. The silver ion binding free energies at 298 K to each of these peptides was calculated to be 60.8, 52.0, 54.3, 71.2 and 63.3 kcal mol⁻¹, respectively. Structural information and relative free energies are presented for several isomers for each of the five complexes. Each of the global minima found for the five complexes is a charge-solvated ion. An important finding is that the Ag⁺-ProGly is the only complex where a salt bridge structure is energetically favored occurring at 4.0 kcal mol⁻¹ higher in free energy than the global minimum. The Ag⁺ ion in this salt bridge structure is attached to the carboxylate anion of zwitterionic ProGly in which the terminal amino nitrogen is protonated. For all the other complexes studied, the salt bridge structure occurs at much higher energies. All the dipeptide complexes with Ag⁺, but one, exhibit a di- or tri-coordinate metal where the sites of attachment are amino and carbonyl groups. However, the highest coordination numbers are not always the global minima due to steric costs. The global minimum of the Ag⁺-GlyProAla complex is the only structure found in this study where the metal is tetra-coordinated, binding to the terminal amino nitrogen and all three carbonyl oxygen atoms. Silver binding to sulphur and imidazole nitrogen atoms of MetGly and GlyHis, respectively, are present in the three most energetically favored species in each of these cases.     

Functional and geographical differentiation of candidate balanced polymorphisms in Arabidopsis thaliana

Molecular population genetic analysis of three chromosomal regions in Arabidopsis thaliana suggested that balancing selection might operate to maintain variation at three novel candidate adaptive trait genes, including SOLUBLE STARCH SYNTHASE I (SSI) , PLASTID TRANSCRIPTIONALLY ACTIVE 7(PTAC7) , and BELL-LIKE HOMEODOMAIN 10 (BLH10). If balanced polymorphisms are indeed maintained at these loci, then we would expect to observe functional variation underlying the previously detected signatures of selection. We observe multiple replacement polymorphisms within and in the 32 amino acids just upstream of the protein–protein interacting BELL domain at the BLH10 locus. While no clear protein sequence differences are found between allele types in SSI and PTAC7, these two genes show evidence for allele-specific variation in expression levels. Geographical patterns of allelic differentiation seem consistent with population stratification in this species and a significant longitudinal cline was observed at all three candidate loci. These data support a hypothesis of balancing selection at all three candidate loci and provide a basis for more detailed functional work by identifying possible functional differences that might be selectively maintained.     

Algorithm independent properties of RNA secondary structure predictions

Algorithms predicting RNA secondary structures based on different folding criteria – minimum free energies (mfe), kinetic folding (kin), maximum matching (mm) – and different parameter sets are studied systematically. Two base pairing alphabets were used: the binary GC and the natural four-letter AUGC alphabet. Computed structures and free energies depend strongly on both the algorithm and the parameter set. Statistical properties, such as mean number of base pairs, mean numbers of stacks, mean loop sizes, etc., are much less sensitive to the choice of parameter set and even of algorithm. Some features of RNA secondary structures, such as structure correlation functions, shape space covering and neutral networks, seem to depend only on the base pairing logic (GC or AUGC alphabet). Received: 16 May 1996 / Accepted: 10 July 1996     

Solvation effects in calculated electrostatic association free energies for the C3d‐CR2 complex and comparison with experimental data

The complement system is an integral part of the innate immune system that participates in the clearance of pathogens from the body. The association between complement protein fragment C3d and B or T cell‐receptor complement receptor (CR) 2 represents a crucial link between innate and adaptive immunities. The goal of this study is to predict association abilities of C3d and CR2 mutants by theoretically calculating electrostatic free energies of association and to assess the importance of solvation effects in the calculations. We demonstrate that calculated solvation free energy differences and Coulombic free energies of association are more sensitive than electrostatic free energies of association in solution and, thus, more accurate in predicting previously published experimental data for the association abilities (relative to the parent proteins) of specific C3d and CR2 mutants. We show that a proportional relationship exists between the predicted solvation free energy differences and the experimental data, while an inversely proportional relationship exists between the predicted Coulombic free energies of association and the experimental data. Our results yield new insights into the physicochemical properties underlying C3d‐CR2 association. We discuss the predictive validity of Coulombic, solvation, and solution electrostatic free energies of association and the generalization of our method for theoretical mutagenesis studies of other systems. This is a basic study, aimed toward improving our understanding of the theoretical basis of immune system regulation at the molecular level. Such insight can serve as the groundwork for the design of regulators with tailored properties, vaccines, and other biotechnology products. © 2010 Wiley Periodicals, Inc. Biopolymers 93:509–519, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com