A proteolytic reaction of papain with a simple peptide model substrate N-methylacetamide has been studied. Our aim was twofold: (i) we proposed a plausible reaction mechanism with the aid of potential energy surface scans and second geometrical derivatives calculated at the stationary points, and (ii) we investigated the applicability of the dispersion corrected density functional methods in comparison with the popular hybrid generalized gradient approximations (GGA) method (B3LYP) without such a correction in the QM/MM calculations for this particular problem. In the resting state of papain the ion pair and neutral forms of the Cys–His catalytic dyad have approximately the same energy and they are separated by only a small barrier. Zero point vibrational energy correction shifted this equilibrium slightly to the neutral form. On the other hand, the electrostatic solvation free energy corrections, calculated using the Poisson–Boltzmann method for the structures sampled from molecular dynamics simulation trajectories, resulted in a more stable ion-pair form. All methods we applied predicted at least a two elementary step acylation process via a zwitterionic tetrahedral intermediate. Using dispersion corrected DFT methods the thioester S–C bond formation and the proton transfer from histidine occur in the same elementary step, although not synchronously. The proton transfer lags behind (or at least does not precede) the S–C bond formation. The predicted transition state corresponds mainly to the S–C bond formation while the proton is still on the histidine Nδ atom. In contrast, the B3LYP method using larger basis sets predicts a transition state in which the S–C bond is almost fully formed and the transition state can be mainly featured by the Nδ(histidine) to N(amid) proton transfer. Considerably lower activation energy was predicted (especially by the B3LYP method) for the next amide bond breaking elementary step of acyl–enzyme formation. Deacylation appeared to be a single elementary step process in all the methods we applied.

Tetrodotoxin (TTX) is believed to be one of the most selective inhibitors of voltage-gated fast Na⁺ channels in excitable tissues. Recently, however, TTX has been shown to block L-type Ca²⁺ current (I_Ca) in canine cardiac cells. In the present study, the TTX-sensitivity of I_Ca was studied in isolated canine ventricular myocytes as a function of (1) channel phosphorylation, (2) extracellular pH and (3) the redox potential of the bathing medium using the whole cell voltage clamp technique. Fifty-five micromoles of TTX (IC₅₀ value obtained under physiological conditions) caused 60% ± 2% inhibition of I_Ca in acidic (pH = 6.4), while only a 26% ± 2% block in alkaline (pH = 8.4) milieu. Similarly, the same concentration of TTX induced 62% ± 6% suppression of I_Ca in a reductant milieu (containing glutathione + ascorbic acid + dithiothreitol, 1 mM each), in contrast to the 31% ± 3% blockade obtained in the presence of a strong oxidant (100 μM H₂O₂). Phosphorylation of the channel protein (induced by 3 μM forskolin) failed to modify the inhibiting potency of TTX; an IC₅₀ value of 50 ± 4 μM was found in forskolin. The results are in a good accordance with the predictions of our model, indicating that TTX binds, in fact, to the selectivity filter of cardiac L-type Ca channels.

While antithrombin (AT) has small basal inhibitory activity, it reaches its full inhibitory potential against activated blood coagulation factors, FXa, FIXa, and FIIa (thrombin), via an allosteric and/or template (bridging) mechanism by the action of heparin, heparan sulfate, or heparin-mimetic pentasaccharides (PS). From the numerous X-ray structures available for different conformational states of AT, only indirect and incomplete conclusions can be drawn on the inherently dynamic properties of AT. As a typical example, the basal inhibitory activity of AT cannot be interpreted on the basis of “non-activated” free antithrombin X-ray structures since the Arg393 side chain, playing crucial role in antithrombin-proteinase interaction, is not exposed. In order to reveal the intrinsic dynamic properties and the reason of basal inhibitory activity of antithrombin, 2 μs molecular dynamics simulations were carried out on its native free-forms. It was shown from the simulation trajectories that the reactive center loop which is functioning as “bait” for proteases, even without any biasing potential can populate conformational state in which the Arg393 side chain is solvent exposed. It is revealed from the trajectory analysis that the peptide sequences correspond to the helix D extension, and new helix P formation can be featured with especially large root-mean-square fluctuations. Mutual information analyses of the trajectory showed remarkable (generalized) correlation between those regions of antithrombin which changed their conformations as the consequence of AT–PS complex formation. This suggests that allosteric information propagation pathways are present even in the non-activated native form of AT.

Von Hippel-Lindau disease is an autosomal dominantly inherited highly penetrant tumor syndrome predisposing to retinal and central nervous system hemangioblastomas, renal cell carcinoma and phaeochromocytoma among other less frequent complications.

A simplified molecular-dynamics-based electronic circular dichroism (ECD) approach was tested on three condensed derivatives with limited conformational flexibility and an isochroman-2H-chromene hybrid, the ECD spectra of which could not be precisely reproduced by the conventional ECD calculation protocol. Application of explicit solvent molecules at the molecular mechanics (MD) level in the dynamics simulations and subsequent TDDFT-ECD calculation for the unoptimized MD structures was able to improve the agreements between experimental and computed spectra. Since enhancements were achieved even for molecules with limited conformational flexibility, deformations caused by the solvent molecules and multitudes of conformers produced with unoptimized geometries seem to be key factors for better agreement. The MD approach could confirm that aggregation of the phenanthrene natural product luzulin A had a significant contribution to a specific wavelength range of the experimental ECD. The MD approach has proved that dimer formation occurred in solution and this was responsible for the anomalous ECD spectrum. The scope and limitations of the method have also been discussed.