Investigation of Indazole Unbinding Pathways in CYP2E1 by Molecular Dynamics Simulations

Human microsomal cytochrome P450 2E1 (CYP2E1) can oxidize not only low molecular weight xenobiotic compounds such as ethanol, but also many endogenous fatty acids. The crystal structure of CYP2E1 in complex with indazole reveals that the active site is deeply buried into the protein center. Thus, the unbinding pathways and associated unbinding mechanisms remain elusive. In this study, random acceleration molecular dynamics simulations combined with steered molecular dynamics and potential of mean force calculations were performed to identify the possible unbinding pathways in CYP2E1. The results show that channel 2c and 2a are most likely the unbinding channels of CYP2E1. The former channel is located between helices G and I and the B-C loop, and the latter resides between the region formed by the F-G loop, the B-C loop and the β1 sheet. Phe298 and Phe478 act as the gate keeper during indazole unbinding along channel 2c and 2a, respectively. Previous site-directed mutagenesis experiments also supported these findings

Molecular Dynamics Analysis Reveals Structural Insights into Mechanism of Nicotine N-Demethylation Catalyzed by Tobacco Cytochrome P450 Mono-Oxygenase

CYP82E4, a cytochrome P450 monooxygenase, has nicotine N-demethylase (NND) activity, which mediates the bioconversion of nicotine into nornicotine in senescing tobacco leaves. Nornicotine is a precursor of the carcinogen, tobacco-specific nitrosamine. CYP82E3 is an ortholog of CYP82E4 with 95% sequence identity, but it lacks NND activity. A recent site-directed mutagenesis study revealed that a single amino acid substitution, i.e., cysteine to tryptophan at the 330 position in the middle of protein, restores the NND activity of CYP82E3 entirely. However, the same amino acid change caused the loss of the NND activity of CYP82E4. To determine the mechanism of the functional turnover of the two molecules, four 3D structures, i.e., the two molecules and their corresponding cys–trp mutants were modeled. The resulting structures exhibited that the mutation site is far from the active site, which suggests that no direct interaction occurs between the two sites. Simulation studies in different biological scenarios revealed that the mutation introduces a conformation drift with the largest change at the F-G loop. The dynamics trajectories analysis using principal component analysis and covariance analysis suggests that the single amino acid change causes the opening and closing of the transfer channels of the substrates, products, and water by altering the motion of the F-G and B-C loops. The motion of helix I is also correlated with the motion of both the F-G loop and the B-C loop and; the single amino acid mutation resulted in the curvature of helix I. These results suggest that the single amino acid mutation outside the active site region may have indirectly mediated the flexibility of the F-G and B-C loops through helix I, causing a functional turnover of the P450 monooxygenase

Single-Molecule Force Spectroscopy: Experiments, Analysis, and Simulations

International audienceThe mechanical properties of cells and of subcellular components are important to obtain a mechanistic molecular understanding of biological processes. The quantification of mechanical resistance of cells and biomolecules using biophysical methods matured thanks to the development of nanotechnologies such as optical and magnetic tweezers, the biomembrane force probe and atomic force microscopy (AFM). The quantitative nature of force spectroscopy measurements has converted AFM into a valuable tool in biophysics. Force spectroscopy allows the determination of the forces required to unfold protein domains and to disrupt individual receptor/ligand bonds. Molecular simulation as a computational microscope allows investigation of similar biological processes with an atomistic detail. In this chapter, we first provide a step-by-step protocol of force spectroscopy including sample preparation, measurement and analysis of force spectroscopy using AFM and its interpretation in terms of available theories. Next, we present the background for molecular dynamics (MD) simulations focusing on steered molecular dynamics (SMD) and the importance of bridging of computational tools with experimental technique

Assessing the Microbial Community and Functional Genes in a Vertical Soil Profile with Long-Term Arsenic Contamination

Conceived and designed the experiments: GW. Performed the experiments: JX GL. Analyzed the data: JX JZ GW. Contributed reagents/materials/analysis tools: ST JZ GW. Wrote the paper: JX ZH JDVN JZ GW.Arsenic (As) contamination in soil and groundwater has become a serious problem to public health. To examine how microbial communities and functional genes respond to long-term arsenic contamination in vertical soil profile, soil samples were collected from the surface to the depth of 4 m (with an interval of 1 m) after 16-year arsenic downward infiltration. Integrating BioLog and functional gene microarray (GeoChip 3.0) technologies, we showed that microbial metabolic potential and diversity substantially decreased, and community structure was markedly distinct along the depth. Variations in microbial community functional genes, including genes responsible for As resistance, carbon and nitrogen cycling, phosphorus utilization and cytochrome c oxidases were detected. In particular, changes in community structures and activities were correlated with the biogeochemical features along the vertical soil profile when using the rbcL and nifH genes as biomarkers, evident for a gradual transition from aerobic to anaerobic lifestyles. The C/N showed marginally significant correlations with arsenic resistance (p = 0.069) and carbon cycling genes (p = 0.073), and significant correlation with nitrogen fixation genes (p = 0.024). The combination of C/N, NO3− and P showed the highest correlation (r = 0.779, p = 0.062) with the microbial community structure. Contradict to our hypotheses, a long-term arsenic downward infiltration was not the primary factor, while the spatial isolation and nutrient availability were the key forces in shaping the community structure. This study provides new insights about the heterogeneity of microbial community metabolic potential and future biodiversity preservation for arsenic bioremediation management.Yeshttp://www.plosone.org/static/editorial#pee

A proportion of proteinase 3 (PR3)-specific anti-neutrophil cytoplasmic antibodies (ANCA) only react with PR3 after cleavage of its N-terminal activation dipeptide

ANCA directed against PR3 are highly specific for Wegener's granulomatosis and microscopic polyangiitis, and have been implicated in the pathogenesis of small vessel vasculitis. Most PR3-ANCA are directed against conformational epitopes on PR3. This study was designed to determine whether the cleavage of the N-terminal activation dipeptide of PR3 is required for the binding of PR3-ANCA. Recombinant PR3 (rPR3) variants were expressed in the epithelial cell line, 293. As confirmed by radiosequencing, the rPR3 secreted into the 293 cell culture supernatant is N-terminally unprocessed. Two enzymatically inactive rPR3 mutants were expressed in 293 cells: rPR3-S176A and δ-rPR3-S176A. rPR3-S176A contains the N-propetide Ala-2-Glu-1, δ-rPR3-S176A does not. Culture supernatants of rPR3-S176A and δ-rPR3-S176A expressing 293 cells were used as sources of target antigen for PR3-ANCA testing by capture ELISA. Forty unselected consecutive PR3-ANCA+ sera were tested. With δ-rPR3-S176A as antigen all 40 were recognized, compared with only 34 of 40 when rPR3-S176A served as target antigen. The majority of the serum samples contained a mixture of antibodies reacting with epitopes accessible on the mature and on the proform of PR3. In conclusion, the cleavage of the N-terminal activation dipeptide of PR3 is not an absolute requirement for recognition by all PR3-ANCA. However, a substantial proportion of PR3-ANCA recognize (a) target antigen(s) exposed only after the conformational change of PR3 associated with the N-terminal processing. In 15% of sera this PR3-ANCA subset occurred exclusively. PR3-ANCA subtypes can be differentiated using specifically designed rPR3 variants as target antigens, and non-haematopoietic mammalian cells without regulated secretory pathway can be used for their expression

A novel 6.14 Mb duplication of chromosome 8p21 in a patient with autism and self mutilation

Autism spectrum disorders (ASDs) are a group of neurodevelopmental disorders with a strong genetic etiology. Cytogenetic abnormalities have been detected in 5–10% of the patients with autism. In this study, we present the clinical, cytogenetic and array-comparative genomic hybridization (array-CGH) evaluation of a 13-year-old male with severe developmental delay, facial dysmorphic features, autism and self mutilation. The patient was found to carry a de novo duplication of chromosome region 8p21 of minimally 6.14 and maximally 6.58 Mb as ascertained by bacterial artificial chromosome (BAC)-based array-CGH. Hitherto, only a few patients with autism with cytogenetically visible duplications involving the chromosome 8p21 region have been described, but the extent of these duplications has not been determined at the molecular level. This represents the smallest rearrangement of chromosomal region 8p21 as yet found in a patient with autism. For 11 of the 36 genes with known functions located within this duplication clear transcription in the brain was found. Of those the STMN4 and DPYSL2 genes are the most likely candidate genes to be involved in neuronal development, and, if altered in gene-dosage, in the autistic phenotype of our patient