Search

We report the results of computational studies of the guanosine triphosphate (GTP) hydrolysis in the active site of the KRas-NF1 protein complex, where KRas stands for the K-isoform of the Ras (ras sarcoma) protein and NF1 (neurofbromin-1) is the activating protein. The model system was constructed using coordinates of heavy atoms from the crystal structure PDB ID 6OB2 with the GTP analog GMPPNP. Large-scale classical molecular dynamics (MD) calculations were performed to analyze conformations of the enzyme-substrate complexes. The Gibbs energy profiles for the hydrolysis reaction were computed using MD simulations with quantum mechanics/molecular mechanics (QM/MM) interaction potentials. The density functional theory DFT(ωB97X-D3/6-31G**) approach was applied in QM and the CHARMM36 force field parameters in MM. The most likely scenario of the chemical step of the GTP hydrolysis in KRas-NF1 corresponds to the water-assisted mechanism of the formation of the inorganic phosphate coupled with the dissociation of GTP to GDP.

Interaction of molecular oxygen 3O2 with the flavin-dependent protein miniSOG after light illumination results in creation of singlet oxygen 1O2 and superoxide O2●−. Despite the recently resolved crystal structures of miniSOG variants, oxygen-binding sites near the flavin chromophore are poorly characterized. We report the results of computational studies of the protein−oxygen systems using molecular dynamics (MD) simulations with force-field interaction potentials and quantum mechanics/molecular mechanics (QM/MM) potentials for the original miniSOG and the mutated protein. We found several oxygen-binding pockets and pointed out possible tunnels bridging the bulk solvent and the isoalloxazine ring of the chromophore. These findings provide an essential step toward understanding photophysical properties of miniSOG—an important singlet oxygen photosensitizer.

Oxygenase activity of the flavin-dependent enzyme RutA is commonly associated with the formation of flavin-oxygen adducts in the enzyme active site. We report the results of quantum mechanics/molecular mechanics (QM/MM) modeling of possible reaction pathways initiated by various triplet state complexes of the molecular oxygen with the reduced flavin mononucleotide (FMN) formed in the protein cavities. According to the calculation results, these triplet-state flavin-oxygen complexes can be located at both re-side and si-side of the isoalloxazine ring of flavin. In both cases, the dioxygen moiety is activated by electron transfer from FMN, stimulating the attack of the arising reactive oxygen species at the C4a, N5, C6, and C8 positions in the isoalloxazine ring after the switch to the singlet state potential energy surface. The reaction pathways lead to the C(4a)-peroxide, N(5)-oxide, or C(6)-hydroperoxide covalent adducts or directly to the oxidized flavin, depending on the initial position of the oxygen molecule in the protein cavities.

We report the results of a computational study of the mechanism of the light-induced chemical reaction of chromophore hydration in the fluorescent protein Dreiklang, responsible for its switching from the fluorescent ON-state to the dark OFF-state. We explore the relief of the charge-transfer excited-state potential energy surface in the ON-state to locate minimum energy conical intersection points with the ground-state energy surface. Simulations of the further evolution of model systems allow us to characterize the ground-state reaction intermediate tentatively suggested in the femtosecond studies of the light-induced dynamics in Dreiklang and finally to arrive at the reaction product. The obtained results clarify the details of the photoswitching mechanism in Dreiklang, which is governed by the chemical modification of its chromophore.

We report the results of a computational study of the hydrolysis reaction mechanism of N-acetyl-l-aspartyl-l-glutamate (NAAG) catalyzed by glutamate carboxypeptidase II. Analysis of both mechanistic and electronic structure aspects of this multistep reaction is in the focus of this work. In these simulations, model systems are constructed using the relevant crystal structure of the mutated inactive enzyme. After selection of reaction coordinates, the Gibbs energy profiles of elementary steps of the reaction are computed using molecular dynamics simulations with ab initio type QM/MM potentials (QM/MM MD). Energies and forces in the large QM subsystem are estimated in the DFT(PBE0-D3/6-31G**) approximation. The established mechanism includes four elementary steps with the activation energy barriers not exceeding 7 kcal/mol. The models explain the role of point mutations in the enzyme observed in the experimental kinetic studies; namely, the Tyr552Ile substitution disturbs the “oxyanion hole”, and the Glu424Gln replacement increases the distance of the nucleophilic attack. Both issues diminish the substrate activation in the enzyme active site. To quantify the substrate activation, we apply the QTAIM-based approaches and the NBO analysis of dynamic features of the corresponding enzyme-substrate complexes. Analysis of the 2D Laplacian of electron density maps allows one to define structures with the electron density deconcentration on the substrate carbon atom, i.e., at the electrophilic site of reactants. The similar electronic structure element in the NBO approach is a lone vacancy on the carbonyl carbon atom in the reactive species. The electronic structure patterns revealed in the NBO and QTAIM-based analyses consistently clarify the reactivity issues in this system.