Interprotein electron transfer in a confined space: uncoupling protein dynamics from electron transfer by sol-gel encapsulation

In this paper, we describe the first observations of photoinitiated interprotein electron transfer (ET) within sol-gels. We have encapsulated three protein-protein complexes, specifically selected because they represent a full range of affinities, are sensitive to different types of dynamic processes, and thus are expected to respond differently to sol-gel encapsulation. The three systems are (i) the [Zn, Fe(3+)L] mixed-metal hemoglobin hybrids, where the alpha(1)-Zn and beta(2)-Fe subunits correspond to a "predocked" protein-protein complex with a crystallographically defined interface (Natan, M. J.; Baxter, W. W.; Kuila, D.; Gingrich, D. J.; Martin, G. S.; Hoffman, B. M. Adv. Chem. Ser. 1991, 228 (Electron-Transfer Inorg., Org., Biol. Syst.), 201-213), (ii) the Zn-cytochrome c peroxidase complex with cytochrome c, [ZnCcP, Fe(3+)Cc], having an intermediate affinity between its partners (Nocek, J. M.; Zhou, J. S.; De Forest, S.; Priyadarshy, S.; Beratan, D. N.; Onuchic, J. N.; Hoffman, B. M. Chem. Rev. 1996, 96, 2459-2489), and (iii) the [Zn-deuteromyoglobin, ferricytochrome b(5)] complex, [ZnDMb, Fe(3+)b(5)], which is loosely bound and highly dynamic (Liang, Z.-X.; Nocek, J.; Huang, K.; Hayes, R. T.; Kurnikov, I. V.; Beratan, D. N.; Hoffman, B. M. J. Am. Chem. Soc. 2002, 124, 6849-6859. Intersubunit ET within the hybrid does not involve second-order processes or subunit rearrangements, and thus is influenced only by perturbations of high-frequency motions coupled to ET. For the latter two complexes, sol-gel encapsulation eliminates second-order processes: protein partners encapsulated as a complex must stay together throughout a photoinitiated ET cycle, while proteins encapsulated alone cannot acquire a partner. It further modulates intracomplex motions of the two partners.     

Beating cancer in multiple ways using nanogold

Gold nanoparticles possess a unique combination of properties which allow them to act as highly multifunctional anti-cancer agents (X. H. Huang, P. K. Jain, I. H. El-Sayed and M. A. El-Sayed, Nanomedicine, 2007, 2, 681-693; P. Ghosh, G. Han, M. De, C. K. Kim and V. M. Rotello, Adv. Drug Delivery Rev., 2008, 60, 1307-1315; S. Lal, S. E. Clare and N. J. Halas, Acc. Chem. Res., 2008, 41, 1842-1851; D. A. Giljohann, D. S. Seferos, W. L. Daniel, M. D. Massich, P. C. Patel and C. A. Mirkin, Angew. Chem., Int. Ed., 2010, 49, 3280-3294). Not only can they be used as targeted contrast agents for photothermal cancer therapy, they can serve as scaffolds for increasingly potent cancer drug delivery, as transfection agents for selective gene therapy, and as intrinsic antineoplastic agents. This tutorial review will highlight some of the many forms and recent applications of these gold nanoparticle conjugates by our lab and others, as well as their rational design and physiologic interactions.     

BOOK REVIEWS

Book reviewed in this article:
Chemical and Biochemical Applications of Lasers. Volume II. Edited by C. B radley M oore
Laser Spectroscopy III, Springer Series in Optical Spectroscopy. Volume 7. Edited by J. L. H all and J. L. C arlsten
Photosynthesis. Second Edition, D. O. H all and K. K. R ao
The Amphibian Visual System. Edited by K atherine V. F ite
The Physiology and Biochemistry of Seed Dormancy and Germination. Edited by A. A. K han
Phytochrome and Plant Growth. Studies in Biology NO. 68. Edited by R. E. K endrick and B. F rankland
The Science of Photobiology. Edited by K endric C. S mith
Photosensitized Reactions, By M. K oizumi , S. K ato , N. M ataga , T. M atsuura and Y. U sui
Bioenergetics of Membranes, Edited by L. P acker , G. C. P apageorgiou and A. T rebst     

Mass selective gas phase study of ClO, OClO, ClOO and ClAr by anion-ZEKE-photoelectron spectroscopy

Anion-ZEKE-photoelectron spectra of ClO-, OClO-, ClOO- and the van der Waals cluster ArCl- have been measured. Refined or new values for the electron affinity of ClO, OClO and ClOO have been found. The peak positions in these spectra are in very good agreement with former ClO- and OClO- anion-photoelectron spectra (K. M. Gilles, M. L. Polak and W. C. Lineberger, J. Chem. Phys., 1992, 96, 8012) and a recent ArCl- anion-ZEKE spectrum (T. Lenzer, I. Yourshaw, M. Furlanetto, G. Reiser and D. Neumark, J. Chem. Phys., 1992, 110, 9578). The higher resolution of our anion-ZEKE-photoelectron spectrum of OClO- led to a refined assignment of the corresponding anion-photoelectron spectrum. In addition, a strong difference in the relative intensities of the vibrational peaks has been found in the anion-ZEKE-spectrum of OClO- in comparison with the anion-photoelectron spectrum. For the first time, mass selective spectroscopic information has been obtained for ClOO. The strong similarity to the ArCl- spectrum indicates a weakly bound van der Waals cluster Cl.O2. Binding energies of the anion, neutral ground and neutral excited state could be deduced. These are in good agreement with the electron affinities of Cl and ClOO, but differ from theoretical values (K. A. Peterson and H. J. Werner, J. Chem. Phys., 1992, 96, 8948) by a factor of 4.5 and from thermochemically determined values (J. M. Nicovich, K. D. Kreutter, C. J. Shackelford and P. H. Wine, Chem. Phys. Lett., 1991, 179, 367 and S. Baer, H. Hippler, R. Rahn, M. Siefke, N. Seitzinger and J. Troe, J. Chem. Phys., 1991, 95, 6463) by a factor of 9.     

Hydrogen adsorption strength and sites in the metal organic framework MOF5: Comparing experiment and model calculations

Hydrogen adsorption in porous, high surface area, and stable metal organic frameworks (MOF’s) appears a novel route towards hydrogen storage materials [N.L. Rosi, J. Eckert, M. Eddaoudi, D.T. Vodak, J. Kim, M. O’Keeffe, O.M. Yaghi, Science 300 (2003) 1127; J.L.C. Rowsell, A.R. Millward, K. Sung Park, O.M. Yaghi, J. Am. Chem. Soc. 126 (2004) 5666; G. Ferey, M. Latroche, C. Serre, F. Millange, T. Loiseau, A. Percheron-Guegan, Chem. Commun. (2003) 2976; T. Loiseau, C. Serre, C. Huguenard, G. Fink, F. Taulelle, M. Henry, T. Bataille, G. Férey, Chem. Eur. J. 10 (2004) 1373]. A prerequisite for such materials is sufficient adsorption interaction strength for hydrogen adsorbed on the adsorption sites of the material because this facilitates successful operation under moderate temperature and pressure conditions. Here we report detailed information on the geometry of the hydrogen adsorption sites, based on the analysis of inelastic neutron spectroscopy (INS). The adsorption energies for the metal organic framework MOF5 equal about 800 K for part of the different sites, which is significantly higher than for nanoporous carbon materials (550 K) [H.G. Schimmel, G.J. Kearley, M.G. Nijkamp, C.T. Visser, K.P. de Jong, F.M. Mulder, Chem. Eur. J. 9 (2003) 4764], and is in agreement with what is found in first principles calculations [T. Sagara, J. Klassen, E. Ganz, J. Chem. Phys. 121 (2004) 12543; F.M. Mulder, T.J. Dingemans, M. Wagemaker, G.J. Kearley, Chem. Phys. 317 (2005) 113]. Assignments of the INS spectra is realized using comparison with independently published model calculations [F.M. Mulder, T.J. Dingemans, M. Wagemaker, G.J. Kearley, Chem. Phys. 317 (2005) 113] and structural data [T. Yildirim, M.R. Hartman, Phys. Rev. Lett. 95 (2005) 215504].     

Excluded volume effect for large and small solutes in water

The cavitation effect, i.e., the process of the creation of a void of excluded volume in bulk solvent (a cavity), is considered. The cavitation free energy is treated in terms of the information theory (IT) approach [Hummer, G.; Garde, S.; Garcia, A. E.; Paulaitis, M. E.; Pratt, L. R. J. Phys. Chem. B 1998, 102, 10469]. The binomial cell model suggested earlier is applied as the IT default distribution p(m) for the number m of solute (water) particles occupying a cavity of given size and shape. In the present work, this model is extended to cover the entire range of cavity size between small ordinary molecular solutes and bulky biomolecular structures. The resulting distribution consists of two binomial peaks responsible for producing the free energy contributions, which are proportional respectively to the volume and to the surface area of a cavity. The surface peak dominates in the large cavity limit, when the two peaks are well separated. The volume effects become decisive in the opposite limit of small cavities, when the two peaks reduce to a single-peak distribution as considered in our earlier work. With a proper interpolation procedure connecting these two regimes, the MC simulation results for model spherical solutes with radii increasing up to R = 10 A [Huang, D. H.; Geissler, P. L.; Chandler, D. J. Phys. Chem. B 2001, 105, 6704] are well reproduced. The large cavity limit conforms to macroscopic properties of bulk water solvent, such as surface tension, isothermal compressibility and Tolman length. The computations are extended to include nonspherical solutes (hydrocarbons C1-C6).     

Comments on “Cofino statistics”

Cofino et al. have put forward a method for the analysis of data from interlaboratory studies. The method uses ideas and notation from quantum chemistry, and is quite sophisticated. A critical examination shows that it has some relationship with both kernel density estimation and robust estimation methods. The quantitative results it provides are not founded on any statistical probability model, and it is not at all clear how the variance in particular is to be interpreted. There is no clear argument for adopting Cofino statistics in preference to simple problem-free methods that are at least as effective in estimating a concensus value.This report was prepared for the Statistical Subcommittee by Professor T. Fearn, with the assistance of Mr T. Robinson, and approved by the AMC on 25/11/03. The Statistical Subcommittee comprised: Dr. S.L.R. Ellison, Professor T. Fearn, Mr M. Gardner, Dr. F. Hollywood, Professor R.J. Howarth, Dr. P. Lowthian, Professor J.N. Miller, Dr. E.J. Newman, Professor B.D. Ripley, Professor M. Thompson (Chair), Dr. R. Wood, Dr. A. Williams and Mr J.J. Wilson (Secretary).     

The reaction OH + C2H4: an example of rotational channel switching

The low-temperature data for the reaction between OH and C(2)H(4) is treated canonically as either a two-well or one-well problem using the "Multiwell" suite of codes, in which a "well" refers to a minimum in the potential energy surface. The former is analogous to the two transition state model of Greenwald et al. [Greenwald, E. E.; North, S. W.; Georgievskii, Y.; Klippenstein, S. J. J. Phys. Chem. A2005, 109, 6031], while the latter reflects the dominance of the so-called "inner transition state". External rotations are treated adiabatically, causing changes in the magnitude of effective barriers as a function of temperature. Extant data are well-described with either model using only the average energy transferred in a downward direction, upon collision, ΔE(d)(T), as a fitting parameter. The best value for the parameters describing the rate coefficient as a function of temperature (200 < T/K < 400) (Data at lower temperature is too sparse to yield a recommendation.) and pressure in the form used in the NASA/JPL format [Sander, S. P.; Abbatt, J.; Barker, J. R.; Burkholder, J. B.; Friedl, R. R.; Golden, D. M.; Huie, R. E.; Kolb, C. E.; Kurylo, M. J.; Moortgat, G. K et al., Chemical Kinetics and Photochemical Data for Use in Atmospheric Studies, Evaluation Number 17, Jet Propulsion Laboratory, 2011] are k(0) = 1.0 × 10(-28)(T/300)(-3.5) cm(6) molecule(-2) s(-1) and k(∞) to 8.0 × 10(-12)(T/300)(-2.3) cm(3) molecule(-1) s(-1).     

Searching for the second oxidant in the catalytic cycle of cytochrome P450: a theoretical investigation of the iron(III)-hydroperoxo species and its epoxidation pathways

Iron(III)-hydroperoxo, [Por(CysS)Fe(III)-OOH](-), a key species in the catalytic cycle of cytochrome P450, was recently identified by EPR/ENDOR spectroscopies (Davydov, R.; Makris, T. M.; Kofman, V.; Werst, D. E.; Sligar, S. G.; Hoffman, B. M. J. Am. Chem. Soc. 2001, 123, 1403-1415). It constitutes the last station of the preparative steps of the enzyme before oxidation of an organic compound and is implicated as the second oxidant capable of olefin epoxidation (Vaz, A. D. N.; McGinnity, D. F.; Coon, M. J. Proc. Natl. Acad. Sci. U.S.A. 1998, 95, 3555-3560), in addition to the penultimate active species, Compound I (Groves, J. T.; Han, Y.-Z. In Cytochrome P450: Structure, Mechanism and Biochemistry, 2nd ed.; Ortiz de Montellano, P. R., Ed.; Plenum Press: New York, 1995; pp 3-48). In response, we present a density functional study of a model species and its ethylene epoxidation pathways. The study characterizes a variety of properties of iron(III)-hydroperoxo, such as the O-O bonding, the Fe-S bonding, Fe-O and Fe-S stretching frequencies, its electron attachment, and ionization energies. Wherever possible these properties are compared with those of Compound I. The proton affinities for protonation on the proximal and distal oxygen atoms of iron(III)-hydroperoxo, and the effect of the thiolate ligand thereof, are determined. In accordance with previous results (Harris, D. L.; Loew, G. H. J. Am. Chem. Soc. 1998, 120, 8941-8948), iron(III)-hydroperoxo is a strong base (as compared with water), and its distal protonation leads to a barrier-free formation of Compound I. The origins of this barrier-free process are discussed using a valence bond approach. It is shown that the presence of the thiolate is essential for this process, in line with the "push effect" deduced by experimentalists (Sono, M.; Roach, M. P.; Coulter, E. D.; Dawson, J. H. Chem. Rev. 1996, 96, 2841-2887). Finally, four epoxidation pathways of iron(III)-hydroxperoxo are located, in which the species transfers oxygen to ethylene either from the proximal or from the distal sites, in both concerted and stepwise manners. The barriers for the four mechanisms are 37-53 kcal/mol, in comparison with 14 kcal/mol for epoxidation by Compound I. It is therefore concluded that iron(III)-hydroperoxo, as such, cannot be a second oxidant, in line with its significant basicity and poor electron-accepting capability. Possible versions of a second oxidant are discussed.     

Effects of a chemical chaperone on genetic mutations in ��-galactosidase A in Korean patients with Fabry disease

Fabry disease is an X-linked inborn error of glycosphingolipid catabolism that results from mutations in the gene encoding the α-galactosidase A (GLA) enzyme. We have identified 15 distinct mutations in the GLA gene in 13 unrelated patients with classic Fabry disease and 2 unrelated patients with atypical Fabry disease. Two of the identified mutations were novel (i.e., the D231G missense mutation and the L268delfsX1 deletion mutation). This study evaluated the effects of the chemical chaperones 1-deoxygalactonojirimycin (DGJ) on the function of GLA in vitro, in cells containing missense mutations in the GLA gene. Nine missense and a nonsense mutations, including one novel mutation were cloned into mammalian expression vectors. After transient expression in COS-7 cells, GLA enzyme activity and protein expression were analyzed using fluorescence spectrophotometry and Western blot analysis, respectively. DGJ enhanced GLA enzyme activity in the M42V, I91T, R112C and F113L mutants. Interestingly, the I91T and F113L mutations are associated with the atypical form of Fabry disease. However, DGJ treatment did not have any significant effect on the GLA enzyme activity and protein expression of other mutants, including C142W, D231G, D266N, and S297F. Of note, GLA enzyme activity was not detected in the novel mutant (i.e., D231G), although protein expression was similar to the wild type. In the absence of DGJ, the E66Q mutant had wild-type levels of GLA protein expression and approximately 40% GLA activity, indicating that E66Q is either a mild mutation or a functional single nucleotide polymorphism (SNP). Thus, the results of this study suggest that the chemical chaperone DGJ enhances GLA enzyme activity and protein expression in milder mutations associated with the atypical form of Fabry disease.     

Marked difference in the electronic structure of cyanide-ligated ferric protoglobins and myoglobin due to heme ruffling

Electron paramagnetic resonance experiments reveal a significant difference between the principal g values (and hence ligand-field parameters) of the ferric cyanide-ligated form of different variants of the protoglobin of Methanosarcina acetivorans (MaPgb) and of horse heart myoglobin (hhMb). The largest principal g value of the ferric cyanide-ligated MaPgb variants is found to be significantly lower than for any of the other globins reported so far. This is at least partially caused by the strong heme distortions as proven by the determination of the hyperfine interaction of the heme nitrogens and mesoprotons. Furthermore, the experiments confirm recent theoretical predictions [Forti, F.; Boechi, L., Bikiel, D., Martí, M.A.; Nardini, M.; Bolognesi, M.; Viappiani, C.; Estrin, D.; Luque, F. J. J. Phys. Chem. B2011, 115, 13771-13780] that Phe(G8)145 plays a crucial role in the ligand modulation in MaPgb. Finally, the influence of the N-terminal 20 amino-acid chain on the heme pocket in these protoglobins is also proven.     

Design principles for Brownian molecular machines: how to swim in molasses and walk in a hurricane

Protein molecular motors-perfected over the course of millions of years of evolution-play an essential role in moving and assembling biological structures. Recently chemists have been able to synthesize molecules that emulate in part the remarkable capabilities of these biomolecular motors (for extensive reviews see the recent papers: E. R. Kay, D. A. Leigh and F. Zerbetto, Angew. Chem., Int. Ed., 2006, 46, 72-191; W. R. Browne and B. L. Feringa, Nat. Nanotechnol., 2006, 1, 25-35; M. N. Chatterjee, E. R. Kay and D. A. Leigh, J. Am. Chem. Soc., 2006, 128, 4058-4073; G. S. Kottas, L. I. Clarke, D. Horinek and J. Michl, Chem. Rev., 2005, 105, 1281-1376; M. A. Garcia-Garibay, Proc. Natl. Acad. Sci., U. S. A., 2005, 102, 10771-10776)). Like their biological counterparts, many of these synthetic machines function in an environment where viscous forces dominate inertia-to move they must "swim in molasses". Further, the thermal noise power exchanged reversibly between the motor and its environment is many orders of magnitude greater than the power provided by the chemical fuel to drive directed motion. One might think that moving in a specific direction would be as difficult as walking in a hurricane. Yet biomolecular motors (and increasingly, synthetic motors) move and accomplish their function with almost deterministic precision. In this Perspective we will investigate the physical principles that govern nanoscale systems at the single molecule level and how these principles can be useful in designing synthetic molecular machines.     

Melting temperature of ice Ih calculated from coexisting solid-liquid phases

We carried out molecular-dynamics simulations by using the two-phase coexistence method with the constant pressure, particle number, and enthalpy ensemble to compute the melting temperature of proton-disordered hexagonal ice I(h) at 1-bar pressure. Four models of water were considered, including the widely used TIP4P [W. L. Jorgensen, J. Chandrasekha, J. D. Madura, R. W. Impey, and M. L. Klein, J. Chem. Phys.79, 926 (1983)] and TIP5P [M. W. Mahoney and W. L. Jorgensen J. Chem. Phys.112, 8910 (2000)] models, as well as recently improved TIP4P and TIP5P models for use with Ewald techniques-the TIP4P-Ew [W. Horn, W. C. Swope, J. W. Pitera, J. C. Madura, T. J. Dick, G. L. Hura, and T. Head-Gordon, J. Chem. Phys.120, 9665 (2004)] and TIP5P-Ew [S. W. Rick, J. Chem. Phys.120, 6085 (2004)] models. The calculated melting temperature at 1 bar is T(m) = 229 +/- 1 K for the TIP4P and T(m) = 272.0 +/- 0.6 K for the TIP5P ice I(h), both are consistent with previous simulations based on free-energy methods. For the TIP4P-Ew and TIP5P-Ew models, the calculated melting temperature is T(m) = 257.0 +/- 1.1 K and T(m) = 253.9 +/- 1.1 K, respectively.     

Geminate carbon monoxide rebinding to a c-type haem

A chemically modified form of cytochrome c(cyt. c), termed carboxymethyl cytochrome c(cm cyt. c), possesses a vacant sixth coordination site to the haem iron that is available to bind external ligands. We present data on the rapid flash photolysis of CO from the ferrous haem iron of cm cyt. c and describe the kinetics and spectral transitions that accompany the recombination. This was achieved using 30-femtosecond laser pulses and a white light continuum to monitor spectral transitions. Whereas the photo-dissociation quantum yield is close to 1, the yield of CO escape from the protein (the apparent quantum yield, varphi) relative to myoglobin (varphi=1) is small due to rapid geminate recombination of CO. On ligand photo-dissociation the haem undergoes a spin-state transition from low-spin ferrous CO bound to penta-coordinate high-spin. Subsequently the system reverts to the CO bound form. The data were fitted with a minimum number of exponentials using global analysis. Recombination of CO with the haem iron of cm cyt. c is multiphasic (tau=16 ps, 120 ps and 1 ns), involving three spectrally distinct components. The fraction of haem (0.11) not recombining with CO within 4 ns is similar to the value of varphi(0.12) measured on the same preparation by the "pulse method" (M. Brunori, G. Giacometti, E. Antonini and J. Wyman, Proc. Natl. Acad. Sci. USA, 1973, 70, 3141-3144, ). This implies that no further geminate recombination occurs at t>4 ns. This unusually efficient CO-haem geminate recombination indicates the sterically hindered ("caged") nature of the distal haem pocket in cm cyt. c from which it is difficult for CO to escape. The large geminate phase may be contrasted with the behaviour of myoglobin in which geminate recombination is small. This is in general agreement with the well-documented extensive structural dynamics in myoglobin that allow ligand passage, and a higher structural rigidity in cyt. c imposed by the restraints of minimising reorganisation energy for electron transfer (M. Brunori, D. Bourgeois and D. Vallone, J. Struct. Biol., 2004, 147, 223-234, ). The high pH ferrous form of cm cyt. c is a low-spin species having a lysine bound to the central iron atom of the haem (M. Brunori, M. Wilson and E. Antonini, J. Biol. Chem., 1972, 247, 6076-6081; G. Silkstone, G. Stanway, P. Brzezinski and M. Wilson, Biophys. Chem., 2002, 98, 65-77, ). This high pH (pH approximately 8) form of deoxy cm cyt. c undergoes photo-dissociation of lysine (although the proximal histidine is possible) after photo-excitation. Recombination occurs with a time constant (tau) of approximately 7 ps. This is similar to that observed for the geminate rebinding of the Met80 residue in native ferrous cyt. c(tau approximately 6 ps) following its photo-dissociation (S. Cianetti, M. Negrerie, M. Vos, J.-L. Martin and S. Kruglik, J. Am. Chem. Soc., 2004, 126, 13 932-13 933; W. Wang, X. Ye, A. Demidov, F. Rosca, T. Sjodin, W. Cao, M. Sheeran and P. Champion, J. Phys. Chem., 2000, 104, 10 789-10 801, ).     

Network regularised Cox regression and multiplex network models to predict disease comorbidities and survival of cancer

In cancer genomics, gene expression levels provide important molecular signatures for all types of cancer, and this could be very useful for predicting the survival of cancer patients. However, the main challenge of gene expression data analysis is high dimensionality, and microarray is characterised by few number of samples with large number of genes. To overcome this problem, a variety of penalised Cox proportional hazard models have been proposed. We introduce a novel network regularised Cox proportional hazard model and a novel multiplex network model to measure the disease comorbidities and to predict survival of the cancer patient. Our methods are applied to analyse seven microarray cancer gene expression datasets: breast cancer, ovarian cancer, lung cancer, liver cancer, renal cancer and osteosarcoma. Firstly, we applied a principal component analysis to reduce the dimensionality of original gene expression data. Secondly, we applied a network regularised Cox regression model on the reduced gene expression datasets. By using normalised mutual information method and multiplex network model, we predict the comorbidities for the liver cancer based on the integration of diverse set of omics and clinical data, and we find the diseasome associations (disease–gene association) among different cancers based on the identified common significant genes. Finally, we evaluated the precision of the approach with respect to the accuracy of survival prediction using ROC curves. We report that colon cancer, liver cancer and renal cancer share the CXCL5 gene, and breast cancer, ovarian cancer and renal cancer share the CCND2 gene. Our methods are useful to predict survival of the patient and disease comorbidities more accurately and helpful for improvement of the care of patients with comorbidity. Software in Matlab and R is available on our GitHub page: https://github.com/ssnhcom/NetworkRegularisedCox.git.     

Global genetic analysis of all single nucleotide polymorphisms in exons of the human deoxyribonuclease I-like 3 gene and their effect on its catalytic activity

Deoxyribonucleases (DNases) have been suggested to be implicated in the pathophysiology of autoimmune diseases. In the DNASE1L3 gene encoding human DNase I‐like 3 (DNase 1L3), a member of the DNase I family, only two non‐synonymous (R178 H and R206C) single nucleotide polymorphisms (SNPs) have been examined [Ueki et al., Clin. Chim. Acta 2009, 407, 20–24]. Three other non‐synonymous (G82R, K96N, and I243M) and four synonymous (S17S, T84T, R92R, and A181A) SNPs, in addition to R206C and R178H, have been identified in DNASE1L3. We investigated the distribution of all these SNPs in exons of the gene in eight Asian, three African, and three Caucasian populations worldwide using newly devised genotyping methods. SNP T84T showed polymorphism in all the populations, and R92R was polymorphic in the three African and three Caucasian populations; R206C was distributed only in Caucasian populations. In contrast, no minor allele was found in five SNPs (S17S, G82R, K96N, A181A, and I243M) in DNASE1L3. Generally, the DNase 1L3 gene shows relatively low genetic diversity with regard to exonic SNPs. When the effect of amino acid/nucleotide substitutions resulting from the SNPs on DNase 1L3 activity was examined, none of the synonymous SNPs had any effect on the DNase 1L3 activity, whereas among non‐synonymous SNPs, SNP G82R diminished the activity of the enzyme, being similar to R206C. These findings permit us to assume that, although only R206 exhibits polymorphisms in a Caucasian‐specific manner, at least SNPs G82R and R206C in DNASE1L3 might be potential risk factors for autoimmune disease.     

Validation studies of the site-directed docking program LibDock

The performance of the site-features docking algorithm LibDock has been evaluated across eight GlaxoSmithKline targets as a follow-up to a broad validation study of docking and scoring software (Warren, G. L.; Andrews, W. C.; Capelli, A.; Clarke, B.; Lalonde, J.; Lambert, M. H.; Lindvall, M.; Nevins, N.; Semus, S. F.; Senger, S.; Tedesco, G.; Walls, I. D.; Woolven, J. M.; Peishoff, C. E.; Head, M. S. J. Med. Chem. 2006, 49, 5912-5931). Docking experiments were performed to assess both the accuracy in reproducing the binding mode of the ligand and the retrieval of active compounds in a virtual screening protocol using both the DJD (Diller, D. J.; Merz, K. M., Jr. Proteins 2001, 43, 113-124) and LigScore2 (Krammer, A. K.; Kirchoff, P. D.; Jiang, X.; Venkatachalam, C. M.; Waldman, M. J. Mol. Graphics Modell. 2005, 23, 395-407) scoring functions. This study was conducted using DJD scoring, and poses were rescored using all available scoring functions in the Accelrys LigandFit module, including LigScore2. For six out of eight targets at least 30% of the ligands were docked within a root-mean-square difference (RMSD) of 2.0 A for the crystallographic poses when the LigScore2 scoring function was used. LibDock retrieved at least 20% of active compounds in the top 10% of screened ligands for four of the eight targets in the virtual screening protocol. In both studies the LigScore2 scoring function enhanced the retrieval of crystallographic poses or active compounds in comparison with the results obtained using the DJD scoring function. The results for LibDock accuracy and ligand retrieval in virtual screening are compared to 10 other docking and scoring programs. These studies demonstrate the utility of the LigScore2 scoring function and that LibDock as a feature directed docking method performs as well as docking programs that use genetic/growing and Monte Carlo driven algorithms.     

Theoretical analysis and computer simulation of fluorescence lifetime measurements. I. Kinetic regimes and experimental time scales

The configuration-controlled regime and the diffusion-controlled regime of conformation-modulated fluorescence emission are systematically studied for Markovian and non-Markovian dynamics of the reaction coordinate. A path integral simulation is used to model fluorescence quenching processes on a semiflexible chain. First-order inhomogeneous cumulant expansion in the configuration-controlled regime defines a lower bound for the survival probability, while the Wilemski-Fixman approximation in the diffusion-controlled regime defines an upper bound. Inclusion of the experimental time window of the fluorescence measurement adds another dimension to the two kinetic regimes and provides a unified perspective for theoretical analysis and experimental investigation. We derive a rigorous generalization of the Wilemski-Fixman approximation [G. Wilemski and M. Fixman, J. Chem. Phys. 60, 866 (1974)] and recover the 1/D expansion of the average lifetime derived by Weiss [G. H. Weiss, J. Chem. Phys. 80, 2880 (1984)].