The electron–nucleus coupling: A breakthrough in the investigation of paramagnetic metalloproteins

Water proton relaxation measurements as a function of magnetic field on solutions containing paramagnetic metalloproteins contain information on (i) number and distance of metal-coordinated water molecules, (ii) dynamics of water exchange, (iii) electron relaxation rates, and (iv) dynamics of the metal coordination sphere influencing electron relaxation. The use of appropriate theoretical tools permits us in many cases to learn about some or all of the above properties.     

Gadolinium (III) ion in liquid water: structure, dynamics, and magnetic interactions from first principles

We applied first principles molecular dynamics (MD) technique to study structure, dynamics, and magnetic interactions of the Gd(3+) aqua ion dissolved in liquid water, a prototypical system for Gd-based complexes used as contrast agents for magnetic resonance imaging. The first coordination sphere contains eight water molecules with an average Gd-O distance of 2.37 A and an average geometric arrangement close to a square antiprism. The mean tilt angle of the electric dipole vector of these water molecules is theta=145 degrees . In our picosecond time scale simulation we observe no exchange event from the first coordination sphere but only fast "wagging" motions. The second coordination sphere is well pronounced though water molecules in this sphere are subjected to large amplitude dynamic motions. The isotropic hyperfine coupling constants for the inner sphere water molecules [A(iso)((17)O(I))=0.65+/-0.03 MHz, A(iso)((1)H(I))=0.085+/-0.005 MHz] are in good agreement with experimental data and with an earlier study using classical MD. Second sphere Fermi contact hyperfine coupling constants calculated are more than one order of magnitude smaller and of opposite sign as those of the first coordination sphere. The effect of spin polarization induced by the paramagnetic Gd(3+) ion on the dipolar hyperfine interaction was found to be sizable only for the (17)O nuclei of inner sphere water molecules and has a screening character.     

Design of Gd(III)-based magnetic resonance imaging contrast agents: static and transient zero-field splitting contributions to the electronic relaxation and their impact on relaxivity

A multiple-frequency (9.4-325 GHz) and variable-temperature (276-320 K) electron paramagnetic resonance (EPR) study on low molecular weight gadolinium(III) complexes for potential use as magnetic resonance imaging (MRI) contrast agents has been performed. Peak-to-peak linewidths Delta Hpp and central magnetic fields have been analyzed within the Redfield approximation taking into account the static zero-field splitting (ZFS) up to the sixth order and the transient ZFS up to the second order. Longitudinal electronic relaxation is dominated by the static ZFS contribution at low magnetic fields (B < 0.3 T) and by the transient ZFS at high magnetic fields (B > 1.5 T). Whereas the static ZFS clearly depends on the nature of the chelating ligand, the transient ZFS does not. For the relatively fast rotating molecules studied water proton relaxivity is mainly limited by the fast rotation and electronic relaxation has only a marked influence at frequencies below 30 MHz. From our EPR results we can conclude that electronic relaxation will have no influence on the efficiency of Gd(III)-based MRI contrast agents designed for studies at very high magnetic fields (B > 3T).     

Two-dimensional (2D) pulsed electron paramagnetic resonance study of VO(2+)-triphosphate interactions: evidence for tridentate triphosphate coordination, and relevance to bone uptake and insulin enhancement by vanadium pharmaceuticals

Two- and four-pulse electron spin echo envelope modulation (ESEEM) and four-pulse two-dimensional hyperfine sublevel correlation (HYSCORE) spectroscopies have been used to determine the solution structure of a 3:1 triphosphate:vanadyl solution at pH 5.0. Limited quantitative data were extracted from the two pulse spectra; however, HYSCORE proved to be more useful in the detection and interpretation of the (31)P and (1)H couplings. Three sets of cross-peaks were observed for each nucleus. For the (31)P couplings, three sets of cross-peaks were observed in the HYSCORE spectrum, and contour line shape analysis yielded coupling constants of approximately 15, 9, and 1 MHz. HYSCORE cross-peaks in the proton region were partially overlapping; however, interpretation of the proton coupling was simplified through the use of one-dimensional four-pulse ESEEM and subsequent analysis of the sum combination peaks. Comparison of the derived isotropic and anisotropic coupling constants with results from earlier ESEEM and electron nuclear double resonance (ENDOR) studies was consistent with the presence of at least one, and most likely two, water molecules coordinated in the equatorial plane of the vanadyl cation. The vanadyl-triphosphate system was shown to be an accurate model of the in vivo vanadyl-phosphate coupling constants determined in an earlier study (Dikanov, S. A.; Liboiron, B. D.; Thompson, K. H.; Vera, E.; Yuen, V. G.; McNeill, J. H.; Orvig, C. J. Am. Chem. Soc. 1999, 121, 11004.) Comparison of these values to those found in previous spectroscopic studies of vanadyl-triphosphate interactions, along with a detailed structural interpretation, are presented. This work represents the first detection of tridentate polyphosphate coordination to the vanadyl ion, and the first observation of an axial phosphate interaction not previously reported in earlier ENDOR and pulsed electron paramagnetic resonance studies.     

Dinuclear complexes formed with the triazacyclononane derivative ENOTA4-: high-pressure 17O NMR evidence of an associative water exchange on [MnII2(ENOTA)(H2O)2

Mn2+ has five unpaired d-electrons, a long electronic relaxation time, and labile water exchange, all of which make it an attractive candidate for contrast agent application in medical magnetic resonance imaging. In the quest for stable and nonlabile Mn2+ complexes, we explored a novel dimeric triazacyclononane-based ligand bearing carboxylate functional groups, H4ENOTA. The protonation constants of the ligand and the stability constants of the complexes formed with some endogenously important metals (Ca2+, Cu2+, Zn2+), as well as with Mn2+ and Ce3+, have been assessed by NMR methods, potentiometry, and UV-vis spectrophotometry. Overall, the thermodynamic stability of the complexes is lower as compared to that of the corresponding NOTA analogues (H3NOTA, 1,4,7-triaazacyclononane-1,4,7-triacetic acid). The crystal structure of Mn2(ENOTA)(H2O) x 5H2O contains two six-coordinated Mn2+, in addition to the three amine nitrogens and the two oxygens from the pendent monodentate carboxylate groups, and one water (Mn2) or one bridging carboxylate oxygen (Mn1) completes the coordination sphere of the metal ion. In an aqueous solution, this bridging carboxylate is replaced by a water molecule, as evidenced by the 17O chemical shifts and proton relaxivity data that point to monohydration for both metal ions in the dinuclear complex. A variable-temperature and -pressure 17O NMR study has been performed on [Mn2(ENOTA)(H2O)2] to assess the rate and, for the first time on a Mn2+ chelate, also the mechanism of the water exchange. The inner sphere water is slightly more labile in [Mn2(ENOTA)(H2O)2] (k298ex = 5.5 x 107 s-1) than in the aqua ion (2.1 x 107 s-1, Merbach, A. E.; et al. Inorg. Chem. 1980, 19, 3696). The water exchange proceeds via an almost limiting associative mechanism, as evidenced by the large negative activation volume (deltaV = -10.7 cm3 mol-1). The proton relaxivities measured on [Mn2(ENOTA)(H2O)2] show a low-field dispersion at approximately 0.1 MHz arising from a contact interaction between the MnII electron spin and the water proton nuclear spins.     

NMR relaxometric study of new Gd(III) macrocyclic complexes and their interaction with human serum albumin

Five novel Gd(iii) complexes based on the structure of the heptadentate macrocyclic 1,4,7,10-tetraazacyclododecane-1,4,7-triacetic acid (DO3A) ligand have been synthesized and their (1)H and (17)O NMR relaxometric properties investigated in detail. The complexes have been functionalised on the secondary nitrogen atom of the macrocyclic ring with different pendant groups for promoting their ability to interact non-covalently with human serum albumin (HSA). The analysis of the proton relaxivity, measured as a function of pH and magnetic field strength, have revealed that the three complexes bearing a poly(ethylene glycol)(PEG) chain possess a single coordinated water molecule, whereas the complexes functionalised with 1-[3-(2-hydroxyphenyl)]-propyl and 1-[3-(2-carboxyphenyloxy)]-propyl pendant groups have two inner sphere water molecules. The water exchange rates, measured by variable temperature (17)O NMR, cover a broad range of values (from 18 to 770 ns) as a function of their charge, the chemical nature of the substituent and its ability to organize a second sphere of hydration near the water(s) binding site. All the complexes have shown some degree of interaction with HSA, with a stronger binding affinity measured for those bearing an aromatic moiety on the pendant group. However, upon binding the expected relaxation enhancement has not been observed and this has been explained with the displacement of the coordinated water molecules by the protein and formation of ternary adducts.     

Structural and dynamical analysis of the hydration of the Alzheimer's beta-amyloid peptide

An analysis of the water molecules in the first solvation shell obtained from the molecular dynamics simulation of the amyloid beta(10-35)NH2 peptide and the amyloid beta(10-35)NH2E22Q "Dutch" mutant peptide is presented. The structure, energetics, and dynamics of water in the hydration shell have been investigated using a variety of measures, including the hydrogen bond network, the water residence times for all the peptide residues, the diffusion constant, experimentally determined HN amide proton exchange, and the transition probabilities for water to move from one residue to another or into the bulk. The results of the study indicate that: (1) the water molecules at the peptide-solvent interface are organized in an ordered structure similar for the two peptide systems but different from that of the bulk, (2) the peptide structure inhibits diffusion perpendicular to the peptide surface by a factor of 3 to 5 relative to diffusion parallel to the peptide surface, which is comparable to diffusion of bulk water, (3) water in the first solvation shell shows dynamical relaxation on fast (1-2 ps) and slow (10-40 ps) time scales, (4) a novel solvent relaxation master equation is shown to capture the details of the fast relaxation of water in the peptide's first solvation shell, (5) the interaction between the peptide and the solvent is stronger in the wild type than in the E22Q mutant peptide, in agreement with earlier results obtained from computer simulations [Massi, F.; Straub, J. E. Biophys J 2001, 81, 697] correlated with the observed enhanced activity of the E22Q mutant peptide.     

Theoretical studies of manganese and iron superoxide dismutases: superoxide binding and superoxide oxidation

Density-functional calculations indicate that the second sphere of coordination around the metal centers of manganese and iron superoxide dismutases (MnSODs and FeSODs) plays an important role in the binding of O2(-). In these systems, O2(-) prefers to bind to Mn or Fe in end-on configurations. For human and E. coli MnSODs, the bound O2(-) forms hydrogen bonds with the tyrosine and glutamine amino acid residues in the second sphere of coordination. In the cases of E. coli and T. elongates FeSODs, hydrogen bonding occurs between the bound O2(-) and the tyrosine amino acid only because the glutamine is too far away for an effective bonding interaction. The manner in which the O2(-) binds to the metal center in MnSODs and FeSODs can affect the rate of subsequent protonation and determine the mechanism for the formation of H2O2. Both Mn- and Fe-containing superoxide dismutases contain a metal-bound solvent molecule that has been suggested to be involved in the uptake of a H+ upon reduction of the metal center [Bull, C.; Fee, J. A. J. Am. Chem. Soc. 1985, 107, 3295; Miller, A.-F.; Padmakumar, K.; Sorkin, D. L.; Karapetian, A.; Vance, C. K. J. Inorg. Biochem. 2003, 93, 71]. Using density-functional theory, we confirm this suggestion and show the involvement of the second sphere of coordination in the process. We show that the oxidation of superoxide by Mn- or Fe-containing superoxide dismutases is facilitated by a cooperative effect between superoxide binding, protonation of the OH- bound to the metal, and electron transfer from the superoxide molecule to the oxidized metal. In particular, proton transfer through tyrosine-34 on the absence of a bound superoxide is uphill while, once superoxide is bound, the energetic barrier is lowered. It is this barrier that likely keeps the resting state (Mn(III)SOD) of the enzyme with a bound hydroxide, instead of a water. This work provides a model for the mechanism of reaction of superoxide with the oxidized form of the metal within Mn- and FeSODs.     

Chelators for radioimmunotherapy: I. NMR and ab initio calculation studies on 1,4,7,10-tetra(carboxyethyl)-1,4,7,10-tetraazacyclododecane (DO4Pr) and 1,4,7-tris(carboxymethyl)-10-(carboxyethyl)-1,4,7,10-tetraazacyclododecane (DO3A1Pr)

This work describes the modification of the chelating agent 1,4,7,10-tetraazacyclododecane-N,N',N' ',N' "-tetraacetic acid (DOTA) to improve the rate of metal loading for radioimmunotherapy applications. Previous ab initio calculations predicted that the compounds 1,4,7,10-tetra(carboxyethyl)-1,4,7,10-tetraazacyclododecane (DO4Pr) and 1,4,7-tris(carboxymethyl)-10-(carboxyethyl)-1,4,7,10-tetraazacyclododecane (DO3A1Pr) have a ca. 2000-fold improvement in yttrium metal loading rates compared to those of DOTA (Jang, Y. H.; Blanco, M.; Dasgupta, S.; Keire, D. A.; Shively, J. E.; Goddard, W. A., III. J. Am. Chem. Soc. 1999, 121, 6142-6151). In this study, we report the synthesis, purification, (1)H-NMR chemical shift assignments, pK(a) values, metal loading rate measurements, and additional ab initio calculations of these two compounds. The yttrium loading rates of DO3A1Pr are approximately twice those of DOTA, at pH 4.6 and 37 degrees C. The NMR data indicates that the DO4Pr analogue forms a stable type I complex but does not form a type II complex. The new ab initio calculations performed on DO4Pr and DO3A1Pr indicate that the rate-determining step is the deprotonation of the first macrocycle amine proton, not the second proton as assumed in the previous calculations. The new calculations predict an improvement in the rate of metal loading that more closely matches the experimentally observed change in the rate.     

Structure of interfacial water on fluorapatite (100) surface

The structure relaxation mechanism of the fluorapatite (100) surface under completely hydrated ambient conditions has been investigated with the grazing incidence X-ray diffraction (GIXRD) technique. Detailed information on lateral as well as perpendicular ordering corresponding to the water molecules and atomic relaxations of the (100) surface of fluorapatite (FAp) crystal was obtained from the experimental analysis of the CTR intensities. Two laterally ordered water layers are present at the water/mineral interface. The first layer consists of four water molecules located at 1.6(1) A above the relaxed fluorapatite (100) surface while the second shows the presence of only two water molecules at a distance of 3.18(10) A from the mineral surface. Thus, the first layer water molecules complete the truncated coordination sites of the topmost surface Ca atoms, while the second water layer molecules remain bonded by means of H-bonding to the first layer molecules and the surface phosphate groups. Molecular mechanics simulations using force field techniques are in good agreement with this general structural behavior determined from the experiment.     

Role for bound water and CH-pi aromatic interactions in photosynthetic electron transfer

Photosystem I (PSI) is one of two photosynthetic reaction centers present in plants, algae, and cyanobacteria and catalyzes the reduction of ferredoxin and the oxidation of cytochrome c or plastocyanin. The PSI primary chlorophyll donor, which is oxidized in the primary electron-transfer events, is a heterodimer of chl a and a' called P700. It has been suggested that protein relaxation accompanies light-induced electron transfer in this reaction center (Dashdorj, N.; Xu, W.; Martinsson, P.; Chitnis, P. R.; Savikhin, S. Biophys. J. 2004, 86, 3121. Kim, S.; Sacksteder, C. A.; Bixby, K. A.; Barry, B. A. Biochemistry 2001, 40, 15384). To investigate the details of electron transfer and relaxation events in PSI, we have employed several experimental approaches. First, we report a pH-dependent viscosity effect on P700+ reduction; this result suggests a role for proton transfer in the PSI electron-transfer reactions. Second, we find that changes in hydration alter the rate of P700+ reduction and the interactions of P700 with the protein environment. This result suggests a role for bound water in electron transfer to P700+. Third, we present evidence that deuteration of the tyrosine aromatic side chain perturbs the vibrational spectrum, associated with P700+ reduction. We attribute this result to a linkage between CH-pi interactions and electron transfer to P700+.     

Broadband dielectric and resistivity spectroscopy of WO(3) x H(2)O in the range of 10(3)-10(10) Hz: particle size effect

The interest in studying the electrical properties of WO(3) x H(2)O powders is made absolutely necessary because their infrared modulation properties depend on their morphologies and electronic populations. Broadband dielectric and resistivity spectra of WO(3) x H(2)O powders were recorded in a frequency range of 10(3)-10(10) Hz at temperatures varying between 200 and 300 K. Complex resistivity and permittivity diagrams have permitted thermal behavior of both dc-conductivity and permittivity to be obtained. A dielectric relaxation is found, attributed to water molecules motions. The role of the powder morphology has been investigated on two types of compounds: the first one being constituted by nanometric particles and the second by micrometric particles. Strong differences are observed in the thermal behaviors of the dc-conductivities (activation energies). Particle size effect is evidenced, giving rise to stronger electron localization on the nanometric particles. The permittivity values and the dynamical behavior of the structural water are also influenced by the particle size effect. A strong interaction between moving polarons and water molecules has been determined.     

CPS K 221 Model of proton longitudinal magnetic relaxation in hydrated crosslinked poly(methacrylic acid)

Proton longitudinal magnetic relaxation time (T₁) measurements have been made at 30 MHz over a wide range of temperature for crosslinked poly(methacrylic acid), PMA, hydrated with H₂O as well as with D₂O. From the point of view of nuclear magnetic relaxation, PMA hydrogel is a multiregion system in which three proton regions (a, b, c) can be distinguished. Region a is regarded as to be formed by the nonexchangeable polymer protons, region b by the protons of -COOH · H₂O combinations, and region c by the protons of remaining water molecules. Cross relaxation between polymer and water protons and a log normal distribution of correlation times have been assumed to take place. Temperature dependences of the T₁ time for the particular regions have been determined, from which the distribution width parameter, the second moment and the intramolecular proton-proton distance for sorbed water have been calculated.     

EuII-cryptate with optimal water exchange and electronic relaxation: a synthon for potential pO2 responsive macromolecular MRI contrast agents

The cryptate [EuII(2.2.2)(H2O)2]2+ displays several interesting features with respect to pO2 responsive MRI contrast agent applications: it is relatively stable against oxidation, it has two inner sphere water molecules, and the water exchange and electron spin relaxation rates are in the optimal range to attain high proton relaxivities, provided the rotation is also optimized.     

The possibility of the existence of (H2O) n − anions withn=5, 6

Structural changes that occur in cyclic and chain-like water pentamers and hexamers when an electron is added were analyzed at the unrestricted Hartree—Fock level. The vertical and adiabatic electron affinity of the oligomers, the energies of vertical detachment of an electron from the anions (VDE), and the stability of the anions against dissociation into individual water molecules and a free electron were estimated taking into account the second order perturbation theory (MP2) corrections to the energy. All water anions considered are stable against dissociation, but theirVDE values are negative, and only the chain-like hexamer anion has a value ofVDE close to zero, which means metastability of the anion. The energy of attachment of an electron to the oligomers is lower in the case of chain-like structures. The process is accompanied by structure relaxation, which is more substantial for cycles, especially for the cyclic hexamer. In the chain-like anions, the excess electron density is localized on the hydrogen nuclei of that terminal water molecule that acts as an acceptor of the H-bond proton, while in the cyclic anions it is distributed on the orbitals of those free hydrogen atoms that significantly approach each other due to structural relaxation. Translated fromIzvestiya Akademii Nauk, Seriya Khimicheskaya, No. 1, pp. 47–53, January 1997     

Quantum chemical investigation of hyperfine coupling constants on first coordination sphere water molecule of gadolinium(III) aqua complexes

Hyperfine interactions (HFI) on the nuclei of the first coordination sphere water molecules in a model [Gd(H(2)O)(8)](3+) aqua complex and in the magnetic resonance imaging contrast agent [Gd(DOTA)(H(2)O)](-) were studied theoretically. Density functional theory (DFT) calculations combined with classical molecular dynamics (MD) simulations have been used in order to take into account dynamic effects in aqueous solution. DFT relativistic calculations show a strong spin-polarization of the first coordination sphere water molecules. This spin-polarization leads to a positive (17)O isotropic hyperfine coupling constant (A(iso)((17)O) = 0.58 +/- 0.11 MHz) and to a significant increase of the effective distance (r(eff)(Gd-O) = 2.72 +/- 0.06 A) of dipolar interaction compared to the mean internuclear distance (r(Gd-O) = 2.56 +/- 0.06 A) obtained from the MD trajectory of [Gd(DOTA)(H(2)O)](-) in aqueous solution. The point-dipole model for anisotropic hyperfine interaction overestimates therefore the longitudinal relaxation rate of the (17)O nucleus by approximately 45%. The (1)H isotropic hyperfine coupling constant of the bound water molecule is predicted to be very small (A(iso)((1)H) = 0.03 +/- 0.02 MHz), and the point-dipole approximation for first coordination sphere water protons holds. The calculated hyperfine parameters are in good agreement with available experimental data.     

Long-range proton transfer in aqueous Acid-base reactions

We study the mechanism of proton transfer (PT) in the aqueous acid-base reaction between the photoacid 8-hydroxy-1,3,6-pyrenetrisulfonic acid (HPTS) and acetate by probing the vibrational resonances of HPTS, acetate, and the hydrated proton with femtosecond mid-infrared laser pulses. We find that PT takes place in a distribution of hydrogen-bound reaction complexes that differ in the number of water molecules separating the acid and the base. The number of intervening water molecules ranges from 0 to 5, which, together with a strongly distance-dependent PT rate, explains the observed highly nonexponential reaction kinetics. The kinetic isotope effect for the reaction is determined to be 1.5, indicating that tunneling does not play a significant role in the transfer of the proton. Rather, the transfer mechanism is best described in terms of the adiabatic PT picture as it has been formulated by Hynes and co-workers [Staib, A.; Borgis, D.; Hynes, J. T. J. Chem. Phys. 1995, 102, 2487. Ando, K.; Hynes, J. T. J. Phys. Chem. B 1997, 101, 10464.], where solvent fluctuations play an essential role in forming the correct hydrogen-bond configuration and solvent polarization to facilitate PT.     

Bimodal system (luminophore and paramagnetic contrastophore) derived from Ln(III) complexes based on a bipyridine-containing macrocyclic ligand

The synthesis of a new 15-membered polyaza-macrocyclic ligand L3H3, which is based on a 2,2'-bipyridine moiety and a diethylenetriaminetriacetic acid core, is reported. The lanthanide chelates of this octadentate ligand were programmed for bimodal probes, luminescent agents (Sm, Eu, Tb, Dy), and magnetic resonance imaging agents (Gd3+). The neutral 1:1 complexes with these Ln3+ ions were prepared and studied in aqueous solution by luminescence and NMR techniques. The main photophysical characteristics of these complexes (i.e., the absorption and luminescence spectra, the metal-centered lifetimes, and the overall luminescence yields, Phi) were measured. In addition, the role played by nonradiative pathways (vibrational energy transfer involving coordinated water molecules, involvement of ligand-to-metal charge-transfer excited states, or metal --> ligand back transfer) is discussed. The L3.Eu and L3.Tb complexes show very bright luminescence when photoexcited from the lowest-energy absorption band of the bipyridine chromophore. The luminescence quantum yields in an air-equilibrated water solution at room temperature are 0.10 and 0.21, respectively, despite the presence of one water molecule in the first coordination sphere of the metal ion. NMR data show that L3.Gd contains also one H2O molecule in the inner sphere. The proton longitudinal relaxivity, r1, of this complex is 3.4 s(-1) mM(-1) (0.47 T, 310 K) and the rotational correlation time, tau(R), is 57 ps (310 K). These values are comparable to those of the clinically used Gd-DTPA. Interestingly, the water exchange rate between the coordination site and the bulk solvent is slow (tau(M) = 3.5 micros at 310 K). The presence of water molecules in the second sphere and in rapid exchange with the solvent is discussed. Finally, it was found by luminescence and NMR experiments that these lanthanide complexes are stable versus transmetalation by several cations (especially Ca2+ and Zn2+) at physiological pH and have no interaction with blood proteins.     

Structure and dynamics of the hydrogen-bond network around (R,R)-pterocarpans with biological activity in aqueous solution

Molecular dynamics simulations were carried out in the presence of 2380 water molecules (TIP3P) to explore the conformational preferences of 3,9-dimethoxy-4-prenylpterocarpan (bitucarpin A) and 3,9-dihydroxy-4,8-diprenylpterocarpan (erybraedin C) and the H-bond network around them, using the empirical general AMBER force field (GAFF). Specific angle and torsional parameters have been improved in order to match the geometries of the minimum energy structures obtained from an earlier DFT/ab initio study in vacuo, taking into account a few configurations [Alagona, G.; Ghio, C.; Monti, S. Phys. Chem. Chem. Phys. 2004, 6, 2849-2857]. RESP partial charges were assigned to reproduce the electrostatic potential determined at the HF/6-31G level of theory. The analysis of trajectories allowed the conformations of bitucarpin and erybraedin as well as the distribution of water molecules around them to be elucidated. During one of the simulations only, the scaffold of erybraedin undergoes interconversion from its most stable Ht conformation to the Ot one and vice versa. Radial distribution functions, coordination numbers, and angular distributions put forward the extent of solvent structure and the hydrogen bonding behavior of their various (methoxy, hydroxyl, or ethereal) oxygen atoms. The distribution of solvent molecules in the first and second solvation shells as well as the residence times for the different solute-solvent interacting sites have been considered.