Effect of ferric citrate on amyloid-beta peptides behavior

Amyloid beta (Aβ) aggregation and oxidative stress are two of the central events in Alzheimer's Disease (AD). Both these phenomena can be caused by the interaction of Aβ with metal ions. In the last years the interaction between Zn^II , Cu^II , and Aβ was much studied, but between iron and Aβ it is still little known. In this work we determine how three Aβ peptides, present in AD, interact with Fe^III -citrate. The three Aβ peptides are: full length Aβ1-42, an isoform truncated at Glutamic acid in position three, Aβ3-42, and its pyroglutamated form AβpE3-42. Conformation and morphology of the three peptides, aggregated with and without Fe^III -citrate were studied. Besides, we have determined the strength of the interactions Aβ/Fe^III -citrate studying the effect of ethylenediaminetetraacetic acid as chelator. Results reported here demonstrate that Fe^III -citrate promotes the aggregation in all the three peptides. Moreover, Aspartic acid 1, Glutamic acid 3, and Tyrosine 10 have an important role in the coordination with iron, generating a more stable complex for Aβ1-42 compared to that for the truncated peptides.

Electrostatic theory of the assembly of PAMAM dendrimers and DNA

The electrostatic interactions mediated by counterions between a cationic PAMAM dendrimer, modelized as a sphere of radius and cationic surface charge highly increasing with generation, and a DNA, modelized as an anionic elastic line, are analytically calculated in the framework of condensation theory. Under these interactions the DNA is wrapped around the sphere. For excess phosphates relative to dendrimer primary amines, the free energy of the DNA-dendrimer complex displays an absolute minimum when the complex is weakly negatively overcharged. This overcharging opposes gene delivery. For a highly positive dendrimer and a DNA fixed by experimental conditions to a number of phosphates less than the number of dendrimer primary amines, excess amine charges, the dendrimer may at the same time bind stably DNA and interact with negative cell membranes to activate cell transfection in fair agreement with molecular simulations and experiments.

A critical concentration of N-terminal pyroglutamylated amyloid beta drives the misfolding of Ab1-42 into more toxic aggregates

A wide consensus based on robust experimental evidence indicates pyroglutamylated amyloid-β isoform (AβpE3-42) as one of the most neurotoxic peptides involved in the onset of Alzheimer's disease. Furthermore, AβpE3-42 co-oligomerized with excess of Aβ1-42, produces oligomers and aggregates that are structurally distinct and far more cytotoxic than those made from Aβ1-42 alone. Here, we investigate quantitatively the influence of AβpE3-42 on biophysical properties and biological activity of Aβ1-42. We tested different ratios of AβpE3-42/Aβ1-42 mixtures finding a correlation between the biological activity and the structural conformation and morphology of the analyzed mixtures. We find that a mixture containing 5% AβpE3-42, induces the highest disruption of intracellular calcium homeostasis and the highest neuronal toxicity. These data correlate to an high content of relaxed antiparallel β-sheet structure and the coexistence of a population of big spheroidal aggregates together with short fibrils. Our experiments provide also evidence that AβpE3-42 causes template-induced misfolding of Aβ1-42 at ratios below 33%. This means that there exists a critical concentration required to have seeding on Aβ1-42 aggregation, above this threshold, the seed effect is not possible anymore and AβpE3-42 controls the total aggregation kinetics.

Lamellar cationic lipid-DNA complexes from lipids with a strong preference for planar geometry: A Minimal Electrostatic Model

We formulate and analyze a minimal model, based on condensation theory, of the lamellar cationic lipid (CL)-DNA complex of alternately charged lipid bilayers and DNA monolayers in a salt solution. Each lipid bilayer, composed by a random mixture of cationic and neutral lipids, is assumed to be a rigid uniformly charged plane. Each DNA monolayer, located between two lipid bilayers, is formed by the same number of parallel DNAs with a uniform separation distance. For the electrostatic calculation, the model lipoplex is collapsed to a single plane with charge density equal to the net lipid and DNA charge. The free energy difference between the lamellar lipoplex and a reference state of the same number of free lipid bilayers and free DNAs, is calculated as a function of the fraction of CLs, of the ratio of the number of CL charges to the number of negative charges of the DNA phosphates, and of the total number of planes. At the isoelectric point the free energy difference is minimal. The complex formation, already favoured by the decrease of the electrostatic charging free energy, is driven further by the free energy gain due to the release of counterions from the DNAs and from the lipid bilayers, if strongly charged. This minimal model compares well with experiment for lipids having a strong preference for planar geometry and with major features of more detailed models of the lipoplex.

Differential toxicity, conformation and morphology of typical initial aggregation states of Aβ1-42 and Aβpy3-42 beta-amyloids

Among the different species of water-soluble β-peptides (Aβ1-42, Aβ1-40 and N-terminal truncated Aβ-peptides), Aβpy3-42 is thought to play a relevant role in Alzheimer's pathogenesis due to its abundance, resistance to proteolysis, fast aggregation kinetics, dynamic structure and high neurotoxicity. To evaluate the specific structural characteristics and neurotoxicity of Aβpy3-42, we separated different aggregation states of Aβ1-42 and Aβpy3-42 using fast protein liquid chromatography, isolating in both cases three peaks that corresponded to sa (small), ma (medium) and la (large) aggregates. Conformational analysis, by circular dichroism showed a prevailing random coil conformation for sa and ma, and typical β-sheet conformation for la. AFM and TEM show differential structural features between the three aggregates of a given β-peptide and among the aggregate of the two β-peptides. The potential toxic effects of the different aggregates were evaluated using human neuroblastoma SH-SY5Y cells in the MTT reduction, in the xCELLigence System, and in the Annexin V binding experiments. In the case of Aβ1-42 the most toxic aggregate is la, while in the case of Aβpy3-42 both sa and la are equally toxic. Aβ aggregates were found to be internalized in the cells, as estimated by confocal immunofluorescence microscopy, with a higher effect observed for Aβpy3-42, showing a good correlation with the toxic effects. Together these experiments allowed the discrimination of the intermediate states more responsible of oligomer toxicity, providing new insights on the correlation between the aggregation process and the toxicity and confirming the peculiar role in the pathogenesis of Alzheimer disease of Aβpy3-42 peptide.

[Biogenic amines in epigeal spontaneous mushrooms: indicators of quality and freshness?]

Purpose of the present study is to determine the qualitative and quantitative composition of nine biogenic amines (BA) in 153 samples of epigeal spontaneous mushrooms (Boletus edulis and Amanita Caesarea), sampled at large and small food retail shops in the Florence area. Each BA has been identified and quantified using a liquid chromatography system with a reversed phase, C18 column and post-column derivatization with o-oftalaldeide. The results have shown a widespread presence of each amine, although with varying frequency and concentration, and an increasing concentration of tyramine and putrescine in samples with poor organoleptic quality than those in good conditions.

Wrapped-around models for the lac operon complex

The protein-DNA complex, involved in the lac operon of enteric bacteria, is paradigmatic in understanding the extent of DNA bending and plasticity due to interactions with protein assemblies acting as DNA regulators. For the lac operon, two classes of structures have been proposed: 1), with the protein tetramer lying away from the DNA loop (wrapped-away model); and 2), with the protein tetramer lying inside the DNA loop (wrapped-around model). A recently developed electrostatic analytical model shows that the size and net charge of the Lac protein tetramer allow the bending of DNA, which is consistent with another wrapped-around model from the literature. Coarse-grained models, designed based on this observation, are extensively investigated and show three kinds of wrapped-around arrangements of DNA and a lower propensity for wrapped-away configurations. Molecular dynamics simulations of an all-atom model, built on the basis of the most tightly collapsed coarse-grained model, show that most of the DNA double-helical architecture is maintained in the region between O3 and O1 DNA operators, that the DNA distortion is concentrated in the chain beyond the O1 operator, and that the protein tetramer can adapt the N-terminal domains to the DNA tension.

N-terminal truncated pyroglutamyl beta amyloid peptide Abetapy3-42 shows a faster aggregation kinetics than the full-length Abeta1-42

We tested directly the differences in the aggregation kinetics of three important beta amyloid peptides, the full-length Abeta1-42, and the two N-terminal truncated and pyroglutamil modified Abetapy3-42 and Abetapy11-42 found in different relative concentrations in the brains in normal aging and in Alzheimer disease. By following the circular dichroism signal and the ThT fluorescence of the solution in phosphate buffer, we found substantially faster aggregation kinetics for Abetapy3-42. This behavior is due to the particular sequence of this peptide, which is also responsible for the specific oligomeric aggregation states, found by TEM, during the fibrillization process, which are very different from those of Abeta1-42, more prone to fibril formation. In addition, Abetapy3-42 is found here to have an inhibitory effect on Abeta1-42 fibrillogenesis, coherently with its known greater infective power. This is an indication of the important role of this peptide in the aggregation process of beta-peptides in Alzheimer disease.

Reference values of urinary trans,trans-muconic acid: Italian Multicentric Study

This article reports the results of a study, conducted in the framework of the scientific activities of the Italian Society for Reference Values, aimed at defining reference values of urinary trans,trans-muconic acid (t,t-MA) in the general population not occupationally exposed to benzene. t,t-MA concentrations detected in 376 subjects of the resident population in three areas of Italy, two in central (Florence and southern Tuscany) and one in northern Italy (Padua), by three laboratories, compared by repeated interlaboratory controls, showed an interval of 14.4-225.0 microg/L (5th-95th percentile) and a geometric mean of 52.5 microg/L. The concentrations measured were influenced by tobacco smoking in a statistically significant way: Geometric mean concentrations were 44.8 microg/L and 76.1 microg/Ll in nonsmokers (264 subjects) and smokers (112 subjects), respectively. In the nonsmoking population, a significant influence of gender was found when concentrations were corrected for urinary creatinine, geometric mean concentrations being 36.7 microg/g creatinine in males (128 subjects) and 44.7 microg/g creatinine in females (136 subjects). The place of residence of subjects did not seem to influence urinary excretion of the metabolite, although personal inhalation exposure to benzene over a 24-h period showed slightly higher concentrations in Padua and Florence (geometric means of 6.5 microg/m(3) and 6.6 microg/m(3), respectively) than in southern Tuscany (geometric mean of 3.9 microg/m(3)). Concentration of t,t-MA in urine samples collected at the end of personal air sampling showed little relationship to personal inhalation exposure to benzene, confirming the importance of other factors in determining excretion of t,t-MA when concentrations in personal air samples are very low.

Intracellular accumulation of a mild-denatured monomer of the human PrP fragment 90-231, as possible mechanism of its neurotoxic effects

Because of high tendency of the prion protein (PrP) to aggregate, the exact PrP isoform responsible for prion diseases as well as the pathological mechanism that it activates remains still controversial. In this study, we show that a pre-fibrillar, monomeric or small oligomeric conformation of the human PrP fragment 90-231 (hPrP90-231), rather than soluble or fibrillar large aggregates, represents the neurotoxic species. In particular, we demonstrate that monomeric mild-denatured hPrP90-231 (incubated for 1 h at 53 degrees C) induces SH-SY5Y neuroblastoma cell death, while, when structured in large aggregates, it is ineffective. Using spectroscopic and cellular techniques we demonstrate that this toxic conformer is characterized by a high exposure of hydrophobic regions that favors the intracellular accumulation of the protein. Inside the cells hPrP90-231 is mainly compartmentalized into the lysosomes where it may trigger pro-apoptotic 'cell death' signals. The PrP toxic conformation, which we have obtained inducing a controlled in vitro conformational change of the protein, might mimic mild-unfolding events occurring in vivo, in the presence of specific mutations, oxidative reactions or proteolysis. Thus, in light of this model, we propose that novel therapeutic strategies, designed to inhibit the interaction of the toxic PrP with the plasmamembrane, could be beneficial to prevent the formation of intracellular neurotoxic aggregates and ultimately the neuronal death.

Generalized electrostatic model of the wrapping of DNA around oppositely charged proteins

Histonelike proteins in prokaryotes and histone octamers in eukaryotes carry large positive charges, which are responsible of strong electrostatic interactions with DNA. As a result, DNA wraps around proteins and genetic information is condensed. We describe a generalized model of these electrostatic interactions mediated by salt that explains the wrapping of DNA around the nucleosome octamer, around remodeling factors in eukaryotes and around histonelike proteins in prokaryotes. It comes out that small changes in protein dimension and charge produce large effects in the supramolecular DNA-protein architecture.

Most sensitive urinary cotinine cut-off level for environmental tobacco smoke exposure assessment: a pilot study

The aim of this pilot study was to determine the most sensitive urinary cotinine level able to assess environmental tobacco smoke (ETS) exposure. 54 Florentine subjects (29 males and 25 females), reporting to be nonsmokers and exposed (E) or not exposed (NE) to ETS at home, at work or in places of recreation, were examined. The urinary cotinine concentration was determined using gaschromatographic analysis in samples collected on three consecutive days. 18 subjects (33.3%) reported to be exposed to ETS had a greater median cotinine concentration than 36 ETS-NE subjects (E = 3.3 pg/L vs NE = 2.2 microg/L, median values), with borderline statistical significance (P = 0.05). The 2.5 microg/L cotinine concentration was the only statistically significant cut-off (P = 0.04) discriminating between ETS-E to ETS-NE subjects, identifying 51.9% of the subjects examined as exposed (E). Considering the expanded uncertainty of measurement of the method used (20%), urinary cotinine concentrations higher than 3.1 microg/L, a value whose confidence interval is higher than our proposed cut-off of 2.5 microg/L, mean that to be sure that a subject is exposed to ETS.

Modeling H3 histone N-terminal tail and linker DNA interactions

Molecular dynamics computer simulations were performed for the 25-residue N-terminal tail of the H3 histone protein in the proximity of a DNA segment of 10 base pairs (bp), representing a model for the linker DNA in chromatin. Several least biased configurations were used as initial configurations. The secondary structure content of the protein was increased by the presence of DNA close to it, but the locations of the secondary motifs were different for different initial orientations of the DNA grooves with respect to the protein. As a common feature to all simulations, the electrostatic attraction between negatively charged DNA and positively charged protein was screened by the water solvent and counterbalanced by the intrinsic compaction of the protein due to hydrophobic effects. The protein secondary structure limited the covering of DNA by the protein to 4-5 bp. The degree of compaction and charge density of the bound protein suggests a possible role of H3 tail in a nonspecific bending and plasticity of the linker DNA when the protein is located in the crowded dense chromatin.

Hyaluronan chain conformation and dynamics

An overview of the present state of research in the field of hyaluronan chain conformational aspects is presented. The relationship between structure and dynamics are illustrated for a series of hyaluronan oligomers. Conformational characteristics of hyaluronan chains are discussed, together with the dynamic chain patterns, evaluated by using a theoretical approach to diffusive polymer dynamics. The dependence of correlation times and NMR relaxation parameters from the chain dimension are investigated. Topological features and dimensional properties are related to the structural determinants by using classical computational methods of molecular mechanics and Monte Carlo simulation.

Modeling the backbone dynamics of reduced and oxidized solvated rat microsomal cytochrome b5

In this article, a description of the statistics and dynamics of cytochrome b(5) in both reduced and oxidized forms is given. Results of molecular dynamics computer simulations in the explicit solvent have been combined with mode-coupling diffusion models including and neglecting the molecule-solvent correlations. R(1) and R(1 rho) nuclear magnetic relaxation parameters of (15)N in the protein backbone have been calculated and compared with experiments. Slight changes in charge density in the heme upon oxidation produces a cascade of changes in charge distributions from heme propionates up to charged residues approximately 1.5 nm from Fe. These changes in charge distributions modify the molecular surface and the water shell surrounding the protein. The statistical changes upon oxidation can be included in diffusive models that physically explain the upper and lower limits of R(1 rho) relaxation parameters at high off-resonance fields.

Electrostatic interactions with histone tails may bend linker DNA in chromatin

Is linker DNA bent in the 30-nm chromatin fiber at physiological conditions? We show here that electrostatic interactions between linker DNA and histone tails including salt condensation and release may bend linker DNA, thus affecting the higher order organization of chromatin.

Designing generalized statistical ensembles for numerical simulations of biopolymers

Conformational properties of polymers, such as average dihedral angles or molecular alpha-helicity, display a rather weak dependence on the detailed arrangement of the elementary constituents (atoms). We propose a computer simulation method to explore the polymer phase space using a variant of the standard multicanonical method, in which the density of states associated to suitably chosen configurational variables is considered in place of the standard energy density of states. This configurational density of states is used in the Metropolis acceptance/rejection test when configurations are generated with the help of a hybrid Monte Carlo algorithm. The resulting configurational probability distribution is then modulated by exponential factors derived from the general principle of the maximal constrained entropy by requiring that certain average configurational quantities take preassigned (possibly temperature dependent) values. Thermal averages of other configurational quantities can be computed by using the probability distributions obtained in this way. Moments of the energy distribution require an extra canonical sampling of the system phase space at the desired temperature, in order to locally thermalize the configurational degrees of freedom. As an application of these ideas we present the study of the structural properties of two simple models: a bead-and-spring model of polyethylene with independent hindered torsions and an all-atom model of alanine and glycine oligomers with 12 amino acids in vacuum.

Determination of N-acetyl-S-(N-methylcarbamoyl)cysteine (AMCC) in urine by high-performance liquid chromatography combined with mass spectrometry

A simple, sensitive and specific method for the quantification of N-acetyl-S-(N-methylcarbamoyl)cysteine (AMCC) in urines of the general population and of exposed workers has been developed. Samples were first diluted with phosphate buffer followed by purification by solid-phase extraction (SPE) with SAX columns prior to analysis by liquid chromatography/mass spectrometry (LC/MS). An external standard was used for the quantification together with selected ion monitoring of m/z 219, [M-H]-. The linear calibration curve showed a good correlation (r2 = 0.98, p < 0.001) between 0.1-110 mg/L of AMCC, and the detection limit was calculated to be 2 microg/L. The within-day and between-day precision, calculated for exposed workers and general population AMCC levels in urine samples, were determined as 1.2 and 3.6%, respectively (as relative standard deviation (RSD)).

[Evaluation of professional exposure to antiblastic chemotherapeutic agents in a Tuscan hospital environment]

The aim of this study, carried out from 1997 to 2000, was to evaluate the exposure to antiblastic drugs of hospital staff (38 nurses of 7 oncological Day Hospital care in Tuscany). To evaluate the internal exposure was used biological monitoring, in particular was determined the concentration of urinary cyclophosfamide (CF); pad was made to evaluate the skin contamination. The contamination of working surfaces was determinate, by wipe-test, to verify the decontamination procedures; were used, as markers, CF and 5-fluorouracil (5-FU). We assessed the permeability thorough gloves usually used by the staff and the degradative activity to these drugs by agents used to decontaminate the working surfaces. The results shows, in urine, value lower than the detection limits, although was found a diffused contamination of the working environment. The results confirmed that NaClO is the best product to decontaminate working surfaces and nitrile gloves the ones with less permeability.

Modeling the dynamics of the solvated SL1 domain of HIV-1 genomic RNA

A mode-coupling solution of the Smoluchowski diffusion equation (MCD theory), designed to describe the dynamics of wobbling macromolecules in water, is applied to a macromolecular bead model including water beads in the nearest layers. The necessary statistical averages are evaluated by time averaging along a molecular dynamics (MD) trajectory where both solute and water are introduced as atomistic models. The cross peaks in (1)H nuclear Overhauser effect spectroscopy (NOESY) NMR spectra that are routinely measured to determine biological structures are here calculated for the mutated 23 nucleotides stem-loop fragment of the SL1 domain in the HIV-1(Lai) genomic RNA. The calculations are in acceptable agreement with experiments without requiring any screening of the hydrodynamic interactions. The screening of hydrodynamics was necessary in previous MCD calculations obtained by using the same full atomistic MD trajectory, but a nonsolvated frictional model.

[Reference values of solvents and metabolites in biological samples]

This article presents a review of reference values for organic solvent biomarkers. Some of these results were obtained in the research activities of the Italian Society Reference Values (SIVR). Most experiences show data obtained from control groups during occupational exposures assessment investigations. We considered only data related to the following biomarkers: immodified solvents in blood and urine, their main urinary metabolites. The reference values of the following solvents are reported: benzene, toluene, xylene, nhexane, cyclohexane, trichloroethylene, perchloroethylene, methanol, acetone, N,N-dimethylformamide, carbon disulphide. In the text also the influence of some confounding factors is discussed.

Assessment of exposure to polycyclic aromatic hydrocarbons in police in Florence, Italy, through personal air sampling and biological monitoring of the urinary metabolite 1-hydroxypyrene

In this study, the authors evaluated exposure to airborne polycyclic aromatic hydrocarbons (PAHs) in workers exposed to exhaust gas from cars, and they assessed the efficiency of urinary 1-hydroxypyrene as an indicator of exposure to pyrene and PAHs. The authors selected 2 groups of police who worked in 2 areas in the city of Florence: 1 group was highly exposed to high-density traffic emissions during the winter and summer of 1997, and the 2nd group experienced low exposure to traffic emissions during the same period. Ambient monitoring was achieved with personal sampling of airborne PAHs during each workshift. Eight hydrocarbons were used as indicators of pollution caused by PAHs (e.g., pyrene, benzo[a]pyrene, benzo[a]anthracene, dibenzo[a,h]anthracene). Biological monitoring was performed through dosing of 1-hydroxypyrene (pyrene metabolite) in urine samples taken at the end of each workshift. The ambient monitoring revealed that PAH concentrations were influenced by both season of sampling and varying intensities of traffic in the different areas. The median concentration of benzo[a]pyrene in winter was twice as high in the high-density traffic area as in the low-density traffic area (i.e., 4.1 ng/m3 versus 1.8 ng/m3). In summer, the high-density traffic area experienced benzo[a]pyrene concentrations that were 6 times higher than in the low-density traffic area (i.e., 1.2 ng/m3 versus 0.2 ng/m3). Benzo[a]pyrene was also correlated highly (r(s) = .92, p < .0001) with the mixture of total PAHs analyzed, thus confirming its function as a good indicator of exposure to PAHs in an urban environment. Levels of urinary 1-hydroxypyrene appeared to be generally influenced by the intensity of traffic, especially during the winter (i.e., median value in winter was 199.2 ng/gm creatinine in the high-density traffic area and 120.5 ng/gm creatinine in the low-density traffic area). An analysis of the general data revealed that 1-hydroxypyrene was, to some degree, related to pyrene, benzo[a]pyrene, and airborne total PAHs, whereas analysis of separate data for the area and the season revealed an emergence of a closer correlation during the winter in the high-traffic area. Therefore, 1 -hydroxypyrene can be considered a good biological indicator of exposure to airborne PAHs in the urban environment, especially in winter and in high-density traffic areas.

Backbone dynamics for the wild type and a double H52R/T56W mutant of the vnd/NK-2 homeodomain from Drosophila melanogaster

The (15)N relaxation behavior and heteronuclear Overhauser effect data for the wild type and an H52R/T56W double mutant protein that encompasses the vnd/NK-2 homeodomain from Drosophila melanogaster were used to characterize and describe the protein backbone dynamics. This investigation, which includes a description of a model structure for the H52R/T56W double mutant vnd/NK-2 homeodomain, was carried out for the two proteins in both the free and DNA-bound states. The double residue replacement at positions 52 and 56 within the DNA recognition helix of vnd/NK-2 has been shown to lead to a significant secondary structural modification resulting in an increase in the length of the recognition helix for the unbound protein. These structural changes are accompanied by corresponding changes in the T(1) and T(1)(rho) relaxation times as well as in the heteronuclear Overhauser effect (XNOE) values that show that the structural stability of the protein is enhanced by the two residue replacements. The values of the rotational anisotropy, D(parallel)/D(perpendicular), derived from analysis of the (15)N T(1) and T(1)(rho) relaxation values are small (1.189 for the unbound homeodomain and 1.110 for the bound homeodomain; both analyzed as prolate ellipsoids of revolution). A comparison of the T(2) values of the wild type and double mutant homeodomain reveals the presence of a low-frequency exchange contribution for the wild type analogue. These relaxation studies show that the motional behavior of the protein primarily reflects the tertiary structure and stability of the homeodomain backbone as well as the respective changes induced upon site-directed residue replacement or DNA binding.

Modeling the dynamics of a mutated stem-loop in the SL1 domain of HIV-1Lai genomic RNA by 1H-NOESY spectra

The cross-peaks of 1H-NOESY spectra at different time delays are compared to a mode-coupling diffusion (MCD) calculation, including the evaluation of the full 1H relaxation matrix, in the case of a 23 nucleotide fragment of the stem-loop SL1 domain of HIV-1Lai genomic RNA mutated in a single position. The MCD theory gives significant agreement with 1H relaxation experiments enabling a thorough understanding of the differential local dynamics along the sequence and particularly of the dynamics of nucleotides in the stem and in the loop. The differential dynamics of this hairpin structure is important in directing the dimerization of the retroviral genome, a fundamental step in the infectious process. The demonstration of a reliable use of time dependent NOE cross-peaks, largely available from NMR solution structure determination, coupled to MCD analysis, to probe the local dynamics of biological macromolecules, is a result of general interest of this paper.

Smoluchowski dynamics of the vnd/NK-2 homeodomain from Drosophila melanogaster: second-order maximum correlation approximation

The mode coupling diffusion theory is applied to the derivation of local dynamics in proteins in solution. The rotational dynamics of the bonds along the protein sequence are calculated and compared to the experimentally measured nmr (15)N spin-lattice relaxation time T(1), at 36.5, 60.8, and 81.1 MHz of the vnd/NK-2 homeodomain from Drosophila melanogaster. The starting point for the calculations is the experimental three-dimensional solution structure of the homeodomain determined by multidimensional nmr spectroscopy. The higher order mode-coupling computations are compared also with the recently published first-order approximation calculations. The more accurate calculations improve substantially the first-order ORZLD calculations and show that the role of the strength of the hydrodynamic interactions becomes crucial to fix the order of magnitude of the rotational dynanics for these very compact molecules characterized by partial screening of the internal atoms to water. However, the relative mobility of the bonds along the sequence and the differential fluctuations depend only weakly on the hydrodynamic strength but strongly on the geometry of the three-dimensional structure and on the statistics incorporated into the theory. Both rigid and fluctuating dynamic models are examined, with fluctuations evaluated using molecular dynamics simulations. The comparison with nmr data shows that mode coupling diffusion accounts for the T(1) relaxation pattern at low frequency where the rotational tumbling dominates. An important contribution of internal motions in the nanosecond time scale is seen at high frequencies and is discussed in terms of diffusive concepts.

Dynamics of a double stranded DNA oligomer: mode-coupling diffusion approach and reduced rigid fragment models

The local dynamics of a double stranded DNA fragment [d(CpGpCpApApApTpTpTpGpCpG)]2 of twelve base pairs is obtained to second order in the mode-coupling expansion of the Smoluchowski diffusion theory. The DNA is considered a fluctuating three-dimensional (3D) structure undergoing rotational diffusion. The starting structure for the calculations is the B canonical structure of the fragment, while the fluctuations are evaluated using molecular dynamics simulations, with the ensemble averages approximated by time averages along a trajectory of length 1.5 ns. The rotational dynamics of the bonds along the double strands are calculated and compared to experimental NMR relaxation rates of different 13C along the sequence: R(Cz), R(Cxy) and R(Hz-->Cz). For a fluctuating 3D structure the mode-coupling diffusion theory is found to be in good agreement with several relative characteristics of the experimental relaxation parameters, while motivations are given for the few differences which are due mainly to poor statistics or to inaccuracies in the diffusion model. With a view to application to larger DNA fragments, discussion is dedicated to the validity of reducing the number of degrees of freedom in the double helix statistics by grouping the atoms in rigid fragments (e.g. the backbone atoms, the sugar atoms and the base atoms of each nucleotide). Consideration is given to the effect on local dynamics properties of reduced descriptions that include only three or four rigid bodies per nucleotide as well as five rigid bodies per base pair. It is found that in general these approximations almost uniformly produce slight increase in the correlation time pattern, which grows as the rigidity in the model increases. The relative effects on the dynamic pattern for the most accurate rigid body models are modest. The errors in C1' and C5' mobilities are more significant if C5' is included in the backbone rigid body. These results offer new tools to analyse NMR relaxation behaviour and new perspectives in studying the role of dynamics in biological macromolecules.