Elucidation of stable intermediates in urea-induced unfolding pathway of human carbonic anhydrase IX

Human carbonic anhydrase IX (CAIX) has evolved as a promising biomarker for cancer prognosis, due to its overexpression in various cancers and restricted expression in normal tissue. However, limited information is available on its biophysical behavior. The unfolding of CAIX in aqueous urea solution was studied using all-atom molecular dynamics simulation approach. The results of this study revealed a stable intermediate state along the unfolding pathway of CAIX. At intermediate concentrations of urea (2.0–4.0 M), the protein displays a native-like structure with a large population of its secondary structure and hydrophobic contacts remaining intact in addition to small confined overall motions. Beyond 4.0 M urea, the unfolding is more gradual and at 8.0 M urea the structure is largely collapsed due to the solvent effect. The hydrophobic contact analysis suggests that the contact in terminal α-helices is separated initially which propagates in the loss of contacts from centrally located β-sheets. The reduction of 60–65% tertiary contacts in 7.0–8.0 M urea suggested the presence of residual structure in unfolded state and is confirmed with structural snap shot. Free energy landscape analysis suggested that unfolding of CAIX exists through the different intermediate states.     

Effects of osmolytes on solvent features of water in aqueous solutions

The solvatochromic solvent features of water (dipolarity/polarizability, π*, hydrogen bond donor acidity, α, and hydrogen bond acceptor basicity, β) of water have been determined in aqueous solutions of erythritol, glucose, inositol, sarcosine, xylitol and urea with concentrations from 0 to ~3 M and higher. The concentration effects of the osmolytes on the solvent features of water were characterized and compared with those reported previously for sorbitol, sucrose, trimethylamine N-oxide (TMAO), and trehalose. The solvent features of water in solutions of all osmolytes except TMAO and sarcosine were established to be linearly interrelated. It is shown that the concentration effects of essentially all nonionic osmolytes depend on osmolytes’ lipophilicity, molecular polarizability, and polar surface area. It is demonstrated that solubility of various compounds in aqueous solutions of glucose, sucrose, sorbitol, and urea of varied concentrations may be described in terms of solvent dipolarity/polarizability of water in these solutions. Surface tension of aqueous solutions of sucrose and sorbitol may also be described in the same terms. The relative permittivity of aqueous solutions of glucose and sucrose may be described in terms of the solvent hydrogen bond donor acidity of water. It is suggested that the effects of nonionic osmolytes on behavior of proteins and nucleic acids in aqueous media may be considered in terms of the altered solvent features of water instead of “nano-molecular crowding” effect.     

Effect of salt additives on protein partition in polyethylene glycol–sodium sulfate aqueous two-phase systems

Partitioning of 15 proteins in polyethylene glycol (PEG)–sodium sulfate aqueous two-phase systems (ATPS) formed by PEG of two different molecular weights, PEG-600 and PEG-8000 in the presence of different buffers at pH 7.4 was studied. The effect of two salt additives (NaCl and NaSCN) on the protein partition behavior was examined. The salt effects on protein partitioning were analyzed by using the Collander solvent regression relationship between the proteins partition coefficients in ATPS with and without salt additives. The results obtained show that the concentration of buffer as well as the presence and concentration of salt additives affects the protein partition behavior. Analysis of ATPS in terms of the differences between the relative hydrophobicity and electrostatic properties of the phases does not explain the protein partition behavior. The differences between protein partitioning in PEG-600–salt and PEG-8000–salt ATPS cannot be explained by the protein size or polymer excluded volume effect. It is suggested that the protein–ion and protein–solvent interactions in the phases of ATPS are primarily important for protein partitioning.     

Exploring the differences and similarities between urea and thermally driven denaturation of bovine serum albumin: intermolecular forces and solvation preferences

The interactions of bovine serum albumin (BSA) with urea/water were investigated by computer simulation. It was revealed that the BSA-hydrophobic residues in urea solutions favored contact with urea more than with water. Energy decomposition analysis showed that van der Waals energy was the dominant driving force behind urea affinity for hydrophobic residues, whereas coulombic attraction was largely responsible for water affinity for these residues. Meanwhile, urea–BSA hydrogen bond energies were found to be weaker than water–BSA hydrogen bond energies. The greater strength of water–BSA hydrogen bonds than urea–BSA hydrogen bonds, and the opposing preferential interaction between the BSA and urea suggest that disruption of hydrophobic interaction predominates urea–protein denaturation. In pure water, hydrophobic residues showed aggregation tendencies at 323 K, suggesting an increase in hydrophobicity, while at 353 K the residues were partly denatured due to loss of hydrogen bonds; thus, disruption of hydrophobic interactions appeared to contribute less to thermal denaturation.     

Effects of urea and sodium hydroxide on the molecular weight and conformation of alpha-(1-->3)-D-glucan from Lentinus edodes in aqueous solution

The weight-average molecular weight (Mw) and intrinsic viscosity ([eta]) of the alpha-(1-->3)-D-glucan (L-FV-II) from Lentinus edodes in 0.5 and 1.0 M NaOH aqueous solution containing urea, were studied by light scattering and viscometry. The Mw value of the glucan decreased with increase of the urea and NaOH concentration. A strong intermolecular hydrogen bonding confers water-insolubility on the glucan, but NaOH and especially urea, broke this hydrogen bonding leading to enhanced water-solubility. Use of 1.0 M urea-1.0 M NaOH as solvent broke not only intermolecular hydrogen bonds but also partial covalent bonds of the alpha-glucan in aqueous solution, resulting in a decrease of Mw and [eta]. The urea and NaOH concentrations, storage time with stirring, and mode of preparation of the polysaccharide in aqueous solution significantly affected the determination of Mw and [eta]. The dependences of specific rotation and fluorescence emission ratio of a probe on urea concentration showed that a change in the molecular conformation of the alpha-glucan in 0.5 M NaOH aqueous solution containing urea occurred in the range 0.4-0.6 M urea. The 0.5 M urea-0.5 M NaOH aqueous solution is a suitable solvent for the glucan, and the Mw and [eta] values obtained were 5.21 x 10(5) and 148 cm3 g(-1), respectively. Degradation of the glucan was obvious after storage for 15 months.     

Urea effects on protein stability: Hydrogen bonding and the hydrophobic effect

The effects of urea on protein stability have been studied using a model system in which we have determined the energetics of dissolution of a homologous series of cyclic dipeptides into aqueous urea solutions of varying concentration at 25°C using calorimetry. The data support a model in which urea denatures proteins by decreasing the hydrophobic effect and by directly binding to the amide units via hydrogen bonds. The data indicate also that the enthalpy of amide hydrogen bond formation in water is considerably higher than previously estimated. Previous estimates included the contribution of hydrophobic transfer of the α-carbon resulting in an overestimate of the binding between urea and the amide unit of the backbone and an underestimate of the binding enthalpy. Proteins 31:107–115, 1998. © 1998 Wiley-Liss, Inc.     

Enhancement of nitrogen release properties of urea–kaolinite fertilizer with chitosan binder

The use of controlled release fertilizer (CRF) has become a new trend to minimize environmental pollution. In this study, urea–kaolinite containing 20 wt% urea after one hour dry grinding was mixed with different concentrations of chitosan as a binder to prepare nitrogen-based CRF. Fourier transform infrared spectroscopy confirmed the hydrogen bonding between urea and kaolinite. Covalent interaction between urea–kaolinite and chitosan make the granules stronger. The nitrogen release was measured in 5 days interval using a diacetylmonoxime calorimetric method at a wavelength of 527 nm. The results illustrated that by increasing the chitosan concentration from 3 to 7.5%, nitrogen release decreased from 41.23 to 25.25% after one day and from 77.31 to 59.27% after 30 days incubation in water. Compressive stress at break tests confirmed that granules with chitosan 6% had the highest resistance and were chosen for ammonia volatilization tests. Ammonia volatilization was carried out using the forced-draft technique for a period of 10 weeks. The results showed that the total amount of ammonia loss for conventional urea fertilizer and urea–kaolinite–chitosan granules was 68.63 and 56.75%, respectively. This controlled release product could be applied in agricultural crop production purpose due to its controlled solubility in the soil, high nutrient use efficiency and potential economic benefits.     

Structural characteristics of thermostable immunogenic outer membrane protein from Salmonella enterica serovar Typhi

In this work, we explored the acid-induced unfolding pathway of non-porin outer membrane protein (OMP), an immunogenic protein from Salmonella Typhi, by monitoring the conformational changes over a pH range of 1.0–7.0 by circular dichroism, intrinsic fluorescence, ANS binding, acrylamide quenching, and dynamic light scattering. The spectroscopic measurements showed that OMP in its native state at pH 7.0 exists in more stable and compact conformation. In contrast, at pH 2.0, OMP retains substantial amount of secondary structure, disrupted side chain interactions, increased hydrodynamic radii, and nearly four-fold increase in ANS fluorescence with respect to the native state, indicating that MG state exists at pH 2.0. Quenching of tryptophan fluorescence by acrylamide further confirmed the accumulation of a partially unfolded state between native and unfolded state. The effect of pH on the conformation and thermostability of OMP points towards its heat resistance at neutral pH (T _m?~?69 °C at pH 7.0, monitored by change in MRE_{222 nm}). Acid unfolded state was also characterized by the lack of a cooperative thermal transition. All these results suggested that acid-induced unfolded state of OMP at pH 2.0 represented the molten globule state. The chemical denaturation studies with GuHCl and urea as denaturants showed dissimilar results. The chemical unfolding experiments showed that in both far-UV CD and fluorescence measurements, GuHCl is more efficient than urea. GuHCl is characterized by low C _m (~1 M), while urea is characterized by high C _m (~3 M). The fully unfolded states were reached at 2 M GuHCl and 4 M urea concentration, respectively. This study adds to several key considerations of importance in the development of therapeutic agents against typhoid fever for clinical purposes.     

Evidence for the residual tertiary structure in the urea-unfolded form of bacteriophage T5 endolysin

Using high-resolution NMR spectroscopy, we studied peculiarities of the unfolding process of the bacteriophage T5 endolysin (EndoT5) by strong denaturants. It was shown that in the absence of zinc ions this protein is mostly unfolded in the solution of 8 M urea or 6 M guanidine hydrochloride. However, in the presence of zinc ions EndoT5 unfolding can be achieved only in acidic solutions (at pH < 4.0), whereas at pH > 4.0 NMR spectra of the metal-bound protein (Zn²⁺–Ca²⁺–EndoT5 or Zn²⁺–EndoT5 complexes) exhibit a few chemical shifts characteristic of the native or native-like proteins. Our data, including the pH–titration curve with the pK of ~5, suggested involvement of the zinc-binding histidines in the stabilization of this protein. Up-field signals that appear in the NMR spectra of apo-EndoT5 in the presence of high concentrations of strong denaturants are probably derived from the amino acid residues included in the formation of structured hydrophobic cluster, which likely corresponds to the 81–93 region of EndoT5 and contains some residual tertiary structure. It is possible also that this hydrophobic fragment serves as a foundation for the formation of structured cluster in the unfolded state.     

Mesoscopic simulation of self-assembly of linear and dendritic copolymer poly(styrene)-b-poly(ethyleneglycol) in polar and non-polar solvents

In this work, we simulate the microphase separation of aqueous and non-aqueous solutions of diblock copolymer poly(styrene)-b-poly(ethyleneglycol) under different architectures (linear and linear–dendritic) by dissipative particle dynamics. The observed morphologies in water where poly(ethyleneglycol) (PEG) block is soluble are as follows: (1) at low concentrations spherical micelles, cylinders and bicontinuous structures are formed in dendritic structures and spheres, cylinders and perforated lamellas in linear structure. The architectures simulated at low–moderate concentrations show an evolution sphere → cylinder → bicontinuous or perforated lamellas as the concentration is increased. (2) At high concentrated solutions rich defect structures of the sponge type are formed. In a non-aqueous non-polar solution such as cyclohexane, which is a good solvent for the polystyrene block, the formation of well-defined aggregates at low concentrations is not observed; however, irregular structures are achieved in concentrated solutions. We compare these results with a polymeric chimera consisting of a mixture of linear poly(styrene) homopolymer and PEG homopolymer in the linear, G1 or G2 dendritic configurations. Our simulations are in agreement with the experimentally observed structures of these polymers.     

Theoretical studies on the structures,intra- and inter-molecular hydrogen bonding interactions in HNF and HNF–H2O clusters in the gaseous,aqueous and solid phases

Hydrazimium nitroformate ([N₂H₅]⁺[C(NO₂)₃]^? , HNF) is an ionic oxidiser used in solid propellants. Its properties are easily affected by H₂O because of its hygroscopicity. In this article, density functional theory (DFT) and molecular dynamics (MD) were employed to study the isolated HNF molecule and the HNF–H₂O cluster in gas phase and in the aqueous solution. Three stable conformations were obtained for HNF in the gas phase and in the aqueous solution, respectively, and each conformation can form several different HNF–H₂O clusters. Irrespective of whether it is in gas phase or in solution, intramolecular hydrogen bond interactions and other interactions (e.g. the binding energy, the dispersion energy, the second-order perturbation energy and the energy gap between frontier orbitals) of HNF are weaker in the clusters than in the isolated state. The initial decomposition energy of the cluster is lower than that of the isolated HNF molecule in both gaseous and aqueous phases, while the dissociation processes are the same. Molecular dynamic simulations showed that the clustered H₂O elongates and weakens the C–NO₂ bond in the solid HNF–H₂O cluster compared with that in the solid HNF. H₂O reduces and weakens intramolecular N–HΛO bonds too, and O–HΛN is the dominant intermolecular hydrogen bond between HNF and H₂O.     

Theoretical study on the structures and properties of mixtures of urea and choline chloride

In this work, we investigated in detail the structural characteristics of mixtures of choline chloride and urea with different urea contents by performing molecular dynamic (MD) simulations, and offer possible explanations for the low melting point of the eutectic mixture of choline chloride and urea with a ratio of 1:2. The insertion of urea molecules was found to change the density distribution of cations and anions around the given cations significantly, disrupting the long-range ordered structure of choline chloride. Moreover, with increasing urea concentration, the hydrogen bond interactions between choline cations and Cl^? anions decreased, while those among urea molecules obviously increased. From the hydrogen bond lifetimes, it was found that a ratio of 1:2 between choline chloride and urea is necessary for a reasonable strength of hydrogen bond interaction to maintain the low melting point of the mixture of choline chloride with urea. In addition, it was also deduced from the interaction energies that a urea content of 67.7 % may make the interactions of cation–anion, cation–urea and anion–urea modest, and thus results in the lower melting point of the eutectic mixture of choline chloride and urea. The present results may offer assistance to some extent for understanding the physicochemical properties of the eutectic mixture of choline chloride and urea, and give valuable information for the further development and application of deep eutectic solvents.     

The role of CS<Subscript>2</Subscript> in CS<Subscript>2</Subscript>/NMP mixed solvent in weakening the hydrogen bond of OH–N in coal: a DFT investigation

The interaction processes of trace amounts of N-methyl-2-pyrrolidinone (NMP), CS₂/NMP (1:1 by volume) and pure NMP solvent with the hydrogen bond of OH?N in coal were constructed and simulated by density functional theory methods. The distances and bond orders between the main related atoms, and the hydrogen bond energy of OH?N were calculated. The calculated results show that pure NMP solvent does not weaken the hydrogen bond of OH?N in coal. However, trace amounts of NMP and CS₂/NMP (1:1 by volume) have a strong capacity to weaken the hydrogen bond of OH?N in coal. The H2–N3 distances are elongated from 1.87 Å to 3.80 Å and 3.44 Å, the bond orders of H2–N3 all disappear, and the corresponding hydrogen bond energies of OH?N in coal decrease from 45.72 kJ mol^?1 to 7.06 and 11.24 kJ mol^?1, respectively. These results show that CS₂ added to pure NMP solvent plays an important role in releasing the original capacity of NMP to weaken the hydrogen bond of OH?N in coal, in agreement with experimental observations.     

Characterization of the denaturation of human alpha-lactalbumin in urea by molecular dynamics simulations

Molecular dynamics (MD) simulations were used to characterize the non-cooperative denaturation of the molten globule A-state of human alpha-lactalbumin by urea. A solvent of explicit urea and water molecules was used, corresponding to a urea concentration of approximately 6M. Three simulations were performed at temperatures of 293K, 360K and 400K, with lengths of 2 ns, 8 ns and 8 ns respectively. The results of the simulations were compared with experimental data from NMR studies of human alpha-lactalbumin and related peptides. During the simulations, hydrogen bonds were formed from the protein to both urea and water molecules as intra-protein hydrogen bonds were lost. Urea was shown to compete efficiently with water as both a hydrogen bond donor and acceptor. Radial distribution functions of water and urea around hydrophobic side chain atoms showed a significant increase in urea molecules in the solvation shell as the side chains became exposed during denaturation. A considerable portion of the native-like secondary structure persisted throughout the simulations. However, in the simulations at 360K and 400K, there were substantial changes in the packing of aromatic and other hydrophobic side chains in the protein, and many native contacts were lost. The results suggest that during the non-cooperative denaturation of the molten globule, secondary structure elements are stabilized by non-specific, non-native interactions.     

Effects of hesperidin,a flavanone glycoside interaction on the conformation,stability, and aggregation of lysozyme: multispectroscopic and molecular dynamic simulation studies?

Hesperidin (HESP), a flavanone glycoside, shows high antioxidant properties and posses ability to go through the blood–brain barrier. Therefore, it could be a potential drug molecule against aggregation based diseases such as Alzheimer’s, Parkinson’s, and systemic amyloidoses. In this work, we investigated the potential of HESP to interact with hen egg-white lysozyme (HEWL) monomer and prevent its aggregation. The HESP–HEWL binding studies were performed using a fluorescence quenching technique, molecular docking and molecular dynamics simulations. We found a strong interaction of HESP with the lysozyme monomer (K_a, ~ 5 × 10⁴ M^?1) mainly through hydrogen bonding, water bridges, and hydrophobic interactions. We showed that HESP molecule spanned the highly aggregation prone region (amino acid residues 48-101) of HEWL and prevented its fibrillar aggregation. Further, we found that HESP binding completely inhibited amorphous aggregation of the protein induced by disulfide-reducing agent tries-(2-carboxyethyl) phosphine. Conformational and stability studies as followed by various tertiary and secondary structure probes revealed that HESP binding only marginally affected the lysozyme monomer conformation and increased both stability and reversibility of the protein against thermal denaturation. Future studies should investigate detail effects of HESP on solvent dynamics, structure, and toxicity of various aggregates. The answers to these questions will not only target the basic sciences, but also have application in biomedical and biotechnological sciences.     

Residual structure in a peptide fragment of the outer membrane protein X under denaturing conditions: a molecular dynamics study

The Escherichia coli outer membrane protein X (OmpX) contains two polypeptide segments that present nonrandom residual structure in 8 M aqueous urea, whereas the remainder of the protein is in a flexibly disordered conformation (Tafer et al. in Biochemistry 43:860–869, 2004). In the present study, the results of two long-timescale (0.4 μs) unrestrained explicit-solvent molecular dynamics (MD) simulations of a tetradecapeptide representative of one of these two segments in 8 M aqueous urea are reported and analyzed. The two simulations were initiated either from the conformation of the corresponding segment in an NMR model structure of the unfolded protein or from an entirely extended configuration. The sampled conformational ensembles agree qualitatively with the experimentally observed NOEs, but not quantitatively, suggesting that a number of relevant configurations were not visited on the 2 × 0.4 μs timescale. Major conformational transitions occur on the 0.1 μs timescale, and the ensembles corresponding to the two independent simulations overlap only to a limited extent. However, both simulations show in multiple events the reversible formation and disruption of α-helical secondary structure (characteristic of the urea-denatured state) and β-turn secondary structure (characteristic of the native state). Events of helix formation are correlated with the appearance of hydrogen bonds between two side chains (Asp75–Ser78) and of a persistent hydrophobic contact (Trp76–Tyr80). They also evidence a peculiar helix stabilization and N-terminal capping role for a negatively charged residue (Asp75). These features are in good qualitative agreement with the NMR model for the structured state of the corresponding segment in the urea-denatured protein. The analysis of the simulations provides a detailed picture of the structural and dynamic features of the considered peptide at atomic resolution that is of high relevance in the understanding of the OmpX folding process.     

分子动力学模拟尿素对水溶液中蛋白质构象转变的影响

采用分子动力学方法和全原子模型研究尿素和水分子对模型蛋白S-肽链结构转化的影响。模拟结果显示S-肽链的变性速率常数k值随着尿素浓度的增加而先降低后升高,在尿素浓度为2.9 mol/L时达到最低值。模拟了不同尿素浓度下尿素-肽链、水-肽链以及肽链分子氢键的形成状况。结果表明:尿素浓度较低时,尿素分子与S-肽链的极性氨基酸侧链形成氢键,但不破坏其分子内的骨架氢键,尿素在S-肽链水化层外形成限制性空间,增强了S-肽链的稳定性。随着尿素的升高,尿素分子进入S-肽链内部并与其内部氨基酸残基形成氢键,导致S-肽链的骨架氢键丧失,S-肽链发生去折叠。上述模拟结果与文献报道的实验结果一致,从分子水平上揭示了尿素对蛋白质分子结构变化的影响机制,对于研究和发展蛋白质折叠及稳定化技术具有指导意义。     

Investigation and optimisation of the drying of reduced glutathione by steered molecular dynamics simulation

The drying of reduced glutathione from a series of aqueous–ethanol binary solutions at 300 K (below human body temperature) and 330 K (above human body temperature) was investigated in detail by steered molecular simulation and an umbrella sampling method with the Gromacs software package and Gromos96(53a6) united atomic force field. The results show that electrostatic interactions between glutathione and solvent represent the main resistance to drying. When the aqueous solution was gradually changed to pure ethanol, the energy of electrostatic interaction between glutathione and solvent molecules increased by 445.088 kJ/mol, and the drying potential of mean force (PMF) free energy also fell by 253.040 kJ/mol. However, an increase in temperature from 300 to 330 K in the aqueous solution only results in an increase of 23.013 kJ/mol in electrostatic interaction energy and a decrease of 34.956 kJ/mol in drying PMF free energy. Furthermore, we show that hydrogen bonding is the major form of electrostatic interaction involved, and directly affects the drying of glutathione. Therefore, choosing water-miscible solvents that minimise hydrogen-bond formation with glutathione will enhance its drying rate, and this is likely to be more efficient than increasing the temperature of the process. Thus, a power-saving technology can be used to produce the high bioactivity medicines.     

Solvent interaction analysis as a proteomic approach to structure-based biomarker discovery and clinical diagnostics

Proteins have several measurable features in biological fluids that may change under pathological conditions. The current disease biomarker discovery is mostly based on protein concentration in the sample as the measurable feature. Changes in protein structures, such as post-translational modifications and in protein–partner interactions are known to accompany pathological processes. Changes in glycosylation profiles are well-established for many plasma proteins in various types of cancer and other diseases. The solvent interaction analysis method is based on protein partitioning in aqueous two-phase systems and is highly sensitive to changes in protein structure and protein–protein- and protein–partner interactions while independent of the protein concentration in the biological sample. It provides quantitative index: partition coefficient representing changes in protein structure and interactions with partners. The fundamentals of the method are presented with multiple examples of applications of the method to discover and monitor structural protein biomarkers as disease-specific diagnostic indicators.