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

Do mammalian cytochrome P450s exhibit multiple ligand access pathways and ligand channelling?

Understanding substrate binding and product release in cytochrome P450 (CYP) enzymes is important for explaining their key role in drug metabolism, toxicity, xenobiotic degradation and biosynthesis. Here, molecular simulations of substrate and product exit from the buried active site of a mammalian P450, the microsomal CYP2C5, identified a dominant exit channel, termed pathway (pw) 2c. Previous simulations with soluble bacterial P450s showed a different dominant egress channel, pw2a. Combining these, we propose two mechanisms in CYP2C5: (i) a one-way route by which lipophilic substrates access the enzyme from the membrane by pw2a and hydroxylated products egress along pw2c; and (ii) a two-way route for access and egress, along pw2c, for soluble compounds. The proposed differences in substrate access and product egress routes between membrane-bound mammalian P450s and soluble bacterial P450s highlight the adaptability of the P450 fold to the requirements of differing cellular locations and substrate specificity profiles

Computational study and peptide inhibitors design for the CDK9 – cyclin T1 complex

Cyclin dependent kinase 9 (CDK9) is a protein that belongs to the cyclin-dependent kinases family, and its main role is in the regulation of the cell transcription processes. Since the increased activity of CDK9 is connected with the development of pathological processes such as tumor growth and survival and HIV-1 replication, inhibition of the CDK9 could be of particular interest for treating such diseases. The activation of CDK9 is initiated by the formation of CDK9/cyclin T1 complex, therefore disruption of its formation could be a promising strategy for the design of CDK9 inhibitors. In order to assist in the design of potential inhibitors of CDK9/cyclin T1 complex formation, a computational study of the CDK9/cyclin T1 interface was conducted. Ten peptides were designed using the information from the analysis of the complex, hot spot residues and fragment based design. The designed peptides were docked to CDK9 structures obtained by molecular dynamics simulations of CDK9/cyclin T1 complex and the CDK9 alone and their binding affinities were evaluated using molecular mechanics Poisson Boltzman surface area (MM-PBSA) method and steered molecular dynamics (SMD). Designed peptide sequences LQTLGF and ESIILQ, both derived from the surface of cyclin T1, as well as the peptide sequence PRWPE, derived from fragment based design, showed the most favorable binding properties and were selected for our further studies