Mechanism of proteolysis in matrix metalloproteinase-2 revealed by QM/MM modeling

The mechanism of enzymatic peptide hydrolysis in matrix metalloproteinase-2 (MMP-2) was studied at atomic resolution through quantum mechanics/molecular mechanics (QM/MM) simulations. An all-atom three-dimensional molecular model was constructed on the basis of a crystal structure from the Protein Data Bank (ID: 1QIB), and the oligopeptide Ace-Gln-Gly∼Ile-Ala-Gly-Nme was considered as the substrate. Two QM/MM software packages and several computational protocols were employed to calculate QM/MM energy profiles for a four-step mechanism involving an initial nucleophilic attack followed by hydrogen bond rearrangement, proton transfer, and C—N bond cleavage. These QM/MM calculations consistently yield rather low overall barriers for the chemical steps, in the range of 5–10 kcal/mol, for diverse QM treatments (PBE0, B3LYP, and BB1K density functionals as well as local coupled cluster treatments) and two MM force fields (CHARMM and AMBER). It, thus, seems likely that product release is the rate-limiting step in MMP-2 catalysis. This is supported by an exploration of various release channels through QM/MM reaction path calculations and steered molecular dynamics simulations

Role of the flat-designed surface in improving the cyclic fatigue resistance of endodontic NiTi rotary instruments

The aim of this study was to investigate the role of the flat-designed surface in improving the resistance to cyclic fatigue by comparing heat-treated F-One (Fanta Dental, Shanghai, China) nickel-titanium (NiTi) rotary instruments and similar prototypes, differing only by the absence of the flat side. The null hypothesis was that there were no differences between the two tested instruments in terms of cyclic fatigue lifespan. A total of 40 new NiTi instruments (20 F-One and 20 prototypes) were tested in the present study. The instruments were rotated with the same speed (500 rpm) and torque (2 N) using an endodontic motor (Elements Motor, Kerr, Orange, CA, USA) in the same stainless steel, artificial canal (90° angle of curvature and 5 mm radius). A Wilcoxon-Mann-Whitney test was performed to assess the differences in terms of time to fracture and the length of the fractured segment between the flat- and non-flat-sided instruments. Significance was set at p = 0.05. The differences in terms of time to fracture between non-flat and flat were statistically significant (p < 0.001). In addition, the differences in terms of fractured segment length were statistically significant (p = 0.034). The results of this study highlight the importance of flat-sided design in increasing the cyclic fatigue lifespan of NiTi rotary instruments

BLUF Domain Function Does Not Require a Metastable Radical Intermediate State

BLUF (blue light using flavin) domain proteins are an important family of blue light-sensing proteins which control a wide variety of functions in cells. The primary light-activated step in the BLUF domain is not yet established. A number of experimental and theoretical studies points to a role for photoinduced electron transfer (PET) between a highly conserved tyrosine and the flavin chromophore to form a radical intermediate state. Here we investigate the role of PET in three different BLUF proteins, using ultrafast broadband transient infrared spectroscopy. We characterize and identify infrared active marker modes for excited and ground state species and use them to record photochemical dynamics in the proteins. We also generate mutants which unambiguously show PET and, through isotope labeling of the protein and the chromophore, are able to assign modes characteristic of both flavin and protein radical states. We find that these radical intermediates are not observed in two of the three BLUF domains studied, casting doubt on the importance of the formation of a population of radical intermediates in the BLUF photocycle. Further, unnatural amino acid mutagenesis is used to replace the conserved tyrosine with fluorotyrosines, thus modifying the driving force for the proposed electron transfer reaction; the rate changes observed are also not consistent with a PET mechanism. Thus, while intermediates of PET reactions can be observed in BLUF proteins they are not correlated with photoactivity, suggesting that radical intermediates are not central to their operation. Alternative nonradical pathways including a keto–enol tautomerization induced by electronic excitation of the flavin ring are considered

Photoinduced electron transfer facilitates tautomerization of the conserved signaling glutamine side chain in BLUF protein light sensors

The BLUF domain (sensor of blue light using flavin adenine dinucleotide) from a bacterial photoreceptor protein AppA undergoes a cascade of chemical transformations, including hydrogen bond rearrangements around the flavin adenine dinucleotide (FAD) chromophore, in response to light illumination. These transformations are initiated by photoinduced electron and proton transfer from a tyrosine residue to the photoexcited flavin which is assisted by a glutamine residue. According to the recent studies, the proton−coupled electron transfer leads to formation of a radical−pair intermediate Tyr•···FADH• and a tautomeric EE form of glutamine in the ground electronic state. This intermediate is a precursor of the light−induced state of the BLUF photoreceptor implicated in biological signaling. In order to describe evolution of the radical pair, we computed reaction pathways on the ground state potential energy surface employing quantum−chemical calculations in the DFT PBE0/cc−pVDZ approximation for a molecular cluster mimicking the chromophore containing pocket of the AppA BLUF protein. We found a minimum−energy pathway comprised of the following consecutive reaction steps: (1) rotation of the imidic group of the EE glutamine side chain around the Cγ−Cδ bond; (2) flip of the OεH group and formation of the ZE form of the glutamine side chain; and (3) biradical recombination via coupled proton and electron transfer, leading to the ZZ form of the glutamine side chain. The potential−energy barriers for stages 1−3 do not exceed 9 kcal/mol. Energy barrier 3 describing the ZE to ZZ glutamine tautomerization is significantly smaller in the BLUF model than in isolated glutamine, since tautomerization in BLUF is facilitated by electron transfer and radical recombination. Thus, our study shows that tautomerization of the conserved glutamine is coupled to the light−induced electron transfer process in BLUF and, thus, is a viable candidate for the photoactivation mechanism which at present is very much debate

Molecular mechanism of the dark-state recovery in BLUF photoreceptors

Photoactivation of the bacterial photoreceptor BLUF is achieved by a unique tautomerization and rotation of the conserved Gln residue. To identify interactions controlling the lifetime of the high-energy ZZ-imide tautomer, we consider thermal dark-state recovery of the BLUF photoreceptor BlrB using MD and QM/MM calculations. We found that protonation/deprotonation of Gln51 by protonated His73 via Tyr9 was energetically favorable. The rate-limiting steps correspond to rearrangements of hydrogen bonds and rotation of the protonated Gln51 side chain. Notably, strengthening the hydrogen bond between protonated Gln51 and flavin O4 plays a central role in extending the lifetime of the photoactivated state

Molecular Mechanism of Light-Induced Acceleration of the Adenosine Triphosphate Conversion to the Cyclic Adenosine Monophosphate Catalyzed by the Photoactivated Adenylyl Cyclase bPAC

We report the first computational characterization of an optogenetic system composed of two photosensing BLUF (blue light sensor using flavin adenine dinucleotide) domains and two catalytic adenylyl cyclase (AC) domains. Conversion of adenosine triphosphate (ATP) to cyclic adenosine monophosphate (cAMP) and pyrophosphate (PPi) catalyzed by ACs coupled with excitation in photosensing domains has emerged in the focus of modern optogenetic applications because of the request in photoregulated enzymes to modulate cellular concentrations of signaling messengers. The photoactivated adenylyl cyclase from the soil bacterium Beggiatoa sp. (bPAC) is an important model showing considerable increase of the ATP to cAMP conversion rate in the catalytic domain after the illumination of the BLUF domain. The 1 μs classical molecular dynamics simulations reveal that the activation of the BLUF domain leading to tautomerization of Gln49 in the chromophore binding pocket results in switching of position of the side chain of Arg278 in the active site of AC. Allosteric signal transmission pathways between Gln49 from BLUF and Arg278 from AC were revealed by the dynamical network analysis. The Gibbs energy profiles of the ATP → cAMP + PPi reaction computed using QM(DFT(ωB97X-D3/6-31G**))/MM(CHARMM) molecular dynamics simulations for both Arg278 conformations in AC clarify the reaction mechanism. In the light-activated system, the corresponding arginine conformation stabilizes the pentacoordinated phosphorus of the α-phosphate group in the transition state, thus lowering the activation energy. Simulations of the bPAC system with the Tyr7Phe replacement in BLUF demonstrate occurrence of both arginine conformations in an equal ratio, explaining the experimentally observed intermediate catalytic activity of the bPAC-Y7F variant as compared with the dark and light states of the wild type bPAC. </p

Methodological aspects of QM/MM calculations: A case study on matrix metalloproteinase-2

We address methodological issues in quantum mechanics/molecular mechanics (QM/MM) calculations on a zinc-dependent enzyme. We focus on the first stage of peptide bond cleavage by matrix metalloproteinase-2 (MMP-2), that is, the nucleophilic attack of the zinc-coordinating water molecule on the carbonyl carbon atom of the scissile fragment of the substrate. This step is accompanied by significant charge redistribution around the zinc cation, bond cleavage, and bond formation. We vary the size and initial geometry of the model system as well as the computational protocol to demonstrate the influence of these choices on the results obtained. We present QM/MM potential energy profiles for a set of snapshots randomly selected from QM/MM-based molecular dynamics simulations and analyze the differences in the computed profiles in structural terms. Since the substrate in MMP-2 is located on the protein surface, we investigate the influence of the thickness of the water layer around the enzyme on the QM/MM energy profile. Thin water layers (0–2 Å) give unrealistic results because of structural reorganizations in the active-site region at the protein surface. A 12 Å water layer appears to be sufficient to capture the effect of the solvent; the corresponding QM/MM energy profile is very close to that obtained from QM/MM/SMBP calculations using the solvent macromolecular boundary potential (SMBP). We apply the optimized computational protocol to explain the origin of the different catalytic activity of the Glu116Asp mutant: the energy barrier for the first step is higher, which is rationalized on structural grounds