Bacterial cytochrome P450s: a bioinformatics odyssey of substrate discovery.

Bacterial P450 cytochromes (BacCYPs) are versatile heme-containing proteins responsible for oxidation reactions on a wide range of substrates, contributing to the production of valuable natural products with limitless biotechnological potential. While the sequencing of microbial genomes has provided a wealth of BacCYP sequences, functional characterization lags behind, hindering our understanding of their roles. This study employs a comprehensive approach to predict BacCYP substrate specificity, bridging the gap between sequence and function. We employed an integrated approach combining sequence and functional data analysis, genomic context exploration, 3D structural modeling with molecular docking, and phylogenetic clustering. The research begins with an in-depth analysis of BacCYP sequence diversity and structural characteristics, revealing conserved motifs and recurrent residues in the active site. Phylogenetic analysis identifies distinct groups within the BacCYP family based on sequence similarity. However, our study reveals that sequence alone does not consistently predict substrate specificity, necessitating additional perspectives. The study delves into the genetic context of BacCYPs, utilizing neighboring gene information to infer potential substrates, a method proven very effective in many cases. Molecular docking is employed to assess BacCYP-substrate interactions, confirming potential substrates and providing insights into selectivity. Finally, a comprehensive strategy is proposed for predicting BacCYP substrates, involving all the evaluated approaches. The effectiveness of this strategy is demonstrated with two case studies, highlighting its potential for substrate discovery.

Efficient Calculation of Enzyme Reaction Free Energy Profiles Using a Hybrid Differential Relaxation Algorithm: Application to Mycobacterial Zinc Hydrolases.

Determination of the free energy profile for an enzyme reaction mechanism is of primordial relevance, paving the way for our understanding of the enzyme's catalytic power at the molecular level. Although hybrid, mostly DFT-based, QM/MM methods have been extensively applied to this type of studies, achieving accurate and statistically converged results at a moderate computational cost is still an open challenge. Recently, we have shown that accurate results can be achieved in less computational time, combining Jarzynski's relationship with a hybrid differential relaxation algorithm (HyDRA), which allows partial relaxation of the solvent during the nonequilibrium steering of the reaction. In this work, we have applied this strategy to study two mycobacterial zinc hydrolases. Mycobacterium tuberculosis infections are still a worldwide problem and thus characterization and validation of new drug targets is an intense field of research. Among possible drug targets, recently two essential zinc hydrolases, MshB (Rv1170) and MA-amidase (Rv3717), have been proposed and structurally characterized. Although possible mechanisms have been proposed by analogy to the widely studied human Zn hydrolases, several key issues, particularly those related to Zn coordination sphere and its role in catalysis, remained unanswered. Our results show that mycobacterial Zn hydrolases share a basic two-step mechanism. First, the attacking water becomes deprotonated by the conserved base and establishes the new C-O bond leading to a tetrahedral intermediate. The intermediate requires moderate reorganization to allow for proton transfer to the amide N and C-N bond breaking to occur in the second step. Zn ion plays a key role in stabilizing the tetrahedral intermediate and balancing the negative charge of the substrate during hydroxide ion attack. Finally, comparative analysis of other Zn hydrolases points to a convergent mechanistic evolution.

QM/MM study of the C-C coupling reaction mechanism of CYP121, an essential cytochrome p450 of Mycobacterium tuberculosis.

Among 20 p450s of Mycobacterium tuberculosis (Mt), CYP121 has received an outstanding interest, not only due to its essentiality for bacterial viability but also because it catalyzes an unusual carbon-carbon coupling reaction. Based on the structure of the substrate bound enzyme, several reaction mechanisms were proposed involving first Tyr radical formation, second Tyr radical formation, and C-C coupling. Key and unknown features, being the nature of the species that generate the first and second radicals, and the role played by the protein scaffold each step. In the present work we have used classical and quantum based computer simulation methods to study in detail its reaction mechanism. Our results show that substrate binding promotes formation of the initial oxy complex, Compound I is the responsible for first Tyr radical formation, and that the second Tyr radical is formed subsequently, through a PCET reaction, promoted by the presence of key residue Arg386. The final C-C coupling reaction possibly occurs in bulk solution, thus yielding the product in one oxygen reduction cycle. Our results thus contribute to a better comprehension of MtCYP121 reaction mechanism, with direct implications for inhibitor design, and also contribute to our general understanding of these type of enzymes.