NMR studies of protein hydration and TEMPOL accessibility

Solvent paramagnetic relaxation enhancement as a versatile method for studying structure and dynamics of biomolecular systems

Solvent paramagnetic relaxation enhancement (sPRE) is a versatile nuclear magnetic resonance (NMR)-based method that allows characterization of the structure and dynamics of biomolecular systems through providing quantitative experimental information on solvent accessibility of NMR-active nuclei. Addition of soluble paramagnetic probes to the solution of a biomolecule leads to paramagnetic relaxation enhancement in a concentration-dependent manner. Here we review recent progress in the sPRE-based characterization of structural and dynamic properties of biomolecules and their complexes, and aim to deliver a comprehensive illustration of a growing number of applications of the method to various biological systems. We discuss the physical principles of sPRE measurements and provide an overview of available co-solute paramagnetic probes. We then explore how sPRE, in combination with complementary biophysical techniques, can further advance biomolecular structure determination, identification of interaction surfaces within protein complexes, and probing of conformational changes and low-population transient states, as well as deliver insights into weak, nonspecific, and transient interactions between proteins and co-solutes. In addition, we present examples of how the incorporation of solvent paramagnetic probes can improve the sensitivity of NMR experiments and discuss the prospects of applying sPRE to NMR metabolomics, drug discovery, and the study of intrinsically disordered proteins.

Paramagnetic Chemical Probes for Studying Biological Macromolecules

Paramagnetic chemical probes have been used in electron paramagnetic resonance (EPR) and nuclear magnetic resonance (NMR) spectroscopy for more than four decades. Recent years witnessed a great increase in the variety of probes for the study of biological macromolecules (proteins, nucleic acids, and oligosaccharides). This Review aims to provide a comprehensive overview of the existing paramagnetic chemical probes, including chemical synthetic approaches, functional properties, and selected applications. Recent developments have seen, in particular, a rapid expansion of the range of lanthanoid probes with anisotropic magnetic susceptibilities for the generation of structural restraints based on residual dipolar couplings and pseudocontact shifts in solution and solid state NMR spectroscopy, mostly for protein studies. Also many new isotropic paramagnetic probes, suitable for NMR measurements of paramagnetic relaxation enhancements, as well as EPR spectroscopic studies (in particular double resonance techniques) have been developed and employed to investigate biological macromolecules. Notwithstanding the large number of reported probes, only few have found broad application and further development of probes for dedicated applications is foreseen.

Deriving Ligand Orientation in Weak Protein-Ligand Complexes by DEEP-STD NMR Spectroscopy in the Absence of Protein Chemical-Shift Assignment

Differential epitope mapping saturation transfer difference (DEEP‐STD) NMR spectroscopy is a recently developed powerful approach for elucidating the structure and pharmacophore of weak protein–ligand interactions, as it reports key information on the orientation of the ligand and the architecture of the binding pocket.1 The method relies on selective saturation of protein residues in the binding site and the generation of a differential epitope map by observing the ligand, which depicts the nature of the protein residues making contact with the ligand in the bound state. Selective saturation requires knowledge of the chemical‐shift assignment of the protein residues, which can be obtained either experimentally by NMR spectroscopy or predicted from 3D structures. Herein, we propose a simple experimental procedure to expand the DEEP‐STD NMR methodology to protein–ligand cases in which the spectral assignment of the protein is not available. This is achieved by experimentally identifying the chemical shifts of the residues present in binding hot‐spots on the surface of the receptor protein by using 2D NMR experiments combined with a paramagnetic probe.

Hot spot mapping of protein surfaces with TEMPOL: Bovine pancreatic RNase A as a model system

TEMPOL spin-label has been used to identify surface exposure of protein nuclei from NMR analysis of the induced paramagnetic relaxation enhancements (PRE). The absence of linear dependence between atom depths and observed PRE reveals that specific mechanisms drive the approach of the paramagnet to the protein surface. RNase A represents a unique protein system to explore the fine details of the information offered by TEMPOL induced PRE, due to the abundance of previous results, obtained in solution and in the crystal, dealing with surface dynamics behavior of this protein. MD simulations in explicit solvent have been performed, also in the presence of TEMPOL, in order to delineate the role of intermolecular hydrogen bonds (HB) on PRE extents. Comparison of our results with the ones obtained from multiple solvent crystal structure (MSCS) studies yields information on the specificities that these two techniques have for characterizing protein-ligand interactions, a fundamental step in the development of reliable surface druggability predictors.

Searching for protein binding sites from Molecular Dynamics simulations and paramagnetic fragment-based NMR studies

Hotspot delineation on protein surfaces represents a fundamental step for targeting protein-protein interfaces. Disruptors of protein-protein interactions can be designed provided that the sterical features of binding pockets, including the transient ones, can be defined. Molecular Dynamics, MD, simulations have been used as a reliable framework for identifying transient pocket openings on the protein surface. Accessible surface area and intramolecular H-bond involvement of protein backbone amides are proposed as descriptors for characterizing binding pocket occurrence and evolution along MD trajectories. TEMPOL induced paramagnetic perturbations on (1)H-(15)N HSQC signals of protein backbone amides have been analyzed as a fragment-based search for surface hotspots, in order to validate MD predicted pockets. This procedure has been applied to CXCL12, a small chemokine responsible for tumor progression and proliferation. From combined analysis of MD data and paramagnetic profiles, two CXCL12 sites suitable for the binding of small molecules were identified. One of these sites is the already well characterized CXCL12 region involved in the binding to CXCR4 receptor. The other one is a transient pocket predicted by Molecular Dynamics simulations, which could not be observed from static analysis of CXCL12 PDB structures. The present results indicate how TEMPOL, instrumental in identifying this transient pocket, can be a powerful tool to delineate minor conformations which can be highly relevant in dynamic discovery of antitumoral drugs.

Studying the structure and dynamics of biomolecules by using soluble paramagnetic probes

Characterisation of the structure and dynamics of large biomolecules and biomolecular complexes by NMR spectroscopy is hampered by increasing overlap and severe broadening of NMR signals. As a consequence, the number of available NMR spectroscopy data is often sparse and new approaches to provide complementary NMR spectroscopy data are needed. Paramagnetic relaxation enhancements (PREs) obtained from inert and soluble paramagnetic probes (solvent PREs) provide detailed quantitative information about the solvent accessibility of NMR-active nuclei. Solvent PREs can be easily measured without modification of the biomolecule; are sensitive to molecular structure and dynamics; and are therefore becoming increasingly powerful for the study of biomolecules, such as proteins, nucleic acids, ligands and their complexes in solution. In this Minireview, we give an overview of the available solvent PRE probes and discuss their applications for structural and dynamic characterisation of biomolecules and biomolecular complexes.

Solution NMR studies on the orientation of membrane-bound peptides and proteins by paramagnetic probes

Many peptides and proteins are attached to or immersed in a biological membrane. In order to understand their function not only the structure but also their topology in the membrane is important. Solution NMR spectroscopy is one of the most often used approaches to determine the orientation and localization of membrane-bound peptides and proteins. Here we give an application-oriented overview on the use of paramagnetic probes for the investigation of membrane-bound peptides and proteins. The examples discussed range from the large pool of antimicrobial peptides, bacterial toxins, cell penetrating peptides to domains of larger proteins or the calcium regulating protein phospholamban. Topological information is obtained in all these examples by the use of either attached or freely mobile paramagnetic tags. For some examples information obtained from the paramagnetic probes was included in the structure determination.

Determining the orientation and localization of membrane-bound peptides

Many naturally occurring bioactive peptides bind to biological membranes. Studying and elucidating the mode of interaction is often an essential step to understand their molecular and biological functions. To obtain the complete orientation and immersion depth of such compounds in the membrane or a membrane-mimetic system, a number of methods are available, which are separated in this review into four main classes: solution NMR, solid-state NMR, EPR and other methods. Solution NMR methods include the Nuclear Overhauser Effect (NOE) between peptide and membrane signals, residual dipolar couplings and the use of paramagnetic probes, either within the membrane-mimetic or in the solvent. The vast array of solid state NMR methods to study membrane-bound peptide orientation and localization includes the anisotropic chemical shift, PISA wheels, dipolar waves, the GALA, MAOS and REDOR methods and again the use of paramagnetic additives on relaxation rates. Paramagnetic additives, with their effect on spectral linewidths, have also been used in EPR spectroscopy. Additionally, the orientation of a peptide within a membrane can be obtained by the anisotropic hyperfine tensor of a rigidly attached nitroxide label. Besides these magnetic resonance techniques a series of other methods to probe the orientation of peptides in membranes has been developed, consisting of fluorescence-, infrared- and oriented circular dichroism spectroscopy, colorimetry, interface-sensitive X-ray and neutron scattering and Quartz crystal microbalance.

Hydration studies on the archaeal protein Sso7d using NMR measurements and MD simulations

Background

How proteins approach surrounding molecules is fundamental to our understanding of the specific interactions that occur at the surface of proteins. The enhanced surface accessibility of small molecules such as organic solvents and paramagnetic probes to protein binding sites has been observed; however, the molecular basis of this finding has not been fully established. Recently, it has been suggested that hydration dynamics play a predominant role in controlling the distribution of hot spots on surface of proteins.

Results

In the present study, the hydration of the archaeal multifunctional protein Sso7d from Solfolobus solfataricus was investigated using a combination of computational and experimental data derived from molecular dynamics simulations and ePHOGSY NMR spectroscopy.

Conclusions

We obtained a convergent protein hydration landscape that indicated how the shape and stability of the Sso7d hydration shell could modulate the function of the protein. The DNA binding domain overlaps with the protein region involved in chaperon activity and this domain is hydrated only in a very small central region. This localized hydration seems to favor intermolecular approaches from a large variety of ligands. Conversely, high water density was found in surface regions of the protein where the ATP binding site is located, suggesting that surface water molecules play a role in protecting the protein from unspecific interactions.

Nitroxides: Versatile Reporters and Reactants

Stable nitroxyl radicals are important tools in chemistry, biophysics, biology, and materials science. Their stability and the sensitivity of their EPR spectra to the local environment make them valuable molecular probes. This review seeks to give an overview of the developments in the field of nitroxide spin probes and their various applications, with the main focus on the pH-sensitive imidazoline nitroxide family.

Approaches for the measurement of solvent exposure in proteins by 19F NMR

Fluorine NMR is a useful tool to probe protein folding, conformation and local topology owing to the sensitivity of the chemical shift to the local electrostatic environment. As an example we make use of (19)F NMR and 3-fluorotyrosine to evaluate the conformation and topology of the tyrosine residues (Tyr-99 and Tyr-138) within the EF-hand motif of the C-terminal domain of calmodulin (CaM) in both the calcium-loaded and calcium-free states. We critically compare approaches to assess topology and solvent exposure via solvent isotope shifts, (19)F spin-lattice relaxation rates, (1)H-(19)F nuclear Overhauser effects, and paramagnetic shifts and relaxation rates from dissolved oxygen. Both the solvent isotope shifts and paramagnetic shifts from dissolved oxygen sensitively reflect solvent exposed surface areas.

An NMR method for the determination of protein binding interfaces using TEMPOL-induced chemical shift perturbations

BACKGROUND

The determination of protein-protein interfaces is of crucial importance to understand protein function and to guide the design of compounds. To identify protein-protein interface by NMR spectroscopy, 13C NMR paramagnetic shifts induced by freely diffusing 4-hydroxy-2, 2, 6, 6-tetramethyl-piperidine-1-oxyl (TEMPOL) are promising, because TEMPOL affects distinct 13C NMR chemical shifts of the solvent accessible nuclei belonging to proteins of interest, while 13C nuclei within the interior of the proteins may be distinguished by a lack of such shifts.

METHOD

We measured the 13C NMR paramagnetic shifts induced by TEMPOL by recording 13C-(13)C TOCSY spectra for ubiquitin in the free state and the complex state with yeast ubiquitin hydrolase1 (YUH1).

RESULTS

Upon complexation of ubiquitin with YUH1, 13C NMR paramagnetic shifts associated with the protein binding interface were reduced by 0.05 ppm or more. The identified interfacial atoms agreed with the prior X-ray crystallographic data.

CONCLUSIONS

The TEMPOL-induced 13C chemical shift perturbation is useful to determine precise protein-protein interfaces.

GENERAL SIGNIFICANCE

The present method is a useful method to determine protein-protein interface by NMR, because it has advantages in easy sample preparations, simple data analyses, and wide applicabilities.

Positioning of micelle-bound peptides by paramagnetic relaxation enhancements

Many peptides, proteins, and drugs interact with biological membranes, and knowing the mode of binding is essential to understanding their biological functions. To obtain the complete orientation and immersion depth of such a compound, the membrane-mimetic system (micelle) is placed in an aqueous buffer containing the soluble and inert paramagnetic contrast agent Gd(DTPA-BMA). Paramagnetic relaxation enhancements (PREs) of a specific nucleus then depend only on its distance from the surface. The positioning of a structurally characterized compound can be obtained by least-squares fitting of experimental PREs to the micelle center position. This liquid-state NMR approach, which does not rely on isotopic labeling or chemical modification, has been applied to determine the location of the presumed transmembrane region 7 of yeast V-ATPase (TM7) and the membrane-bound antimicrobial peptide CM15 in micelles. TM7 binds in a trans-micelle orientation with the N-terminus being slightly closer to the surface than the C-terminus. CM15 is immersed unexpectedly deep into the micelle with the more hydrophilic side of the helix being closer to the surface than the hydrophobic one.

Guanidinoneomycin B recognition of an HIV-1 RNA helix

Aminoglycoside antibiotics are small-molecule drugs that bind RNA. The affinity and specificity of aminoglycoside binding to RNA can be increased through chemical modification, such as guanidinylation. Here, we report the binding of guanidinoneomycin B (GNB) to an RNA helix from the HIV-1 frameshift site. The binding of GNB increases the melting temperature (T(m)) of the frameshift-site RNA by at least 10 degrees C, to a point at which a melting transition is not even observed in 2 M urea. A structure of the complex was obtained by using multidimensional heteronuclear NMR spectroscopic methods. We also used a novel paramagnetic-probe assay to identify the site of GNB binding to the surface of the RNA. GNB makes major-groove contacts to two sets of Watson-Crick bases and is in van der Waals contact with a highly structured ACAA tetraloop. Rings I and II of GNB fit into the major groove and form the binding interface with the RNA, whereas rings III and IV are exposed to the solvent and disordered. The binding of GNB causes a broadening of the major groove across the binding site.

Measuring the dynamic surface accessibility of RNA with the small paramagnetic molecule TEMPOL

The surface accessibility of macromolecules plays a key role in modulating molecular recognition events. RNA is a complex and dynamic molecule involved in many aspects of gene expression. However, there are few experimental methods available to measure the accessible surface of RNA. Here, we investigate the accessible surface of RNA using NMR and the small paramagnetic molecule TEMPOL. We investigated two RNAs with known structures, one that is extremely stable and one that is dynamic. For helical regions, the TEMPOL probing data correlate well with the predicted RNA surface, and the method is able to distinguish subtle variations in atom depths, such as the relative accessibility of pyrimidine versus purine aromatic carbon atoms. Dynamic motions are also detected by TEMPOL probing, and the method accurately reports a previously characterized pH-dependent conformational transition involving formation of a protonated C-A pair and base flipping. Some loop regions are observed to exhibit anomalously high accessibility, reflective of motions that are not evident within the ensemble of NMR structures. We conclude that TEMPOL probing can provide valuable insights into the surface accessibility and dynamics of RNA, and can also be used as an independent means of validating RNA structure and dynamics in solution.

NMR studies of lysozyme surface accessibility by using different paramagnetic relaxation probes

Paramagnetic probes, whose approach to proteins can be monitored by nuclear magnetic resonance (NMR) studies, have been found to be of primary relevance for investigating protein surfaces' accessibility. Here, a Gd(III) neutral complex which contains two metal ions, [Gd2(L7)(H2O)2], is suggested as a paramagnetic probe particularly suited for systematic NMR investigation of protein surface accessibility, due to an expected high relaxivity and to the lack of electric charge which could favor specific interactions. Hen egg white lysozyme has been used as a model system to verify the absence of preferential approaches of this paramagnetic probe to specific protein moieties by comparing paramagnetic perturbation profiles of 1H-13C HSQC signals obtained in the presence of TEMPOL and [Gd2(L7)(H2O)2]. From the similarity of the measured paramagnetic perturbation profiles induced by the two different probes, specific interactions of [Gd2(L7)(H2O)2] with the enzyme could be ruled out. The large size of the latter probe is suggested to be responsible for the strong paramagnetic perturbations observed for CalphaH groups which are located in convex surface-exposed regions. The combined use of the two probes reveals fine details of the dynamics controlling their approach toward the protein surface.

MD and NMR studies of alpha-bungarotoxin surface accessibility

Protein surface accessibility represents a dimension of structural biology which has not been discussed in details so far, in spite of its fundamental role in controlling the molecular recognition process. In the present report the surface accessibility of alpha-bungarotoxin, a small and well characterized protein, has been investigated by analyzing its interaction with solvent and paramagnetic molecules in an integrated way. The presence of strong hydration sites, identified by a combined analysis of MD simulation and NMR results, seems to prevent the access of Gd(III)DTPA-BMA to the protein surface. On the contrary, the limited hydration of the alpha-bungarotoxin active site favors frequent encounters between the paramagnetic probe and the protein in the latter region. All the data obtained here for alpha-bungarotoxin suggest that shape and stability of the solvation shell control its surface accessibility and, hence, intermolecular interactions in a way which could be common to many other proteins.

Long-Range-Distance NMR Effects in a Protein Labeled with a Lanthanide–DOTA Chelate

A two-thiol reactive lanthanide-DOTA (1,4,7,10-tetraazacyclododecane-N,N',N'',N'''-tetraacetic acid) chelate, CLaNP-3 (CLaNP=caged lanthanide NMR probe), was synthesized for the rigid attachment to cysteine groups on a protein surface, and used to obtain long-range-distance information from the {15N,1H} HSQC spectra of the protein-lanthanide complex. The DOTA ring exhibits several isomers that are in exchange; however, single resonances were observed for most amide groups in the protein, allowing determination of a single, apparent magnetic-susceptibility tensor. Pseudocontact shifts caused by Yb-containing CLaNP-3 were observed for atoms at 15-35 A from the metal. By using Gd-containing CLaNP-3, relaxation effects were observed, allowing distances up to 30 A from the paramagnetic center to be determined accurately. Similar results were obtained with a Gd-DTPA (diethylene-triaminepentaacetic acid) chelate, CLaNP-1, bound in the same bidentate manner to the protein. This study demonstrates that bidentate attachment of a paramagnetic probe enables determination of long-range distances.

The dynamics of water evaporation from partially solvated cytochrome c in the gas phase

The study of evaporation of water from biological macromolecules is important for the understanding of electrospray mass spectrometry experiments. In electrospray ionization (ESI), electrically charged nanoscale droplets are formed from solutions of, for example, proteins. Then evaporation of the solvent leads to dry protein ions that can be analyzed in the mass spectrometer. In this work the dynamics of water evaporation from native cytochrome c covered by a monolayer of water is studied by molecular dynamics (MD) simulations at constant energy. A model of the initial conditions of the process is introduced. The temperature of the protein drops by about 100 K during the 400 picoseconds of the simulations. This sharp drop in temperature causes the water evaporation rate to decrease by about an order of magnitude, leaving the protein with 50% to 90% of the original water molecules, depending on the initial temperature of the simulation. The structural changes of the protein upon desolvation were considered through calculations of the radius of gyration and the root mean square (RMS) of the protein. A variation of 0.4 A in the radius of gyration, together with an RMS value of less than 3 A, indicates only minor changes in the overall shape of the protein structure. The water coordination number of the solvation shell is much smaller than that for bulk water. The mobility of water is high at the beginning of the simulations and drops as the simulation progresses and the temperature decreases. Incomplete desolvation of protein ions was also observed in recent experiments.

19F NMR studies of solvent exposure and peptide binding to an SH3 domain

(19)F NMR was used to study topological features of the SH3 domain of Fyn tyrosine kinase for both the free protein and a complex formed with a binding peptide. Metafluorinated tyrosine was biosynthetically incorporated into each of 5 residues of the G48M mutant of the SH3 domain (i.e. residues 8, 10, 49 and 54 in addition to a single residue in the linker region to the C-terminal polyhistidine tag). Distinct (19)F NMR resonances were observed and subsequently assigned after separately introducing single phenylalanine mutations. (19)F NMR chemical shifts were dependent on protein concentration above 0.6 mM, suggestive of dimerization via the binding site in the vicinity of the tyrosine side chains. (19)F NMR spectra of Fyn SH3 were also obtained as a function of concentration of a small peptide (2-hydroxynicotinic-NH)-Arg-Ala-Leu-Pro-Pro-Leu-Pro-diaminopropionic acid -NH(2), known to interact with the canonical polyproline II (PPII) helix binding site of the SH3 domain. Based on the (19)F chemical shifts of Tyr8, Tyr49, and Tyr54, as a function of peptide concentration, an equilibrium dissociation constant of 18 +/- 4 microM was obtained. Analysis of the line widths suggested an average exchange rate, k(ex), associated with the peptide-protein two-site exchange, of 5200 +/- 600 s(-1) at a peptide concentration where 96% of the FynSH3 protein was assumed to be bound. The extent of solvent exposure of the fluorine labels was studied by a combination of solvent isotope shifts and paramagnetic effects from dissolved oxygen. Tyr54, Tyr49, Tyr10, and Tyr8, in addition to the Tyr on the C-terminal tag, appear to be fully exposed to the solvent at the metafluoro position in the absence of binding peptide. Tyr54 and, to some extent, Tyr10 become protected from the solvent in the peptide bound state, consistent with known structural data on SH3-domain peptide complexes. These results show the potential utility of (19)F-metafluorotyrosine to probe protein-protein interactions in conjunction with paramagnetic contrast agents.

NMR studies of BPTI aggregation by using paramagnetic relaxation reagents

Paramagnetic probes, whose approach to proteins can be monitored by nuclear magnetic resonance (NMR) studies, have been found of primary relevance for investigating protein surfaces accessibility. Here, paramagnetic probes are also suggested for a systematic investigation on protein aggregation. Bovine pancreatic trypsin inhibitor (BPTI) was used as a model system for aggregation by analyzing its interaction with TEMPOL and Gd(III)DTPA-BMA. Some of the measured paramagnetic relaxation rates of BPTI protons exhibited a reverse dependence on protein concentration, which can be attributed to the formation of transient BPTI aggregates.

High-resolution protein hydration NMR experiments: Probing how protein surfaces interact with water and other non-covalent ligands

High-resolution solution NMR experiments are extremely useful to characterize the location and the dynamics of hydrating water molecules at atomic resolution. However, these methods are severely limited by undesired incoherent transfer pathways such as those arising from exchange-relayed intra-molecular cross-relaxation. Here, we review several complementary exchange network editing methods that can be used in conjunction with other types of NMR hydration experiments such as magnetic relaxation dispersion and 1J(NC') measurements to circumvent these limitations. We also review several recent contributions illustrating how the original solution hydration NMR pulse sequence architecture has inspired new approaches to map other types of non-covalent interactions going well beyond the initial scope of hydration. Specifically, we will show how hydration NMR methods have evolved and have been adapted to binding site mapping, ligand screening, protein-peptide and peptide-lipid interaction profiling.

Toward the understanding of MNEI sweetness from hydration map surfaces

The binding mechanism of sweet proteins to their receptor, a G-protein-coupled receptor, is not supported by direct structural information. In principle, the key groups responsible for biological activity (glucophores) can be localized on a small structural unit (sweet finger) or spread on a larger surface area. A recently proposed model, called "wedge model", implies a large surface of interaction with the receptor. To explore this model in greater detail, it is necessary to examine the physicochemical features of the surfaces of sweet proteins, since their interaction with the receptor, with respect to that of small sweeteners, is more dependent on general physicochemical properties of the interface, such as electrostatic potential and hydration. In this study, we performed exhaustive molecular dynamics simulations in explicit water of the sweet protein MNEI and of its structural mutant G-16A, whose sweetness is one order of magnitude lower than that of MNEI. Solvent density and self-diffusion calculated from molecular dynamics simulations suggest a likely area of interaction delimited by four stretches arranged as a tetrahedron whose shape is complementary to that of a cavity on the surface of the receptor, in agreement with the wedge model. The suggested area of interaction is amazingly consistent with known mutagenesis data. In addition, the asymmetric hydration of the only helix in both proteins hints at a specific role for this secondary structure element in orienting the protein during the binding process.

Molecular dynamics of a protein surface: ion-residues interactions

Time-resolved measurements indicated that protons could propagate on the surface of a protein or a membrane by a special mechanism that enhanced the shuttle of the proton toward a specific site. It was proposed that a suitable location of residues on the surface contributes to the proton shuttling function. In this study, this notion was further investigated by the use of molecular dynamics simulations, where Na(+) and Cl(-) are the ions under study, thus avoiding the necessity for quantum mechanical calculations. Molecular dynamics simulations were carried out using as a model a few Na(+) and Cl(-) ions enclosed in a fully hydrated simulation box with a small globular protein (the S6 of the bacterial ribosome). Three independent 10-ns-long simulations indicated that the ions and the protein's surface were in equilibrium, with rapid passage of the ions between the protein's surface and the bulk. However, it was noted that close to some domains the ions extended their duration near the surface, thus suggesting that the local electrostatic potential hindered their diffusion to the bulk. During the time frame in which the ions were detained next to the surface, they could rapidly shuttle between various attractor sites located under the electrostatic umbrella. Statistical analysis of the molecular dynamics and electrostatic potential/entropy consideration indicated that the detainment state is an energetic compromise between attractive forces and entropy of dilution. The similarity between the motion of free ions next to a protein and the proton transfer on the protein's surface are discussed.

Three-dimensional computation of atom depth in complex molecular structures

MOTIVATION

For a complex molecular system the delineation of atom-atom contacts, exposed surface and binding sites represents a fundamental step to predict its interaction with solvent, ligands and other molecules. Recently, atom depth has been also considered as an additional structural descriptor to correlate protein structure with folding and functional properties. The distance between an atom and the nearest water molecule or the closest surface dot has been proposed as a measure of the atom depth, but, in both cases, the 3D character of depth is largely lost. In the present study, a new approach is proposed to calculate atom depths in a way that the molecular shape can be taken into account.

RESULTS

An algorithm has been developed to calculate intersections between the molecular volume and spheres centered on the atoms whose depth has to be quantified. Many proteins with different size and shape have been chosen to compare the results obtained from distance-based and volume-based depth calculations. From the wealth of experimental data available for hen egg white lysozyme, H/D exchange rates and TEMPOL induced paramagnetic perturbations have been analyzed both in terms of depth indexes and of atom distances to the solvent accessible surface. The algorithm here proposed yields better correlations between experimental data and atom depth, particularly for those atoms which are located near to the protein surface.

AVAILABILITY

Instructions to obtain source code and the executable program are available either from http://sienabiografix.com or http://sadic.sourceforge.net

CONTACT

niccolai@unisi.it

SUPPLEMENTARY INFORMATION

http://www.Sienabiogzefix.com/publication.

An NMR Method for the Determination of Protein-Binding Interfaces Using Dioxygen-Induced Spin−Lattice Relaxation Enhancement

Using oxygen as a paramagnetic probe, researchers can routinely study topologies and protein-binding interfaces by NMR. The paramagnetic contribution to the amide (1)H spin-lattice relaxation rates (R(1)(P)) have been studied for uniformly (2)H,(15)N-labeled FB protein, a 60-residue three-helix bundle, constituting the B domain of protein A. Through TROSY versions of inversion-recovery experiments, R(1)(P) could be determined. R(1)(P) was then measured in the presence of a stoichiometric equivalent of an unlabeled Fc fragment of immunoglobulin (Ig) G, and the ratio of R(1)(P) of the FB-Fc complex to that of free FB [i.e., R(1)(P)(complex)/R(1)(P)(free)] was determined for each observable residue. Regions of helix I and helix II, which were previously known to interact with Fc, were readily identified as belonging to the binding interface by their characteristically reduced values of R(1)(P)(complex)/R(1)(P)(free). The method of comparing oxygen-induced spin-lattice relaxation rates of free protein and protein-protein complexes, to detect binding interfaces, offers greater sensitivity than chemical shift perturbation, while it is not necessary to heavily deuterate the labeled protein, as is the case in cross saturation experiments.

On the dynamics of water molecules at the protein solute interfaces

Proteins, with the large variety of chemical groups they present at their molecular surface, are a class of molecules which can be very informative on most of the possible solute-solvent interactions. Hen egg white lysozyme has been used as a probe to investigate the complex solvent dynamics occurring at the protein surface, by analysing the results obtained from Nuclear Magnetic Resonance, X-ray diffractometry and Molecular Dynamics simulations. A consistent overall picture for the dynamics of water molecules close to the protein is obtained, suggesting that a rapid exchange occurs, in a picosecond timescale, among all the possible hydration surface sites both in solution and the solid state, excluding the possibility that solvent molecules can form liquid-crystal-like supramolecular adducts, which have been proposed as a molecular basis of 'memory of water'.

High sensitivity 1H-NMR spectroscopy of homeopathic remedies made in water

Background

The efficacy of homeopathy is controversial. Homeopathic remedies are made via iterated shaking and dilution, in ethanol or in water, from a starting substance. Remedies of potency 12 C or higher are ultra-dilute (UD), i.e. contain zero molecules of the starting material. Various hypotheses have been advanced to explain how a UD remedy might be different from unprepared solvent. One such hypothesis posits that a remedy contains stable clusters, i.e. localized regions where one or more hydrogen bonds remain fixed on a long time scale. High sensitivity proton nuclear magnetic resonance spectroscopy has not previously been used to look for evidence of differences between UD remedies and controls.

Methods

Homeopathic remedies made in water were studied via high sensitivity proton nuclear magnetic resonance spectroscopy. A total of 57 remedy samples representing six starting materials and spanning a variety of potencies from 6 C to 10 M were tested along with 46 controls.

Results

By presaturating on the water peak, signals could be reliably detected that represented H-containing species at concentrations as low as 5 μM. There were 35 positions where a discrete signal was seen in one or more of the 103 spectra, which should theoretically have been absent from the spectrum of pure water. Of these 35, fifteen were identified as machine-generated artifacts, eight were identified as trace levels of organic contaminants, and twelve were unexplained. Of the unexplained signals, six were seen in just one spectrum each. None of the artifacts or unexplained signals occurred more frequently in remedies than in controls, using a p < .05 cutoff. Some commercially prepared samples were found to contain traces of one or more of these small organic molecules: ethanol, acetate, formate, methanol, and acetone.

Conclusion

No discrete signals suggesting a difference between remedies and controls were seen, via high sensitivity ¹H-NMR spectroscopy. The results failed to support a hypothesis that remedies made in water contain long-lived non-dynamic alterations of the H-bonding pattern of the solvent.

Refolding of lysozyme by quasistatic and direct dilution reaction paths: a first-order-like state transition

A first-order-like state transition model is considered to be a global reaction mechanism to directly folded proteins from an unfolded state to its native form. In order to verify the general applicability of this mechanism, we used lysozyme as a model protein. It was fully unfolded by 4.5 M urea, 0.1 M dithiothreitol (DTT) in pH 3 and refolded to its native form by way of an overcritical reaction path (a quasistatic process) or directly crossing transition boundary path (a directly dilution process). In addition to the two states coexisting in the direct folding path, lysozyme might be trapped in a glassy state. However, it can escape from the glassy state by concentration twice. This indicates the existence of a state transition line or boundary in the direct folding reaction. However, lysozyme can continuously fold from unfolded to native by an overcritical reaction path. During the overcritical path, four stable folding intermediates and native lysozyme were obtained. The secondary structures, particle size distributions, thermal stabilities, and oxidation state of disulfide bonds of folding intermediates were analyzed by circular dichroism spectra, dynamic light scattering, differential scanning calorimetry, and Raman spectra, respectively. According to the data, the intermediates of both the overcritical reaction and the direct crossing transition boundary paths can be described by a common concept pertaining to a model that undergoes collapse, sequential, and first-order-like state transition. This indicated that protein folding by way of different reaction paths might follow a similar folding mechanism-i.e., a mechanism of overcritical folding of intermediates. A protein folding reaction diagram is postulated and discussed. In spite of a global interaction mechanism the alpha -helix is formed prior to the beta -sheet, which may indicate that protein folding is initiated by local interactions.