Evaluation of structural features in fungal cytochromes P450 predicted to rule catalytic diversification

Fungi belong to the large kingdom of lower eukaryotic organisms encompassing yeasts along with filamentous and dimorphic members. Microbial P450 enzymes have contributed to exploration of and adaptation to diverse ecological niches such as conversion of lipophilic compounds to more hydrophilic derivatives or degradation of a vast array of environmental toxicants. To better understand diversification of the catalytic behavior of fungal P450s, detailed insight into the molecular machinery steering oxidative attack on the distinctly structured endogenous and xenobiotic substrates is of preeminent interest. Based on a general, CYP102A1-related template the bulk of predicted substrate/inhibitor-binding determinants were shown to cluster near the distal heme face within the six known substrate recognition sites (SRSs) made up by the α-helical B′/F/G/I tetrad, the B′–C interhelical loop and strands of the β6-sheet, population density being highest in the structurally flexible SRS-1 and SRS-4 domains, showing a low degree of conservation. Reactivity toward ligands favorably coincides with the lipophilicity/hydrophilicity profile and bulkiness of critical amino acids acting as selective filters. Some decisive elements may also serve in maintenance of catalytic competence via their action as gatekeepers directing substrate access/positioning or stabilizers of the heme environment enabling dioxygen activation. Non-SRS residues seem to control spin state equilibria and attract redox partners by electrostatic forces. Of note, the inhibitory potency of azole-type fungicides is likely to arise from perturbation of the complex interplay of the mechanistic principles addressed above. Knowledge-supported exploitation of the topological data will be helpful in the manufacture of commodity/specialty chemicals as well as therapeutic agents. Also, engineered fungal P450s may be used to improve pollutant-specific bioremediation of contaminated soils.     

Molecular Mechanics Calculations of Proteins. Comparison of Different Energy Minimization Strategies

Abstract

A general strategy for performing energy minimization of proteins using the SYBYL molecular modelling program has been developed. The influence of several variables including energy minimization procedure, solvation, dielectric function and dielectric constant have been investigated in order to develop a general method, which is capable of producing high quality protein structures. Avian pancreatic polypeptide (APP) and bovine pancreatic phospholipase A2 (BP PLA2) were selected for the calculations, because high quality X-ray structures exist and because all classes of secondary structure are represented in the structures. The energy minimized structures were evaluated relative to the corresponding X-ray structures. The overall similarity was checked by calculating RMS distances for all atom positions. Backbone conformation was checked by Ramachandran plots and secondary structure elements evaluated by the length on hydrogen bonds. The dimensions of active site in BP PLA2 is very dependent on electrostatic interactions, due to the presence of the positively charged calcium ion. Thus, the distances between calcium and the calcium-coordinating groups were used as a quality index for this protein. Energy minimized structures of the trimeric PLA2 from Indian cobra (N.n.n. PLA2) were used for assessing the impact of protein-protein interactions. Based on the above mentioned criteria, it could be concluded that using the following conditions: Dielectric constant ? = 4 or 20; a distance dependent dielectric function and stepwise energy minimization, it is possible to reproduce X-ray structures very accurately without including explicit solvent molecules.