Evolution of P450 Monooxygenases toward Formation of Transient Channels and Exclusion of Nonproductive Gases

Unusual catalytic strategy by non-heme Fe(ii)/2-oxoglutarate-dependent aspartyl hydroxylase AspH

Biocatalytic C-H oxidation reactions are of important synthetic utility, provide a sustainable route for selective synthesis of important organic molecules, and are an integral part of fundamental cell processes. The multidomain non-heme Fe(ii)/2-oxoglutarate (2OG) dependent oxygenase AspH catalyzes stereoselective (3R)-hydroxylation of aspartyl- and asparaginyl-residues. Unusually, compared to other 2OG hydroxylases, crystallography has shown that AspH lacks the carboxylate residue of the characteristic two-His-one-Asp/Glu Fe-binding triad. Instead, AspH has a water molecule that coordinates Fe(ii) in the coordination position usually occupied by the Asp/Glu carboxylate. Molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) studies reveal that the iron coordinating water is stabilized by hydrogen bonding with a second coordination sphere (SCS) carboxylate residue Asp721, an arrangement that helps maintain the six coordinated Fe(ii) distorted octahedral coordination geometry and enable catalysis. AspH catalysis follows a dioxygen activation-hydrogen atom transfer (HAT)-rebound hydroxylation mechanism, unusually exhibiting higher activation energy for rebound hydroxylation than for HAT, indicating that the rebound step may be rate-limiting. The HAT step, along with substrate positioning modulated by the non-covalent interactions with SCS residues (Arg688, Arg686, Lys666, Asp721, and Gln664), are essential in determining stereoselectivity, which likely proceeds with retention of configuration. The tetratricopeptide repeat (TPR) domain of AspH influences substrate binding and manifests dynamic motions during catalysis, an observation of interest with respect to other 2OG oxygenases with TPR domains. The results provide unique insights into how non-heme Fe(ii) oxygenases can effectively catalyze stereoselective hydroxylation using only two enzyme-derived Fe-ligating residues, potentially guiding enzyme engineering for stereoselective biocatalysis, thus advancing the development of non-heme Fe(ii) based biomimetic C-H oxidation catalysts, and supporting the proposal that the 2OG oxygenase superfamily may be larger than once perceived.

X-ray Crystal Structures of Methane Monooxygenase Hydroxylase Complexes with Variants of Its Regulatory Component: Correlations with Altered Reaction Cycle Dynamics

Full activity of soluble methane monooxygenase (sMMO) depends upon the formation of a 1:1 complex of the regulatory protein MMOB with each alpha subunit of the (αβγ)₂ hydroxylase, sMMOH. Previous studies have shown that mutations in the core region of MMOB and in the N- and C-termini cause dramatic changes in the rate constants for steps in the sMMOH reaction cycle. Here, X-ray crystal structures are reported for the sMMOH complex with two double variants within the core region of MMOB, DBL1 (N107G/S110A), and DBL2 (S109A/T111A), as well as two variants in the MMOB N-terminal region, H33A and H5A. DBL1 causes a 150-fold decrease in the formation rate constant of the reaction cycle intermediate P, whereas DBL2 accelerates the reaction of the dinuclear Fe(IV) intermediate Q with substrates larger than methane by three- to fourfold. H33A also greatly slows P formation, while H5A modestly slows both formation of Q and its reactions with substrates. Complexation with DBL1 or H33A alters the position of sMMOH residue R245, which is part of a conserved hydrogen-bonding network encompassing the active site diiron cluster where P is formed. Accordingly, electron paramagnetic resonance spectra of sMMOH:DBL1 and sMMOH:H33A complexes differ markedly from that of sMMOH:MMOB, showing an altered electronic environment. In the sMMOH:DBL2 complex, the position of M247 in sMMOH is altered such that it enlarges a molecular tunnel associated with substrate entry into the active site. The H5A variant causes only subtle structural changes despite its kinetic effects, emphasizing the precise alignment of sMMOH and MMOB required for efficient catalysis.

Tunnel Architectures in Enzyme Systems that Transport Gaseous Substrates

Molecular tunnels regulate delivery of substrates/intermediates in enzymes which either harbor deep-seated reaction centers or are for transport of reactive/toxic intermediates that need to be specifically delivered. Here, we focus on the importance of structural diversity in tunnel architectures, especially for the gaseous substrate translocation, in rendering differential substrate preferences and directionality. Two major types of tunnels have been discussed, one that transports stable gases from the environment to the active site, namely, external gaseous (EG) tunnels, and the other that transports molecules between active sites, namely, internal gaseous (IG) tunnels. Aspects as to how the gaseous tunnels have shaped during the course of evolution and their potential to modulate the substrate flow and enzymatic function are examined. In conclusion, the review highlights our perspective on the pulsation mechanism that could facilitate unidirectional translocation of the gaseous molecules through buried tunnels.

Small-Molecule Tunnels in Metalloenzymes Viewed as Extensions of the Active Site

Rigorous substrate selectivity is a hallmark of enzyme catalysis. This selectivity is generally ascribed to a thermodynamically favorable process of substrate binding to the enzyme active site based upon complementary physiochemical characteristics, which allows both acquisition and orientation. However, this chemical selectivity is more difficult to rationalize for diminutive molecules that possess too narrow a range of physical characteristics to allow either precise positioning or discrimination between a substrate and an inhibitor. Foremost among these small molecules are dissolved gases such as H₂, N₂, O₂, CO, CO₂, NO, N₂O, NH₃, and CH₄ so often encountered in metalloenzyme catalysis. Nevertheless, metalloenzymes have evolved to metabolize these small-molecule substrates with high selectivity and efficiency.The soluble methane monooxygenase enzyme (sMMO) acts upon two of these small molecules, O₂ and CH₄, to generate methanol as part of the C1 metabolic pathway of methanotrophic organisms. sMMO is capable of oxidizing many alternative hydrocarbon substrates. Remarkably, however, it will preferentially oxidize methane, the substrate with the fewest discriminating physical characteristics and the strongest C-H bond. Early studies led us to broadly attribute this specificity to the formation of a "molecular sieve" in which a methane- and oxygen-sized tunnel provides a size-selective route from bulk solvent to the completely buried sMMO active site. Indeed, recent cryogenic and serial femtosecond ambient temperature crystallographic studies have revealed such a route in sMMO. A detailed study of the sMMO tunnel considered here in the context of small-molecule tunnels identified in other metalloenzymes reveals three discrete characteristics that contribute to substrate selectivity and positioning beyond that which can be provided by the active site itself. Moreover, the dynamic nature of many tunnels allows an exquisite coordination of substrate binding and reaction phases of the catalytic cycle. Here we differentiate between the highly selective molecular tunnel, which allows only the one-dimensional transit of small molecules, and the larger, less-selective channels found in typical enzymes. Methods are described to identify and characterize tunnels as well as to differentiate them from channels. In metalloenzymes which metabolize dissolved gases, we posit that the contribution of tunnels is so great that they should be considered to be extensions of the active site itself. A full understanding of catalysis by these enzymes requires an appreciation of the roles played by tunnels. Such an understanding will also facilitate the use of the enzymes or their synthetic mimics in industrial or pharmaceutical applications.

Multi-scale simulation reveals that an amino acid substitution increases photosensitizing reaction inputs in Rhodopsins

Evaluating the availability of molecular oxygen (O₂ ) and energy of excited states in the retinal binding site of rhodopsin is a crucial challenging first step to understand photosensitizing reactions in wild-type (WT) and mutant rhodopsins by absorbing visible light. In the present work, energies of the ground and excited states related to 11-cis-retinal and the O₂ accessibility to the β-ionone ring are evaluated inside WT and human M207R mutant rhodopsins. Putative O₂ pathways within rhodopsins are identified by using molecular dynamics simulations, Voronoi-diagram analysis, and implicit ligand sampling while retinal energetic properties are investigated through density functional theory, and quantum mechanical/molecular mechanical methods. Here, the predictions reveal that an amino acid substitution can lead to enough energy and O₂ accessibility in the core hosting retinal of mutant rhodopsins to favor the photosensitized singlet oxygen generation, which can be useful in understanding retinal degeneration mechanisms and in designing blue-lighting-absorbing proteic photosensitizers.

Role of Water Hydrogen Bonding on Transport of Small Molecules inside Hydrophobic Channels

We conducted a systematic analysis of water networking inside smooth hyperboloid hydrophobic structures (cylindrical, barrel, and hourglass shapes) to elucidate the role of water hydrogen bonding on the transport of small hydrophobic molecules (ligands). Through a series of molecular dynamics simulations, we established that a hydrogen-bonded network forming along the centerline resulted in a water exclusion zone adjacent to the walls. The size of the exclusion zone is a function of the geometry and the nonbonded interaction strength, defining the effective hydrophobicity of the structure. Exclusion of water molecules from this zone results in lower apparent viscosity, leading to acceleration of ligand transport up to 7 times faster than that measured in the bulk. Transport of ligands into and out of the hydrophobic structures was shown to be controlled by a single water molecule that capped the narrow regions in the structure. This mechanism provides physical insights into the behavior and role of water in the bottleneck regions of hydrophobic enzyme channels. These findings were then used in a sister publication [ Escalante , D. E. , Comput. Struct. Biotechnol. J. 2019 17 757 760 ] to develop a model that can accurately predict the transport of ligands along nanochannels of broad-substrate specificity enzymes.

Engineering of lysine cyclodeaminase conformational dynamics for relieving substrate and product inhibitions in the biosynthesis of l-pipecolic acid

Efficient biocatalytic process construction by relieving substrate and product inhibitions via identification and engineering of enzyme conformational dynamics.

How phosphorylation influences E1 subunit pyruvate dehydrogenase: A computational study

Pyruvate (PYR) dehydrogenase complex (PDC) is an enzymatic system that plays a crucial role in cellular metabolism as it controls the entry of carbon into the Krebs cycle. From a structural point of view, PDC is formed by three different subunits (E1, E2 and E3) capable of catalyzing the three reaction steps necessary for the full conversion of pyruvate to acetyl-CoA. Recent investigations pointed out the crucial role of this enzyme in the replication and survival of specific cancer cell lines, renewing the interest of the scientific community. Here, we report the results of our molecular dynamics studies on the mechanism by which posttranslational modifications, in particular the phosphorylation of three serine residues (Ser-264-α, Ser-271-α, and Ser-203-α), influence the enzymatic function of the protein. Our results support the hypothesis that the phosphorylation of Ser-264-α and Ser-271-α leads to (1) a perturbation of the catalytic site structure and dynamics and, especially in the case of Ser-264-α, to (2) a reduction in the affinity of E1 for the substrate. Additionally, an analysis of the channels connecting the external environment with the catalytic site indicates that the inhibitory effect should not be due to the occlusion of the access/egress pathways to/from the active site.

ChannelsDB: database of biomacromolecular tunnels and pores

ChannelsDB (http://ncbr.muni.cz/ChannelsDB) is a database providing information about the positions, geometry and physicochemical properties of channels (pores and tunnels) found within biomacromolecular structures deposited in the Protein Data Bank. Channels were deposited from two sources; from literature using manual deposition and from a software tool automatically detecting tunnels leading to the enzymatic active sites and selected cofactors, and transmembrane pores. The database stores information about geometrical features (e.g. length and radius profile along a channel) and physicochemical properties involving polarity, hydrophobicity, hydropathy, charge and mutability. The stored data are interlinked with available UniProt annotation data mapping known mutation effects to channel-lining residues. All structures with channels are displayed in a clear interactive manner, further facilitating data manipulation and interpretation. As such, ChannelsDB provides an invaluable resource for research related to deciphering the biological function of biomacromolecular channels.

Ligand diffusion in proteins via enhanced sampling in molecular dynamics

Computational simulations in biophysics describe the dynamics and functions of biological macromolecules at the atomic level. Among motions particularly important for life are the transport processes in heterogeneous media. The process of ligand diffusion inside proteins is an example of a complex rare event that can be modeled using molecular dynamics simulations. The study of physical interactions between a ligand and its biological target is of paramount importance for the design of novel drugs and enzymes. Unfortunately, the process of ligand diffusion is difficult to study experimentally. The need for identifying the ligand egress pathways and understanding how ligands migrate through protein tunnels has spurred the development of several methodological approaches to this problem. The complex topology of protein channels and the transient nature of the ligand passage pose difficulties in the modeling of the ligand entry/escape pathways by canonical molecular dynamics simulations. In this review, we report a methodology involving a reconstruction of the ligand diffusion reaction coordinates and the free-energy profiles along these reaction coordinates using enhanced sampling of conformational space. We illustrate the above methods on several ligand-protein systems, including cytochromes and G-protein-coupled receptors. The methods are general and may be adopted to other transport processes in living matter.