Comparison between computational alanine scanning and per-residue binding free energy decomposition for protein-protein association using MM-GBSA: application to the TCR-p-MHC complex.

Recognition by the T-cell receptor (TCR) of immunogenic peptides (p) presented by Class I major histocompatibility complexes (MHC) is the key event in the immune response against virus-infected cells or tumor cells. A study of the 2C TCR/SIYR/H-2K(b) system using a computational alanine scanning and a much faster binding free energy decomposition based on the Molecular Mechanics-Generalized Born Surface Area (MM-GBSA) method is presented. The results show that the TCR-p-MHC binding free energy decomposition using this approach and including entropic terms provides a detailed and reliable description of the interactions between the molecules at an atomistic level. Comparison of the decomposition results with experimentally determined activity differences for alanine mutants yields a correlation of 0.67 when the entropy is neglected and 0.72 when the entropy is taken into account. Similarly, comparison of experimental activities with variations in binding free energies determined by computational alanine scanning yields correlations of 0.72 and 0.74 when the entropy is neglected or taken into account, respectively. Some key interactions for the TCR-p-MHC binding are analyzed and some possible side chains replacements are proposed in the context of TCR protein engineering. In addition, a comparison of the two theoretical approaches for estimating the role of each side chain in the complexation is given, and a new ad hoc approach to decompose the vibrational entropy term into atomic contributions, the linear decomposition of the vibrational entropy (LDVE), is introduced. The latter allows the rapid calculation of the entropic contribution of interesting side chains to the binding. This new method is based on the idea that the most important contributions to the vibrational entropy of a molecule originate from residues that contribute most to the vibrational amplitude of the normal modes. The LDVE approach is shown to provide results very similar to those of the exact but highly computationally demanding method.

Plutonium from above-ground nuclear tests in milk teeth: investigation of placental transfer in children born between 1951 and 1995 in Switzerland.

BACKGROUND: Occupational risks, the present nuclear threat, and the potential danger associated with nuclear power have raised concerns regarding the metabolism of plutonium in pregnant women. OBJECTIVE: We measured plutonium levels in the milk teeth of children born between 1951 and 1995 to assess the potential risk that plutonium incorporated by pregnant women might pose to the radiosensitive tissues of the fetus through placenta transfer. METHODS: We used milk teeth, whose enamel is formed during pregnancy, to investigate the transfer of plutonium from the mother's blood plasma to the fetus. We measured plutonium using sensitive sector field inductively coupled plasma mass spectrometry techniques. We compared our results with those of a previous study on strontium-90 ((90)Sr) released into the atmosphere after nuclear bomb tests. RESULTS: Results show that plutonium activity peaks in the milk teeth of children born about 10 years before the highest recorded levels of plutonium fallout. By contrast, (90)Sr, which is known to cross the placenta barrier, manifests differently in milk teeth, in accordance with (90)Sr fallout deposition as a function of time. CONCLUSIONS: These findings demonstrate that plutonium found in milk teeth is caused by fallout that was inhaled around the time the milk teeth were shed and not from any accumulation during pregnancy through placenta transfer. Thus, plutonium may not represent a radiologic risk for the radiosensitive tissues of the fetus.

Host and invader impact of transfer of the clc genomic island into Pseudomonas aeruginosa PAO1

Genomic islands, large potentially mobile regions of bacterial chromosomes, are a major contributor to bacteria evolution. Here, we investigated the fitness cost and phenotypic differences between the bacterium Pseudomonas aeruginosa PAO1 and a derivative carrying one integrated copy of the clc element, a 103-kb genomic island [and integrative and conjugative element (ICE)] originating in Pseudomonas sp. strain B13 and a close relative of genomic islands found in clinical and environmental isolates of P. aeruginosa. By using a combination of whole genome transcriptome profiling, phenotypic arrays, competition experiments, and biofilm formation studies, only few differences became apparent, such as reduced biofilm growth and fourfold stationary phase repression of genes involved in acetoin metabolism in PAO1 containing the clc element. In contrast, PAO1 carrying the clc element acquired the capacity to grow on 3-chlorobenzoate and 2-aminophenol as sole carbon and energy substrates. No fitness loss >1% was detectable in competition experiments between PAO1 and PAO1 carrying the clc element. The genes from the clc element were not silent in PAO1, and excision was observed, although transfer of clc from PAO1 to other recipient bacteria was reduced by two orders of magnitude. Our results indicate that newly acquired mobile DNA not necessarily invoke an important fitness cost on their host. Absence of immediate detriment to the host may have contributed to the wide distribution of genomic islands like clc in bacterial genomes

Postischemic recovery of heart metabolism and function: role of mitochondrial fatty acid transfer.

Postischemic recovery of contractile function is better in hearts from fasted rats than in hearts from fed rats. In this study, we examined whether feeding-induced inhibition of palmitate oxidation at the level of carnitine palmitoyl transferase I is involved in the mechanism underlying impaired recovery of contractile function. Hearts isolated from fasted or fed rats were submitted to no-flow ischemia followed by reperfusion with buffer containing 8 mM glucose and either 0.4 mM palmitate or 0.8 mM octanoate. During reperfusion, oxidation of palmitate was higher after fasting than after feeding, whereas oxidation of octanoate was not influenced by the nutritional state. In the presence of palmitate, recovery of left ventricular developed pressure was better in hearts from fasted rats. Substitution of octanoate for palmitate during reperfusion enhanced recovery of left ventricular developed pressure in hearts from fed rats. However, the chain length of the fatty acid did not influence diastolic contracture. The results suggest that nutritional variation of mitochondrial fatty acid transfer may influence postischemic recovery of contractile function.

Theoretical study of the second-order vibrational Stark effect

The behaviour of the harmonic infrared frequency of diatomic molecules subjected to moderate static uniform electric fields is analysed. The potential energy expression has been developed as a function of a static uniform electric field, which brings about a formulation describing the frequency versus field strength curve. With the help of the first and second derivatives of the expressions obtained, which correspond to the first- and second-order Stark effects, it was possible to find the maxima of the frequency versus field strength curves for a series of molecules using a Newton-Raphson search. A method is proposed which requires only the calculation of a few energy derivatives at a particular value of the field strength. At the same time, the expression for the dependence of the interatomic distance on the electric field strength is derived and the minimum of this curve is found for the same species. Derived expressions and numerical results are discussed and compared with other studi

Activity-dependent regulation of energy metabolism by astrocytes: an update.

Astrocytes play a critical role in the regulation of brain metabolic responses to activity. One detailed mechanism proposed to describe the role of astrocytes in some of these responses has come to be known as the astrocyte-neuron lactate shuttle hypothesis (ANLSH). Although controversial, the original concept of a coupling mechanism between neuronal activity and glucose utilization that involves an activation of aerobic glycolysis in astrocytes and lactate consumption by neurons provides a heuristically valid framework for experimental studies. In this context, it is necessary to provide a survey of recent developments and data pertaining to this model. Thus, here, we review very recent experimental evidence as well as theoretical arguments strongly supporting the original model and in some cases extending it. Aspects revisited include the existence of glutamate-induced glycolysis in astrocytes in vitro, ex vivo, and in vivo, lactate as a preferential oxidative substrate for neurons, and the notion of net lactate transfer between astrocytes and neurons in vivo. Inclusion of a role for glycogen in the ANLSH is discussed in the light of a possible extension of the astrocyte-neuron lactate shuttle (ANLS) concept rather than as a competing hypothesis. New perspectives offered by the application of this concept include a better understanding of the basis of signals used in functional brain imaging, a role for neuron-glia metabolic interactions in glucose sensing and diabetes, as well as novel strategies to develop therapies against neurodegenerative diseases based upon improving astrocyte-neuron coupled energetics.

Basis set superposition error-counterpoise corrected potential energy surfaces : application to hydrogen peroxide ... X (X=F-, Cl-, Br-, Li+, Na+) complexes

Møller-Plesset (MP2) and Becke-3-Lee-Yang-Parr (B3LYP) calculations have been used to compare the geometrical parameters, hydrogen-bonding properties, vibrational frequencies and relative energies for several X- and X+ hydrogen peroxide complexes. The geometries and interaction energies were corrected for the basis set superposition error (BSSE) in all the complexes (1-5), using the full counterpoise method, yielding small BSSE values for the 6-311 + G(3df,2p) basis set used. The interaction energies calculated ranged from medium to strong hydrogen-bonding systems (1-3) and strong electrostatic interactions (4 and 5). The molecular interactions have been characterized using the atoms in molecules theory (AIM), and by the analysis of the vibrational frequencies. The minima on the BSSE-counterpoise corrected potential-energy surface (PES) have been determined as described by S. Simón, M. Duran, and J. J. Dannenberg, and the results were compared with the uncorrected PES

Simple finite field method for calculation of static and dynamic vibrational hyperpolarizabilities: curvature contributions

In the static field limit, the vibrational hyperpolarizability consists of two contributions due to: (1) the shift in the equilibrium geometry (known as nuclear relaxation), and (2) the change in the shape of the potential energy surface (known as curvature). Simple finite field methods have previously been developed for evaluating these static field contributions and also for determining the effect of nuclear relaxation on dynamic vibrational hyperpolarizabilities in the infinite frequency approximation. In this paper the finite field approach is extended to include, within the infinite frequency approximation, the effect of curvature on the major dynamic nonlinear optical processes

Field-induced coordinates for the determination of dynamic vibrational nonlinear optical properties

Three conjugated organic molecules that span a range of polarity and valence-bond/charge transfer characteristics were studied. It was found that dispersion can be insignificant, and that adequate treatment can be achieved with frequency-dependent field-induced vibrational coordinates (FD-FICs)

Variational calculation of static and dynamic vibrational nonlinear optical properties

The vibrational configuration interaction method used to obtain static vibrational (hyper)polarizabilities is extended to dynamic nonlinear optical properties in the infinite optical frequency approximation. Illustrative calculations are carried out on H2 O and N H3. The former molecule is weakly anharmonic while the latter contains a strongly anharmonic umbrella mode. The effect on vibrational (hyper)polarizabilities due to various truncations of the potential energy and property surfaces involved in the calculation are examined

Variational calculation of vibrational linear and nonlinear optical properties

A variational approach for reliably calculating vibrational linear and nonlinear optical properties of molecules with large electrical and/or mechanical anharmonicity is introduced. This approach utilizes a self-consistent solution of the vibrational Schrödinger equation for the complete field-dependent potential-energy surface and, then, adds higher-level vibrational correlation corrections as desired. An initial application is made to static properties for three molecules of widely varying anharmonicity using the lowest-level vibrational correlation treatment (i.e., vibrational Møller-Plesset perturbation theory). Our results indicate when the conventional Bishop-Kirtman perturbation method can be expected to break down and when high-level vibrational correlation methods are likely to be required. Future improvements and extensions are discussed

Charge transfer in DNA: hole charge is confined to a single base pair due to solvation effects

We include solvation effects in tight-binding Hamiltonians for hole states in DNA. The corresponding linear-response parameters are derived from accurate estimates of solvation energy calculated for several hole charge distributions in DNA stacks. Two models are considered: (A) the correction to a diagonal Hamiltonian matrix element depends only on the charge localized on the corresponding site and (B) in addition to this term, the reaction field due to adjacent base pairs is accounted for. We show that both schemes give very similar results. The effects of the polar medium on the hole distribution in DNA are studied. We conclude that the effects of polar surroundings essentially suppress charge delocalization in DNA, and hole states in (GC)n sequences are localized on individual guanines

Systematic study of the static electrical properties of the CO molecule: influence of the basis set size and correlation energy

The influence of the basis set size and the correlation energy in the static electrical properties of the CO molecule is assessed. In particular, we have studied both the nuclear relaxation and the vibrational contributions to the static molecular electrical properties, the vibrational Stark effect (VSE) and the vibrational intensity effect (VIE). From a mathematical point of view, when a static and uniform electric field is applied to a molecule, the energy of this system can be expressed in terms of a double power series with respect to the bond length and to the field strength. From the power series expansion of the potential energy, field-dependent expressions for the equilibrium geometry, for the potential energy and for the force constant are obtained. The nuclear relaxation and vibrational contributions to the molecular electrical properties are analyzed in terms of the derivatives of the electronic molecular properties. In general, the results presented show that accurate inclusion of the correlation energy and large basis sets are needed to calculate the molecular electrical properties and their derivatives with respect to either nuclear displacements or/and field strength. With respect to experimental data, the calculated power series coefficients are overestimated by the SCF, CISD, and QCISD methods. On the contrary, perturbation methods (MP2 and MP4) tend to underestimate them. In average and using the 6-311 + G(3df) basis set and for the CO molecule, the nuclear relaxation and the vibrational contributions to the molecular electrical properties amount to 11.7%, 3.3%, and 69.7% of the purely electronic μ, α, and β values, respectively

Laser-induced forward transfer of liquids: Study of the droplet ejection process

Laser-induced forward transfer (LIFT) is a laser direct-write technique that offers the possibility of printing patterns with a high spatial resolution from a wide range of materials in a solid or liquid state, such as conductors, dielectrics, and biomolecules in solution. This versatility has made LIFT a very promising alternative to lithography-based processes for the rapid prototyping of biomolecule microarrays. Here, we study the transfer process through the LIFT of droplets of a solution suitable for microarray preparation. The laser pulse energy and beam size were systematically varied, and the effect on the transferred droplets was evaluated. Controlled transfers in which the deposited droplets displayed optimal features could be obtained by varying these parameters. In addition, the transferred droplet volume displayed a linear dependence on the laser pulse energy. This dependence allowed determining a threshold energy density value, independent of the laser focusing conditions, which acted as necessary conditions for the transfer to occur. The corresponding sufficient condition was given by a different total energy threshold for each laser beam dimension. The threshold energy density was found to be the dimensional parameter that determined the amount of the transferred liquid per laser pulse, and there was no substantial loss of material due to liquid vaporization during the transfer.

Time-resolved imaging of the laser forward transfer of liquids

Time-resolved imaging is carried out to study the dynamics of the laser-induced forward transfer of an aqueous solution at different laser fluences. The transfer mechanisms are elucidated, and directly correlated with the material deposited at the analyzed irradiation conditions. It is found that there exists a fluence range in which regular and well-defined droplets are deposited. In this case, laser pulse energy absorption results in the formation of a plasma, which expansion originates a cavitation bubble in the liquid. After the further expansion and collapse of the bubble, a long and uniform jet is developed, which advances at a constant velocity until it reaches the receptor substrate. On the other hand, for lower fluences no material is deposited. In this case, although a jet can be also generated, it recoils before reaching the substrate. For higher fluences, splashing is observed on the receptor substrate due to the bursting of the cavitation bubble. Finally, a discussion of the possible mechanisms which lead to such singular dynamics is also provided.