Implicit solvent simulations of DPC micelle formation

The formation of micelles by dodecylphosphocholine (DPC) is modeled by treating the surfactants in atomic detail and the solvent implicitly, in the spirit of the EEF1 solvation model for proteins. The solvation parameters of the DPC atoms are carried over from those of similar atoms in proteins. A slight adjustment of the parameters for the headgroup was found necessary for obtaining an aggregation number consistent with experiment. Molecular dynamics simulations of 960 DPC molecules at different concentrations are used to obtain the aggregation number, the micelle size distribution, and the CMC. At 20 mM concentration we obtain an aggregation number of 53-56 and a CMC of 1.25 mM, values close to the experimental ones. At 100 mM the aggregation number increases to 90. Simulations of individual micelles of varying size show that the effective energy per surfactant molecule is initially a decreasing function of aggregation number but stabilizes at about 60 molecules. The van der Waals term and the desolvation of nonpolar groups contribute to micellization, whereas the desolvation of polar groups opposes it. From the difference between the effective energy and the free energy (calculated from the CMC), the translational and rotational entropy contributions to the free energy are estimated at about 7 kcal/mol per monomer. The micelles obtained here are more irregular than those obtained in explicit water simulations. This modeling approach allows the study of larger surfactant aggregates for longer times and the extraction of thermodynamic in addition to structural information.     

Membrane protein native state discrimination by implicit membrane models

Four implicit membrane models [IMM1, generalized Born (GB)‐surface area‐implicit membrane (GBSAIM), GB with a simple switching (GBSW), and heterogeneous dielectric GB (HDGB)] were tested for their ability to discriminate the native conformation of five membrane proteins from 450 decoys generated by the Rosetta‐Membrane program. The energy ranking of the native state and Z‐scores were used to assess the performance of the models. The effect of membrane thickness was examined and was found to be substantial. Quite satisfactory discrimination was achieved with the all‐atom IMM1 and GBSW models at 25.4 Å thickness and with the HDGB model at 28.5 Å thickness. The energy components by themselves were not discriminative. Both van der Waals and electrostatic interactions contributed to native state discrimination, to a different extent in each model. Computational efficiency of the models decreased in the order: extended‐atom IMM1 > all‐atom IMM1 > GBSAIM > GBSW > HDGB. These results encourage the further development and use of implicit membrane models for membrane protein structure prediction. © 2012 Wiley Periodicals, Inc.     

Comparison of a Rapid Light-Induced and Forced Test to Study the Oxidative Stability of White Wines

The oxidation processes of white wines can occur during storage and commercialization due to several factors, and these can negatively affect the color, aroma, and quality of the wine. Wineries should have faster and simpler methods that provide valuable information on oxidation stability of wines and allow fast decision-making procedures, able to trigger suitable technological interventions. Using a portable prototype instrument for light irradiations at different wavelengths and times was considered and evaluated on sensorial, spectrophotometric, and colorimetric parameters of white wines. The sensorial analysis revealed that white and light blue were the most significant, after only 1 h of irradiation. The experimental results showed that hydrogen peroxide could enhance the effect of light treatment, allowing a contemporary evaluation of the oxidation stability of wine against light and chemical stresses. As expected, a good correlation (R² > 0.89) between optical density at 420 nm and b* parameter was highlighted. The synergic effect of light and H₂O₂ was also studied on the hydrolyzable and condensed tannins’ additions to white wine. The proposed methodology could be used to evaluate the oxidative stability of white wines, but also to evaluate the effect of some oenological adjuvants on wine stability.     

Calculations of pH-dependent binding of proteins to biological membranes

Binding of proteins to membranes is often accompanied by titration of ionizable residues and is, therefore, dependent on pH. We present a theoretical treatment and computational approach for predicting absolute, pH-dependent membrane binding free energies. The standard free energy of binding, DeltaG, is defined as -RTln(P(b)/P(f)), where P(b) and P(f) are the amounts of bound and free protein. The apparent pK(a) of binding is the pH value at which P(b) and P(f) are equal. Proteins bind to the membrane in the pH range where DeltaG is negative. The components of the binding free energy are (a) the free energy cost of ionization state changes (DeltaG(ion)), (b) the effective energy of transfer from solvent to the membrane surface, (c) the translational/rotational entropy cost of binding, and (d) an ideal entropy term that depends on the relative volume of the bound and free state and therefore depends on lipid concentration. Calculation of the first term requires determination of pK(a) values in solvent and on the membrane surface. All energies required by the method are obtained from molecular dynamics trajectories on an implicit membrane (IMM1-GC). The method is tested on pentalysine and the helical peptide VEEKS, derived from the membrane-binding domain of phosphocholine cytidylyltransferase. The agreement between the measured and the calculated free energies of binding of pentalysine is good. The extent of membrane binding of VEEKS is, however, underestimated compared to experiment. Calculations of the interaction energy between two VEEKS helices on the membrane suggest that the discrepancy is mainly due to the neglect of protein-protein interactions on the membrane surface.     

Natural time analysis of critical phenomena: The case of acoustic emissions in triaxially deformed Etna basalt

Acoustic emissions exhibit complex correlations between space, time, and magnitude. As such, they present a unique example for a complex time series. We apply the recently introduced method of natural time analysis, which enables the detection of long-range temporal correlations even in the presence of heavy tails and find that the acoustic emissions exhibits features similar to that of other equilibrium or non-equilibrium critical systems such as worldwide seismicity as presented in the Centennial earthquake catalogue which includes global seismicity events with magnitude _Mw>7.0$M_w>7.0$. These results support the universal behaviour of the fracturing processes from the laboratory to the global scale.     

Exact solutions for linear propagation of chirped pulses using a chirped Gauss Hermite orthogonal basis

A generalized solution of the linear propagation equation is proposed in terms of chirped Gauss-Hermite orthogonal functions. Some well-known special cases are pointed out, and the usefulness of this approach in analyzing arbitrarily shaped chirped pulses in rapidly converging series is discussed.     

Water at biomolecular binding interfaces

Water molecules are often found at the binding interface of biomolecular complexes mediating the interaction between polar groups via hydrogen bonds, or simply filling space providing van der Waals interactions. Recent studies have demonstrated the importance of taking such water molecules into account in docking and binding affinity prediction. Here, we review the recent experimental and theoretical work aimed at quantifying the influence of interfacial water on the thermodynamic properties of binding. We highlight especially our recent results obtained by inhomogeneous fluid solvation theory in several systems and the prediction of the thermodynamic consequences of displacement of the bound water molecule by ligand modification. Finally, we discuss possible directions for further progress in this field.     

On the non-extensive nature of the isothermal depolarization relaxation currents in cement mortars

The study of isothermal depolarization current relaxation during the cement hardening process follows a generalized exponential function, which explicitly introduces hierarchically constrained dynamics and macroscopic interactions. The interactions are associated to the non-extensive entropy parameter q, which varies during hardening evolution. Combining ideas of Levy and Tsallis statistics we argue that a Levy walk type mechanism can organize the geometry of the cement mortar heterogeneous system. The estimation of non-extensive q parameter safe concludes that cement's hardening is a subextensive process.     

Potentials of mean force between ionizable amino acid side chains in water

Potentials of mean force (PMF) between all possible ionizable amino acid side chain pairs in various protonation states were calculated using explicit solvent molecular dynamics simulations with umbrella sampling and the weighted histogram analysis method. The side chains were constrained in various orientations inside a spherical cluster of 200 water molecules. Beglov and Roux's Spherical Solvent Boundary Potential was used to account for the solvent outside this sphere. This approach was first validated by calculating PMFs between monatomic ions (K(+), Na(+), Cl(-)) and comparing them to results from the literature and results obtained using Ewald summation. The strongest interaction (-4.5 kcal/mol) was found for the coaxial Arg(+).Glu(-) pair. Many like-charged side chains display a remarkable lack of repulsion, and occasionally a weak attraction. The PMFs are compared to effective energy curves obtained with common implicit solvation models, namely Generalized Born (GB), EEF1, and uniform dielectric of 80. Overall, the EEF1 curves are too attractive, whereas the GB curves in most cases match the minima of the PMF curves quite well. The uniform dielectric model, despite some fortuitous successes, is grossly inadequate.