Assessing the Performance of Peptide Force Fields for Modeling the Solution Structural Ensembles of Cyclic Peptides.

Molecular dynamics simulation is a powerful tool for characterizing the solution structural ensembles of cyclic peptides. However, the ability of simulation to recapitulate experimental results and make accurate predictions largely depends on the force fields used. In our work here, we evaluate the performance of seven state-of-the-art force fields in recapitulating the experimental NMR results in water of 12 benchmark cyclic peptides, consisting of 6 cyclic pentapeptides, 4 cyclic hexapeptides, and 2 cyclic heptapeptides. The results show that RSFF2+TIP3P, RSFF2C+TIP3P, and Amber14SB+TIP3P exhibit similar and the best performance, all recapitulating the NMR-derived structure information on 10 cyclic peptides. Amber19SB+OPC successfully recapitulates the NMR-derived structure information on 8 cyclic peptides. In contrast, OPLS-AA/M+TIP4P, Amber03+TIP3P, and Amber14SBonlysc+GB-neck2 could only recapitulate the NMR-derived structure information on 5 cyclic peptides, the majority of which are not well-structured.

Photolytic properties of B12-dependent enzymes: A theoretical perspective.

The biologically active vitamin B12 derivates, methylcobalamin (MeCbl) and adenosylcobalamin (AdoCbl), are ubiquitous organometallic cofactors. In addition to their key roles in enzymatic catalysis, B12 cofactors have complex photolytic properties which have been the target of experimental and theoretical studies. With the recent discovery of B12-dependent photoreceptors, there is an increased need to elucidate the underlying photochemical mechanisms of these systems. This book chapter summarizes the photolytic properties of MeCbl- and AdoCbl-dependent enzymes with particular emphasis on the effect of the environment of the cofactor on the excited state processes. These systems include isolated MeCbl and AdoCbl as well as the enzymes, ethanolamine ammonia-lyase (EAL), glutamate mutase (GLM), methionine synthase (MetH), and photoreceptor CarH. Central to determining the photodissociation mechanism of each system is the analysis of the lowest singlet excited state (S1) potential energy surface (PES). Time-dependent density functional theory (TD-DFT), employing BP86/TZVPP, is widely used to construct such PESs. Regardless of the environment, the topology of the S1 PES of AdoCbl or MeCbl is marked by characteristic features, namely the metal-to-ligand charge transfer (MLCT) and ligand field (LF) regions. Conversely, the relative energetics of these electronic states are affected by the environment. Applications and outlooks for Cbl photochemistry are also discussed.

Aerobic photolysis of methylcobalamin: unraveling the photoreaction mechanism.

The photo-reactivity of cobalamins (Cbls) is influenced by the nature of axial ligands and the cofactor's environment. While the biologically active forms of Cbls with alkyl axial ligands, such as methylcobalamin (MeCbl) and adenosylcobalamin (AdoCbl), are considered to be photolytically active, in contrast, the non-alkyl Cbls are photostable. In addition to these, the photolytic properties of Cbls can also be modulated in the presence of molecular oxygen, i.e., under aerobic conditions. Herein, the photoreaction of the MeCbl in the presence of oxygen has been explored using density functional theory (DFT) and time-dependent DFT (TD-DFT). The first stage of the aerobic photoreaction is the activation of the Co-C bond and the formation of the ligand field (LF) electronic state through the displacement of axial bonds. Once the photoreaction reaches the LF excited state, three processes can occur: namely the formation of OO-CH3 through the reaction of CH3 with molecular oxygen, de-activation of the {Im⋯[CoII(corrin)]⋯CH3}+ sub-system from the LF electronic state by changing the electronic configuration from (dyz)1(dz2)2 to (dyz)2(dz2)1 and the formation of the deactivation complex (DC) complex via the recombination of OO-CH3 species with the de-excited [CoII(corrin)] system. In the proposed mechanism, the deactivation of the [CoII(corrin)] subsystem may coexist with the formation of OO-CH3, followed by immediate relaxation of the subsystems in the ground state. Moreover, the formation of the OO-CH3 species followed by the formation of the {[CoIII(corrin)]-OO-CH3}+ complex stabilizes the system compared to the reactant complex.

Methyl transfer reactions catalyzed by cobalamin-dependent enzymes: Insight from molecular docking.

Methyl transfer reactions, mediated by methyltransferases (MeTrs), such as methionine synthase (MetH) or monomethylamine: CoM (MtmBC), constitute one of the most important classes of vitamin B12-dependent reactions. The challenge in exploring the catalytic function of MeTrs is related to their modular structure. From the crystallographic point of view, the structure of each subunit has been determined, but there is a lack of understanding of how each subunit interacts with each other. So far, theoretical studies of methyl group transfer were carried out for the structural models of the active site of each subunit. However, those studies do not include the effect of the enzymatic environment, which is crucial for a comprehensive understanding of enzyme-mediated methyl transfer reactions. Herein, to explore how two subunits interact with each other and how the methyl transfer reaction is catalyzed by MeTrs, molecular docking of the functional units of MetH and MtmBC was carried out. Along with the interactions of the functional units, the reaction coordinates, including the Co-C bond distance for methylation of cob(I)alamin (CoICbl) and the C-S bond distance in demethylation reaction of cob(III)alamin (CoIIICbl), were considered. The functional groups should be arranged so that there is an appropriate distance to transfer a methyl group and present results indicate that steric interactions can limit the number of potential arrangements. This calls into question the possibility of SN2-type mechanism previously proposed for MeTrs. Further, it leads to the conclusion that the methyl transfer reaction involves some spatial changes of modules suggesting an alternate radical-based pathway for MeTrs-mediated methyl transfer reactions. The calculations also showed that changes in torsion angles induce a change in reaction coordinates, namely Co-C and C-S bond distances, for the methylation and demethylation reactions catalyzed both by MetH and MtmBC.

Aerobic photolysis of methylcobalamin: structural and electronic properties of the Cbl-O-O-CH3 intermediate.

Photolysis of methylcobalamin (MeCbl) in the presence of molecular oxygen (O2) has been investigated using density functional theory (DFT) and time-dependent DFT (TD-DFT). The key step involves the formation of the Cbl-O-O-CH3 intermediate as a result of triplet O2 insertion in the Co-C bond in the presence of light. Analysis of low-lying excited states shows that the presence of light is only needed to activate the Co-C bond via the formation of the ligand field (LF) state. The insertion of O2, as well as the change in the spin state, takes place in the ground state. The analysis of the structural and electronic properties of the Cbl-O-O-CH3 intermediate is presented and possible decomposition also discussed.

How does the mutation in the cap domain of methylcobalamin-dependent methionine synthase influence the photoactivation of the Co-C bond?

Methionine synthase (MetH) is a methylcobalamin (MeCbl)-dependent mammalian enzyme which plays a critical role in carrying out the transfer of a methyl group from methyl tetrahydrofolate to homocysteine to generate methionine and tetrahydrofolate. This catalytic cycle proceeds via cleavage of a Co-C bond which is formally heterolytic. This cleavage results in a structural change in the MeCbl cofactor bound to an enzyme. Unlike the native catalysis, upon photoexcitation, the Co-C bond in MeCbl-bound MetH generates the Co(ii)/CH3 radical pairs (RPs). Protein residues of the cap domain, particularly phenylalanine708 (F708) and leucine 715 (L715), which surrounds the upper face of the MeCbl cofactor, inhibit the photolysis of MeCbl by caging the CH3 radical and inducing the geminate recombination of the Co(ii)/CH3 RP. A molecular-level understanding of these effects requires a detailed investigation of the low-lying electronic states. Toward this, we have mutated the F708 residue with alanine (A708) and constructed the potential energy surfaces (PESs) for the low-lying S1 electronic state using a combined quantum mechanics/molecular mechanics (QM/MM) approach. The S1 PESs for the wild-type (WT) and mutant enzymes are the result of crossing of two electronic states, namely metal-to-ligand charge transfer (MLCT) and ligand field (LF) states, indicated by a seam. It is shown that the topologies of the S1 PESs are significantly modulated by introducing a mutation at the F708 position. Specifically, for the WT enzyme, the energy barrier of photoreaction and the energy difference between MLCT and LF minima are markedly higher than those of its mutant counterpart. Moreover, mutation influences the photoactivation of the Co-C bond in enzyme-bound MeCbl by decreasing the rate of geminate recombination and altering the rate of radical pair formation. This theoretical insight was also compared with transient absorption spectroscopic (TAS) studies which are in good agreement with the present findings.

Mechanism of the photo-induced activation of CoC bond in methylcobalamin-dependent methionine synthase.

Methylcobalamin (MeCbl)-dependent enzyme methionine synthase (MetH), plays a critical role in the catalysis of methyl group transfer from methyltetrahydrofolate (CH3-H4folate) to homocysteine. It often performs a side reaction to generate cob(II)alamin through photolysis of the organometallic CoC σ bond. A hybrid QM/MM method has been applied to explore the photochemistry of MeCbl-bound MetH. The photolytic properties of MeCbl inside MetH are mediated by its manifold of low-lying excited states. The corresponding potential energy surfaces (PESs) of the electronically excited S1 state has been constructed as a function of axial bond lengths to elucidate the mechanism of photo-induced activation of CoC bond inside the enzyme. The analysis of the S1 PES has revealed that the two different electronic states of the S1 PES, namely metal-to-ligand charge transfer (MLCT) and the ligand field (LF), are relevant to the photodissociation of the CoC bond. There are two possible pathways identified, Path A and Path B, that connect the MLCT to LF state that represent possible photodissociation mechanisms. In the case of MetH, one possible photodissociation pathway (Path B) was identified based on the energetics of the MLCT and LF states. The energetically accessible Path B involves the initial detachment of the Co-NIm bond followed by a subsequent displacement of the CoC bond prior to the formation of cob(II)alamin / CH3 radical pair (RP). The photochemical data of base-on MeCbl in solution was compared with the computed result of MeCbl-bound MetH to understand the effect of the enzymatic environment on the photolytic properties of MeCbl.