Biomedical subjects

D Asthagiri

Publications and source records attributed to D Asthagiri.

14 recordsLinked to original sources

PubMed2006

Role of fluctuations in a snug-fit mechanism of KcsA channel selectivity.

The thermodynamic exclusion of Na+ relative to K+ in potassium channels is examined by calculating the distribution of binding energies for Na+ and K+ in a model of the selectivity filter of the KcsA potassium channel. These distributions are observed to take a surprisingly simple form: Gaussian with a slight positive skewness that is insignificant in the present context. Complications that might be anticipated from these distributions are not problematic here. Na+ occupies the filter with a mean binding energy substantially lower than that of K+. The difference is comparable to the difference in hydration free energies of Na+ and K+ in bulk aqueous solution. Thus, the average energies of binding to the filter do not discriminate Na+ from K+ when measured from a baseline of the difference in bulk hydration free energies. The strong binding of Na+ constricts the filter, consistent with a negative partial molar volume of Na+ in water in contrast with a positive partial molar volume of K+ in water. Discrimination in favor of K+)can be attributed to the scarcity of favorable binding configurations for Na+ compared to K+. That relative scarcity is quantified as enhanced binding energy fluctuations, which reflects both the energetically stronger binding of Na+ and the constriction of the filter induced by Na+ binding.

Journal Article↗

PubMed2006

An analysis of molecular packing and chemical association in liquid water using quasichemical theory.

We calculate the hydration free energy of liquid TIP3P water at 298 K and 1 bar using a quasi-chemical theory framework in which interactions between a distinguished water molecule and the surrounding water molecules are partitioned into chemical associations with proximal (inner-shell) waters and classical electrostatic-dispersion interactions with the remaining (outer-shell) waters. The calculated free energy is found to be independent of this partitioning, as expected, and in excellent agreement with values derived from the literature. An analysis of the spatial distribution of inner-shell water molecules as a function of the inner-shell volume reveals that water molecules are preferentially excluded from the interior of large volumes as the occupancy number decreases. The driving force for water exclusion is formulated in terms of a free energy for rearranging inner-shell water molecules under the influence of the field exerted by outer-shell waters in order to accommodate one water molecule at the center. The results indicate a balance between chemical association and molecular packing in liquid water that becomes increasingly important as the inner-shell volume grows in size.

Chemistry, Physical↗

PubMed2005

Light-scattering studies of protein solutions: role of hydration in weak protein-protein interactions.

We model the hydration contribution to short-range electrostatic/dispersion protein interactions embodied in the osmotic second virial coefficient, B(2), by adopting a quasi-chemical description in which water molecules associated with the protein are identified through explicit molecular dynamics simulations. These water molecules reduce the surface complementarity of highly favorable short-range interactions, and therefore can play an important role in mediating protein-protein interactions. Here we examine this quasi-chemical view of hydration by predicting the interaction part of B(2) and comparing our results with those derived from light-scattering measurements of B(2) for staphylococcal nuclease, lysozyme, and chymotrypsinogen at 25 degrees C as a function of solution pH and ionic strength. We find that short-range protein interactions are influenced by water molecules strongly associated with a relatively small fraction of the protein surface. However, the effect of these strongly associated water molecules on the surface complementarity of short-range protein interactions is significant, and must be taken into account for an accurate description of B(2). We also observe remarkably similar hydration behavior for these proteins despite substantial differences in their three-dimensional structures and spatial charge distributions, suggesting a general characterization of protein hydration.

Biophysics↗

PubMed2005

Ab initio molecular dynamics and quasichemical study of H+(aq).

The excess proton in water, H(+)(aq), plays a fundamental role in aqueous solution chemistry. Its solution thermodynamic properties are essential to molecular descriptions of that chemistry and for validation of dynamical calculations. Within the quasichemical theory of solutions those thermodynamic properties are conditional on recognizing underlying solution structures. The quasichemical treatment identifies H(3)O(+) and H(2)O(5)(+) as natural inner-shell complexes, corresponding to the cases of n = 1, 2 water molecule ligands, respectively, of a distinguished H(+) ion. A quantum-mechanical treatment of the inner-shell complex with both a dielectric continuum and a classical molecular dynamics treatment of the outer-shell contribution identifies the latter case (the Zundel complex) as the more numerous species. Ab initio molecular dynamics simulations, with two different electron density functionals, suggest a preponderance of Zundel-like structures, but a symmetrical ideal Zundel cation is not observed.

Biophysical Phenomena↗

PubMed2005

A consistent experimental and modeling approach to light-scattering studies of protein-protein interactions in solution.

The osmotic second virial coefficient, B(2), obtained by light scattering from protein solutions has two principal components: the Donnan contribution and a contribution due to protein-protein interactions in the limit of infinite dilution. The Donnan contribution accounts for electroneutrality in a multicomponent solution of (poly)electrolytes. The importance of distinguishing this ideal contribution to B(2) is emphasized, thereby allowing us to model the interaction part of B(2) by molecular computations. The model for protein-protein interactions that we use here extends earlier work (Neal et al., 1998) by accounting for long-range electrostatic interactions and the specific hydration of the protein by strongly associated water molecules. Our model predictions are compared with measurements of B(2) for lysozyme at 25 degrees C over pH from 5.0 to 9.0, and 7-60 mM ionic strength. We find that B(2) is positive at all solution conditions and decreases with increasing ionic strength, as expected, whereas the interaction part of B(2) is negative at all conditions and becomes progressively less negative with increasing ionic strength. Although long-range electrostatic interactions dominate this contribution, particularly at low ionic strength, short-range electrostatic/dispersion interactions with specific hydration are essential for an accurate description of B(2) derived from experiment.

Animals↗

PubMed2004

On the role of the conserved aspartate in the hydrolysis of the phosphocysteine intermediate of the low molecular weight tyrosine phosphatase.

The usual rate-determining step in the catalytic mechanism of the low molecular weight tyrosine phosphatases involves the hydrolysis of a phosphocysteine intermediate. To explain this hydrolysis, general base-catalyzed attack of water by the anion of a conserved aspartic acid has sometimes been invoked. However, experimental measurements of solvent deuterium kinetic isotope effects for this enzyme do not reveal a rate-limiting proton transfer accompanying dephosphorylation. Moreover, base activation of water is difficult to reconcile with the known gas-phase proton affinities and solution phase pK(a)'s of aspartic acid and water. Alternatively, hydrolysis could proceed by a direct nucleophilic attack by a water molecule. To understand the hydrolysis mechanism, we have used high-level density functional methods of quantum chemistry combined with continuum electrostatics models of the protein and the solvent. Our calculations do not support a catalytic activation of water by the aspartate. Instead, they indicate that the water oxygen directly attacks the phosphorus, with the aspartate residue acting as a H-bond acceptor. In the transition state, the water protons are still bound to the oxygen. Beyond the transition state, the barrier to proton transfer to the base is greatly diminished; the aspartate can abstract a proton only after the transition state, a result consistent with experimental solvent isotope effects for this enzyme and with established precedents for phosphomonoester hydrolysis.

Animals↗

PubMed2004

Hydration and mobility of HO-(aq).

The hydroxide anion plays an essential role in many chemical and biochemical reactions. But a molecular-scale description of its hydration state, and hence also its transport, in water is currently controversial. The statistical mechanical quasichemical theory of solutions suggests that HO.[H2O]3(-) is the predominant species in the aqueous phase under standard conditions. This result agrees with recent spectroscopic studies on hydroxide water clusters and with the available thermodynamic hydration free energies. In contrast, a recent ab initio molecular dynamics simulation has suggested that HO.[H2O]4(-) is the only dominant aqueous solution species. We apply adiabatic ab initio molecular dynamics simulations and find good agreement with both the quasichemical theoretical predictions and experimental results. The present results suggest a picture that is simpler, more traditional, but with additional subtlety. These coordination structures are labile but the tricoordinate species is the prominent case. This conclusion is unaltered with changes in the electronic density functional. No evidence is found for rate-determining activated interconversion of a HO.[H2O]4(-) trap structure to HO.[H2O]3(-) mediating hydroxide transport. The view of HO- diffusion as the hopping of a proton hole has substantial validity, the rate depending largely on the dynamic disorder of the water hydrogen-bond network.

Journal Article↗

PubMed2004

Hydration structure and free energy of biomolecularly specific aqueous dications, including Zn2+ and first transition row metals.

The hydration of some of the alkaline earth divalent metal cations and first row transition metal cations is considered within the quasi-chemical theory of solutions. Quantum chemical calculations provide information on the chemically important interactions between the ion and its first-shell water molecules. A dielectric continuum model supplies the outer-shell contribution. The theory then provides the framework to mesh these quantities together. The agreement between the calculated and experimental quantities is good. For the transition metal cations, it is seen that the ligand field contributions play an important role in the physics of hydration. Removing these bonding contributions from the computed hydration free energy results in a linear decrease in the hydration free energy along the period. It is precisely such effects that molecular mechanics force fields have not captured. The implications and extensions of this study to metal atoms in proteins are suggested.

Calcium↗

PubMed2003

Free energy of liquid water on the basis of quasichemical theory and ab initio molecular dynamics.

We use ab initio molecular dynamics as a basis for quasichemical theory evaluation of the free energy of water near conventional liquid thermodynamic states. The Perdew-Wang-91 (PW91), Perdew-Burke-Ernzerhof (PBE), and revised PBE (rPBE) functionals are employed. The oxygen radial density distribution using the rPBE functional is in reasonable agreement with current experiments, whereas the PW91 and PBE functionals predict a more structured oxygen radial density distribution. The diffusion coefficient with the rPBE functional is in reasonable accord with experiments. Using a maximum entropy procedure, we obtain x0 from the coordination number distribution xn for oxygen atoms having n neighbors. Likewise, we obtain p0 from pn, the probability of observing cavities of specified radius containing n water molecules. The probability x0 is a measure of the local chemical interactions and is central to the quasichemical theory of solutions. The probability p0, central to the theory of liquids, is a measure of the free energy required to open cavities of defined sizes in the solvent. Using these values and a reasonable model for electrostatic and dispersion effects, the hydration free energy of water in water at 314 K is calculated to be -5.1 kcal/mole with the rPBE functional, in encouraging agreement with the experimental value of -6.1 kcal/mole.

Computer Simulation↗

PubMed2003

Hydration of krypton and consideration of clathrate models of hydrophobic effects from the perspective of quasi-chemical theory.

Ab initio molecular dynamics (AIMD) results on a krypton-water liquid solution are presented and compared to recent XAFS results for the radial hydration structure for a Kr atom in liquid water solution. Though these AIMD calculations have important limitations of scale, the comparisons with the liquid solution results are satisfactory and significantly different from the radial distributions extracted from the data on the solid Kr/H(2)O clathrate hydrate phase. The calculations also produce the coordination number distribution that can be examined for metastable coordination structures suggesting possibilities for clathrate-like organization; none are seen in these results. Clathrate pictures of hydrophobic hydration are discussed, as is the quasi-chemical theory that should provide a basis for clathrate pictures. Outer shell contributions are discussed and estimated; they are positive and larger than the positive experimental hydration free energy of Kr(aq), implying that inner shell contributions must be negative and of comparable size. Clathrate-like inner shell hydration structures on a Kr atom solute are obtained for some, but not all, of the coordination number cases observed in the simulation. The structures found have a delicate stability. Inner shell coordination structures extracted from the simulation of the liquid, and then subjected to quantum chemical optimization, always decomposed. Interactions with the outer shell material are decisive in stabilizing coordination structures observed in liquid solution and in clathrate phases. The primitive quasi-chemical estimate that uses a dielectric model for the influence of the outer shell material on the inner shell equilibria gives a contribution to hydration free energy that is positive and larger than the experimental hydration free energy. The 'what are we to tell students' question about hydrophobic hydration, often answered with structural clathrate pictures, is then considered; we propose an alternative answer that is consistent with successful molecular theories of hydrophobic effects and based upon distinctive observable properties of liquid water. Considerations of parsimony, for instance Ockham's razor, then suggest that additional structural hypotheses in response to 'what are we to tell students' are not required at this stage.

Hydrophobic and Hydrophilic Interactions↗

PubMed2002

Continuum and atomistic modeling of ion partitioning into a peptide nanotube.

Continuum and atomistic descriptions of the partitioning of ions into a self-assembled (D,L)-octapeptide nanotube, cyclo[-(L-Ala-D-Ala)(4)-], are presented. Perturbation free energy calculations, including Ewald electrostatics, are used to estimate the electrostatic component of the excess free energy of charging Li(+), Na(+), Rb(+), and Cl(minus sign) ions inside the nanotube. The radial density and orientational distribution of water around the ion is calculated for the ion at two different positions inside the tube; it is seen that the calculated distributions are sensitive to the location of the ions. Two different continuum electrostatic models are formulated to describe the ion solvation inside the nanotube. When enhanced orientational structuring of water dipoles is evidenced, explicitly including the first solvation shell as part of the low dielectric nanotube environment provides good agreement with molecular dynamics simulations. When water orientational structuring is as in the reference bulk solvent, we find that treating the first shell water explicitly or as a high dielectric continuum leads to similar results. These results are discussed, and their importance for continuum electrostatic modeling of ion channels are highlighted.

Biophysical Phenomena↗

PubMed2000

Measured and calculated effects of mutations in bacteriophage T4 lysozyme on interactions in solution.

Understanding the molecular determinants of protein interactions in solution has fundamental implications for understanding protein solution thermodynamics and, hence, processes as diverse as separations performance and cellular self-organization. Our earlier theoretical calculations indicate that the protein-protein interactions are dominated by a small number of configurations in which highly complementary surface regions are apposed, rather than by the overall colloidal interactions. To examine this paradigm more explicitly, we investigated the effects of protein structural modifications on protein-protein interactions. Experimental measurements are presented of B(22)(') values of a set of mutants of Ser44 in bacteriophage T4 lysozyme. Effects are seen with both charged and uncharged substitutions. The results with the charged substitutions follow the expected trends, whereas those with the uncharged substitutions may be explained by the impact of the mutations on the local protein geometry, which directly affects the complementarity of protein interactions. These effects are also captured well by molecular calculations that account for the mutations. The interaction energetics between protein pairs could provide information on the propensity for adventitious interactions, which can have important implications for separations and for normal and pathological self-assembly. Thus, protein structural data implicit in genomic information, coupled with appropriate calculational and experimental tools, can ultimately provide insights into protein interactions in vivo and in bioprocessing.

Bacteriophage T4↗

PubMed1999

Calculation of short-range interactions between proteins.

Macromolecular association is an integral component of numerous cellular and technologically relevant processes. Most molecular theories of such association neglect the explicit solvent structure and rely on continuum concepts such as surface energies for calculating short-range interactions. We present a new such method for calculating the non-electrostatic component of the interaction-free energy, based on formalisms for calculating dispersion interactions between macromolecules. The interactions are separated into a short-ranged component that is treated atomistically, and a longer range component that is treated within the continuum Lifshitz-Hamaker approach. This description avoids the singularities inherent in the continuum dispersion formulation, and its effectiveness in characterizing the shape complementarity between interacting surfaces is shown to be comparable to that of surface area-based methods of similar parametric complexity. An advantage of the new method is that it allows facile calculation of the interaction free energy as a function of intermolecular separation, including steric effects; this makes it suitable for use in simulations of protein solutions.

Journal Article↗

PubMed1998

Molecular origins of osmotic second virial coefficients of proteins.

The thermodynamic properties of protein solutions are determined by the molecular interactions involving both solvent and solute molecules. A quantitative understanding of the relationship would facilitate more systematic procedures for manipulating the properties in a process environment. In this work the molecular basis for the osmotic second virial coefficient, B22, is studied; osmotic effects are critical in membrane transport, and the value of B22 has also been shown to correlate with protein crystallization behavior. The calculations here account for steric, electrostatic, and short-range interactions, with the structural and functional anisotropy of the protein molecules explicitly accounted for. The orientational dependence of the protein interactions is seen to have a pronounced effect on the calculations; in particular, the relatively few protein-protein configurations in which the apposing surfaces display geometric complementarity contribute disproportionately strongly to B22. The importance of electrostatic interactions is also amplified in these high-complementarity configurations. The significance of molecular recognition in determining B22 can explain the correlation with crystallization behavior, and it suggests that alteration of local molecular geometry can help in manipulating protein solution behavior. The results also have implications for the role of protein interactions in biological self-organization.

Animals↗