Biomedical subjects

K A Sharp

Publications and source records attributed to K A Sharp.

At least 19 recordsLinked to original sources

PubMed2001

Hydration heat capacity of nucleic acid constituents determined from the random network model.

The heat capacities of hydration (dCp) of the five nucleic acid bases A, G, C, T, and U, the sugars ribose and deoxyribose, and the phosphate backbone were determined using Monte Carlo simulations and the random network model. Solute-induced changes in the mean length and root mean square angle of hydrogen bonds between hydration shell waters were used to compute dCp for these solutes. For all solutes the dCp is significantly more positive than predicted from accessible surface area (ASA) models of heat capacity. In ASA models, nitrogen, oxygen, and phosphorus atoms are considered as uniformly polar, therefore making a negative contribution to dCp. However, the simulations show that many of these polar atoms are hydrated by water whose hydrogen bonds are less distorted than in bulk, leading to a positive dCp. This is in contrast to the effect of polar groups seen previously in small molecules and amino acids, which increase the water H-bond distortion, giving negative dCp contributions. Our results imply that dCp accompanying DNA dehydration in DNA-ligand and DNA-protein binding reactions may be significantly more negative than previously believed and that dehydration is a significant contributor to the large decrease in heat capacity seen in experiments.

Adenine↗

PubMed1999

The maximal affinity of ligands.

We explore the question of what are the best ligands for macromolecular targets. A survey of experimental data on a large number of the strongest-binding ligands indicates that the free energy of binding increases with the number of nonhydrogen atoms with an initial slope of approximately -1.5 kcal/mol (1 cal = 4.18 J) per atom. For ligands that contain more than 15 nonhydrogen atoms, the free energy of binding increases very little with relative molecular mass. This nonlinearity is largely ascribed to nonthermodynamic factors. An analysis of the dominant interactions suggests that van der Waals interactions and hydrophobic effects provide a reasonable basis for understanding binding affinities across the entire set of ligands. Interesting outliers that bind unusually strongly on a per atom basis include metal ions, covalently attached ligands, and a few well known complexes such as biotin-avidin.

Binding Sites↗

PubMed1999

Exploration of the structural features defining the conduction properties of a synthetic ion channel.

The finite-difference Poisson-Boltzmann methodology was applied to a series of parallel, alpha-helical bundle models of the designed ion channel peptide Ac-(LSSLLSL)3-CONH2. This method is able to fully describe the current-voltage curves for this channel and quantitatively explains their cation selectivity and rectification. We examined a series of energy-minimized models representing different aggregation states, side-chain rotamers, and helical rotations, as well as an ensemble of structures from a molecular dynamics trajectory. Potential energies were computed for single, permeating K+ and Cl- ions at a series of positions along a central pathway through the models. A variable-electric-field Nernst-Planck electrodiffusion model was used, with two adjustable parameters representing the diffusion coefficients of K+ and Cl- to scale the individual ion current magnitudes. The ability of a given DelPhi potential profile to fit the experimental data depended strongly on the magnitude of the desolvation of the permeating ion. Below a pore radius of 3.8 A, the predicted profiles showed large energy barriers, and the experimental data could be fit only with unrealistically high values for the K+ and Cl- diffusion coefficients. For pore radii above 3.8 A, the desolvation energies were 2kT or less. The electrostatic calculations were sensitive to positioning of the Ser side chains, with the best fits associated with maximum exposure of the Ser side-chain hydroxyls to the pore. The backbone component was shown to be the major source of asymmetry in the DelPhi potential profiles. Only two of the energy-minimized structures were able to explain the experimental data, whereas an average of the dynamics structures gave excellent agreement with experimental results. Thus this method provides a promising approach to prediction of current-voltage curves from three-dimensional structures of ion channel proteins.

Cations↗

PubMed1999

Calculating the electrostatic properties of RNA provides new insights into molecular interactions and function.

Solutions to the nonlinear Poisson-Boltzmann equation were used to obtain the electrostatic potentials of RNA molecules that have known three-dimensional structures. The results are described in terms of isopotential contours and surface electrostatic potential maps. Both representations have unexpected features: 'cavities' within isopotential contours and areas of enhanced negative potential on molecular surfaces. Intriguingly, the sites of unusual electrostatic features correspond to functionally important regions, suggesting that electrostatic properties play a key role in RNA recognition and stabilization. These calculations reveal that the electrostatic potentials generated by RNA molecules have a variety of functionally important characteristics that cannot be discerned by simple visual inspection of the molecular structure.

Anticodon↗

PubMed1998

Calculation of HyHel10-lysozyme binding free energy changes: effect of ten point mutations.

The change in free energy of binding of hen egg white lysozyme (HEL) to the antibody HyHel-10 arising from ten point mutations in HEL (D101K, D101G, K96M, K97D, K97G, K97G, R21E, R21K, W62Y, and W63Y) was calculated using a combination of the finite difference Poisson-Boltzmann method for the electrostatic contribution, a solvent accessible surface area term for the non-polar contribution, and rotamer counting for the sidechain entropy contribution. Comparison of experimental and calculated results indicate that because of pKa shifts in some of the mutated residues, primarily those involving Aspartate or Glutamate, proton uptake or release occurs in binding. When this effect was incorporated into the binding free energy calculations, the agreement with experiment improved significantly, and resulted in a mean error of about 1.9 kcal/mole. Thus these calculations predict that there should be a significant pH dependence to the change in binding caused by these mutations. The other major contributions to binding energy changes comes from solvation and charge charge interactions, which tend to oppose each other. Smaller contributions come from nonpolar interactions and sidechain entropy changes. The structures of the HyHel-10-HEL complexes with mutant HEL were obtained by modeling, and the effect of the modeled structure on the calculations was also examined. "Knowledge based" modeling and automatic generation of models using molecular mechanics produced comparable results.

Animals↗

PubMed1998

The role of protonation and metal chelation preferences in defining the properties of mercury-binding coiled coils.

To define the delicate interplay between metal chelation, protein folding and function in metalloproteins, a family of de novo-designed peptides was synthesized that self-assemble in aqueous solution to form two and three-stranded alpha-helical coiled coils. Each peptide contains a single Cys residue at an a or d position of the heptad repeat. Peptide association thus produces a Cys-rich coordination environment that has been used to bind Hg(II) ions. These peptides display a pH-dependent association, with trimers observed above the pKa of Glu side-chains and dimers below this value. Finite-difference Poisson-Boltzmann calculations suggest that the dimeric state decreases the unfavorable electrostatic interactions between positively charged Lys side-chains (relative to the trimer). The Cys-containing peptides bind Hg(II) in a position-dependent fashion. Cys at a positions form three-coordinate Hg complexes at high pH where the trimeric aggregation state predominates, and two-coordinate complexes at lower pH. A d position Cys, however, is only able to generate the two-coordinate complex, illustrating the difference in coordination geometry between the two positions in the coiled coil. The binding of Hg(II) was also shown to substantially increase the stability of the helical aggregates.

Amino Acid Sequence↗

PubMed1998

Effect of charge interactions on the carboxylate vibrational stretching frequency in c-type cytochromes investigated by continuum electrostatic calculations and FTIR spectroscopy.

The FTIR spectra of the asymmetric carboxylate absorption region of three c-type cytochromes--namely horse heart, yeast and bonito cytochromes c--as well as continuum electrostatic calculations performed on their respective protein matrices, show that these combined methods can target specific protein regions and yield pertinent protein charge information that correlates with the observed spectral data. Deconvolution of the IR carboxylate stretch frequency region (1525-1675 cm-1) in the three cytochromes yield different v(oco)a distributions. In the case of the bonito cytochrome c carboxylates, two v(oco)a populations are clearly distinguishable in the deconvoluted spectra--which is not the case for the more complex v(oco)a deconvolutions of the other two cytochromes. The frequency distributions of the calculated potentials are consistent with the experimental observations and we conclude that the IR carboxylate absorption in proteins can be modified by the electrostatic environment.

Animals↗

PubMed1998

Calculation of electron transfer reorganization energies using the finite difference Poisson-Boltzmann model.

A description is given of a method to calculate the electron transfer reorganization energy (lambda) in proteins using the linear or nonlinear Poisson-Boltzmann (PB) equation. Finite difference solutions to the linear PB equation are then used to calculate lambda for intramolecular electron transfer reactions in the photosynthetic reaction center from Rhodopseudomonas viridis and the ruthenated heme proteins cytochrome c, myoglobin, and cytochrome b and for intermolecular electron transfer between two cytochrome c molecules. The overall agreement with experiment is good considering both the experimental and computational difficulties in estimating lambda. The calculations show that acceptor/donor separation and position of the cofactors with respect to the protein/solvent boundary are equally important and, along with the overall polarizability of the protein, are the major determinants of lambda. In agreement with previous studies, the calculations show that the protein provides a low reorganization environment for electron transfer. Agreement with experiment is best if the protein polarizability is modeled with a low (<8) average effective dielectric constant. The effect of buried waters on the reorganization energy of the photosynthetic reaction center was examined and found to make a contribution ranging from 0.05 eV to 0.27 eV, depending on the donor/acceptor pair.

2,2'-Dipyridyl↗

PubMed1997

Empirical free energy calculations: a blind test and further improvements to the method.

Empirical Gibbs functions estimate free energies of non-covalent reactions (deltaG) from atomic coordinates of reaction products (e.g. antibody-antigen complexes). The function previously developed by us has four terms that quantify the effects of hydrophobic, electrostatic and entropy changes (conformational, association) upon complexation. The function was used to calculate delta deltaG of ten lysozyme mutants affecting the stability of the HyHEL-10 antibody-lysozyme complex. The mutants were computer-modeled from the X-ray structure of the wild-type, and free energy calculations produced a correlation coefficient of 0.5 with the experimental delta deltaG data (average absolute error +/-3 kcal). The following changes were then introduced into the Gibbs function: (1) the hydrophobic force was made proportional to the molecular surface, as calculated by the GEPOL93 algorithm, with the scaling constant of 70 cal/mol/A2; (2) calculation of the electrostatics of binding was carried out by the finite difference Poisson-Boltzmann algorithm, which employed uniform grid charging, dielectric boundary smoothing and charge anti-aliasing; and (3) side-chain conformational entropy was estimated from the CONGEN sampling of torsional degrees of freedom. In the new calculations, correlation with experimental data improved to 0.6 or 0.8 if a single outlying mutant, K96M, was neglected. Analysis of the errors remaining in our calculations indicated that molecular mechanics-based modeling of the mutants, rather than the form of our amended Gibbs function, was the main factor limiting the accuracy of the free energy estimates.

Animals↗

PubMed1997

Energetics of cyclic dipeptide crystal packing and solvation.

Calculations of the thermodynamics of transfer of the cyclic alanine-alanine (cAA) and glycine-glycine (cGG) dipeptides between the gas, water, and crystal phases were carried out using a combination of molecular mechanics, normal mode analysis, and continuum electrostatics. The experimental gas-to-water solvation free energy and the enthalpy of gas-to-crystal transfer of cGG are accurately reproduced by the calculations. The enthalpies of cGG and cAA crystal-to-water transfer are close to the experimental values. A combination of experimental data and normal mode analysis of cGG provides an accurate estimate of the association entropy penalty (loss of rational and translational entropy and gain in vibrational entropy) for "binding" in the crystalline phase of -14.1 cal/mol/k. This is a smaller number than most previous theoretical estimates, but it is similar to previous experimental estimates. Calculated entropies of the crystal phase underestimate the experimental entropy by about 15 cal/mol/k because of neglect of long-range lattice motions. Comparison of the intermolecular interactions in the crystals of cGG and cAA provides a possible explanation of the puzzling decrease in enthalpy, with increasing hydrophobicity seen previously for both cyclic dipeptide dissolution and protein unfolding. This decrease arises from a favorable long-range electrostatic interaction between dipeptide molecules in the crystals, which is attenuated by the more hydrophobic side chains.

Crystallization↗

PubMed1997

Entropy in protein folding and in protein-protein interactions.

The reduction of conformational entropy is a major barrier that has to be overcome in protein folding and binding. Changes in solvent entropy are also a major factor. Recent advances include clarification of the fundamental issues concerning the separation of entropy into components, the treatment of association entropy in binding, and the role of size and shape effects in solvation entropy. Advances in the application of entropy calculations include an emerging consensus for estimates of backbone and sidechain entropy loss in protein folding via use of numerically intensive methods for sampling, and use of the expanding protein-structure database.

Algorithms↗

PubMed1996

On the decomposition of free energies.

The decomposition of free energies and entropies into components has recently been discussed within the framework of the free energy perturbation (FEP) and thermodynamic integration (TI) methods. In FEP, the cumulant expansion of the excess free energy contains coupling terms in second and higher orders. It is shown here that this expansion can be expressed in terms of temperature derivatives of the mean energy, suggesting a natural decomposition of the free energy into components corresponding to each term in the Hamiltonian. This result is derived in such a way that it establishes the equivalence to a particular form of component analysis based on TI in which all terms in the interaction energy are turned on simultaneously using 1/kT as the coupling parameter.

Computer Simulation↗

PubMed1996

Electrostatic interactions in hirudin-thrombin binding.

Hirudin is a good anticoagulant owing to potent inhibition of the serine protease thrombin. An aspartate- and glutamate-rich portion of hirudin plays an important part in its tight binding to thrombin through a ladder of salt bridges, and these residues have previously been mutated to asparagine or glutamine. Detailed calculations of the electrostatic contribution to changes in binding from these mutations have been performed using the finite-difference Poisson-Boltzmann method which include charge--charge interactions, solvation interactions, the residual electrostatic interaction of mutant residues, pKa shifts, and ionic strength. Single mutant effects on binding energy were close to experimental values, except for the D55N mutant whose effect is overestimated, perhaps because of displacement of a bound chloride ion from the site where it binds. Multiple mutation values were generally overestimated. The effect of pKa shifts upon the binding is significant for one hirudin residue E58, but this appears to be due to a poor salt bridge with thrombin caused by crystal contacts. Electrostatic interaction between the acidic residues is unfavorable. However, analysis of experimental multiple mutation/single mutation data shows apparently negative interactions between these residues, from which it is concluded that structural changes can occur in the complex to relieve an unfavorable interaction when more than one acidic residue is mutated. In all cases, there is a loss in stability of the complex from mutations due to loss of favorable charge--charge interactions with thrombin, but this is largely compensated for by reduced unfavorable desolvation interactions, and by residual polar interactions in the Asn/Gln mutants.

Anticoagulants↗

PubMed1995

Decomposition of interaction free energies in proteins and other complex systems.

A recent analysis of Mark and van Gunsteren has questioned the validity of separating different free energy components in proteins, or indeed in any complex system. The separability of free energy terms is re-examined from both a theoretical and a numerical perspective. Using a power series expansion of the free energy, it is found that the leading terms are free energy components that arise from individual contributions to the Hamiltonian ("in situ" free energies). The energetic part of an in situ free energy component is given by the ensemble average of the corresponding Hamiltonian component, while the leading term in the entropic part, which was missing in the analysis of Mark and van Gunsteren, is given by the mean square fluctuation. In addition there are correlations between fluctuations in each Hamiltonian component, which give rise to a coupling, or correlation entropy. A simple system, whose configurational degrees of freedom can be completely sampled, was examined in order to determine the relative sizes of these different contributions to the free energy. Under certain conditions, the change in system free energy observed when a particular component of the Hamiltonian is removed or altered is well approximated by the change in the in situ free energy of that component. In practical terms, this allows one in these cases to separate out different free energy contributions.

Energy Transfer↗

PubMed1995

Salt effects on polyelectrolyte-ligand binding: comparison of Poisson-Boltzmann, and limiting law/counterion binding models.

The theory for salt dependence of the free energy, entropy, and enthalpy of a polyelectrolyte in the PB (PB) model is extended to treat the nonspecific salt dependence of polyelectrolyte-ligand binding reactions. The salt dependence of the binding constant (K) is given by the difference in osmotic pressure terms between the reactants and products. For simple 1-1 salts it is shown that this treatment is equivalent to the general preferential interaction model for the salt dependence of binding [C. Anderson and M. Record (1993) Journal of Physical Chemistry, Vol. 97, pp. 7116-7126]. The salt dependence, entropy, and enthalpy are compared for the PB model and one specific form of the preferential interaction coefficient model that uses counterion condensation/limiting law (LL) behavior. The PB and LL models are applied to three ligand-polyelectrolyte systems with the same net ligand charge: a model sphere-cylinder binding reaction, a drug-DNA binding reaction, and a protein-DNA binding reaction. For the small ligands both the PB and limiting law models give (In K vs. In[salt]) slopes close in magnitude to the net ligand charge. However, the enthalpy/entropy breakdown of the salt dependence is quite different. In the PB model there are considerable contributions from electrostatic enthalpy and dielectric (water reorientation) entropy, compared to the predominant ion cratic (release) entropy in the limiting law model. The relative contributions of these three terms in the PB model depends on the ligand: For the protein, ion release entropy is the smallest contribution to the salt dependence of binding. The effect of three approximations made in the LL model is examined: These approximations are (1) the ligand behaves ideally, (2) the preferential interaction coefficient of the polyelectrolyte is unchanged upon ligand binding, and (3) the polyelectrolyte preferential interaction coefficient is given by the limiting law/counterion-condensation value. Analysis of the PB model shows that assumptions 2 and 3 break down at finite salt concentrations. For the small ligands the effects on the slope cancel, however, giving net slopes that are similar in the PB and LL models, but with a different entropy/enthalpy breakdown. For the protein ligand the errors from assumptions 2 and 3 in the LL models do not cancel.

Computer Graphics↗

PubMed1995

Salt effects on nucleic acids.

Salt-dependent electrostatic effects are a major factor in determining the stability, structure, reactivity, and binding behavior of nucleic acids. Increasingly detailed theoretical methods, especially those based on Monte Carlo and Poisson-Boltzmann methodologies, combined with powerful computational algorithms are being used to examine how the shape, charge distribution and dielectric properties of the molecules affect the ion distribution in the surrounding aqueous solution, and how they play a role in ligand binding, structural transitions and other biologically important reactions. These studies indicate that inclusion of detailed structural information about the nucleic acid and its ligands is crucial for improving models of nucleic acid electrostatics, and that better treatment of the ion atmosphere and dielectric effects is also of major importance.

Algorithms↗

PubMed1995

Calculations of the electrostatic potential adjacent to model phospholipid bilayers.

We used the nonlinear Poisson-Boltzmann equation to calculate electrostatic potentials in the aqueous phase adjacent to model phospholipid bilayers containing mixtures of zwitterionic lipids (phosphatidylcholine) and acidic lipids (phosphatidylserine or phosphatidylglycerol). The aqueous phase (relative permittivity, epsilon r = 80) contains 0.1 M monovalent salt. When the bilayers contain < 11% acidic lipid, the -25 mV equipotential surfaces are discrete domes centered over the negatively charged lipids and are approximately twice the value calculated using Debye-Hückel theory. When the bilayers contain > 25% acidic lipid, the -25 mV equipotential profiles are essentially flat and agree well with the values calculated using Gouy-Chapman theory. When the bilayers contain 100% acidic lipid, all of the equipotential surfaces are flat and agree with Gouy-Chapman predictions (including the -100 mV surface, which is located only 1 A from the outermost atoms). Even our model bilayers are not simple systems: the charge on each lipid is distributed over several atoms, these partial charges are non-coplanar, there is a 2 A ion-exclusion region (epsilon r = 80) adjacent to the polar headgroups, and the molecular surface is rough. We investigated the effect of these four factors using smooth (or bumpy) epsilon r = 2 slabs with embedded point charges: these factors had only minor effects on the potential in the aqueous phase.

Biophysical Phenomena↗

PubMed1994

Salt effects on ligand-DNA binding. Minor groove binding antibiotics.

Salt dependent electrostatic effects play a central role in intermolecular interactions involving nucleic acids. In this paper, the finite-difference solution to the nonlinear Poisson-Boltzmann (NLPB) equation is used to evaluate the salt dependent contribution to the electrostatic binding free energy of the minor groove binding antibiotics DAPI, Hoechst 33258 and netropsin to DNA using detailed molecular structures of the complexes. For each of these systems, a treatment based on the NLPB equation accurately describes the variation of the experimentally observed binding constant with bulk salt concentration. A solvation formalism is developed in which salt effects are described in terms of three free energy contributions: the electrostatic ion-molecule interaction free energy, delta delta G degrees im; the electrostatic ion-ion interaction free energy, delta delta G degrees ii; and the entropic ion organization free energy, delta delta G degrees org. The electrostatic terms, delta delta G degrees im and delta delta G degrees ii, have both enthalpic and entropic components, while the term delta delta G degrees org is purely a cratic entropy. Each of these terms depends significantly on salt dependent changes in the counterion and coion concentrations around the DNA. In each of the systems studied, univalent ions substantially destabilize charged ligand-DNA complexes at physiological salt concentrations. This effect involves a salt dependent redistribution of counterions near the DNA. The free energy associated with the redistribution of counterions upon binding is dominated by the unfavorable change in the electrostatic ion-molecule interactions, delta delta G degrees im, rather than the change in the cratic entropy of ion organization, delta delta G degrees org. In addition, the observed slope of the salt dependence of the free energy is determined by electrostatic ion-molecule and ion-ion interactions as well as the cratic entropy of ion release. These findings are in contrast to models in which the cratic entropy of counterion release drives binding.

Anti-Bacterial Agents↗