We propose a quantitative structure–property
relationship
(QSPR) model for prediction of spectral tuning in cyan, green, orange,
and red fluorescent proteins, which are engineered by motifs of the
green fluorescent protein. Protein variants, in which their chromophores
are involved in the π-stacking interaction with amino acid residues
tyrosine, phenylalanine, and histidine, are prospective markers useful
in bioimaging and super-resolution microscopy. In this work, we constructed
training sets of the π-stacked complexes of four fluorescent
protein chromophores (of the green, orange, red, and cyan series)
with various substituted benzenes and imidazoles and tested the use
of dipole moment variation upon excitation (DMV) as a descriptor to
evaluate the vertical excitation energies in these systems. To validate
this approach, we computed and analyzed electron density distributions
of the π-stacked complexes and correlated the QSPR predictions
with the reference values of the transition energies obtained using
the high-level ab initio quantum chemistry methods.
According to our results, the use of the DMV descriptor allows one
to predict excitation energies in the π-stacked complexes with
errors not exceeding 0.1 eV, which makes this model a practically
useful tool in the development of efficient fluorescent markers for in vivo imaging.
We report the results of calculations of the Gibbs energy profiles of the guanosine triphosphate (GTP) hydrolysis by the Arl3-RP2 protein complex using molecular dynamics (MD) simulations with ab initio type QM/MM potentials. The chemical reaction of GTP hydrolysis to guanosine diphosphate (GDP) and inorganic phosphate (Pi) is catalyzed by GTPases, the enzymes, which are responsible for signal transduction in live cells. A small GTPase Arl3, catalyzing the GTP → GDP reaction in complex with the activating protein RP2, constitute an essential part of the human vision cycle. To simulate the reaction mechanism, a model system is constructed by motifs of the crystal structure of the Arl3-RP2 complexed with a substrate analog. After selection of reaction coordinates, energy profiles for elementary steps along the reaction pathway GTP + H2O → GDP + Pi are computed using the umbrella sampling and umbrella integration procedures. QM/MM MD calculations are carried out, interfacing the molecular dynamics program NAMD and the quantum chemistry program TeraChem. Ab initio type QM(DFT)/MM potentials are computed with atom-centered basis sets 6-31G** and two hybrid functionals (PBE0-D3 and ωB97x-D3) of the density functional theory, describing a large QM subsystem. Results of these simulations of the reaction mechanism are compared to those obtained with QM/MM calculations on the potential energy surface using a similar description of the QM part. We find that both approaches, QM/MM and QM/MM MD, support the mechanism of GTP hydrolysis by GTPases, according to which the catalytic glutamine side chain (Gln71, in this system) actively participates in the reaction. Both approaches distinguish two parts of the reaction: the cleavage of the phosphorus-oxygen bond in GTP coupled with the formation of Pi, and the enzyme regeneration. Newly performed QM/MM MD simulations confirmed the profile predicted in the QM/MM minimum energy calculations, called here the pathway-I, and corrected its relief at the first elementary step from the enzyme–substrate complex. The QM/MM MD simulations also revealed another mechanism at the part of enzyme regeneration leading to pathway-II. Pathway-II is more consistent with the experimental kinetic data of the wild-type complex Arl3-RP2, whereas pathway-I explains the role of the mutation Glu138Gly in RP2 slowing down the hydrolysis rate.
Tautomerizarion between keto, cis‐ and trans‐enol forms of the 4,5‐dimethyl‐2‐(2′‐hydroxyphenyl)imidazole (DMHI) is studied using molecular dynamics simulations with the quantum mechanics/molecular mechanics (QM/MM) potentials. We applied the live solvent selection approach that presumes division of water molecules to either QM or MM subsystems on the fly depending on their distances from the chromophore. This allowed us to treat a water shell around the chromophore at the QM level and to describe properly all chromophore‐water interactions. According to the calculated Gibbs energy profiles the most stable is a keto form and a cis‐enol form is ~1 kcal/mol higher. The trans‐enol form is much higher in energy and should not be observed in experiment. Vertical S0,min − S1 excitation energies of all neutral tautomers were calculated at the XMCQDPT2 level: both enol forms have similar values ~4 eV and keto form is red‐shifted, with the excitation energy of 3.3 eV. Thus, we conclude that two absorption bands in experimental absorption spectrum of DMHI in water solution are due to the presence of cis‐enol and keto forms of the chromophore.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.