No Arabic abstract
In the past couple of years several studies have shown that hybridization in Affymetrix DNA microarrays can be rather well understood on the basis of simple models of physical chemistry. In the majority of the cases a Langmuir isotherm was used to fit experimental data. Although there is a general consensus about this approach, some discrepancies between different studies are evident. For instance, some authors have fitted the hybridization affinities from the microarray fluorescent intensities, while others used affinities obtained from melting experiments in solution. The former approach yields fitted affinities that at first sight are only partially consistent with solution values. In this paper we show that this discrepancy exists only superficially: a sufficiently complete model provides effective affinities which are fully consistent with those fitted to experimental data. This link provides new insight on the relevant processes underlying the functioning of DNA microarrays.
We analyze publicly available data on Affymetrix microarrays spike-in experiments on the human HGU133 chipset in which sequences are added in solution at known concentrations. The spike-in set contains sequences of bacterial, human and artificial origin. Our analysis is based on a recently introduced molecular-based model [E. Carlon and T. Heim, Physica A 362, 433 (2006)] which takes into account both probe-target hybridization and target-target partial hybridization in solution. The hybridization free energies are obtained from the nearest-neighbor model with experimentally determined parameters. The molecular-based model suggests a rescaling that should result in a collapse of the data at different concentrations into a single universal curve. We indeed find such a collapse, with the same parameters as obtained before for the older HGU95 chip set. The quality of the collapse varies according to the probe set considered. Artificial sequences, chosen by Affymetrix to be as different as possible from any other human genome sequence, generally show a much better collapse and thus a better agreement with the model than all other sequences. This suggests that the observed deviations from the predicted collapse are related to the choice of probes or have a biological origin, rather than being a problem with the proposed model.
In a recent paper [Phys. Rev. E 68, 011906 (2003)], Naef and Magnasco suggested that the bright mismatches observed in Affymetrix microarray experiments are caused by the fluorescent molecules used to label RNA target sequences, which would impede target-probe hybridization. Their conclusion is based on the observation of unexpected asymmetries in the affinities obtained by fitting microarray data from publicly available experiments. We point out here that the observed asymmetry is due to the inequivalence of RNA and DNA, and that the reported affinities are consistent with stacking free energies obtained from melting experiments of unlabeled nucleic acids in solution. The conclusion of Naef and Magnasco is therefore based on an unjustified assumption.
We have developed a global analysis model for randomly oriented, fully hydrated inverted hexagonal (H$_text{II}$) phases formed by many amphiphiles in aqueous solution, including membrane lipids. The model is based on a structure factor for hexagonally packed rods and a compositional model for the scattering length density (SLD) enabling also the analysis of positionally weakly correlated H$_text{II}$ phases. For optimization of the adjustable parameters we used Bayesian probability theory, which allows to retrieve parameter correlations in much more detail than standard analysis techniques, and thereby enables a realistic error analysis. The model was applied to different phosphatidylethanolamines including previously not reported H$_text{II}$ data for diC14:0 and diC16:1 phosphatidylethanolamine. The extracted structural features include intrinsic lipid curvature, hydrocarbon chain length and area per lipid at the position of the neutral plane.
Many biological electron transfer (ET) reactions are mediated by metal centres in proteins. NADH:ubiquinone oxidoreductase (complex I) contains an intramolecular chain of seven iron-sulphur (FeS) clusters, one of the longest chains of metal centres in biology and a test case for physical models of intramolecular ET. In biology, intramolecular ET is commonly described as a diffusive hopping process, according to the semi-classical theories of Marcus and Hopfield. However, recent studies have raised the possibility that non-trivial quantum mechanical effects play a functioning role in certain biomolecular processes. Here, we extend the semi-classical model for biological ET to incorporate both semi-classical and coherent quantum phenomena using a quantum master equation based on the Holstein Hamiltonian. We test our model on the structurally-defined chain of FeS clusters in complex I. By exploring a wide range of realistic parameters we find that, when the energy profile for ET along the chain is relatively flat, just a small coherent contribution can provide a robust and significant increase in ET rate (above the semi-classical diffusive-hopping rate), even at physiologically-relevant temperatures. Conversely, when the on-site energies vary significantly along the chain the coherent contribution is negligible. For complex I, a crucial respiratory enzyme that is linked to many neuromuscular and degenerative diseases, our results suggest a new contribution towards ensuring that intramolecular ET does not limit the rate of catalysis. For the emerging field of quantum biology, our model is intended as a basis for elucidating the general role of coherent ET in biological ET reactions.
We present a structural data set of the 20 proteinogenic amino acids and their amino-methylated and acetylated (capped) dipeptides. Different protonation states of the backbone (uncharged and zwitterionic) were considered for the amino acids as well as varied side chain protonation states. Furthermore, we studied amino acids and dipeptides in complex with divalent cations (Ca2+, Ba2+, Sr2+, Cd2+, Pb2+, and Hg2+). The database covers the conformational hierarchies of 280 systems in a wide relative energy range of up to 4 eV (390 kJ/mol), summing up to an overall of 45,892 stationary points on the respective potential-energy surfaces. All systems were calculated on equal first-principles footing, applying density-functional theory in the generalized gradient approximation corrected for long-range van der Waals interactions. We show good agreement to available experimental data for gas-phase ion affinities. Our curated data can be utilized, for example, for a wide comparison across chemical space of the building blocks of life, for the parametrization of protein force fields, and for the calculation of reference spectra for biophysical applications.