Structural details of urea binding to barnase: a molecular dynamics analysis

Molecular mechanism of the common and opposing cosolvent effects of fluorinated alcohol and urea on a coiled coil protein

Alcohols and urea are widely used as effective protein denaturants. Among monohydric alcohols, 2,2,2-trifluoroethanol (TFE) has large cosolvent effects as a helix stabilizer in proteins. In contrast, urea efficiently denatures ordered native structures, including helices, into coils. These opposing cosolvent effects of TFE and urea are well known, even though both preferentially bind to proteins; however, the underlying molecular mechanism remains controversial. Cosolvent-dependent relative stability between native and denatured states is rigorously related to the difference in preferential binding parameters (PBPs) between these states. In this study, GCN4-p1 with two-stranded coiled coil helices was employed as a model protein, and molecular dynamics simulations for the helix dimer and isolated coil were conducted in aqueous solutions with 2 M TFE and urea. As 2 M cosolvent aqueous solutions did not exhibit clustering of cosolvent molecules, we were able to directly investigate the molecular origin of the excess PBP without considering the enhancement effect of PBPs arising from the concentration fluctuations. The calculated excess PBPs of TFE for the helices and those of urea for the coils were consistent with experimentally observed stabilization of helix by TFE and that of coil by urea. The former was caused by electrostatic interactions between TFE and side chains of the helices, while the latter was attributed to both electrostatic and dispersion interactions between urea and the main chains. Unexpectedly, reverse-micelle-like orientations of TFE molecules strengthened the electrostatic interactions between TFE and the side chains, resulting in strengthening of TFE solvation.

NADP/H binding nearly doubles the stability of a Mycobacterium drug target: an unfolding study

Mycobacterium Aspartate beta semialdehyde dehydrogenase (ASADH) was studied using various spectroscopic techniques and size exclusion chromatography to examine the unfolding of free (apo) and NADP/H-bound (holo) forms of ASADH. Non-cooperative guanidinium chloride (GdnHCl)-induced unfolding of the apo ASADH was discovered, and no partially folded intermediate structures were stabilized. On the other hand, it was observed that GdnHCl's unfolding of holoenzyme was a cooperative process without any stable intermediate structure. The native form of holoenzyme is found to be stable against the lower concentration of GdnHCl only (namely up to 1.25 M GdnHCl). The tryptophan environment appears to unfold cooperatively in case of the holoenzyme and is in well coordination with the overall unfolding of the holoenzyme. The presence of NADP/H shows a stabilizing effect on the tryptophan environment as well as on the native NADP/H-bound enzyme. Δ G Solvent o values reveal nearly two-fold (∼1.9) conformationally more stable folded holoenzyme compared to its native apo state. The Cm for the apo and holo forms of ASADH are 1.3 and 1.9 M, respectively. Novel drug leads targeting the NADP/H binding domain of ASADH could offer promising drugs against extremely infective Mycobacterium tuberculosis.Communicated by Ramaswamy H. Sarma.

Using Metadynamics To Explore the Free Energy of Dewetting in Biologically Relevant Nanopores

Water confined within hydrophobic spaces can undergo cooperative dewetting transitions due to slight changes in water density and pressure that push water toward the vapor phase. Many transmembrane protein ion channels contain nanoscale hydrophobic pores that could undergo dewetting transitions, sometimes blocking the flow of ions without physical blockages. Standard molecular dynamics simulations have been extensively applied to study the behavior of water in nanoscale pores, but the large free energy barriers of dewetting often prevent direct sampling of both wet and dry states and quantitative studies of the hydration thermodynamics of biologically relevant pores. Here, we describe a metadynamics protocol that uses the number of waters within the pore as the collective variable to drive many reversible transitions between relevant hydration states and calculate well-converged free energy profiles of pore hydration. By creating model nanopore systems and changing their radius and morphology and including various cosolvents, we quantify how these pore properties and cosolvents affect the dewetting transition. The results reveal that the dewetting free energy of nanoscale pores is determined by two key thermodynamic parameters, namely, the effective surface tension coefficients of water-air and water-pore interfaces. Importantly, while the effect of salt can be fully captured in the water activity dependence, amphipathic cosolvents such as alcohols modify both dry and wet states of the pore and dramatically shift the wet-dry equilibrium. The metadynamics approach could be applied to studies of dewetting transitions within nanoscale pores of proteins and provide new insights into why different pore properties evolved in biological systems.

Unraveling the Influence of Osmolytes on Water Hydrogen-Bond Network: From Local Structure to Graph Theory Analysis

Water structure in aqueous osmolyte solutions, deduced from the slight alteration in the water-water radial distribution function, the decrease in water-water hydrogen bonding, and tetrahedral ordering based only on the orientation of nearest water molecules derived from the molecular dynamics simulations, appears to have been perturbed. A careful analysis, however, reveals that the hydrogen bonding and the tetrahedral ordering around a water molecule in binary solutions remain intact as in neat water when the contribution of osmolyte-water interactions is appropriately incorporated. Furthermore, the distribution of the water binding energies and the water excess chemical potential of solvation in solutions are also pretty much the same as in neat water. Osmolytes are, therefore, well integrated into the hydrogen-bond network of water. Indeed, osmolytes tend to preferentially hydrogen bond with water molecules and their interaction energies are strongly correlated to their hydrogen-bonding capability. The graph network analysis, further, illustrates that osmolytes act as hubs in the percolated hydrogen-bond network of solutions. The degree of hydrogen bonding of osmolytes predominantly determines all of the network properties. Osmolytes like ethanol that form fewer hydrogen bonds than a water molecule disrupt the water hydrogen-bond network, while other osmolytes that form more hydrogen bonds effectively increase the connectivity among water molecules. Our observation of minimal variation in the local structure and the vitality of osmolyte-water hydrogen bonds on the solution network properties clearly imply that the direct interaction between protein and osmolytes is solely responsible for the protein stability. Further, the relevance of hydrogen bonds on solution properties suggests that the hydrogen-bonding interaction among protein, water, and osmolyte could be the key determinant of the protein conformation in solutions.

Urea changes oocyte competence and gene expression in resultant bovine embryo in vitro

SummaryNutrition influences the microenvironment in the proximity of oocyte and affects early embryonic development. Elevated blood urea nitrogen, even in healthy dairy cows, is associated with reduced fertility and there is high correlation between blood urea levels and follicular fluid urea levels. Using a docking calculation (in silico), urea showed a favorable binding activity towards the ZP-N domain of ZP3, that of ZP2, and towards the predicted full-length sperm receptor ZP3. Supplementation of oocyte maturation medium with nutrition-related levels of urea (20 or 40 mg/dl as seen in healthy dairy cows fed on low or high dietary protein, respectively) dose-dependently increased: (i) the proportion of oocytes that remained uncleaved; and (ii) oocyte degeneration; and reduced cleavage, blastocyst and hatching rates. High levels of urea induced shrinkage in oocytes, visualised using scanning electron microscopy. Urea downregulated NANOG while dose-dependently upregulating OCT4, DNMT1, and BCL2 expression. Urea at 20 mg/dl induced BAX expression. Using mathematical modelling, the rate of oocyte degeneration was sensitive to urea levels; while cleavage, blastocyst and hatching rates exhibited negative sensitivity. The present data imply a novel role for urea in reducing oocyte competence and changing gene expression in the resultant embryos.

Time-dependent X-ray diffraction studies on urea/hen egg white lysozyme complexes reveal structural changes that indicate onset of denaturation

Temporal binding of urea to lysozyme was examined using X-ray diffraction of single crystals of urea/lysozyme complexes prepared by soaking native lysozyme crystals in solutions containing 9 M urea. Four different soak times of 2, 4, 7 and 10 hours were used. The five crystal structures (including the native lysozyme), refined to 1.6 Å resolution, reveal that as the soaking time increased, more and more first-shell water molecules are replaced by urea. The number of hydrogen bonds between urea and the protein is similar to that between protein and water molecules replaced by urea. However, the number of van der Waals contacts to protein from urea is almost double that between the protein and the replaced water. The hydrogen bonding and van der Waals interactions are initially greater with the backbone and later with side chains of charged residues. Urea altered the water-water hydrogen bond network both by replacing water solvating hydrophobic residues and by shortening the first-shell intra-water hydrogen bonds by 0.2 Å. These interaction data suggest that urea uses both 'direct' and 'indirect' mechanisms to unfold lysozyme. Specific structural changes constitute the first steps in lysozyme unfolding by urea.

Wool deconstruction using a benign eutectic melt

Wool fibre is deconstructed in a facile ‘top down’ fabrication process into functional, nano-dimensional α-keratin chains using a benign choline chloride-urea deep eutectic solvent (DES) melt.

Protein denaturants at aqueous-hydrophobic interfaces: self-consistent correlation between induced interfacial fluctuations and denaturant stability at the interface

The notion of direct interaction between denaturing cosolvent and protein residues has been proposed in dialogue relevant to molecular mechanisms of protein denaturation. Here we consider the correlation between free energetic stability and induced fluctuations of an aqueous-hydrophobic interface between a model hydrophobically associating protein, HFBII, and two common protein denaturants, guanidinium cation (Gdm(+)) and urea. We compute potentials of mean force along an order parameter that brings the solute molecule close to the known hydrophobic region of the protein. We assess potentials of mean force for different relative orientations between the protein and denaturant molecule. We find that in both cases of guanidinium cation and urea relative orientations of the denaturant molecule that are parallel to the local protein-water interface exhibit greater stability compared to edge-on or perpendicular orientations. This behavior has been observed for guanidinium/methylguanidinium cations at the liquid-vapor interface of water, and thus the present results further corroborate earlier findings. Further analysis of the induced fluctuations of the aqueous-hydrophobic interface upon approach of the denaturant molecule indicates that the parallel orientation, displaying a greater stability at the interface, also induces larger fluctuations of the interface compared to the perpendicular orientations. The correlation of interfacial stability and induced interface fluctuation is a recurring theme for interface-stable solutes at hydrophobic interfaces. Moreover, observed correlations between interface stability and induced fluctuations recapitulate connections to local hydration structure and patterns around solutes as evidenced by experiment (Cooper et al., J. Phys. Chem. A 2014, 118, 5657.) and high-level ab initio/DFT calculations (Baer et al., Faraday Discuss 2013, 160, 89).

Water-mediated potassium acetate intercalation in kaolinite as revealed by molecular simulation

Molecular simulations are suitable tools to study the adsorption and intercalation of molecules in clays. In this work, a recently proposed thermodynamically consistent force field for inorganic compounds (INTERFACE, Heinz H, Lin TJ, Mishra RK, Emami FS (2013) Langmuir 29:1754-1765), which enables accurate simulations of inorganic-organic interfaces, was tested for a two-sheet type clay mineral. All-atom NpT molecular dynamics simulations were used to describe the characteristics (basal spacing, loading, molecular orientation) of some intercalate complexes of kaolinite with potassium acetate and the results were compared with the available experimental data. The most probable structural configurations of the kaolinite/potassium acetate intercalate complexes were determined from the simulations. Our examinations confirmed some supposed (single- or double-layered) arrangements of guest molecules. The need of interlayer water in the intercalate complex, which can be produced by the basic synthesis procedure in air atmosphere, was verified.

Microsecond molecular dynamics simulation of guanidinium chloride induced unfolding of ubiquitin

All atom molecular dynamics simulations have been used to explore the atomic detail mechanism of guanidinium induced unfolding of the protein ubiquitin.

Effect of urea concentration on aggregation of amyloidogenic hexapeptides (NFGAIL)

We have performed large-scale all-atom molecular dynamics (MD) simulations to study the aggregation behavior of four NFGAIL hexapeptides in the aqueous urea solution, with a urea concentration ranging from 0 to 5 M. We find that urea in general suppresses the peptide aggregation, but suppression slows down in the intermediation concentration regime around 3 M. Two competing mechanisms of urea are determined: urea molecules accumulated near the first solvation shell (FSS) tend to unfold the hexapeptide, which favors aggregation; on the other hand, the tight hydrogen bonds formed between urea and peptide mainchains hinder the association of peptides which disfavors the formation of the β-sheet. Furthermore, the different nonlinear urea concentration dependences of the urea-peptide and peptide-peptide hydrogen bonds lead to a nonmonotonic behavior, with a weak enhancement in the peptide aggregation around 3 M.

Urea-induced denaturation of preQ1-riboswitch

Urea, a polar molecule with a large dipole moment, not only destabilizes folded RNA structures but can also enhance the folding rates of large ribozymes. Unlike the mechanism of urea-induced unfolding of proteins, which is well understood, the action of urea on RNA has barely been explored. We performed extensive all-atom molecular dynamics simulations to determine the molecular underpinnings of urea-induced RNA denaturation. Urea displays its denaturing power in both secondary and tertiary motifs of the riboswitch structure. Our simulations reveal that the denaturation of RNA structures is mainly driven by the hydrogen-bonding and stacking interactions of urea with the bases. Through detailed studies of the simulation trajectories, we found that geminate pairs between urea and bases due to hydrogen bonds and stacks persist only ~0.1-1 ns, which suggests that the urea-base interaction is highly dynamic. Most importantly, the early stage of base-pair disruption is triggered by penetration of water molecules into the hydrophobic domain between the RNA bases. The infiltration of water into the narrow space between base pairs is critical in increasing the accessibility of urea to transiently disrupted bases, thus allowing urea to displace inter-base hydrogen bonds. This mechanism--water-induced disruption of base pairs resulting in the formation of a "wet" destabilized RNA followed by solvation by urea--is the exact opposite of the two-stage denaturation of proteins by urea. In the latter case, initial urea penetration creates a dry globule, which is subsequently solvated by water, leading to global protein unfolding. Our work shows that the ability to interact with both water and polar or nonpolar components of nucleotides makes urea a powerful chemical denaturant for nucleic acids.

Collapse of unfolded proteins in a mixture of denaturants

Both urea and guanidinium chloride (GdmCl) are frequently used as protein denaturants. Given that proteins generally adopt extended or unfolded conformations in either aqueous urea or GdmCl, one might expect that the unfolded protein chains will remain or become further extended due to the addition of another denaturant. However, a collapse of denatured proteins is revealed using atomistic molecular dynamics simulations when a mixture of denaturants is used. Both hen egg-white lysozyme and protein L are found to undergo collapse in the denaturant mixture. The collapse of the protein conformational ensembles is accompanied by a decreased solubility and increased non-native self-interactions of hydrophobic residues in the urea/GdmCl mixture. The increase of non-native interactions rather than the native contacts indicates that the proteins experience a simple collapse transition from the fully denatured states. During the protein collapse, the relatively stronger denaturant GdmCl displays a higher tendency to be absorbed onto the protein surface due to their stronger electrostatic interactions with proteins. At the same time, urea molecules also accumulate near the protein surface, resulting in an enhanced "local crowding" for the protein near its first solvation shell. This rearrangement of denaturants near the protein surface and crowded local environment induce the protein collapse, mainly by burying their hydrophobic residues. These findings from molecular simulations are then further explained by a simple analytical model based on statistical mechanics.

Urea-Water Solvation Forces on Prion Structures

Solvation forces are crucial determinants in the equilibrium between the folded and unfolded state of proteins. Particularly interesting are the solvent forces of denaturing solvent mixtures on folded and misfolded states of proteins involved in neurodegeneration. The C-terminal globular domain of the ovine prion protein (1UW3) and its analogue H2H3 in the α-rich and β-rich conformation were used as model structures to study the solvation forces in 4 M aqueous urea using molecular dynamics. The model structures display very different secondary structures and solvent exposures. Most protein atoms favor interactions with urea over interactions with water. The force difference between protein-urea and protein-water interactions correlates with hydrophobicity; i.e., urea interacts preferentially with hydrophobic atoms, in agreement with results from solvent transfer experiments. Solvent Shannon entropy maps illustrate the mobility gradient of the urea-water mixture from the first solvation shell to the bulk. Single urea molecules replace water in the first solvation shell preferably at locations of relatively high solvent entropy.

Coherent microscopic picture for urea-induced denaturation of proteins

In a previous study, we explored the mechanism of urea-induced denaturation of proteins by performing molecular dynamics (MD) simulations of hen lysozyme in 8 M urea and supported the "direct interaction mechanism" whereby urea denatures protein via dispersion interaction (Hua, L.; Zhou, R. H.; Thirumalai, D.; Berne, B. J. Proc. Natl. Acad. Sci. U.S.A. 2008, 105, 16928). Here we perform large scale MD simulations of five representative protein/peptide systems in aqueous urea to investigate if the above mechanism is common to other proteins. In all cases, accumulations of urea around proteins/peptide are observed, suggesting that urea denatures proteins by directly attacking protein backbones and side chains rather than indirectly disrupting water structure as a "water breaker". Consistent with our previous case study of lysozyme, the current energetic analyses with five protein/peptide systems reveal that urea's preferential binding to proteins mainly comes from urea's stronger dispersion interactions with proteins than with bulk solution, whereas the electrostatic (hydrogen-bonded) interactions only play a relatively minor (even negative) role during this denaturation process. Furthermore, the simulations of the peptide system at different urea concentrations (8 and 4.5 M), and with different force fields (CHARMM and OPLSAA) suggest that the above mechanism is robust, independent of the urea concentration and force field used. Last, we emphasize the importance of periodic boundary conditions in pairwise energetic analyses. This article provides a comprehensive study on the physical mechanism of urea-induced protein denaturation and suggests that the "dispersion-interaction-driven" mechanism should be general.

Structure and permeation mechanism of a mammalian urea transporter

As an adaptation to infrequent access to water, terrestrial mammals produce urine that is hyperosmotic to plasma. To prevent osmotic diuresis by the large quantity of urea generated by protein catabolism, the kidney epithelia contain facilitative urea transporters (UTs) that allow rapid equilibration between the urinary space and the hyperosmotic interstitium. Here we report the first X-ray crystal structure of a mammalian UT, UT-B, at a resolution of 2.36 Å. UT-B is a homotrimer and each protomer contains a urea conduction pore with a narrow selectivity filter. Structural analyses and molecular dynamics simulations showed that the selectivity filter has two urea binding sites separated by an approximately 5.0 kcal/mol energy barrier. Functional studies showed that the rate of urea conduction in UT-B is increased by hypoosmotic stress, and that the site of osmoregulation coincides with the location of the energy barrier.

Comparison of chemical and thermal protein denaturation by combination of computational and experimental approaches. II

Chemical and thermal denaturation methods have been widely used to investigate folding processes of proteins in vitro. However, a molecular understanding of the relationship between these two perturbation methods is lacking. Here, we combined computational and experimental approaches to investigate denaturing effects on three structurally different proteins. We derived a linear relationship between thermal denaturation at temperature T(b) and chemical denaturation at another temperature T(u) using the stability change of a protein (ΔG). For this, we related the dependence of ΔG on temperature, in the Gibbs-Helmholtz equation, to that of ΔG on urea concentration in the linear extrapolation method, assuming that there is a temperature pair from the urea (T(u)) and the aqueous (T(b)) ensembles that produces the same protein structures. We tested this relationship on apoazurin, cytochrome c, and apoflavodoxin using coarse-grained molecular simulations. We found a linear correlation between the temperature for a particular structural ensemble in the absence of urea, T(b), and the temperature of the same structural ensemble at a specific urea concentration, T(u). The in silico results agreed with in vitro far-UV circular dichroism data on apoazurin and cytochrome c. We conclude that chemical and thermal unfolding processes correlate in terms of thermodynamics and structural ensembles at most conditions; however, deviations were found at high concentrations of denaturant.

Sequence-specific mapping of the interaction between urea and unfolded ubiquitin from ensemble analysis of NMR and small angle scattering data

The molecular details of how urea interacts with, and eventually denatures proteins, remain largely unknown. In this study we have used extensive experimental NMR data, in combination with statistical coil ensemble modeling and small-angle scattering, to analyze the conformational behavior of the protein ubiquitin in the presence of urea. In order to develop an atomic resolution understanding of the denatured state, conformational ensembles of full-atom descriptions of unfolded proteins, including side chain conformations derived from rotamer libraries, are combined with random sampling of explicit urea molecules in interaction with the protein. Using this description of the conformational equilibrium, we demonstrate that the direct-binding model of urea to the protein backbone is compatible with available experimental data. We find that, in the presence of 8 M urea, between 30 and 40% of the backbone peptide groups bind a urea molecule, independently reproducing results from a model-free analysis of small-angle neutron and X-ray scattering data. Crucially, this analysis also provides sequence specific details of the interaction between urea and the protein backbone. The pattern of urea-binding along the amino-acid sequence reveals a higher level of binding in the central part of the protein, a trend which resembles independent results derived from chemical shift mapping of the urea-protein interaction. Together these results substantiate the direct-binding model and provide a framework for studying the physical basis of interactions between proteins and solvent molecules.

Dissociation of hydrophobic and charged nano particles in aqueous guanidinium chloride and urea solutions: a molecular dynamics study

It has been a long history that urea and guanidinium chloride (GdmCl) are used as agents for denaturing proteins. The underlying mechanism has been extensively studied in the past several decades. However, the question regarding why GdmCl is much stronger than urea has seldom been touched. Here, through molecular dynamics simulations, we show that a 4 M GdmCl solution is more able than 7 M urea solution to dissociate both hydrophobic and charged nano-particles (NP). Both urea and GdmCl affect the NPs' aggregation through direct binding to the NP surface. The advantages of GdmCl originate from the net charge of bound guanidinium ions which can generate a local positively charged environment around hydrophobic and negatively charged NPs. This effective coating can introduce Coulombic repulsion between all the NPs. Urea shows certain ability to dissociate hydrophobic NPs. However, in the case of charged NPs, urea molecules located between two opposite-charged NPs will form ordered hydrogen bonds, acting like "glue" which prevents separation of the NPs. Although urea can form hydrogen bonds with either hydrophilic amino acids or the protein backbone, which are believed to contribute to protein denaturation, our findings strongly suggest that this property does not always contribute positively to urea's denaturation power.

Interconnection of salt-induced hydrophobic compaction and secondary structure formation depends on solution conditions: revisiting early events of protein folding at single molecule resolution

What happens in the early stage of protein folding remains an interesting unsolved problem. Rapid kinetics measurements with cytochrome c using submillisecond continuous flow mixing devices suggest simultaneous formation of a compact collapsed state and secondary structure. These data seem to indicate that collapse formation is guided by specific short and long range interactions (heteropolymer collapse). A contrasting interpretation also has been proposed, which suggests that the collapse formation is rapid, nonspecific, and a trivial solvent related compaction, which could as well be observed by a homopolymer (homopolymer collapse). We address this controversy using fluorescence correlation spectroscopy (FCS), which enables us to monitor the salt-induced compaction accompanying collapse formation and the associated time constant directly at single molecule resolution. In addition, we follow the formation of secondary structure using far UV CD. The data presented here suggest that both these models (homopolymer and heteropolymer) could be applicable depending on the solution conditions. For example, the formation of secondary structure and compact state is not simultaneous in aqueous buffer. In aqueous buffer, formation of the compact state occurs through a two-state co-operative transition following heteropolymer formalism, whereas secondary structure formation takes place gradually. In contrast, in the presence of urea, a compaction of the protein radius occurs gradually over an extended range of salt concentration following homopolymer formalism. The salt-induced compaction and the formation of secondary structure take place simultaneously in the presence of urea.

Motion of a disordered polypeptide chain as studied by paramagnetic relaxation enhancements, 15N relaxation, and molecular dynamics simulations: how fast is segmental diffusion in denatured ubiquitin?

Molecular dynamics (MD) simulations have been widely used to analyze dynamic conformational equilibria of folded proteins, especially in relation to NMR observables. However, this approach found little use in the studies of disordered proteins, where the sampling of vast conformational space presents a serious problem. In this paper, we demonstrate that the latest advances in computation technology make it possible to overcome this limitation. The experimentally validated (calibrated) MD models allow for new insights into structure/dynamics of disordered proteins. As a test system, we have chosen denatured ubiquitin in solution with 8 M urea at pH 2. High-temperature MD simulations in implicit solvent have been carried out for the wild-type ubiquitin as well as MTSL-tagged Q2C, D32C, and R74C mutants. To recalibrate the MD data (500 K) in relation to the experimental conditions (278 K, 8 M urea), the time axes of the MD trajectories were rescaled. The scaling factor was adjusted such as to maximize the agreement between the simulated and experimental (15)N relaxation rates. The resulting effective length of the trajectories, 311 μs, ensures good convergence properties of the MD model. The constructed MD model was validated against the array of experimental data, including additional (15)N relaxation parameters, multiple sets of paramagnetic relaxation enhancements (PREs), and the radius of gyration. In each case, a near-quantitative agreement has been obtained, suggesting that the model is successful. Of note, the MD-based approach rigorously predicts the quantities that are inherently dynamic, i.e., dependent on the motional correlation times. This cannot be accomplished, other than in empirical fashion, on the basis of static structural models (conformational ensembles). The MD model was further used to investigate the relative translational motion of the MTSL label and the individual H(N) atoms. The derived segmental diffusion coefficients proved to be nearly uniform along the peptide chain, averaging to D = 0.49-0.55 × 10(-6) cm(2)/s. This result was verified by direct analysis of the experimental PRE data using the recently proposed Ullman-Podkorytov model. In this model, MTSL and H(N) moieties are treated as two tethered spheres undergoing mutual diffusion in a harmonic potential. The fitting of the experimental data involving D as a single adjustable parameter leads to D = 0.45 × 10(-6) cm(2)/s, in good agreement with the MD-based analyses. This result can be compared with the range of estimates obtained from the resonance energy transfer experiments, D = 0.2-6.0 × 10(-6) cm(2)/s.

Glucose interactions with a model peptide

Molecular dynamics simulations have been conducted of the helical polypeptide melittin, in concentrated aqueous solutions of the alpha and beta anomers of D-glucopyranose. Glucose is an osmolyte, and it is expected to be preferentially excluded from the surfaces of proteins. This was indeed found to be the case in the simulations. The results indicate that the observed exclusion may have a contribution from an under-representation of hydrogen bonding interactions between glucose groups and exposed side chains, compared to water. However, glucose was found to bind quite specifically to melittin by stacking its hydrophobic face, consisting of aliphatic protons, against the flat hydrophobic face of the indole group of the tryptophan-19 side chain. Although the binding site for this interaction is localized, the binding is weak for both anomers, with a binding free energy estimated as only ∼0.5 kcal/mol (i.e. near k(B)T). The face of the sugar stacked against the Trp indole ring is different for the two anomers of glucose, due to the disruption of the H1-H3-H5 hydrophobic triad of the beta anomer by the axial C1 hydroxyl group in the alpha anomer. The measurable affinity of the sugar for the Trp side chain is consistent with the very frequent occurrence of this group in the binding sites of proteins that complex with sugars.

Binding of nonsteroidal anti-inflammatory drugs to Abeta fibril

Nonsteroidal anti-inflammatory drugs are considered as potential therapeutic agents against Alzheimer's disease. Using replica exchange molecular dynamics and atomistic implicit solvent model, we studied the mechanisms of binding of naproxen and ibuprofen to the Abeta fibril derived from solid-state NMR measurements. The binding temperature of naproxen is found to be almost 40 K higher than of ibuprofen implicating higher binding affinity of naproxen. The key factor, which enhances naproxen binding, is strong interactions between ligands bound to the surface of the fibril. The naphthalene ring in naproxen appears to provide a dominant contribution to ligand-ligand interactions. In contrast, ligand-fibril interactions cannot explain differences in the binding affinities of naproxen and ibuprofen. The concave fibril edge with the groove is identified as the primary binding location for both ligands. We show that confinement of the ligands to the groove facilitates ligand-ligand interactions that lowers the energy of the ligands bound to the concave edge compared with those bound to the convex edge. Our simulations appear to provide microscopic rationale for the differing binding affinities of naproxen and ibuprofen observed experimentally.

Chitosan molecular structure as a function of N-acetylation

Molecular dynamics simulations have been carried out to characterize the structure and solubility of chitosan nanoparticle-like structures as a function of the deacetylation level (0, 40, 60, and 100%) and the spatial distribution of the N-acetyl groups in the particles. The polysaccharide chains of highly N-deacetylated particles where the N-acetyl groups are uniformly distributed present a high flexibility and preference for the relaxed two-fold helix and five-fold helix motifs. When these groups are confined to a given region of the particle, the chains adopt preferentially a two-fold helix with ϕ and ψ values close to crystalline chitin. Nanoparticles with up to 40% acetylation are moderately soluble, forming stable aggregates when the N-acetyl groups are unevenly distributed. Systems with 60% or higher N-acetylation levels are insoluble and present similar degrees of swelling regardless the distribution of their N-acetyl groups. Overall particle solvation is highly affected by electrostatic forces resulting from the degree of acetylation. The water mobility and orientation around the polysaccharide chains affects the stability of the intramolecular O3-HO3((n)) ···O5((n +) (1)) hydrogen bond, which in turn controls particle aggregation.

Peptide conformational preferences in osmolyte solutions: transfer free energies of decaalanine

The nature in which the protecting osmolyte trimethylamine N-oxide (TMAO) and the denaturing osmolyte urea affect protein stability is investigated, simulating a decaalanine peptide model in multiple conformations of the denatured ensemble. Binary solutions of both osmolytes and mixed osmolyte solutions at physiologically relevant concentrations of 2:1 (urea:TMAO) are studied using standard molecular dynamics simulations and solvation free energy calculations. Component analysis reveals the differences in the importance of the van der Waals (vdW) and electrostatic interactions for protecting and denaturing osmolytes. We find that urea denaturation governed by transfer free energy differences is dominated by vdW attractions, whereas TMAO exerts its effect by causing unfavorable electrostatic interactions both in the binary solution and mixed osmolyte solution. Analysis of the results showed no evidence in the ternary solution of disruption of the correlations among the peptide and osmolytes, nor of significant changes in the strength of the water hydrogen bond network.

Influence of urea and guanidine hydrochloride on lysozyme stability and thermal denaturation; a correlation between activity, protein dynamics and conformational changes

The effect of urea and guanidine hydrochloride (GuHCl) on lysozyme stability has been investigated using activity measurements, microcalorimetry and Raman spectroscopy in the low-frequency and amide I regions. Raman investigations on lysozyme dissolved in H(2)O and D(2)O in the presence of up to 10 M denaturants have revealed direct binding between the protein and both denaturants. The analysis of isotopic exchanges in the amide I region allows the identification of binding sites as hydrophilic and hydrophobic groups, respectively, for urea and GuHCl. The weak loss of activity of lysozyme in the presence of urea (∼15% maximum) is mainly assigned to a transformation of the tertiary structure corresponding to a molten globule state without unfolding of α-helix structures, in contrast to GuHCl which clearly induces conformational changes, associated with a larger loss of activity (40% maximum). The denaturing power of urea and guanidine hydrochloride on lysozyme has been related to the solvent and protein dynamics, reflecting direct interaction between denaturants and protein. It clearly appears that solvent dynamics control protein dynamics, and the significant hardening of the dynamics of GuHCl aqueous solutions is considered responsible for its important denaturing power. The comparison between the low-frequency spectra of solvents and lysozyme aqueous solutions in the absence and presence of different types of additives (urea, GuHCl, trehalose) reveals the Raman signature of the hydration water dynamics. This comparison points out the exclusion of trehalose around the protein surface.

Backbone additivity in the transfer model of protein solvation

The transfer model implying additivity of the peptide backbone free energy of transfer is computationally tested. Molecular dynamics simulations are used to determine the extent of change in transfer free energy (DeltaG(tr)) with increase in chain length of oligoglycine with capped end groups. Solvation free energies of oligoglycine models of varying lengths in pure water and in the osmolyte solutions, 2M urea and 2M trimethylamine N-oxide (TMAO), were calculated from simulations of all atom models, and DeltaG(tr) values for peptide backbone transfer from water to the osmolyte solutions were determined. The results show that the transfer free energies change linearly with increasing chain length, demonstrating the principle of additivity, and provide values in reasonable agreement with experiment. The peptide backbone transfer free energy contributions arise from van der Waals interactions in the case of transfer to urea, but from electrostatics on transfer to TMAO solution. The simulations used here allow for the calculation of the solvation and transfer free energy of longer oligoglycine models to be evaluated than is currently possible through experiment. The peptide backbone unit computed transfer free energy of -54 cal/mol/M compares quite favorably with -43 cal/mol/M determined experimentally.

Preferential interactions between small solutes and the protein backbone: a computational analysis

To improve our understanding of the effects of small solutes on protein stability, we conducted atomistic simulations to quantitatively characterize the interactions between two broadly used small solutes, urea and glycine betaine (GB), and a triglycine peptide, which is a good model for a protein backbone. Multiple solute concentrations were analyzed, and each solute-peptide-water ternary system was studied with approximately 200-300 ns of molecular dynamics simulations with the CHARMM force field. The comparison between calculated preferential interaction coefficients (Gamma(23)) and experimentally measured values suggests that semiquantitative agreement with experiments can be obtained if care is exercised to balance interactions among the solute, protein, and water. On the other hand, qualitatively incorrect (i.e., wrong sign in Gamma(23)) results can be obtained if a solute model is constructed by directly taking parameters for chemically similar groups from an existing force field. Such sensitivity suggests that small solute thermodynamic data can be valuable in the development of accurate force field models of biomolecules. Further decomposition of Gamma(23) into group contributions leads to additional insights regarding the effects of small solutes on protein stability. For example, use of the CHARMM force field predicts that urea preferentially interacts with not only amide groups in the peptide backbone but also aliphatic groups, suggesting a role for these interactions in urea-induced protein denaturation; quantitatively, however, it is likely that the CHARMM force field overestimates the interaction between urea and aliphatic groups. The results with GB support a simple thermodynamic model that assumes additivity of preferential interaction between GB and various biomolecular surfaces.