Comparison of intrinsic dynamics of cytochrome p450 proteins using normal mode analysis

Extended Sampling of Macromolecular Conformations from Uniformly Distributed Points on Multidimensional Normal Mode Hyperspheres

Proteins are dynamic entities that adopt diverse conformations, which play a pivotal role in their function. Understanding these conformations is essential, and protein collective motions, particularly those captured by normal mode (NM) and their linear combinations, provide a robust means for conformational sampling. This work introduces a novel approach to obtaining a uniformly oriented set of a given number of lowest frequency NM combined vectors and generating harmonically equidistant restrained structures along them. They are all thus uniformly located on concentric hyperspheres, systematically covering the defined NM space fully. The generated structures are further relaxed with standard molecular dynamics (MD) simulations to explore the conformational space. The efficiency of the approach we termed "distributed points Molecular Dynamics using Normal Modes" (dpMDNM) was assessed by applying it to hen egg-white lysozyme and human cytochrome P450 3A4 (CYP3A4). To this purpose, we compared our new approach with other methods and analyzed the sampling of existing experimental structures. Our results demonstrate the efficacy of dpMDNM in extensive conformational sampling, particularly as more NMs are considered. Ensembles generated by dpMDNM exhibited a broad coverage of the experimental structures, providing valuable insights into the functional aspects of lysozyme and CYP3A4. Furthermore, dpMDNM also covered lysozyme structures with relatively elevated energies corresponding to transient states not easily obtained by standard MD simulations, in conformity with nuclear magnetic resonance structural indications. This method offers an efficient and rational framework for comprehensive protein conformational sampling, contributing significantly to our understanding of protein dynamics and function.

Cytochrome P450 monooxygenase systems: Diversity and plasticity for adaptive stress response

Superfamily of cytochromes P450 (CYPs) is composed of heme-thiolate-containing monooxygenase enzymes, which play crucial roles in the biosynthesis, bioactivation, and detoxification of a variety of organic compounds, both endogenic and exogenic. Majority of CYP monooxygenase systems are multi-component and contain various redox partners, cofactors and auxiliary proteins, which contribute to their diversity in both prokaryotes and eukaryotes. Recent progress in bioinformatics and computational biology approaches make it possible to undertake whole-genome and phylogenetic analyses of CYPomes of a variety of organisms. Considerable variations in sequences within and between CYP families and high similarity in secondary and tertiary structures between all CYPs along with dramatic conformational changes in secondary structure elements of a substrate binding site during catalysis have been reported. This provides structural plasticity and substrate promiscuity, which underlie functional diversity of CYPs. Gene duplication and mutation events underlie CYP evolutionary diversity and emergence of novel selectable functions, which provide the involvement of CYPs in high adaptability to changing environmental conditions and dietary restrictions. In our review, we discuss the recent advancements and challenges in the elucidating the evolutionary origin and mechanisms underlying the CYP monooxygenase system diversity and plasticity. Our review is in the view of hypothesis that diversity of CYP monooxygenase systems is translated into the broad metabolic profiles, and this has been acquired during the long evolutionary time to provide structural plasticity leading to high adaptative capabilities to environmental stress conditions.

Writing a literature review as a class project in an upper-level undergraduate biochemistry course

A literature review is an important part of conducting academic research. Knowing how to conduct a literature search and write a high-quality literature review is a valuable skill. Herein, the authors describe the method of introducing a literature review writing exercise in an upper-level biochemistry course. Since 2020, authors have collaborated with numerous undergraduates writing literature reviews on topics in biochemistry that resulted in peer-reviewed publications. Authors believe that this unique idea of providing a course-based undergraduate research experience (CURE) to many undergraduates, especially those who otherwise do not receive collaborative research experience through traditional research paths, must be shared with other instructors.

Exploring the intrinsic dynamics of SARS-CoV-2, SARS-CoV and MERS-CoV spike glycoprotein through normal mode analysis using anisotropic network model

COVID-19 caused by SARS-CoV-2 have become a global pandemic with serious rate of fatalities. SARS-CoV and MERS-CoV have also caused serious outbreak previously but the intensity was much lower than the ongoing SARS-CoV-2. The main infectivity factor of all the three viruses is the spike glycoprotein. In this study we have examined the intrinsic dynamics of the prefusion, lying state of trimeric S protein of these viruses through Normal Mode Analysis using Anisotropic Network Model. The dynamic modes of the S proteins of the aforementioned viruses were compared by root mean square inner product (RMSIP), spectral overlap and cosine correlation matrix. S proteins show homogenous correlated or anticorrelated motions among their domains but direction of Cα atom among the spike proteins show less similarity. SARS-CoV-2 spike shows high vertically upward motion of the receptor binding motif implying its propensity for binding with the receptor even in the lying state. MERS-CoV spike shows unique dynamical motion compared to the other two S protein indicated by low RMSIP, spectral overlap and cosine correlation value. This study will guide in developing common potential inhibitor molecules against closed state of spike protein of these viruses to prevent conformational switching from lying to standing state.

Cloning and characterization of estrogen hydroxylase (cyp1a1 and cyp1b1) genes in the stinging catfish Heteropneustes fossilis and induction of mRNA expression during final oocyte maturation

Estrogen hydroxylases (EHs) are cytochrome P450 Family 1 (Cyp1, Clan 2) proteins involved in estrogen hydroxylations at 2-, 4- or 16- carbon positions to form catecholestrogens. EHs are encoded by CYP1A1, CYP1A2 and CYP1B1 in mammals. In the catfish Heteropneustes fossilis, cyp1a1 and cyp1b1 cDNAs were cloned and characterized from liver and ovary. The cyp1a1 cDNA is 2071 bp long and codes for a 518 amino acids (aa) long protein. The cloned cyp1b1 cDNA is 1927 bp long and codes for a 509 residue protein. The deduced proteins clustered distinctly into teleost Cyp1a1 and Cyp1b1 clades, distinct from the tetrapod clusters and featured common function domains and homology with other teleost proteins. In the qPCR assay, the transcripts were the most abundant in the liver, followed by brain and ovary, and moderate in gill, kidney and muscle. Evidence was presented to show the involvement of the genes in reproduction. Expression of brain and ovarian transcripts showed significant seasonal variations with the highest abundance in the spawning phase. In situ hybridization showed the transcripts in the follicular layer (theca and granulosa) of the ovarian follicles. Periovulatory changes in the expression cyp1a1 and cyp1b1 were obtained during final oocyte maturation (FOM) and ovulation induced by human chorionic gonadotropin (hCG), both in vivo and in vitro, and by 2-hydroxyestradiol-17β (catecholestrogen) in vitro. In the brain, the transcript levels increased with time but in the ovary, the increase was maximal at 16 h and decreased at 24 h. The periovulatory activation of the cyp1 genes was reported in this study and discussed on the basis of complex regulation of FOM and ovulation.

Some Surprising Implications of NMR-directed Simulations of Substrate Recognition and Binding by Cytochrome P450cam (CYP101A1)

Cytochrome P450_cam (CYP101A1) catalyzes the stereospecific 5-exo hydroxylation of d-camphor by molecular oxygen. Previously, residual dipolar couplings measured for backbone amide ¹H-¹⁵N correlations in both substrate-free and bound forms of CYP101A1 were used as restraints in soft annealing molecular dynamic simulations in order to identify average conformations of the enzyme with and without substrate bound. Multiple substrate-dependent conformational changes remote from the enzyme active site were identified, and site-directed mutagenesis and activity assays confirmed the importance of these changes in substrate recognition. The current work makes use of perturbation response scanning (PRS) and umbrella sampling molecular dynamic of the residual dipolar coupling-derived CYP101A1 structures to probe the roles of remote structural features in enforcing the regio- and stereospecific nature of the hydroxylation reaction catalyzed by CYP101A1. An improper dihedral angle Ψ was defined and used to maintain substrate orientation in the CYP101A1 active site, and it was observed that different values of Ψ result in different PRS response maps. Umbrella sampling methods show that the free energy of the system is sensitive to Ψ, and bound substrate forms an important mechanical link in the transmission of mechanical coupling through the enzyme structure. Finally, a qualitative approach to interpreting PRS maps in terms of the roles of secondary structural features is proposed.

Molecular elucidation of a new allelic variation at the Sg-5 gene associated with the absence of group A saponins in wild soybean

In soybean, triterpenoid saponin is one of the major secondary metabolites and is further classified into group A and DDMP saponins. Although they have known health benefits for humans and animals, acetylation of group A saponins causes bitterness and gives an astringent taste to soy products. Therefore, several studies are being conducted to eliminate acetylated group A saponins. Previous studies have isolated and characterized the Sg-5 (Glyma.15g243300) gene, which encodes the cytochrome P450 72A69 enzyme and is responsible for soyasapogenol A biosynthesis. In this study, we elucidated the molecular identity of a novel mutant of Glycine soja, 'CWS5095'. Phenotypic analysis using TLC and LC-PDA/MS/MS showed that the mutant 'CWS5095' did not produce any group A saponins. Segregation analysis showed that the absence of group A saponins is controlled by a single recessive allele. The locus was mapped on chromosome 15 (4.3 Mb) between Affx-89193969 and Affx-89134397 where the previously identified Glyma.15g243300 gene is positioned. Sequence analysis of the coding region for the Glyma.15g243300 gene revealed the presence of four SNPs in 'CWS5095' compared to the control lines. One of these four SNPs (G1127A) leads to the amino acid change Arg376Lys in the EXXR motif, which is invariably conserved among the CYP450 superfamily proteins. Co-segregation analysis showed that the missense mutation (Arg376Lys) was tightly linked with the absence of group A saponins in 'CWS5095'. Even though Arg and Lys have similar chemical features, the 3D modelled protein structure indicates that the replacement of Arg with Lys may cause a loss-of-function of the Sg-5 protein by inhibiting the stable binding of a heme cofactor to the CYP72A69 apoenzyme.

Innovations in Undergraduate Chemical Biology Education

Chemical biology derives intellectual vitality from its scientific interface: applying chemical strategies and perspectives to biological questions. There is a growing need for chemical biologists to synergistically integrate their research programs with their educational activities to become holistic teacher-scholars. This review examines how course-based undergraduate research experiences (CUREs) are an innovative method to achieve this integration. Because CUREs are course-based, the review first offers strategies for creating a student-centered learning environment, which can improve students' outcomes. Exemplars of CUREs in chemical biology are then presented and organized to illustrate the five defining characteristics of CUREs: significance, scientific practices, discovery, collaboration, and iteration. Finally, strategies to overcome common barriers in CUREs are considered as well as future innovations in chemical biology education.

Incorporating modeling and simulations in undergraduate biophysical chemistry course to promote understanding of structure-dynamics-function relationships in proteins

A project-based biophysical chemistry laboratory course, which is offered to the biochemistry and molecular biology majors in their senior year, is described. In this course, the classroom study of the structure-function of biomolecules is integrated with the discovery-guided laboratory study of these molecules using computer modeling and simulations. In particular, modern computational tools are employed to elucidate the relationship between structure, dynamics, and function in proteins. Computer-based laboratory protocols that we introduced in three modules allow students to visualize the secondary, super-secondary, and tertiary structures of proteins, analyze non-covalent interactions in protein-ligand complexes, develop three-dimensional structural models (homology model) for new protein sequences and evaluate their structural qualities, and study proteins' intrinsic dynamics to understand their functions. In the fourth module, students are assigned to an authentic research problem, where they apply their laboratory skills (acquired in modules 1-3) to answer conceptual biophysical questions. Through this process, students gain in-depth understanding of protein dynamics-the missing link between structure and function. Additionally, the requirement of term papers sharpens students' writing and communication skills. Finally, these projects result in new findings that are communicated in peer-reviewed journals.