How T cell receptors interact with peptide-MHCs: A multiple steered molecular dynamics study.

The T-cell receptor (TCR) interaction with antigenic peptides (p) presented by the major histocompatibility complex (MHC) molecule is a key determinant of immune response. In addition, TCR-pMHC interactions offer examples of features more generally pertaining to protein-protein recognition: subtle specificity and cross-reactivity. Despite their importance, molecular details determining the TCR-pMHC binding remain unsolved. However, molecular simulation provides the opportunity to investigate some of these aspects. In this study, we perform extensive equilibrium and steered molecular dynamics simulations to study the unbinding of three TCR-pMHC complexes. As a function of the dissociation reaction coordinate, we are able to obtain converged H-bond counts and energy decompositions at different levels of detail, ranging from the full proteins, to separate residues and water molecules, down to single atoms at the interface. Many observed features do not support a previously proposed two-step model for TCR recognition. Our results also provide keys to interpret experimental point-mutation results. We highlight the role of water both in terms of interface resolvation and of water molecules trapped in the bound complex. Importantly, we illustrate how two TCRs with similar reactivity and structures can have essentially different binding strategies. Proteins 2011; © 2011 Wiley-Liss, Inc.

Cumulative frequency-dependent selective episodes allow for rapid morph cycles and rock-paper-scissors dynamics in species with overlapping generations.

Rock-paper-scissors (RPS) dynamics, which maintain genetic polymorphisms over time through negative frequency-dependent (FD) selection, can evolve in short-lived species with no generational overlap, where they produce rapid morph frequency cycles. However, most species have overlapping generations and thus, rapid RPS dynamics are thought to require stronger FD selection, the existence of which yet needs to be proved. Here, we experimentally demonstrate that two cumulative selective episodes, FD sexual selection reinforced by FD selection on offspring survival, generate sufficiently strong selection to generate rapid morph frequency cycles in the European common lizard Zootoca vivipara, a multi-annual species with major generational overlap. These findings show that the conditions required for the evolution of RPS games are fulfilled by almost all species exhibiting genetic polymorphisms and suggest that RPS games may be responsible for the maintenance of genetic diversity in a wide range of species.

Equine herpesvirus protein E10 induces membrane recruitment and phosphorylation of its cellular homologue, bcl-10.

v-E10, a caspase recruitment domain (CARD)-containing gene product of equine herpesvirus 2, is the viral homologue of the bcl-10 protein whose gene was found to be translocated in mucosa-associated lymphoid tissue (MALT) lymphomas. v-E10 efficiently activates the c-jun NH(2)-terminal kinase (JNK), p38 stress kinase, and the nuclear factor (NF)-kappaB transcriptional pathway and interacts with its cellular homologue, bcl-10, via a CARD-mediated interaction. Here we demonstrate that v-E10 contains a COOH-terminal geranylgeranylation consensus site which is responsible for its plasma membrane localization. Expression of v-E10 induces hyperphosphorylation and redistribution of bcl-10 from the cytoplasm to the plasma membrane, a process which is dependent on the intactness of the v-E10 CARD motif. Both membrane localization and a functional CARD motif are important for v-E10-mediated NF-kappaB induction, but not for JNK activation, which instead requires a functional v-E10 binding site for tumor necrosis factor receptor-associated factor (TRAF)6. Moreover, v-E10-induced NF-kappaB activation is inhibited by a dominant negative version of the bcl-10 binding protein TRAF1, suggesting that v-E10-induced membrane recruitment of cellular bcl-10 induces constitutive TRAF-mediated NF-kappaB activation.

Habitat-quality effects on metapopulation dynamics in greater white-toothed shrews, Crocidura russula.

The effects of patch size and isolation on metapopulation dynamics have received wide empirical support and theoretical formalization. By contrast, the effects of patch quality seem largely underinvestigated, partly due to technical difficulties in properly assessing quality. Here we combine habitat-quality modeling with four years of demographic monitoring in a metapopulation of greater white-toothed shrews (Crocidura russula) to investigate the role of patch quality on metapopulation processes. Together, local patch quality and connectivity significantly enhanced local population sizes and occupancy rates (R2 = 14% and 19%, respectively). Accounting for the quality of patches connected to the focal one and acting as potential sources improved slightly the model explanatory power for local population sizes, pointing to significant source-sink dynamics. Local habitat quality, in interaction with connectivity, also increased colonization rate (R2 = 28%), suggesting the ability of immigrants to target high-quality patches. Overall, patterns were best explained when assuming a mean dispersal distance of 800 m, a realistic value for the species under study. Our results thus provide evidence that patch quality, in interaction with connectivity, may affect major demographic processes.

Direct simulation of interface dynamics: linking capillary pressure, interfacial area and surface energy

We perform direct numerical simulations of drainage by solving Navier- Stokes equations in the pore space and employing the Volume Of Fluid (VOF) method to track the evolution of the fluid-fluid interface. After demonstrating that the method is able to deal with large viscosity contrasts and to model the transition from stable flow to viscous fingering, we focus on the definition of macroscopic capillary pressure. When the fluids are at rest, the difference between inlet and outlet pressures and the difference between the intrinsic phase average pressure coincide with the capillary pressure. However, when the fluids are in motion these quantities are dominated by viscous forces. In this case, only a definition based on the variation of the interfacial energy provides an accurate measure of the macroscopic capillary pressure and allows separating the viscous from the capillary pressure components.

Orosomucoid (alpha1-acid glycoprotein) plasma concentration and genetic variants: effects on human immunodeficiency virus protease inhibitor clearance and cellular accumulation.

BACKGROUND AND OBJECTIVE: Protease inhibitors are highly bound to orosomucoid (ORM) (alpha1-acid glycoprotein), an acute-phase plasma protein encoded by 2 polymorphic genes, which may modulate their disposition. Our objective was to determine the influence of ORM concentration and phenotype on indinavir, lopinavir, and nelfinavir apparent clearance (CL(app)) and cellular accumulation. Efavirenz, mainly bound to albumin, was included as a control drug. METHODS: Plasma and cells samples were collected from 434 human immunodeficiency virus-infected patients. Total plasma and cellular drug concentrations and ORM concentrations and phenotypes were determined. RESULTS: Indinavir CL(app) was strongly influenced by ORM concentration (n = 36) (r2 = 0.47 [P = .00004]), particularly in the presence of ritonavir (r2 = 0.54 [P = .004]). Lopinavir CL(app) was weakly influenced by ORM concentration (n = 81) (r2 = 0.18 [P = .0001]). For both drugs, the ORM1 S variant concentration mainly explained this influence (r2 = 0.55 [P = .00004] and r2 = 0.23 [P = .0002], respectively). Indinavir CL(app) was significantly higher in F1F1 individuals than in F1S and SS patients (41.3, 23.4, and 10.3 L/h [P = .0004] without ritonavir and 21.1, 13.2, and 10.1 L/h [P = .05] with ritonavir, respectively). Lopinavir cellular exposure was not influenced by ORM abundance and phenotype. Finally, ORM concentration or phenotype did not influence nelfinavir (n = 153) or efavirenz (n = 198) pharmacokinetics. CONCLUSION: ORM concentration and phenotype modulate indinavir pharmacokinetics and, to a lesser extent, lopinavir pharmacokinetics but without influencing their cellular exposure. This confounding influence of ORM should be taken into account for appropriate interpretation of therapeutic drug monitoring results. Further studies are needed to investigate whether the measure of unbound drug plasma concentration gives more meaningful information than total drug concentration for indinavir and lopinavir.

Evolutionary dynamics of coding and non-coding transcriptomes.

Gene expression changes may underlie much of phenotypic evolution. The development of high-throughput RNA sequencing protocols has opened the door to unprecedented large-scale and cross-species transcriptome comparisons by allowing accurate and sensitive assessments of transcript sequences and expression levels. Here, we review the initial wave of the new generation of comparative transcriptomic studies in mammals and vertebrate outgroup species in the context of earlier work. Together with various large-scale genomic and epigenomic data, these studies have unveiled commonalities and differences in the dynamics of gene expression evolution for various types of coding and non-coding genes across mammalian lineages, organs, developmental stages, chromosomes and sexes. They have also provided intriguing new clues to the regulatory basis and phenotypic implications of evolutionary gene expression changes.

Brownian dynamics simulation of DNA condensation.

DNA condensation observed in vitro with the addition of polyvalent counterions is due to intermolecular attractive forces. We introduce a quantitative model of these forces in a Brownian dynamics simulation in addition to a standard mean-field Poisson-Boltzmann repulsion. The comparison of a theoretical value of the effective diameter calculated from the second virial coefficient in cylindrical geometry with some experimental results allows a quantitative evaluation of the one-parameter attractive potential. We show afterward that with a sufficient concentration of divalent salt (typically approximately 20 mM MgCl(2)), supercoiled DNA adopts a collapsed form where opposing segments of interwound regions present zones of lateral contact. However, under the same conditions the same plasmid without torsional stress does not collapse. The condensed molecules present coexisting open and collapsed plectonemic regions. Furthermore, simulations show that circular DNA in 50% methanol solutions with 20 mM MgCl(2) aggregates without the requirement of torsional energy. This confirms known experimental results. Finally, a simulated DNA molecule confined in a box of variable size also presents some local collapsed zones in 20 mM MgCl(2) above a critical concentration of the DNA. Conformational entropy reduction obtained either by supercoiling or by confinement seems thus to play a crucial role in all forms of condensation of DNA.

The effects of physiologically plausible connectivity structure on local and global dynamics in large scale brain models.

Functionally relevant large scale brain dynamics operates within the framework imposed by anatomical connectivity and time delays due to finite transmission speeds. To gain insight on the reliability and comparability of large scale brain network simulations, we investigate the effects of variations in the anatomical connectivity. Two different sets of detailed global connectivity structures are explored, the first extracted from the CoCoMac database and rescaled to the spatial extent of the human brain, the second derived from white-matter tractography applied to diffusion spectrum imaging (DSI) for a human subject. We use the combination of graph theoretical measures of the connection matrices and numerical simulations to explicate the importance of both connectivity strength and delays in shaping dynamic behaviour. Our results demonstrate that the brain dynamics derived from the CoCoMac database are more complex and biologically more realistic than the one based on the DSI database. We propose that the reason for this difference is the absence of directed weights in the DSI connectivity matrix.

Brain dynamics underlying training-induced improvement in suppressing inappropriate action.

Inhibitory control, a core component of executive functions, refers to our ability to suppress intended or ongoing cognitive or motor processes. Mostly based on Go/NoGo paradigms, a considerable amount of literature reports that inhibitory control of responses to "NoGo" stimuli is mediated by top-down mechanisms manifesting ∼200 ms after stimulus onset within frontoparietal networks. However, whether inhibitory functions in humans can be trained and the supporting neurophysiological mechanisms remain unresolved. We addressed these issues by contrasting auditory evoked potentials (AEPs) to left-lateralized "Go" and right NoGo stimuli recorded at the beginning versus the end of 30 min of active auditory spatial Go/NoGo training, as well as during passive listening of the same stimuli before versus after the training session, generating two separate 2 × 2 within-subject designs. Training improved Go/NoGo proficiency. Response times to Go stimuli decreased. During active training, AEPs to NoGo, but not Go, stimuli modulated topographically with training 61-104 ms after stimulus onset, indicative of changes in the underlying brain network. Source estimations revealed that this modulation followed from decreased activity within left parietal cortices, which in turn predicted the extent of behavioral improvement. During passive listening, in contrast, effects were limited to topographic modulations of AEPs in response to Go stimuli over the 31-81 ms interval, mediated by decreased right anterior temporoparietal activity. We discuss our results in terms of the development of an automatic and bottom-up form of inhibitory control with training and a differential effect of Go/NoGo training during active executive control versus passive listening conditions.