Biomedical subjects

Shekhar Garde

Publications and source records attributed to Shekhar Garde.

16 recordsLinked to original sources

PubMed2007

Effects of lengthscales and attractions on the collapse of hydrophobic polymers in water.

We present results from extensive molecular dynamics simulations of collapse transitions of hydrophobic polymers in explicit water focused on understanding effects of lengthscale of the hydrophobic surface and of attractive interactions on folding. Hydrophobic polymers display parabolic, protein-like, temperature-dependent free energy of unfolding. Folded states of small attractive polymers are marginally stable at 300 K and can be unfolded by heating or cooling. Increasing the lengthscale or decreasing the polymer-water attractions stabilizes folded states significantly, the former dominated by the hydration contribution. That hydration contribution can be described by the surface tension model, DeltaG = gamma(T)DeltaA, where the surface tension, gamma, is lengthscale-dependent and decreases monotonically with temperature. The resulting variation of the hydration entropy with polymer lengthscale is consistent with theoretical predictions of Huang and Chandler [Huang DM, Chandler D (2000) Proc Natl Acad Sci USA 97:] that explain the blurring of entropy convergence observed in protein folding thermodynamics. Analysis of water structure shows that the polymer-water hydrophobic interface is soft and weakly dewetted, and is characterized by enhanced interfacial density fluctuations. Formation of this interface, which induces polymer folding, is strongly opposed by enthalpy and favored by entropy, similar to the vapor-liquid interface.

Hydrophobic and Hydrophilic Interactions↗

PubMed2006

Mechanism for intein C-terminal cleavage: a proposal from quantum mechanical calculations.

Inteins are autocatalytic protein cleavage and splicing elements. A cysteine to alanine mutation at the N-terminal of inteins inhibits splicing and isolates the C-terminal cleavage reaction. Experiments indicate an enhanced C-terminal cleavage reaction rate upon decreasing the solution pH for the cleavage mutant, which cannot be explained by the existing mechanistic framework. We use intein crystal structure data and the information about conserved amino acids to perform semiempirical PM3 calculations followed by high-level density functional theory calculations in both gas phase and implicit solvent environments. Based on these calculations, we propose a detailed "low pH" mechanism for intein C-terminal cleavage. Water plays an important role in the proposed reaction mechanism, acting as an acid as well as a base. The protonation of the scissile peptide bond nitrogen by a hydronium ion is an important first step in the reaction. That step is followed by the attack of the C-terminal asparagine side chain on its carbonyl carbon, causing succinimide formation and simultaneous peptide bond cleavage. The computed reaction energy barrier in the gas phase is approximately 33 kcal/mol and reduces to approximately 25 kcal/mol in solution, close to the 21 kcal/mol experimentally observed at pH 6.0. This mechanism is consistent with the observed increase in C-terminal cleavage activity at low pH for the cleavage mutant of the Mycobacterium tuberculosis RecA mini-intein.

Binding Sites↗

PubMed2006

Tail ordering due to headgroup hydrogen bonding interactions in surfactant monolayers at the water-oil interface.

Interactions between surfactants, and the resultant ordering of surfactant assemblies, can be tuned by the appropriate choice of head- and tailgroups. Detailed studies of the ordering of monolayers of long-chain n-alkanoic and n-alkanol monolayers at the water-vapor interface have demonstrated that rigid-rod all-trans ordering of the tailgroups is maintained upon replacing the alcohol with a carboxylic acid headgroup. In contrast, at the water-hexane liquid-liquid interface, we demonstrate that substitution of the -CH(2)OH with the -COOH headgroup produces a major conformational change of the tailgroup from disordered to ordered. This is demonstrated by the electron density profiles of triacontanol (CH(3)(CH(2))(29)OH) and triacontanoic acid (CH(3)(CH(2))(28)COOH) monolayers at the water-hexane interface, as determined by X-ray reflectivity measurements. Molecular dynamics simulations illustrate the presence of hydrogen bonding between the triacontanoic acid headgroups that is likely responsible for the tail ordering. A simple free energy illustrates the interplay between the attractive hydrogen bonding and the ordering of the tailgroup.

Letter↗

PubMed2006

Quantifying the protein core flexibility through analysis of cavity formation.

We present an extensive analysis of cavity statistics in the interior of three different proteins, in liquid n-hexane, and in water performed using molecular-dynamics simulations. The heterogeneity of packing density over atomic length scales in different parts of proteins is evident in the wide range of values observed for the average cavity size, the probability of cavity formation, and the corresponding free energy of hard-sphere insertion. More interestingly, however, the distribution of cavity sizes observed at various points in the protein interior is surprisingly homogeneous in width. That width is significantly smaller than that measured for similar distributions in liquid n-hexane or water, indicating that protein interior is much less flexible than liquid hexane. The width of the cavity size distribution correlates well with the experimental isothermal compressibility data for liquids and proteins. An analysis of cavity statistics thus provides an efficient method to quantify local properties, such as packing, stiffness, or compressibility in heterogeneous condensed media.

Amyloid beta-Protein Precursor↗

PubMed2005

Direct determination of phase behavior of square-well fluids.

We have combined Gibbs ensemble Monte Carlo simulations with the aggregation volume-biased method in conjunction with the Gibbs-Duhem method to provide the first direct estimates for the vapor-solid, vapor-liquid, and liquid-solid phase coexistences of square-well fluids with three different ranges of attraction. Our results are consistent with the previous simulations and verify the notion that the vapor-liquid coexistence behavior becomes metastable for cases where the attraction well becomes smaller than 1.25 times the particle diameter. In these cases no triple point is found.

Journal Article↗

PubMed2005

Hydrophobic hydration from small to large lengthscales: Understanding and manipulating the crossover.

Small and large hydrophobic solutes exhibit remarkably different hydration thermodynamics. Small solutes are accommodated in water with minor perturbations to water structure, and their hydration is captured accurately by theories that describe density fluctuations in pure water. In contrast, hydration of large solutes is accompanied by dewetting of their surfaces and requires a macroscopic thermodynamic description. A unified theoretical description of these lengthscale dependencies was presented by Lum, Chandler, and Weeks [(1999) J. Phys. Chem. B 103, 4570-4577]. Here, we use molecular simulations to study lengthscale-dependent hydrophobic hydration under various thermodynamic conditions. We show that the hydration of small and large solutes displays disparate dependencies on thermodynamic variables, including pressure, temperature, and additive concentration. Understanding these dependencies allows manipulation of the small-to-large crossover lengthscale, which is nanoscopic under ambient conditions. Specifically, applying hydrostatic tension or adding ethanol decreases the crossover length to molecular sizes, making it accessible to atomistic simulations. With detailed temperature-dependent studies, we further demonstrate that hydration thermodynamics changes gradually from entropic to enthalpic near the crossover. The nanoscopic lengthscale of the crossover and its sensitivity to thermodynamic variables imply that quantitative modeling of biomolecular self-assembly in aqueous solutions requires elements of both molecular and macroscopic hydration physics. We also show that the small-to-large crossover is directly related to the Egelstaff-Widom lengthscale, the product of surface tension and isothermal compressibility, which is another fundamental lengthscale in liquids.

Biophysical Phenomena↗

PubMed2005

Osmolyte trimethylamine-N-oxide does not affect the strength of hydrophobic interactions: origin of osmolyte compatibility.

Osmolytes are small organic solutes accumulated at high concentrations by cells/tissues in response to osmotic stress. Osmolytes increase thermodynamic stability of folded proteins and provide protection against denaturing stresses. The mechanism of osmolyte compatibility and osmolyte-induced stability has, therefore, attracted considerable attention in recent years. However, to our knowledge, no quantitative study of osmolyte effects on the strength of hydrophobic interactions has been reported. Here, we present a detailed molecular dynamics simulation study of the effect of the osmolyte trimethylamine-N-oxide (TMAO) on hydrophobic phenomena at molecular and nanoscopic length scales. Specifically, we investigate the effects of TMAO on the thermodynamics of hydrophobic hydration and interactions of small solutes as well as on the folding-unfolding conformational equilibrium of a hydrophobic polymer in water. The major conclusion of our study is that TMAO has almost no effect either on the thermodynamics of hydration of small nonpolar solutes or on the hydrophobic interactions at the pair and many-body level. We propose that this neutrality of TMAO toward hydrophobic interactions-one of the primary driving forces in protein folding-is at least partially responsible for making TMAO a "compatible" osmolyte. That is, TMAO can be tolerated at high concentrations in organisms without affecting nonspecific hydrophobic effects. Our study implies that protein stabilization by TMAO occurs through other mechanisms, such as unfavorable water-mediated interaction of TMAO with the protein backbone, as suggested by recent experimental studies. We complement the above calculations with analysis of TMAO hydration and changes in water structure in the presence of TMAO molecules. TMAO is an amphiphilic molecule containing both hydrophobic and hydrophilic parts. The precise balance of the effects of hydrophobic and hydrophilic segments of the molecule appears to explain the virtual noneffect of TMAO on the strength of hydrophobic interactions.

Computer Simulation↗

PubMed2005

Mesoscale model of polymer melt structure: self-consistent mapping of molecular correlations to coarse-grained potentials.

Development and application of coarse-graining methods to condensed phases of macromolecules is an active area of research. Multiscale modeling of polymeric systems using coarse-graining methods presents unique challenges. Here we apply a coarse-graining method that self-consistently maps structural correlations from detailed molecular dynamics (MD) simulations of alkane oligomers onto coarse-grained potentials using a combination of MD and inverse Monte Carlo methods. Once derived, the coarse-grained potentials allow computationally efficient sampling of ensemble of conformations of significantly longer polyethylene chains. Conformational properties derived from coarse-grained simulations are in excellent agreement with experiments. The level of coarse graining provides a control over the balance of computational efficiency and retention of chemical identity of the underlying polymeric system. Challenges to extension and application of this and similar structure-based coarse-graining methods to model dynamics and phase behavior in polymeric systems are briefly discussed.

Journal Article↗

PubMed2005

On the salt-induced stabilization of pair and many-body hydrophobic interactions.

Salting-out of hydrophobic solutes in aqueous salt solutions and their relevance to salt effects on biophysical phenomena are now well appreciated. Although salt effects on hydrophobic transfer have been well studied, to our knowledge, no quantitative molecular simulation study of salt-induced strengthening of hydrophobic interactions has yet been reported. Here we present quantitative characterization of salt-induced strengthening of hydrophobic interactions at the molecular and nanoscopic length scales through molecular dynamics simulations. Specifically, we quantify the effect of NaCl on the potential of mean force between molecular hydrophobic solutes (methanes) and on conformational equilibria of a 25-mer hydrophobic polymer that efficiently samples ensembles of compact and extended states in water. In both cases, we observe relative stabilization of compact conformations that is accompanied by a clear depletion of salt density (preferential exclusion) and a slight enhancement of water density (preferential hydration) in the solute vicinity. We show that the structural details of salt exclusion can be related to the salt-induced free energy changes using preferential interaction coefficients. We also test the applicability of surface-area-based models to describe the salt-induced free energy changes. These models provide a useful empirical description that can be used to predict the effects of salt on conformational equilibria of hydrophobic solutes. However, we find that the effective increase in the surface tension of the solute-aqueous solution interface depends on the type and concentration of salt as well as the length-scale (i.e., molecular vs nanoscopic) of the conformational change. These calculations underscore the utility of simulation studies to connect quantitatively structural details at the molecular level (described by preferential hydration/exclusion) to macroscopic solvation thermodynamics. The hydrophobic polymer also provides a useful model for studies of effect of thermodynamic variables (P, T, salt/additives) on many-body hydrophobic interactions at nanometer length scales.

Computer Simulation↗

PubMed2005

Thermal resistance of nanoscopic liquid-liquid interfaces: dependence on chemistry and molecular architecture.

Systems with nanoscopic features contain a high density of interfaces. Thermal transport in such systems can be governed by the resistance to heat transfer, the Kapitza resistance (RK), at the interface. Although soft interfaces, such as those between immiscible liquids or between a biomolecule and solvent, are ubiquitous, few studies of thermal transport at such interfaces have been reported. Here we characterize the interfacial conductance, 1/RK, of soft interfaces as a function of molecular architecture, chemistry, and the strength of cross-interfacial intermolecular interactions through detailed molecular dynamics simulations. The conductance of various interfaces studied here, for example, water-organic liquid, water-surfactant, surfactant-organic liquid, is relatively high (in the range of 65-370 MW/m2 K) compared to that for solid-liquid interfaces ( approximately 10 MW/m2 K). Interestingly, the dependence of interfacial conductance on the chemistry and molecular architecture cannot be explained solely in terms of either bulk property mismatch or the strength of intermolecular attraction between the two phases. The observed trends can be attributed to a combination of strong cross-interface intermolecular interactions and good thermal coupling via soft vibration modes present at liquid-liquid interfaces.

Journal Article↗

PubMed2004

Evaluation of selectivity changes in HIC systems using a preferential interaction based analysis.

It is well established that salt enhances the interaction between solutes (e.g., proteins, displacers) and the weak hydrophobic ligands in hydrophobic interaction chromatography (HIC) and that various salts (e.g., kosmotropes, chaotropes, and neutral) have different effects on protein retention. In this article, the solute affinity in kosmotropic, chaotropic, and neutral mobile phases are compared and the selectivity of solutes in the presence of these salts is examined. Since solute binding in HIC systems is driven by the release of water molecules, the total number of released water molecules in the presence of various types of salts was calculated using the preferential interaction theory. Chromatographic retention times and selectivity reversals of both proteins and displacers were found to be consistent with the total number of released water molecules. Finally, the solute surface hydrophobicity was also found to have a significant effect on its retention in HIC systems.

Algorithms↗

PubMed2004

Size dependent ion hydration, its asymmetry, and convergence to macroscopic behavior.

The packing and orientation of water molecules in the vicinity of solutes strongly influence the solute hydration thermodynamics in aqueous solutions. Here we study the charge density dependent hydration of a broad range of spherical monovalent ionic solutes (with solute diameters from approximately 0.4 nm to 1.7 nm) through molecular dynamics simulations in the simple point charge model of water. Consistent with previous experimental and theoretical studies, we observe a distinct asymmetry in the structure and thermodynamics of hydration of ions. In particular, the free energy of hydration of negative ions is more favorable than that of positive ions of the same size. This asymmetry persists over the entire range of solute sizes and cannot be captured by a continuum description of the solvent. The favorable hydration of negative ions arises primarily from the asymmetric charge distribution in the water molecule itself, and is reflected in (i) a small positive electrostatic potential at the center of a neutral solute, and (ii) clear structural (packing and orientation) differences in the hydration shell of positive and negative ions. While the asymmetry arising from the positive potential can be quantified in a straightforward manner, that arising from the structural differences in the fully charged states is difficult to quantify. The structural differences are highest for the small ions and diminish with increasing ion size, converging to hydrophobiclike hydration structure for the largest ions studied here. We discuss semiempirical measures following Latimer, Pitzer, and Slansky [J. Chem. Phys. 7, 108 (1939)] that account for these structural differences through a shift in the ion radius. We find that these two contributions account completely for the asymmetry of hydration of positive and negative ions over the entire range of ion sizes studied here. We also present preliminary calculations of the dependence of ion hydration asymmetry on the choice of water model that demonstrate its sensitivity to the details of ion-water interactions.

Journal Article↗

PubMed2004

Hydration of enzyme in nonaqueous media is consistent with solvent dependence of its activity.

Water plays an important role in enzyme structure and function in aqueous media. That role becomes even more important when one focuses on enzymes in low water media. Here we present results from molecular dynamics simulations of surfactant-solubilized subtilisin BPN' in three organic solvents (octane, tetrahydrofuran, and acetonitrile) and in pure water. Trajectories from simulations are analyzed with a focus on enzyme structure, flexibility, and the details of enzyme hydration. The overall enzyme and backbone structures, as well as individual residue flexibility, do not show significant differences between water and the three organic solvents over a timescale of several nanoseconds currently accessible to large-scale molecular dynamics simulations. The key factor that distinguishes molecular-level details in different media is the partitioning of hydration water between the enzyme and the bulk solvent. The enzyme surface and the active site region are well hydrated in aqueous medium, whereas with increasing polarity of the organic solvent (octane --> tetrahydrofuran --> acetonitrile) the hydration water is stripped from the enzyme surface. Water stripping is accompanied by the penetration of tetrahydrofuran and acetonitrile molecules into crevices on the enzyme surface and especially into the active site. More polar organic solvents (tetrahydrofuran and acetonitrile) replace mobile and weakly bound water molecules in the active site and leave primarily the tightly bound water in that region. In contrast, the lack of water stripping in octane allows efficient hydration of the active site uniformly by mobile and weakly bound water and some structural water similar to that in aqueous solution. These differences in the active site hydration are consistent with the inverse dependence of enzymatic activity on organic solvent polarity and indicate that the behavior of hydration water on the enzyme surface and in the active site is an important determinant of biological function especially in low water media.

Acetonitriles↗

PubMed2003

Osmotic water transport through carbon nanotube membranes.

We use molecular dynamics simulations to study osmotically driven transport of water molecules through hexagonally packed carbon nanotube membranes. Our simulation setup comprises two such semipermeable membranes separating compartments of pure water and salt solution. The osmotic force drives water flow from the pure-water to the salt-solution compartment. Monitoring the flow at molecular resolution reveals several distinct features of nanoscale flows. In particular, thermal fluctuations become significant at the nanoscopic length scales, and as a result, the flow is stochastic in nature. Further, the flow appears frictionless and is limited primarily by the barriers at the entry and exit of the nanotube pore. The observed flow rates are high (5.8 water molecules per nanosecond and nanotube), comparable to those through the transmembrane protein aquaporin-1, and are practically independent of the length of the nanotube, in contrast to predictions of macroscopic hydrodynamics. All of these distinct characteristics of nanoscopic water flow can be modeled quantitatively by a 1D continuous-time random walk. At long times, the pure-water compartment is drained, and the net flow of water is interrupted by the formation of structured solvation layers of water sandwiched between two nanotube membranes. Structural and thermodynamic aspects of confined water monolayers are studied.

Membranes, Artificial↗

PubMed2003

Helix propensities of short peptides: molecular dynamics versus bioinformatics.

Knowledge-based potential functions for protein structure prediction assume that the frequency of occurrence of a given structure or a contact in the protein database is a measure of its free energy. Here, we put this assumption to test by comparing the results obtained from sequence-structure cluster analysis with those obtained from long all-atom molecular dynamics simulations. Sixty-four eight-residue peptide sequences with varying degrees of similarity to the canonical sequence pattern for amphipathic helix were drawn from known protein structures, regardless of whether they were helical in the protein. Each was simulated using AMBER6.0 for at least 10 ns using explicit waters. The total simulation time was 1176 ns. The resulting trajectories were tested for reproducibility, and the helical content was measured. Natural peptides whose sequences matched the amphipathic helix motif with greater than 50% confidence were significantly more likely to form helix during the course of the simulation than peptides with lower confidence scores. The sequence pattern derived from the simulation data closely resembles the motif pattern derived from the database cluster analysis. The difficulties encountered in sampling conformational space and sequence space simultaneously are discussed.

Amino Acid Motifs↗

PubMed2003

Role of backbone hydration and salt-bridge formation in stability of alpha-helix in solution.

We test molecular level hypotheses for the high thermal stability of alpha-helical conformations of alanine-based peptides by performing detailed atomistic simulations of a 20-amino-acid peptide with explicit treatment of water. To assess the contribution of large side chains to alpha-helix stability through backbone desolvation and salt-bridge formation, we simulate the alanine-rich peptide, Ac-YAEAAKAAEAAKAAEAAKAF-Nme, referred to as the EK peptide, that has three pairs of "i, i + 3" glutamic acid(-) and lysine(+) substitutions. Efficient configurational sampling of the EK peptide over a wide temperature range enabled by the replica exchange molecular dynamics technique allows characterization of the stability of alpha-helix with respect to heat-induced unfolding. We find that near ambient temperatures, the EK peptide predominately samples alpha-helical configurations with 80% fractional helicity at 300 K. The helix melts over a broad range of temperatures with melting temperature, T(m), equal to 350 K, that is significantly higher than the T(m) of a 21-residue polyalanine peptide, A(21). Salt-bridges between oppositely charged Glu(-) and Lys(+) side chains can, in principle, provide thermal stability to alpha-helical conformers. For the specific EK peptide sequence, we observe infrequent formation of Glu-Lys salt-bridges (with approximately 10-20% probability) and therefore we conclude that salt-bridge formation does not contribute significantly to the EK peptide's helical stability. However, lysine side chains are found to shield specific "i, i + 4" backbone hydrogen bonds from water, indicating that large side-chain substituents can play an important role in stabilizing alpha-helical configurations of short peptides in aqueous solution through mediation of water access to backbone hydrogen bonds. These observations have implications on molecular engineering of peptides and biomolecules in the design of their thermostable variants where the shielding mechanism can act in concert with other factors such as salt-bridge formation, thereby increasing thermal stability considerably.

Alanine↗