MM-GBSA binding free energy decomposition and T cell receptor engineering.

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. The major determinant of T cell activation is the affinity of the TCR for the peptide-MHC complex, though kinetic parameters are also important. A study of the 2C TCR/SIYR/H-2Kb system using a binding free energy decomposition (BFED) based on the MM-GBSA approach had been performed to assess the performance of the approach on this system. The results showed that the TCR-p-MHC BFED including entropic terms provides a detailed and reliable description of the energetics of the interaction (Zoete and Michielin, 2007). Based on these results, we have developed a new approach to design sequence modifications for a TCR recognizing the human leukocyte antigen (HLA)-A2 restricted tumor epitope NY-ESO-1. NY-ESO-1 is a cancer testis antigen expressed not only in melanoma, but also on several other types of cancers. It has been observed at high frequencies in melanoma patients with unusually positive clinical outcome and, therefore, represents an interesting target for adoptive transfer with modified TCR. Sequence modifications of TCR potentially increasing the affinity for this epitope have been proposed and tested in vitro. T cells expressing some of the proposed TCR mutants showed better T cell functionality, with improved killing of peptide-loaded T2 cells and better proliferative capacity compared to the wild type TCR expressing cells. These results open the door of rational TCR design for adoptive transfer cancer therapy.

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.

Arrested kinetic Li isotope fractionation at the margin of the llimaussaq complex, South Greenland: Evidence for open-system processes during final cooling of peralkaline igneous rocks

Li contents [Li] and isotopic composition (delta Li-7) of mafic minerals (mainly amphibole and clinopyroxene) from the alkaline to peralkaline Ilimaussaq plutonic complex, South Greenland, track the behavior of Li and its isotopes during magmatic differentiation and final cooling of an alkaline igneous system. [Li] in amphibole increase from < 10 ppm in Caamphiboles of the least differentiated unit to >3000 ppm in Na-amphiboles of the highly evolved units. In contrast, [Li] in clinopyroxene are comparatively low (<85 ppm) and do not vary systematically with differentiation. The distribution of Li between amphibole and pyroxene is controlled by the major element composition of the minerals (Ca-rich and Na-rich, respectively) and changes in oxygen fugacity (due to Li incorporation via coupled substitution with ferric iron) during magmatic differentiation. delta(7) Li values of all minerals span a wide range from + 17 to - 8 parts per thousand, with the different intrusive units of the complex having distinct Li isotopic systematics. Amphiboles, which dominate the Li budget of whole-rocks from the inner part of the complex, have constant delta Li-7 of + 1.8 +/- 2.2 parts per thousand (2 sigma, n = 15). This value reflects a homogeneous melt reservoir and is consistent with their mantle derivation, in agreement with published O and Nd isotopic data. Clinopyroxenes of these samples are consistently lighter, with Delta Li-7(amph-cpx). as large as 8 parts per thousand and are thus not in Li isotope equilibrium. These low values probably reflect late-stage diffusion of Li into clinopyroxene during final cooling of the rocks, thus enriching the clinopyroxene in 6 Li. At the margin of the complex delta(7) Li in the syenites increases systematically, from +2 to high values of + 14 parts per thousand. This, coupled with the observed Li isotope systematics of the granitic country rocks, reflects post-magmatic open-system processes occurring during final cooling of the intrusion. Although the shape and magnitude of the Li isotope and elemental profiles through syenite and country rock are suggestive of diffusion-driven isotope fractionation, they cannot be modeled by one-dimensional diffusive transport and point to circulation of a fluid having a high 67 Li value (possibly seawater) along the chilled contact. In all, this study demonstrates that Li isotopes can be used to identify complex fluid- and diffusion-governed processes taking place during the final cooling of such rocks. (c) 2007 Elsevier B.V All rights reserved.

The influence of cationic lipid type on in-vitro release kinetic profiles of antisense oligonucleotide from cationic nanoemulsions.

Novel formulations of cationic nanoemulsions based on three different lipids were developed to strengthen the attraction of the polyanionic oligonucleotide (ODN) macromolecules to the cationic moieties on the oil nanodroplets. These formulations were developed to prolong the release of the ODN from the nanoemulsion under appropriate physiological dilutions as encountered in the eye following topical application. Increasing the concentration of the new cationic lipid exhibiting two cationic amine groups (AOA) in the emulsion from 0.05% to 0.4% did not alter markedly the particle size or zeta potential value of the blank cationic nanoemulsion. The extent of ODN association did not vary significantly when the initial concentration of ODN remained constant at 10 microM irrespective of the cationic lipid nature. However, the zeta potential value dropped consistently with the low concentrations of 0.05% and 0.1% of AOA in the emulsions suggesting that an electrostatic attraction occurred between the cationic lipids and the polyanionic ODN molecules at the o/w interface. Only the nanoemulsion prepared with N-[1-(2,3-dioleoyloxy)propyl]-N,N,N-trimethylammonium salts (DOTAP) remained physically stable over time. DOTAP cationic lipid nanoemulsion was the most efficient formulation capable of retaining the ODN despite the high dilution of 1:100 with simulated tear solution (STS). Less than 10% of the ODN was exchanged in contrast to 40-50% with the other cationic nanoemulsions. The in-vitro release kinetic behavior of ODN exchange with physiological anions present in the STS appears to be complex and difficult to characterize using mathematical fitting model equations. Further pharmacokinetic studies are needed to verify our kinetic assumptions and confirm the in-vitro ODN release profile from DOTAP cationic nanoemulsions.

Spatially-resolved eigenmode decomposition of red blood cells membrane fluctuations questions the role of ATP in flickering.

Red blood cells (RBCs) present unique reversible shape deformability, essential for both function and survival, resulting notably in cell membrane fluctuations (CMF). These CMF have been subject of many studies in order to obtain a better understanding of these remarkable biomechanical membrane properties altered in some pathological states including blood diseases. In particular the discussion over the thermal or metabolic origin of the CMF has led in the past to a large number of investigations and modeling. However, the origin of the CMF is still debated. In this article, we present an analysis of the CMF of RBCs by combining digital holographic microscopy (DHM) with an orthogonal subspace decomposition of the imaging data. These subspace components can be reliably identified and quantified as the eigenmode basis of CMF that minimizes the deformation energy of the RBC structure. By fitting the observed fluctuation modes with a theoretical dynamic model, we find that the CMF are mainly governed by the bending elasticity of the membrane and that shear and tension elasticities have only a marginal influence on the membrane fluctations of the discocyte RBC. Further, our experiments show that the role of ATP as a driving force of CMF is questionable. ATP, however, seems to be required to maintain the unique biomechanical properties of the RBC membrane that lead to thermally excited CMF.

An attempt to refine rockfall hazard zoning based on the kinetic energy, frequency and fragmentation degree

Rockfall hazard zoning is usually achieved using a qualitative estimate of hazard, and not an absolute scale. In Switzerland, danger maps, which correspond to a hazard zoning depending on the intensity of the considered phenomenon (e.g. kinetic energy for rockfalls), are replacing hazard maps. Basically, the danger grows with the mean frequency and with the intensity of the rockfall. This principle based on intensity thresholds may also be applied to other intensity threshold values than those used in Switzerland for rockfall hazard zoning method, i.e. danger mapping. In this paper, we explore the effect of slope geometry and rockfall frequency on the rockfall hazard zoning. First, the transition from 2D zoning to 3D zoning based on rockfall trajectory simulation is examined; then, its dependency on slope geometry is emphasized. The spatial extent of hazard zones is examined, showing that limits may vary widely depending on the rockfall frequency. This approach is especially dedicated to highly populated regions, because the hazard zoning has to be very fine in order to delineate the greatest possible territory containing acceptable risks.

A pseudo-spectral method for the simulation of poro-elastic seismic wave propagation in 2D polar coordinates using domain decomposition

We present a novel numerical approach for the comprehensive, flexible, and accurate simulation of poro-elastic wave propagation in 2D polar coordinates. An important application of this method and its extensions will be the modeling of complex seismic wave phenomena in fluid-filled boreholes, which represents a major, and as of yet largely unresolved, computational problem in exploration geophysics. In view of this, we consider a numerical mesh, which can be arbitrarily heterogeneous, consisting of two or more concentric rings representing the fluid in the center and the surrounding porous medium. The spatial discretization is based on a Chebyshev expansion in the radial direction and a Fourier expansion in the azimuthal direction and a Runge-Kutta integration scheme for the time evolution. A domain decomposition method is used to match the fluid-solid boundary conditions based on the method of characteristics. This multi-domain approach allows for significant reductions of the number of grid points in the azimuthal direction for the inner grid domain and thus for corresponding increases of the time step and enhancements of computational efficiency. The viability and accuracy of the proposed method has been rigorously tested and verified through comparisons with analytical solutions as well as with the results obtained with a corresponding, previously published, and independently bench-marked solution for 2D Cartesian coordinates. Finally, the proposed numerical solution also satisfies the reciprocity theorem, which indicates that the inherent singularity associated with the origin of the polar coordinate system is adequately handled.

Cerebrospinal fluid PCR analysis and biochemistry in bodies with severe decomposition.

The aim of this study was to assess whether Neisseria meningitidis, Listeria monocytogenes, Streptococcus pneumoniae and Haemophilus influenzae can be identified using the polymerase chain reaction technique in the cerebrospinal fluid of severely decomposed bodies with known, noninfectious causes of death or whether postmortem changes can lead to false positive results and thus erroneous diagnostic information. Biochemical investigations, postmortem bacteriology and real-time polymerase chain reaction analysis in cerebrospinal fluid were performed in a series of medico-legal autopsies that included noninfectious causes of death with decomposition, bacterial meningitis without decomposition, bacterial meningitis with decomposition, low respiratory tract infections with decomposition and abdominal infections with decomposition. In noninfectious causes of death with decomposition, postmortem investigations failed to reveal results consistent with generalized inflammation or bacterial infections at the time of death. Real-time polymerase chain reaction analysis in cerebrospinal fluid did not identify the studied bacteria in any of these cases. The results of this study highlight the usefulness of molecular approaches in bacteriology as well as the use of alternative biological samples in postmortem biochemistry in order to obtain suitable information even in corpses with severe decompositional changes.

Kinetic study of neuropeptide Y (NPY) proteolysis in blood and identification of NPY3-35: a new peptide generated by plasma kallikrein.

There is little information on how neuropeptide Y (NPY) proteolysis by peptidases occurs in serum, in part because reliable techniques are lacking to distinguish different NPY immunoreactive forms and also because the factors affecting the expression of these enzymes have been poorly studied. In the present study, LC-MS/MS was used to identify and quantify NPY fragments resulting from peptidolytic cleavage of NPY(1-36) upon incubation with human serum. Kinetic studies indicated that NPY(1-36) is rapidly cleaved in serum into 3 main fragments with the following order of efficacy: NPY(3-36) &gt;&gt; NPY(3-35) &gt; NPY(2-36). Trace amounts of additional NPY forms were identified by accurate mass spectrometry. Specific inhibitors of dipeptidyl peptidase IV, kallikrein, and aminopeptidase P prevented the production of NPY(3-36), NPY(3-35), and NPY(2-36), respectively. Plasma kallikrein at physiological concentrations converted NPY(3-36) into NPY(3-35). Receptor binding assays revealed that NPY(3-35) is unable to bind to NPY Y1, Y2, and Y5 receptors; thus NPY(3-35) may represent the major metabolic clearance product of the Y2/Y5 agonist, NPY(3-36).

Role of prohormone convertases in pro-neuropeptide Y processing: coexpression and in vitro kinetic investigations.

Proneuropeptide Y (ProNPY) undergoes cleavage at a single dibasic site Lys38-Arg39 resulting in the formation of 1-39 amino acid NPY which is further processed successively by carboxypeptidase-like and peptidylglycine alpha-amidating monooxygenase enzymes. To investigate whether prohormone convertases are involved in ProNPY processing, a vaccinia virus derived expression system was used to coexpress recombinant ProNPY with each of the prohormone convertases PC1/3, PC2, furin, and PACE4 in Neuro2A and NIH 3T3 cell lines as regulated neuroendocrine and constitutive prototype cell lines, respectively. The analysis of processed products shows that only PC1/3 generates NPY in NIH 3T3 cells while both PC1/3 and PC2 are able to generate NPY in Neuro2A cells. The convertases furin and PACE4 are unable to process ProNPY in either cell line. Moreover, comparative in vitro cleavage of recombinant NPY precursor by the enzymes PC1/3, PC2 and furin shows that only PC1/3 and PC2 are involved in specific cleavage of the dibasic site. Kinetic studies demonstrate that PC1/3 cleaves ProNPY more efficiently than PC2. The main difference between the cleavage efficiency is observed in the Vmax values whereas no major difference is observed in Km values. In addition the cleavage by PC1/3 and PC2 of two peptides reproducing the dibasic cleavage site with different amino acid sequence lengths namely (20-49)-ProNPY and (28-43)-ProNPY was studied. These shortened ProNPY substrates, when recognized by the enzymes, are more efficiently cleaved than ProNPY itself. The shortest peptide is not cleaved by PC2 while it is by PC1/3. On the basis of these observations it is proposed, first, that the constitutive secreted NPY does not result from the cleavage carried out by ubiquitously expressed enzymes furin and PACE4; second, that PC1/3 and PC2 are not equipotent in the cleavage of ProNPY; and third, substrate peptide length might discriminate PC1/3 and PC2 processing activity.

The kinetic parameters and energy cost of the Hsp70 chaperone as a polypeptide unfoldase.

Hsp70-Hsp40-NEF and possibly Hsp100 are the only known molecular chaperones that can use the energy of ATP to convert stably pre-aggregated polypeptides into natively refolded proteins. However, the kinetic parameters and ATP costs have remained elusive because refolding reactions have only been successful with a molar excess of chaperones over their polypeptide substrates. Here we describe a stable, misfolded luciferase species that can be efficiently renatured by substoichiometric amounts of bacterial Hsp70-Hsp40-NEF. The reactivation rates increased with substrate concentration and followed saturation kinetics, thus allowing the determination of apparent V(max)' and K(m)' values for a chaperone-mediated renaturation reaction for the first time. Under the in vitro conditions used, one Hsp70 molecule consumed five ATPs to effectively unfold a single misfolded protein into an intermediate that, upon chaperone dissociation, spontaneously refolded to the native state, a process with an ATP cost a thousand times lower than expected for protein degradation and resynthesis.

Gait disturbances as specific predictive markers of the first fall onset in elderly people: a two-year prospective observational study

Falls are common in the elderly, and potentially result in injury and disability. Thus, preventing falls as soon as possible in older adults is a public health priority, yet there is no specific marker that is predictive of the first fall onset. We hypothesized that gait features should be the most relevant variables for predicting the first fall. Clinical baseline characteristics (e.g., gender, cognitive function) were assessed in 259 home-dwelling people aged 66 to 75 that had never fallen. Likewise, global kinetic behavior of gait was recorded from 22 variables in 1036 walking tests with an accelerometric gait analysis system. Afterward, monthly telephone monitoring reported the date of the first fall over 24 months. A principal components analysis was used to assess the relationship between gait variables and fall status in four groups: non-fallers, fallers from 0 to 6 months, fallers from 6 to 12 months and fallers from 12 to 24 months. The association of significant principal components (PC) with an increased risk of first fall was then evaluated using the area under the Receiver Operator Characteristic Curve (ROC). No effect of clinical confounding variables was shown as a function of groups. An eigenvalue decomposition of the correlation matrix identified a large statistical PC1 (termed "Global kinetics of gait pattern"), which accounted for 36.7% of total variance. Principal component loadings also revealed a PC2 (12.6% of total variance), related to the "Global gait regularity." Subsequent ANOVAs showed that only PC1 discriminated the fall status during the first 6 months, while PC2 discriminated the first fall onset between 6 and 12 months. After one year, any PC was associated with falls. These results were bolstered by the ROC analyses, showing good predictive models of the first fall during the first six months or from 6 to 12 months. Overall, these findings suggest that the performance of a standardized walking test at least once a year is essential for fall prevention.

Structure-Based, Rational Design of T Cell Receptors.

Adoptive cell transfer using engineered T cells is emerging as a promising treatment for metastatic melanoma. Such an approach allows one to introduce T cell receptor (TCR) modifications that, while maintaining the specificity for the targeted antigen, can enhance the binding and kinetic parameters for the interaction with peptides (p) bound to major histocompatibility complexes (MHC). Using the well-characterized 2C TCR/SIYR/H-2K(b) structure as a model system, we demonstrated that a binding free energy decomposition based on the MM-GBSA approach provides a detailed and reliable description of the TCR/pMHC interactions at the structural and thermodynamic levels. Starting from this result, we developed a new structure-based approach, to rationally design new TCR sequences, and applied it to the BC1 TCR targeting the HLA-A2 restricted NY-ESO-1157-165 cancer-testis epitope. Fifty-four percent of the designed sequence replacements exhibited improved pMHC binding as compared to the native TCR, with up to 150-fold increase in affinity, while preserving specificity. Genetically engineered CD8(+) T cells expressing these modified TCRs showed an improved functional activity compared to those expressing BC1 TCR. We measured maximum levels of activities for TCRs within the upper limit of natural affinity, K D = ∼1 - 5 μM. Beyond the affinity threshold at K D < 1 μM we observed an attenuation in cellular function, in line with the "half-life" model of T cell activation. Our computer-aided protein-engineering approach requires the 3D-structure of the TCR-pMHC complex of interest, which can be obtained from X-ray crystallography. We have also developed a homology modeling-based approach, TCRep 3D, to obtain accurate structural models of any TCR-pMHC complexes when experimental data is not available. Since the accuracy of the models depends on the prediction of the TCR orientation over pMHC, we have complemented the approach with a simplified rigid method to predict this orientation and successfully assessed it using all non-redundant TCR-pMHC crystal structures available. These methods potentially extend the use of our TCR engineering method to entire TCR repertoires for which no X-ray structure is available. We have also performed a steered molecular dynamics study of the unbinding of the TCR-pMHC complex to get a better understanding of how TCRs interact with pMHCs. This entire rational TCR design pipeline is now being used to produce rationally optimized TCRs for adoptive cell therapies of stage IV melanoma.