Structural and dynamic properties of cytochrome P450 BM-3 in pure water and in a dimethylsulfoxide/water mixture

Structure-based mechanisms: On the way to apply alcohol dehydrogenases/reductases to organic-aqueous systems

Alcohol dehydrogenases/reductases catalyze enantioselective syntheses of versatile chiral compounds relying on direct hydride transfer from cofactor to substrates, or to an intermediate and then to substrates. Since most of the substrates catalyzed by alcohol dehydrogenases/reductases are insoluble in aqueous solutions, increasing interest has been turning to organic-aqueous systems. However, alcohol dehydrogenases/reductases are normally instable in organic solvents, leading to the unsatisfied enantioselective synthesis efficiency. The behaviors of these enzymes in organic solvents at an atomic level are unclear, thus it is of great importance to understand its structure-based mechanisms in organic-aqueous systems to improve their relative stability. Here, we summarized the accessible structures of alcohol dehydrogenases/reductases in Protein Data Bank crystallized in organic-aqueous systems, and compared the structures of alcohol dehydrogenases/reductases which have different tolerance towards organic solvents. By understanding the catalytic behaviors and mechanisms of these enzymes in organic-aqueous systems, the efficient enantioselective syntheses mediated by alcohol dehydrogenases/reductases and further challenges are also discussed through solvent engineering and enzyme-immobilization in the last decade.

Organic solvent stability and long-term storage of myoglobin-based carbene transfer biocatalysts

Recent years have witnessed a rapid increase in the application of enzymes for chemical synthesis and manufacturing, including the industrial-scale synthesis of pharmaceuticals using multienzyme processes. From an operational standpoint, these bioprocesses often require robust biocatalysts capable of tolerating high concentrations of organic solvents and possessing long shelflife stability. In this work, we investigated the activity and stability of myoglobin (Mb)-based carbene transfer biocatalysts in the presence of organic solvents and after lyophilization. Our studies demonstrate that Mb-based cyclopropanases possess remarkable organic solvent stability, maintaining high levels of activity and stereoselectivity in the presence of up to 30%-50% (v/v) concentrations of various organic solvents, including ethanol, methanol, N,N-dimethylformamide, acetonitrile, and dimethyl sulfoxide. Furthermore, they tolerate long-term storage in lyophilized form, both as purified protein and as whole cells, without significant loss in activity and stereoselectivity. These stability properties are shared by Mb-based carbene transferases optimized for other type of asymmetric carbene transfer reactions. Finally, we report on simple protocols for catalyst recycling as whole-cell system and for obviating the need for strictly anaerobic conditions to perform these transformations. These findings demonstrate the robustness of Mb-based carbene transferases under operationally relevant conditions and should help guide the application of these biocatalysts for synthetic applications.

In Silico Studies of Small Molecule Interactions with Enzymes Reveal Aspects of Catalytic Function

Small molecules, such as solvent, substrate, and cofactor molecules, are key players in enzyme catalysis. Computational methods are powerful tools for exploring the dynamics and thermodynamics of these small molecules as they participate in or contribute to enzymatic processes. In-depth knowledge of how small molecule interactions and dynamics influence protein conformational dynamics and function is critical for progress in the field of enzyme catalysis. Although numerous computational studies have focused on enzyme-substrate complexes to gain insight into catalytic mechanisms, transition states and reaction rates, the dynamics of solvents, substrates, and cofactors are generally less well studied. Also, solvent dynamics within the biomolecular solvation layer play an important part in enzyme catalysis, but a full understanding of its role is hampered by its complexity. Moreover, passive substrate transport has been identified in certain enzymes, and the underlying principles of molecular recognition are an area of active investigation. Enzymes are highly dynamic entities that undergo different conformational changes, which range from side chain rearrangement of a residue to larger-scale conformational dynamics involving domains. These events may happen nearby or far away from the catalytic site, and may occur on different time scales, yet many are related to biological and catalytic function. Computational studies, primarily molecular dynamics (MD) simulations, provide atomistic-level insight and site-specific information on small molecule interactions, and their role in conformational pre-reorganization and dynamics in enzyme catalysis. The review is focused on MD simulation studies of small molecule interactions and dynamics to characterize and comprehend protein dynamics and function in catalyzed reactions. Experimental and theoretical methods available to complement and expand insight from MD simulations are discussed briefly.

Evaluation of the effects of hydrophilic organic solvents on CYP3A-mediated drug-drug interaction in vitro

This study evaluated the effects of the commonly used hydrophilic organic solvents, acetonitrile, methanol, ethanol, 1-propanol, dimethyl sulfoxide (DMSO), N,N-dimethylformamide, polyethylene glycol and propylene glycol, on CYP3A in pooled human liver microsomes, using testosterone and midazolam as substrates. Furthermore, we examined the modulation effect of organic solvents on CYP3A inhibition by ketoconazole. Testosterone 6b-hydroxylation activity was potently inhibited in the presence of DMSO and 1-propanol in a concentration-dependent manner. Midazolam 1'-hydroxylation activity, however, was weakly inhibited only by 1% of DMSO, the highest concentration used in this study. Moreover, the potency of ketoconazole to inhibit CYP3A activities was variable, depending on the organic solvent used as a dissolving solvent for ketoconazole. Our data indicate that each organic solvent had an effect on CYP3A4 activity, evaluated by both substrates with different magnitudes. Furthermore, it was shown that the effects of organic solvents on CYP3A activity are substrate-dependent. The present study also shows that methanol had little effect on either substrate.

Insight into the redox partner interaction mechanism in cytochrome P450BM-3 using molecular dynamics simulations

Flavocytochrome P450BM-3 is a soluble bacterial reductase composed of two flavin (FAD/FMN) and one HEME domains. In this article, we have performed molecular dynamics simulations on both the isolated FMN and HEME domains and their crystallographic complex, with the aim to study their binding modes and to garner insight into the interdomain electron transfer (ET) mechanism. The results evidenced an interdomain conformational rearrangement that reduces the average distance between the FMN and HEME cofactors from 1.81 nm, in the crystal structure, to an average value of 1.41±0.09 nm along the simulation. This modification is in agreement with previously proposed hypotheses suggesting that the crystallographic FMN/HEME complex is not in the optimal arrangement for favorable ET rate under physiological conditions. The calculation of the transfer rate along the simulation, using the Pathways Path method, demonstrated the occurrence of seven ET pathways between the two redox centers, with three of them providing ET rates (KET ) comparable with the experimental one. The sampled ET pathways comprise the amino acids N319, L322, F390, K391, P392, F393, A399, C400, and Q403 of the HEME domain and M490 of the FMN domain. The values of KET closer to the experiment were found along the pathways FMN(C7)→F390→K391→P392→HEME(Fe) and FMN(C8)→M490→F393→HEME(Fe). Finally, the analysis of the collective modes of the protein complex evidences a clear correlation of the first two essential modes with the activation of the most effective ET pathways along the trajectory.

Structure, dynamics, and function of the monooxygenase P450 BM-3: insights from computer simulations studies

The monooxygenase P450 BM-3 is a NADPH-dependent fatty acid hydroxylase enzyme isolated from soil bacterium Bacillus megaterium. As a pivotal member of cytochrome P450 superfamily, it has been intensely studied for the comprehension of structure-dynamics-function relationships in this class of enzymes. In addition, due to its peculiar properties, it is also a promising enzyme for biochemical and biomedical applications. However, despite the efforts, the full understanding of the enzyme structure and dynamics is not yet achieved. Computational studies, particularly molecular dynamics (MD) simulations, have importantly contributed to this endeavor by providing new insights at an atomic level regarding the correlations between structure, dynamics, and function of the protein. This topical review summarizes computational studies based on MD simulations of the cytochrome P450 BM-3 and gives an outlook on future directions.

Computational prediction of metabolism: sites, products, SAR, P450 enzyme dynamics, and mechanisms

Metabolism of xenobiotics remains a central challenge for the discovery and development of drugs, cosmetics, nutritional supplements, and agrochemicals. Metabolic transformations are frequently related to the incidence of toxic effects that may result from the emergence of reactive species, the systemic accumulation of metabolites, or by induction of metabolic pathways. Experimental investigation of the metabolism of small organic molecules is particularly resource demanding; hence, computational methods are of considerable interest to complement experimental approaches. This review provides a broad overview of structure- and ligand-based computational methods for the prediction of xenobiotic metabolism. Current computational approaches to address xenobiotic metabolism are discussed from three major perspectives: (i) prediction of sites of metabolism (SOMs), (ii) elucidation of potential metabolites and their chemical structures, and (iii) prediction of direct and indirect effects of xenobiotics on metabolizing enzymes, where the focus is on the cytochrome P450 (CYP) superfamily of enzymes, the cardinal xenobiotics metabolizing enzymes. For each of these domains, a variety of approaches and their applications are systematically reviewed, including expert systems, data mining approaches, quantitative structure–activity relationships (QSARs), and machine learning-based methods, pharmacophore-based algorithms, shape-focused techniques, molecular interaction fields (MIFs), reactivity-focused techniques, protein–ligand docking, molecular dynamics (MD) simulations, and combinations of methods. Predictive metabolism is a developing area, and there is still enormous potential for improvement. However, it is clear that the combination of rapidly increasing amounts of available ligand- and structure-related experimental data (in particular, quantitative data) with novel and diverse simulation and modeling approaches is accelerating the development of effective tools for prediction of in vivo metabolism, which is reflected by the diverse and comprehensive data sources and methods for metabolism prediction reviewed here. This review attempts to survey the range and scope of computational methods applied to metabolism prediction and also to compare and contrast their applicability and performance.

Active site substitution A82W improves the regioselectivity of steroid hydroxylation by cytochrome P450 BM3 mutants as rationalized by spin relaxation nuclear magnetic resonance studies

Cytochrome P450 BM3 from Bacillus megaterium is a monooxygenase with great potential for biotechnological applications. In this paper, we present engineered drug-metabolizing P450 BM3 mutants as a novel tool for regioselective hydroxylation of steroids at position 16β. In particular, we show that by replacing alanine at position 82 with a tryptophan in P450 BM3 mutants M01 and M11, the selectivity toward 16β-hydroxylation for both testosterone and norethisterone was strongly increased. The A82W mutation led to a ≤42-fold increase in V(max) for 16β-hydroxylation of these steroids. Moreover, this mutation improves the coupling efficiency of the enzyme, which might be explained by a more efficient exclusion of water from the active site. The substrate affinity for testosterone increased at least 9-fold in M11 with tryptophan at position 82. A change in the orientation of testosterone in the M11 A82W mutant as compared to the orientation in M11 was observed by T(1) paramagnetic relaxation nuclear magnetic resonance. Testosterone is oriented in M11 with both the A- and D-ring protons closest to the heme iron. Substituting alanine at position 82 with tryptophan results in increased A-ring proton-iron distances, consistent with the relative decrease in the level of A-ring hydroxylation at position 2β.

P450(BM3) (CYP102A1): connecting the dots

P450(BM3) (CYP102A1), a fatty acid hydroxylase from Bacillus megaterium, has been extensively studied over a period of almost forty years. The enzyme has been redesigned to catalyse the oxidation of non-natural substrates as diverse as pharmaceuticals, terpenes and gaseous alkanes using a variety of engineering strategies. Crystal structures have provided a basis for several of the catalytic effects brought about by mutagenesis, while changes to reduction potentials, inter-domain electron transfer rates and catalytic parameters have yielded functional insights. Areas of active research interest include drug metabolite production, the development of process-scale techniques, unravelling general mechanistic aspects of P450 chemistry, methane oxidation, and improving selectivity control to allow the synthesis of fine chemicals. This review draws together the disparate research themes and places them in a historical context with the aim of creating a resource that can be used as a gateway to the field.

Methionine sulfoxides on prion protein Helix-3 switch on the alpha-fold destabilization required for conversion

Background

The conversion of the cellular prion protein (PrP^C) into the infectious form (PrP^Sc) is the key event in prion induced neurodegenerations. This process is believed to involve a multi-step conformational transition from an α-helical (PrP^C) form to a β-sheet-rich (PrP^Sc) state. In addition to the conformational difference, PrP^Sc exhibits as covalent signature the sulfoxidation of M213. To investigate whether such modification may play a role in the misfolding process we have studied the impact of methionine oxidation on the dynamics and energetics of the HuPrP(125–229) α-fold.

Methodology/Principal Findings

Using molecular dynamics simulation, essential dynamics, correlated motions and signal propagation analysis, we have found that substitution of the sulfur atom of M213 by a sulfoxide group impacts on the stability of the native state increasing the flexibility of regions preceding the site of the modification and perturbing the network of stabilizing interactions. Together, these changes favor the population of alternative states which maybe essential in the productive pathway of the pathogenic conversion. These changes are also observed when the sulfoxidation is placed at M206 and at both, M206 and M213.

Conclusions/Significance

Our results suggest that the sulfoxidation of Helix-3 methionines might be the switch for triggering the initial α-fold destabilization required for the productive pathogenic conversion.

Activity of human P450 2D6 in biphasic solvent systems

Several limitations have restricted the use of P450 enzymes in synthesis, including the narrow substrate specificity of some P450 isoforms, the need for a redox partner and an expensive cofactor, incompatibility with organic solvents, and poor stability. We previously demonstrated that the natural redox partner and cofactor of the promiscuous P450s 3A4 and 2D6 can be efficiently substituted with some cheap hydrogen peroxide donors or organic peroxides. We report here that P450 2D6 maintains as much as 76% of its activity when used in buffer/organic emulsions. Product formation in biphasic solvent systems is comparable whether the natural redox partner and cofactor are used, or a surrogate. As reported for other enzymes, a correlation is observed between the logP and the suitability of a solvent for enzymatic activity. Moreover, the utility of our system was established by demonstrating the transformation of a novel hydrophobic substrate, not modified by P450 2D6 in the absence of organic solvent.

In Vitro Evolution of a Fungal Laccase in High Concentrations of Organic Cosolvents

Fungal laccases are remarkable green catalysts that have a broad substrate specificity and many potential applications in bioremediation, lignocellulose processing, organic synthesis, and more. However, most of these transformations must be carried out at high concentrations of organic cosolvents in which laccases undergo unfolding, thereby losing their activity. We have tailored a thermostable laccase that tolerates high concentrations of cosolvents, the genetic product of five rounds of directed evolution expressed in Saccharomyces cerevisiae. This evolved laccase--R2 variant--was capable of resisting a wide array of cosolvents at concentrations as high as 50% (v/v). Intrinsic laccase features such as the redox potential and the geometry of catalytic copper varied slightly during the course of the molecular evolution. Some mutations at the protein surface stabilized the laccase by allowing additional electrostatic and hydrogen bonding to occur.

Toward understanding the inactivation mechanism of monooxygenase P450 BM-3 by organic cosolvents: A molecular dynamics simulation study

Cytochrome P450 BM-3 from Bacillus megaterium is an extensively studied enzyme for industrial applications. A major focus of current protein engineering research is directed to improving the catalytic performance of P450 BM-3 toward nonnatural substrates of industrial importance in the presence of organic solvents or cosolvents. For the latter reason, it is important to study the effect of organic cosolvent molecules on the structure and dynamics of the enzyme, in particular, the effect of cosolvent molecules on the active site's structure and dynamics. In this paper, we have studied, using molecular dynamics (MD) simulations, the F87A mutant of P450 BM-3 in the presence of DMSO as cosolvent, to understand the role of the F87A substitution for its catalytic activity. This mutant exhibits an altered regioselectivity and substrate specificity compared with wild-type; however, it has lower tolerance toward DMSO. The simulation results offer an explanation for the DMSO sensitivity of the F87A mutant. Our simulation results show that the F87 side chain prevents the disturbance of the water molecule bound to the heme iron by DMSO molecules. The absence of the phenyl ring in F87A mutant promotes interactions of the DMSO molecule with the heme iron resulting in water displacement by DMSO at the catalytic heme center.