Model structure of human APOBEC3G.

BACKGROUND: APOBEC3G (apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like 3G) has antiretroviral activity associated with the hypermutation of viral DNA through cytosine deamination. APOBEC3G has two cytosine deaminase (CDA) domains; the catalytically inactive amino-terminal domain of APOBEC3G (N-CDA) carries the Vif interaction domain. There is no 3-D structure of APOBEC3G solved by X-ray or nuclear magnetic resonance. METHODOLOGY/PRINCIPAL FINDINGS: We predicted the structure of human APOBEC3G based on the crystal structure of APOBEC2. To assess the model structure, we evaluated 48 mutants of APOBEC3G N-CDA that identify novel variants altering DeltaVif HIV-1 infectivity and packaging of APOBEC3G. Results indicated that the key residue D128 is exposed at the surface of the model, with a negative local electrostatic potential. Mutation D128K changes the sign of that local potential. In addition, two novel functionally relevant residues that result in defective APOBEC3G encapsidation, R122 and W127, cluster at the surface. CONCLUSIONS/SIGNIFICANCE: The structure model identifies a cluster of residues important for packaging of APOBEC3G into virions, and may serve to guide functional analysis of APOBEC3G.

Modeling of the TCR-MHC-peptide complex.

The methodology for generating a homology model of the T1 TCR-PbCS-K(d) class I major histocompatibility complex (MHC) class I complex is presented. The resulting model provides a qualitative explanation of the effect of over 50 different mutations in the region of the complementarity determining region (CDR) loops of the T cell receptor (TCR), the peptide and the MHC's alpha(1)/alpha(2) helices. The peptide is modified by an azido benzoic acid photoreactive group, which is part of the epitope recognized by the TCR. The construction of the model makes use of closely related homologs (the A6 TCR-Tax-HLA A2 complex, the 2C TCR, the 14.3.d TCR Vbeta chain, the 1934.4 TCR Valpha chain, and the H-2 K(b)-ovalbumine peptide), ab initio sampling of CDR loops conformations and experimental data to select from the set of possibilities. The model shows a complex arrangement of the CDR3alpha, CDR1beta, CDR2beta and CDR3beta loops that leads to the highly specific recognition of the photoreactive group. The protocol can be applied systematically to a series of related sequences, permitting the analysis at the structural level of the large TCR repertoire specific for a given peptide-MHC complex.

Fracture network characterization using hydrological and geophysical data

Les décisions de gestion des eaux souterraines doivent souvent être justiffées par des modèles quantitatifs d'aquifères qui tiennent compte de l'hétérogénéité des propriétés hydrauliques. Les aquifères fracturés sont parmi les plus hétérogènes et très difficiles à étudier. Dans ceux-ci, les fractures connectées, d'ouverture millimètrique, peuvent agir comme conducteurs hydrauliques et donc créer des écoulements très localisés. Le manque général d'informations sur la distribution spatiale des fractures limite la possibilité de construire des modèles quantitatifs de flux et de transport. Les données qui conditionnent les modèles sont généralement spatialement limitées, bruitées et elles ne représentent que des mesures indirectes de propriétés physiques. Ces limitations aux données peuvent être en partie surmontées en combinant différents types de données, telles que les données hydrologiques et de radar à pénétration de sol plus commun ément appelé géoradar. L'utilisation du géoradar en forage est un outil prometteur pour identiffer les fractures individuelles jusqu'à quelques dizaines de mètres dans la formation. Dans cette thèse, je développe des approches pour combiner le géoradar avec les données hydrologiques affn d'améliorer la caractérisation des aquifères fracturés. Des investigations hydrologiques intensives ont déjà été réalisées à partir de trois forage adjacents dans un aquifère cristallin en Bretagne (France). Néanmoins, la dimension des fractures et la géométrie 3-D des fractures conductives restaient mal connue. Affn d'améliorer la caractérisation du réseau de fractures je propose dans un premier temps un traitement géoradar avancé qui permet l'imagerie des fractures individuellement. Les résultats montrent que les fractures perméables précédemment identiffées dans les forages peuvent être caractérisées géométriquement loin du forage et que les fractures qui ne croisent pas les forages peuvent aussi être identiffées. Les résultats d'une deuxième étude montrent que les données géoradar peuvent suivre le transport d'un traceur salin. Ainsi, les fractures qui font partie du réseau conductif et connecté qui dominent l'écoulement et le transport local sont identiffées. C'est la première fois que le transport d'un traceur salin a pu être imagé sur une dizaines de mètres dans des fractures individuelles. Une troisième étude conffrme ces résultats par des expériences répétées et des essais de traçage supplémentaires dans différentes parties du réseau local. En outre, la combinaison des données de surveillance hydrologique et géoradar fournit la preuve que les variations temporelles d'amplitude des signaux géoradar peuvent nous informer sur les changements relatifs de concentrations de traceurs dans la formation. Par conséquent, les données géoradar et hydrologiques sont complémentaires. Je propose ensuite une approche d'inversion stochastique pour générer des modèles 3-D de fractures discrètes qui sont conditionnés à toutes les données disponibles en respectant leurs incertitudes. La génération stochastique des modèles conditionnés par géoradar est capable de reproduire les connexions hydrauliques observées et leur contribution aux écoulements. L'ensemble des modèles conditionnés fournit des estimations quantitatives des dimensions et de l'organisation spatiale des fractures hydrauliquement importantes. Cette thèse montre clairement que l'imagerie géoradar est un outil utile pour caractériser les fractures. La combinaison de mesures géoradar avec des données hydrologiques permet de conditionner avec succès le réseau de fractures et de fournir des modèles quantitatifs. Les approches présentées peuvent être appliquées dans d'autres types de formations rocheuses fracturées où la roche est électriquement résistive.

Structural joint inversion of time-lapse crosshole ERT and GPR traveltime data

Time-lapse geophysical monitoring and inversion are valuable tools in hydrogeology for monitoring changes in the subsurface due to natural and forced (tracer) dynamics. However, the resulting models may suffer from insufficient resolution, which leads to underestimated variability and poor mass recovery. Structural joint inversion using cross-gradient constraints can provide higher-resolution models compared with individual inversions and we present the first application to time-lapse data. The results from a synthetic and field vadose zone water tracer injection experiment show that joint 3-D time-lapse inversion of crosshole electrical resistance tomography (ERT) and ground penetrating radar (GPR) traveltime data significantly improve the imaged characteristics of the point injected plume, such as lateral spreading and center of mass, as well as the overall consistency between models. The joint inversion method appears to work well for cases when one hydrological state variable (in this case moisture content) controls the time-lapse response of both geophysical methods. Citation: Doetsch, J., N. Linde, and A. Binley (2010), Structural joint inversion of time-lapse crosshole ERT and GPR traveltime data, Geophys. Res. Lett., 37, L24404, doi: 10.1029/2010GL045482.

Integration of high-resolution geophysical data for the characterization of alluvial aquifers : new approaches and implications for hydrological predictions

Abstract Accurate characterization of the spatial distribution of hydrological properties in heterogeneous aquifers at a range of scales is a key prerequisite for reliable modeling of subsurface contaminant transport, and is essential for designing effective and cost-efficient groundwater management and remediation strategies. To this end, high-resolution geophysical methods have shown significant potential to bridge a critical gap in subsurface resolution and coverage between traditional hydrological measurement techniques such as borehole log/core analyses and tracer or pumping tests. An important and still largely unresolved issue, however, is how to best quantitatively integrate geophysical data into a characterization study in order to estimate the spatial distribution of one or more pertinent hydrological parameters, thus improving hydrological predictions. Recognizing the importance of this issue, the aim of the research presented in this thesis was to first develop a strategy for the assimilation of several types of hydrogeophysical data having varying degrees of resolution, subsurface coverage, and sensitivity to the hydrologic parameter of interest. In this regard a novel simulated annealing (SA)-based conditional simulation approach was developed and then tested in its ability to generate realizations of porosity given crosshole ground-penetrating radar (GPR) and neutron porosity log data. This was done successfully for both synthetic and field data sets. A subsequent issue that needed to be addressed involved assessing the potential benefits and implications of the resulting porosity realizations in terms of groundwater flow and contaminant transport. This was investigated synthetically assuming first that the relationship between porosity and hydraulic conductivity was well-defined. Then, the relationship was itself investigated in the context of a calibration procedure using hypothetical tracer test data. Essentially, the relationship best predicting the observed tracer test measurements was determined given the geophysically derived porosity structure. Both of these investigations showed that the SA-based approach, in general, allows much more reliable hydrological predictions than other more elementary techniques considered. Further, the developed calibration procedure was seen to be very effective, even at the scale of tomographic resolution, for predictions of transport. This also held true at locations within the aquifer where only geophysical data were available. This is significant because the acquisition of hydrological tracer test measurements is clearly more complicated and expensive than the acquisition of geophysical measurements. Although the above methodologies were tested using porosity logs and GPR data, the findings are expected to remain valid for a large number of pertinent combinations of geophysical and borehole log data of comparable resolution and sensitivity to the hydrological target parameter. Moreover, the obtained results allow us to have confidence for future developments in integration methodologies for geophysical and hydrological data to improve the 3-D estimation of hydrological properties.

Functional analysis of human POT1 and RPA : insights into the role of single stranded DNA binding proteins at telomeres

Abstract Telomeres, the natural ends of chromosomes, need to be protected from chromosome end fusions, aberrant homologous recombination and degradation. In humans, chromosome ends are specified through arrays of tandemly repeated 5'-TTAGGG-3' hexamers, ending in a 3' overhang. A complex formed by the six proteins TRF1, TRF2, hRap1, TIN2, TPP1 and POT1 specifically assocìates with and protects telomeres. Telomeres are maintained by semiconservative DNA replication and by a specialized reverse transcriptase, telomerase, that carries an RNA subunit which templates new telomeric repeat synthesis. The telomeric single stranded (ss) DNA binding protein POT1 protects the telomeric 3' overhang and modulates telomerase-mediated telomere elongation. It is possible that POT1 also influences DNA synthesis during semiconservative DNA replication, which is initiated by the DNA polymerase alpha-primase complex. The heterotrimeric ss DNA-binding protein RPA plays essential roles during DNA replication. RPA binds to ss DNA with high affinity in order to stabilize ss DNA and facilitate nascent strand synthesis at the replication fork. Here we investigate how the two proteins RPA and POT1 contribute to telomere maintenance by regulating semi-conservative DNA replication and telomerase. Using chromatin immunoprecipitation experiments, we show that RPA associates with telomeres during S-phase. Analysis of telomere structure in cells shRNA-depleted for RPA and POT1 reveals that loss of RPA and POT1 causes exposure of single-stranded DNA at telomeres, suggestive of incomplete DNA replication. Biochemical experiments using purified recombinant POT1 and RPA show that saturating telomeric oligonucleotides with POT1 or RPA reduces the primase activity of the DNA polymerase alpha-primase complex and the overall activity of telomerase. POT1 and RPA also increase the primer extension by DNA polymerase alpha-primase complex and the processivity of telomerase under certain conditions, although POT1 increases the activities to a greater extent than RPA. We propose that POT1 is required for proper replication of the lagging strand of telomeres and that some phenotypes observed in POT1-depleted cells may stern from incomplete DNA replication rather than de-protection of the single-stranded overhang. Résumé Les télomères, les extrémités normales des chromosomes linéaires, doivent être protégés des fusions chromosomiques, d'événements de recombinaison homologue aberrants et de phénomènes de dégradation. Chez l'Homme, les extrémités des chromosomes sont constitués d'ADN double brin répétitif de séquence 5'-TTAGGG-3', d'une extension simple brin 3' sortante et d'un complexe protéique formé des six facteurs TRF1, TRF2, hRap1, TIN2, TPP1 et POT1 qui, s'associant à cette séquence, protègent l'ADN télomèrique. Les télomères sont maintenus par la télomérase, une transcriptase inverse capable d'allonger l'extension 3' sortante télomérique. POT1 lie l'ADN simple brin télomérique et module l'élongation des télomères par la télomérase. POT1 pourrait en théorie également influencer la réplication semi-conservative de l'ADN. L'ADN-polymérase Pal alpha-primase amorce et initie la synthèse d'ADN. Pendant la réplication, l'ADN simple brin est stabilisé par RPA, un complexe hétérotrimèrique qui lie l'ADN simple brin. RPA facilite la synthèse du brin naissant à la fourche de réplication. Ici nous avons étudié comment ces deux protéines qui lient l'ADN simple brin, RPA et POT1, régulent la réplication des télomères par la télomérase et la machinerie classique de réplication de l'ADN. Par immunoprécipitation de chromatine (ChIP), nous montrons que RPA est localisé aux télomères lors de la phase S du cycle cellulaire. De plus, l'analyse de la structure des télomeres indique que !a perte de RPA ou de POT1 conduit à l'apparition d'ADN simple brin télomérique, suggérant une réplication incomplète de l'ADN télomérique in vivo. Par une approche complémentaire biochimique utilisant les protéines POT1 et RPA recombinantes purifiées, nous montrons également que la liaison de POT1 ou de RPA à des oligonucléotides télomériques bloque l'activité primase du complexe polymérase alpha/primase et réduit l'activité télomérase sur ces substrats. En revanche, leur liaison augmente l'activité ADN-polymérase du complexe polymérase alpha/primase, ainsi que fa processivité de la télomérase dans certaines conditions, POT1 étant le plus efficace des deux facteurs. Nous proposons que POT1 est nécessaire à la réplication du brin retardé au niveau des télomères, ce qui suggère que certains phénotypes des cellules déplétés en POT1 puissent résulter d'une réplication incomplète de l'ADN télémétrique plutôt que d'une déprotection de l'extrémité sortante des télomères.

Radiolabeled fragments of monoclonal antibodies against carcinoembryonic antigen for localization of human colon carcinoma grafted into nude mice.

Four monoclonal antibodies against carcinoembryonic antigen (CEA) have been selected from 32 hybrids that produce antibodies against this antigen, by the criteria of high affinity for CEA and low cross-reactivity with granulocyte glycoprotein(s). The specificity of tumor localization in vivo of the four MAb, and their F(ab')2 and Fab fragments was compared in nude mice bearing grafts of a serially transplanted, CEA-producing, human colon carcinoma. The distribution of radiolabeled MAb and their fragments after intravenous injection was analyzed by direct measurement of radioactivity in tumor and normal organs, as well as by whole-body scanning and by autoradiography of tumor sections. Paired labeling experiments, in which 131I-labeled antibody or fragments and 125I-labeled control IgG are injected simultaneously, were undertaken to determine the relative tumor uptakes of each labeled protein. The tumor antibody uptake divided by that of control IgG defines the specificity index of localization. Tumor antibody uptakes (as compared with the whole mouse), ranging between 7 and 15, and specificity indices ranging between 3.4 and 6.8, were obtained with the four intact MAb at day 4-5 after injection. With F(ab')2 fragments of the four MAb, at day 3, the tumor antibody uptakes ranged between 12 and 24 and the specificity indices between 5.3 and 8.2. With the Fab fragments prepared from the two most promising MAb, the antibody uptakes reached values of 34 and 82 at day 2-3 and the specificity indices were as high as 12 and 19. The scanning results paralleled those obtained by direct measurement of radioactivity. With intact MAb, tumor grafts of 0.5-1 g gave very contrasted positive scans 3 d after injection. Using MAb fragments, tumors of smaller size were detectable earlier. The best results were obtained with Fab fragments of MAb 35, which gave clear detections of tumors weighing only 0.1 g as early as 48 h after injection. Autoradiographs of tumor sections from mice injected with 125I-labeled MAb demonstrated that the radioactivity was localized in the tumor tissues and not in the stromal connective tissue of mouse origin. The highest radioactivity concentration was localized in areas known to contain CEA such as the pseudolumen of glands and the apical side of carcinoma cells. The penetration of radioactivity in the central part of tumor nodules and the pseudolumen appeared to be increased with the use of MAb fragments.

A surface-based approach to quantify local cortical gyrification.

The high complexity of cortical convolutions in humans is very challenging both for engineers to measure and compare it, and for biologists and physicians to understand it. In this paper, we propose a surface-based method for the quantification of cortical gyrification. Our method uses accurate 3-D cortical reconstruction and computes local measurements of gyrification at thousands of points over the whole cortical surface. The potential of our method to identify and localize precisely gyral abnormalities is illustrated by a clinical study on a group of children affected by 22q11 Deletion Syndrome, compared to control individuals.

Robotic Gait Training Is not Superior to Conventional Treadmill Training in Parkinson Disease: A Single-Blind Randomized Controlled Trial.

BACKGROUND: The use of robots for gait training in Parkinson disease (PD) is growing, but no evidence points to an advantage over the standard treadmill. METHODS: In this randomized, single-blind controlled trial, participants aged <75 years with early-stage PD (Hoehn-Yahr <3) were randomly allocated to 2 groups: either 30 minutes of gait training on a treadmill or in the Lokomat for 3 d/wk for 4 weeks. Patients were evaluated by a physical therapist blinded to allocation before and at the end of treatment and then at the 3- and 6-month follow-up. The primary outcome measure was the 6-minute walk test. RESULTS: Of 334 screened patients, the authors randomly allocated 30 to receive gait training with treadmill or the Lokomat. At baseline, the 2 groups did not differ. At the 6-month follow-up, both groups had improved significantly in the primary outcome measure (treadmill: mean = 490.95 m, 95% confidence interval [CI] = 448.56-533.34, P = .0006; Lokomat: 458.6 m, 95% CI = 417.23-499.96, P = .01), but no significant differences were found between the 2 groups (P = .53). DISCUSSION: Robotic gait training with the Lokomat is not superior to treadmill training in improving gait performance in patients with PD. Both approaches are safe, with results maintained for up to 6 months.

Heel and toe clearance estimation for gait analysis using wireless inertial sensors.

Tripping is considered a major cause of fall in older people. Therefore, foot clearance (i.e., height of the foot above ground during swing phase) could be a key factor to better understand the complex relationship between gait and falls. This paper presents a new method to estimate clearance using a foot-worn and wireless inertial sensor system. The method relies on the computation of foot orientation and trajectory from sensors signal data fusion, combined with the temporal detection of toe-off and heel-strike events. Based on a kinematic model that automatically estimates sensor position relative to the foot, heel and toe trajectories are estimated. 2-D and 3-D models are presented with different solving approaches, and validated against an optical motion capture system on 12 healthy adults performing short walking trials at self-selected, slow, and fast speed. Parameters corresponding to local minimum and maximum of heel and toe clearance were extracted and showed accuracy ± precision of 4.1 ± 2.3 cm for maximal heel clearance and 1.3 ± 0.9 cm for minimal toe clearance compared to the reference. The system is lightweight, wireless, easy to wear and to use, and provide a new and useful tool for routine clinical assessment of gait outside a dedicated laboratory.

Transformed and nontransformed human T lymphocytes migrate to skin in a chimeric human skin/SCID mouse model.

To study human T cell migration to human skin in vivo, we grafted severe combined immunodeficient mice with 500-microm thick human skin. Two weeks after grafting, epidermal and dermal structures in the grafts were of human origin. When we intraperitoneally injected grafted mice with clones of the human HUT-78 T cell line derived from a patient with cutaneous T cell lymphoma and Sézary syndrome, we detected in the grafts the rare Vbeta23-Jbeta1.2 T cell receptor transcripts characteristic for the HUT-78 clones. These signals were found 2-6 d after cell injection in about 40% of the grafted and HUT-78 cell injected mice but not in grafts from mice that received no exogenous T cells. In contrast to HUT-78 cells, which only accumulate in low number, grafts topically challenged with nickel sufate in vaseline from mice that were injected with autologous nickel-reactive T cell lines led to massive accumulation of T cells within 3 d. Only scattered T cells accumulated in the skin when grafted mice received vaseline plus T cells, nickel sulfate alone, T cells alone, or nickel sulfate plus an allogeneic nickel-nonreactive T cell clone. When the T cell lines were labeled with the fluorochrome PKH-26 before cell injection, spots of fluorescent label in the size and shape of cells were found in the grafts challenged with nickel. Together, these results clearly demonstrate that human T cells can migrate to human skin in this chimeric human/mouse model.

Measurements of seismic attenuation and transient fluid pressure in partially saturated Berea sandstone: Evidence of fluid flow on the mesoscopic scale

A novel laboratory technique is proposed to investigate wave-induced fluid flow on the mesoscopic scale as a mechanism for seismic attenuation in partially saturated rocks. This technique combines measurements of seismic attenuation in the frequency range from 1 to 100?Hz with measurements of transient fluid pressure as a response of a step stress applied on top of the sample. We used a Berea sandstone sample partially saturated with water. The laboratory results suggest that wave-induced fluid flow on the mesoscopic scale is dominant in partially saturated samples. A 3-D numerical model representing the sample was used to verify the experimental results. Biot's equations of consolidation were solved with the finite-element method. Wave-induced fluid flow on the mesoscopic scale was the only attenuation mechanism accounted for in the numerical solution. The numerically calculated transient fluid pressure reproduced the laboratory data. Moreover, the numerically calculated attenuation, superposed to the frequency-independent matrix anelasticity, reproduced the attenuation measured in the laboratory in the partially saturated sample. This experimental?numerical fit demonstrates that wave-induced fluid flow on the mesoscopic scale and matrix anelasticity are the dominant mechanisms for seismic attenuation in partially saturated Berea sandstone.

Transient Ischemic Dilation (TID) ratio measurement during Rb-82 cardiac PET: Comparison between three different software packages.

Background: TIDratio indirectly reflects myocardial ischemia and is correlated with cardiacprognosis. We aimed at comparing the influence of three different softwarepackages for the assessment of TID using Rb-82 cardiac PET/CT. Methods: Intotal, data of 30 patients were used based on normal myocardial perfusion(SSS<3 and SRS<3) and stress myocardial blood flow 2mL/min/g)assessed by Rb-82 cardiac PET/CT. After reconstruction using 2D OSEM (2Iterations, 28 subsets), 3-D filtering (Butterworth, order=10, ωc=0.5), data were automatically processed, and then manually processed fordefining identical basal and apical limits on both stress and rest images.TIDratio were determined with Myometrix®, ECToolbox® and QGS®software packages. Comparisons used ANOVA, Student t-tests and Lin concordancetest (ρc). Results: All of the 90 processings were successfullyperformed. TID ratio were not statistically different between software packageswhen data were processed automatically (P=0.2) or manually (P=0.17). There was a slight, butsignificant relative overestimation of TID with automatic processing incomparison to manual processing using ECToolbox® (1.07 ± 0.13 vs 1.0± 0.13, P=0.001)and Myometrix® (1.07 ± 0.15 vs 1.01 ± 0.11, P=0.003) but not using QGS®(1.02 ±0.12 vs 1.05 ± 0.11, P=0.16). The best concordance was achieved between ECToolbox®and Myometrix® manual (ρc=0.67) processing.Conclusion: Using automatic or manual mode TID estimation was not significantlyinfluenced by software type. Using Myometrix® or ECToolbox®TID was significantly different between automatic and manual processing, butnot using QGS®. Software package should be account for when definingTID normal reference limits, as well as when used in multicenter studies. QGS®software seemed to be the most operator-independent software package, whileECToolbox® and Myometrix® produced the closest results.

Exporting Contours to DICOM-RT Structure Set

This paper presents an ITK implementation for exportingthe contours of the automated segmentation results toDICOM-RT Structure Set format. The âeurooeradiotherapystructure setâeuro (RTSTRUCT) object of the DICOM standard isused for the transfer of patient structures and relateddata, between the devices found within and outside theradiotherapy department. It mainly contains theinformation of regions of interest (ROIs) and points ofinterest (E.g. dose reference points). In many cases,rather than manually drawing these ROIs on the CT images,one can indeed benefit from the automated segmentationalgorithms already implemented in ITK. But at present, itis not possible to export the ROIs obtained from ITK toRTSTRUCT format. In order to bridge this gap, we havedeveloped a framework for exporting contour data toRTSTRUCT. We provide here the complete implementation ofRTSTRUCT exporter and present the details of the pipelineused. Results on a 3-D CT image of the Head and Neck(H&N) region are presented.

ZHARP: three-dimensional motion tracking from a single image plane.

Three-dimensional imaging and quantification of myocardial function are essential steps in the evaluation of cardiac disease. We propose a tagged magnetic resonance imaging methodology called zHARP that encodes and automatically tracks myocardial displacement in three dimensions. Unlike other motion encoding techniques, zHARP encodes both in-plane and through-plane motion in a single image plane without affecting the acquisition speed. Postprocessing unravels this encoding in order to directly track the 3-D displacement of every point within the image plane throughout an entire image sequence. Experimental results include a phantom validation experiment, which compares zHARP to phase contrast imaging, and an in vivo study of a normal human volunteer. Results demonstrate that the simultaneous extraction of in-plane and through-plane displacements from tagged images is feasible.