Ensemble modeling of substrate binding to cytochromes P450: analysis of catalytic differences between CYP1A orthologs

A novel application of modeling and docking approaches involving ensembles of homology models is used to understand structural bases underlying subtle catalytic differences between related cytochromes P450 (CYPs). Mammalian CYP1A1s and fish CYP1As are orthologous enzymes with similar substrate preferences. With some substrates (3,3',4,4'-tetrachlorobiphenyl, TCB) oxidation rates differ by orders of magnitude, while others (e.g., benzo[a]pyrene; B[a]P) are oxidized at similar rates but with somewhat differing regiospecificity. These two environmental chemical substrates (TCB and B[a]P) as well as 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) were docked to multiple models of rat, human, scup, and/or killifish CYP1As, based on multiple templates, retaining multiple poses from each model, giving ensembles of docked poses for each species. With TCB, more poses were observed closer to the heme in ensembles of rat or human CYP1A1 than of killifish CYP1A. Analysis of interacting residues suggested that differences in TCB pose distributions are due primarily to Leu387 and Val230 in killifish CYP1A. In silico mutations L387V and V230G enabled TCB to dock closer to the heme in killifish CYP1A. Mutating additional interacting residues (Ala127, Thr233, Asn317, and Tyr386) of killifish CYP1A to the corresponding residues of human CYP1A1 resulted in TCB pose distributions nearly identical with those of human CYP1A1. Docking of TCDD to sets of consensus models of killifish, rat, and human CYP1As showed species differences similar to those with TCB, but with further structural constraints possibly contributing to slower oxidation of TCDD. Docking B[a]P to sets of consensus models of the human and fish CYP1As yielded frequencies of substrate orientations correlating with known regiospecificities for metabolism of B[a]P by these enzymes. The results demonstrate the utility of this ensemble modeling method, which can account for uncertainty inherent in homology modeling and docking by producing statistical distributions of ligand positions.     

Cyclic Game Dynamics Driven by Iterated Reasoning

Recent theories from complexity science argue that complex dynamics are ubiquitous in social and economic systems. These claims emerge from the analysis of individually simple agents whose collective behavior is surprisingly complicated. However, economists have argued that iterated reasoning–what you think I think you think–will suppress complex dynamics by stabilizing or accelerating convergence to Nash equilibrium. We report stable and efficient periodic behavior in human groups playing the Mod Game, a multi-player game similar to Rock-Paper-Scissors. The game rewards subjects for thinking exactly one step ahead of others in their group. Groups that play this game exhibit cycles that are inconsistent with any fixed-point solution concept. These cycles are driven by a “hopping” behavior that is consistent with other accounts of iterated reasoning: agents are constrained to about two steps of iterated reasoning and learn an additional one-half step with each session. If higher-order reasoning can be complicit in complex emergent dynamics, then cyclic and chaotic patterns may be endogenous features of real-world social and economic systems.     

An In Vivo Study of Self-Regulated Study Sequencing in Introductory Psychology Courses

Study sequence can have a profound influence on learning. In this study we investigated how students decide to sequence their study in a naturalistic context and whether their choices result in improved learning. In the study reported here, 2061 undergraduate students enrolled in an Introductory Psychology course completed an online homework tutorial on measures of central tendency, a topic relevant to an exam that counted towards their grades. One group of students was enabled to choose their own study sequence during the tutorial (Self-Regulated group), while the other group of students studied the same materials in sequences chosen by other students (Yoked group). Students who chose their sequence of study showed a clear tendency to block their study by concept, and this tendency was positively associated with subsequent exam performance. In the Yoked group, study sequence had no effect on exam performance. These results suggest that despite findings that blocked study is maladaptive when assigned by an experimenter, it may actually be adaptive when chosen by the learner in a naturalistic context.     

Novel escape mutants suggest an extensive TRIM5α binding site spanning the entire outer surface of the murine leukemia virus capsid protein

After entry into target cells, retroviruses encounter the host restriction factors such as Fv1 and TRIM5α. While it is clear that these factors target retrovirus capsid proteins (CA), recognition remains poorly defined in the absence of structural information. To better understand the binding interaction between TRIM5α and CA, we selected a panel of novel N-tropic murine leukaemia virus (N-MLV) escape mutants by a serial passage of replication competent N-MLV in rhesus macaque TRIM5α (rhTRIM5α)-positive cells using a small percentage of unrestricted cells to allow multiple rounds of virus replication. The newly identified mutations, many of which involve changes in charge, are distributed over the outer 'top' surface of N-MLV CA, including the N-terminal β-hairpin, and map up to 29 A(o) apart. Biological characterisation with a number of restriction factors revealed that only one of the new mutations affects restriction by human TRIM5α, indicating significant differences in the binding interaction between N-MLV and the two TRIM5αs, whereas three of the mutations result in dual sensitivity to Fv1(n) and Fv1(b). Structural studies of two mutants show that no major changes in the overall CA conformation are associated with escape from restriction. We conclude that interactions involving much, if not all, of the surface of CA are vital for TRIM5α binding.     

Identification,modeling and ligand affinity of early deuterostome CYP51s,and functional characterization of recombinant zebrafish sterol 14α-demethylase

Background

Sterol 14α-demethylase (cytochrome P450 51, CYP51, P450_14DM) is a microsomal enzyme that in eukaryotes catalyzes formation of sterols essential for cell membrane function and as precursors in biosynthesis of steroid hormones. Functional properties of CYP51s are unknown in non-mammalian deuterostomes.

Methods

PCR-cloning and sequencing and computational analyses (homology modeling and docking) addressed CYP51 in zebrafish Danio rerio, the reef fish sergeant major Abudefduf saxatilis, and the sea urchin Strongylocentrotus purpuratus. Following N-terminal amino acid modification, zebrafish CYP51 was expressed in Escherichia coli, and lanosterol 14α-demethylase activity and azole inhibition of CYP51 activity were characterized using GC-MS.

Results

Molecular phylogeny positioned S. purpuratus CYP51 at the base of the deuterostome clade. In zebrafish, CYP51 is expressed in all organs examined, most strongly in intestine. The recombinant protein bound lanosterol and catalyzed 14α-demethylase activity, at 3.2 nmol/min/nmol CYP51. The binding of azoles to zebrafish CYP51 gave K_S (dissociation constant) values of 0.26 μM for ketoconazole and 0.64 μM for propiconazole. Displacement of carbon monoxide also indicated zebrafish CYP51 has greater affinity for ketoconazole. Docking to homology models showed that lanosterol docks in fish and sea urchin CYP51s with an orientation essentially the same as in mammalian CYP51s. Docking of ketoconazole indicates it would inhibit fish and sea urchin CYP51s.

Conclusions

Biochemical and computational analyses are consistent with lanosterol being a substrate for early deuterostome CYP51s.

General significance

The results expand the phylogenetic view of animal CYP51, with evolutionary, environmental and therapeutic implications.