Refining the treatment of membrane proteins by coarse-grained models.

Obtaining a quantitative description of the membrane proteins stability is crucial for understanding many biological processes. However the advance in this direction has remained a major challenge for both experimental studies and molecular modeling. One of the possible directions is the use of coarse-grained models but such models must be carefully calibrated and validated. Here we use a recent progress in benchmark studies on the energetics of amino acid residue and peptide membrane insertion and membrane protein stability in refining our previously developed coarse-grained model (Vicatos et al., Proteins 2014;82:1168). Our refined model parameters were fitted and/or tested to reproduce water/membrane partitioning energetics of amino acid side chains and a couple of model peptides. This new model provides a reasonable agreement with experiment for absolute folding free energies of several ß-barrel membrane proteins as well as effects of point mutations on a relative stability for one of those proteins, OmpLA. The consideration and ranking of different rotameric states for a mutated residue was found to be essential to achieve satisfactory agreement with the reference data.

Methyltransferases do not work by compression, cratic, or desolvation effects, but by electrostatic preorganization.

The enzyme catechol O-methyltransferase (COMT) catalyzes the transfer of a methyl group from S-adenosylmethionine to dopamine and related catechols. The search for the origin of COMT catalysis has led to different proposals and hypothesis, including the entropic, the NAC, and the compression proposals as well as the more reasonable electrostatic idea. Thus, it is important to understand the catalytic power of this enzyme and to examine the validity of different proposals and in particular the repeated recent implication of the compression idea. The corresponding analysis should be done by well-defined physically-based considerations that involve computations rather than circular interpretations of experimental results. Thus, we explore here the origin of the catalytic efficiency of COMT by using the empirical valence bond and the linear response approximation approaches. The results demonstrate that the catalytic effect of COMT is mainly due to electrostatic preorganization effects. It is also shown that the compression, NAC and entropic proposals do not account for the catalytic effect.

Exploring challenges in rational enzyme design by simulating the catalysis in artificial kemp eliminase.

One of the fundamental challenges in biotechnology and in biochemistry is the ability to design effective enzymes. Doing so would be a convincing manifestation of a full understanding of the origin of enzyme catalysis. Despite an impressive progress, most of the advances on this front have been made by placing the reacting fragments in the proper places, rather than by optimizing the environment preorganization, which is the key factor in enzyme catalysis. Rational improvement of the preorganization would require approaches capable of evaluating reliably the actual catalytic effect. This work takes previously designed kemp eliminases as a benchmark for a computer aided enzyme design, using the empirical valence bond as the main screening tool. The observed absolute catalytic effect and the effect of directed evolution are reproduced and analyzed (assuming that the substrate is in the designed site). It is found that, in the case of kemp eliminases, the transition state charge distribution makes it hard to exploit the active site polarity, even with the ability to quantify the effect of different mutations. Unexpectedly, it is found that the directed evolution mutants lead to the reduction of solvation of the reactant state by water molecules rather that to the more common mode of transition state stabilization used by naturally evolved enzymes. Finally it is pointed out that our difficulties in improving Kemp eliminase are not due to overlooking exotic effect, but to the challenge in designing a preorganized environment that would exploit the small change it charge distribution during the formation of the transition state.

Ketosteroid isomerase provides further support for the idea that enzymes work by electrostatic preorganization.

One of the best systems for exploring the origin of enzyme catalysis has been the reaction of ketosteroid isomerase (KSI). Studies of the binding of phenolates to KSI have been taken as proof that the electrostatic preorganization effect only makes a minor contribution to the binding of the real, multiring, transition state (TS). However, our simulation study has determined that the difference between the phenolates and the TS arises from the fact that the nonpolar state of the phenolate can rotate freely relative to the oxyanion hole and thus loses the preorganization contribution. A recent study explored the reactivity of both small and multiring systems and concluded that their similar reactivity contradicts our preorganization idea. Herein, we establish that the available experiments in fact provide what is perhaps the best proof and clarification of the preorganization idea and its crucial role in enzyme catalysis. First, we analyze the binding energy and the pK(a) of equilenin and identify direct experimental evidence for our prediction about the differential electrostatic stabilization of the large TS and the small phenolates. Subsequently, we show that the similar reactivity of the small and large systems is also due to an electrostatic preorganization effect but that this effect only appears in the intermediate state because the TS is not free to rotate. This establishes the electrostatic origin of enzyme catalysis. We also clarify the crucial importance of having a well-defined physical concept when examining catalytic effects and the need for quantitative tools for analyzing such effects.

Electrostatic basis for the unidirectionality of the primary proton transfer in cytochrome c oxidase.

Gaining detailed understanding of the energetics of the proton-pumping process in cytochrome c oxidase (CcO) is one of the challenges of modern biophysics. Despite promising mechanistic proposals, most works have not related the activation barriers of the different assumed steps to the protein structure, and there has not been a physically consistent model that reproduced the barriers needed to create a working pump. This work reevaluates the activation barriers for the primary proton transfer (PT) steps by calculations that reflect all relevant free energy contributions, including the electrostatic energies of the generated charges, the energies of water insertion, and large structural rearrangements of the donor and acceptor. The calculations have reproduced barriers that account for the directionality and sequence of events in the primary PT in CcO. It has also been found that the PT from Glu-286 (E) to the propionate of heme a(3) (Prd) provides a gate for an initial back leakage from the high pH side of the membrane. Interestingly, the rotation of E that brings it closer to Prd appears to provide a way for blocking competing pathways in the primary PT. Our study elucidates and quantifies the nature of the control of the directionality in the primary PT in CcO and provides instructive insight into the role of the water molecules in biological PT, showing that "bridges" of several water molecules in hydrophobic regions present a problem (rather than a solution) that is minimized in the primary PT.

Multiscale simulations of protein landscapes: using coarse-grained models as reference potentials to full explicit models.

Evaluating the free-energy landscape of proteins and the corresponding functional aspects presents a major challenge for computer simulation approaches. This challenge is due to the complexity of the landscape and the enormous computer time needed for converging simulations. The use of simplified coarse-grained (CG) folding models offers an effective way of sampling the landscape but such a treatment, however, may not give the correct description of the effect of the actual protein residues. A general way around this problem that has been put forward in our early work (Fan et al., Theor Chem Acc 1999;103:77-80) uses the CG model as a reference potential for free-energy calculations of different properties of the explicit model. This method is refined and extended here, focusing on improving the electrostatic treatment and on demonstrating key applications. These applications include: evaluation of changes of folding energy upon mutations, calculations of transition-states binding free energies (which are crucial for rational enzyme design), evaluations of catalytic landscape, and evaluations of the time-dependent responses to pH changes. Furthermore, the general potential of our approach in overcoming major challenges in studies of structure function correlation in proteins is discussed.

On catalytic preorganization in oxyanion holes: highlighting the problems with the gas-phase modeling of oxyanion holes and illustrating the need for complete enzyme models.

Oxyanion holes play a major role in catalyzing enzymatic reactions, yet the corresponding energetics is frequently misunderstood. The main problem may be associated with the nontrivial nature of the electrostatic preorganization effect, without following the relevant formulation. That is, although the energetics of oxyanion holes have been fully quantified in early studies (which include both the enzymatic and reference solution reactions), the findings of these studies are sometimes overlooked, and, in some cases, it is assumed that gas-phase calculations with a fixed model of an oxyanion hole are sufficient for assessing the corresponding effect in the protein. Herein, we present a systematic analysis of this issue, clarifying the problems associated with modeling oxyanions by means of two fixed water molecules (or related constructs). We then re-emphasize the point that the effect of the oxyanion hole is mainly due to the fact that the relevant dipoles are already set in an orientation that stabilizes the TS charges, whereas the corresponding dipoles in solution are randomly oriented, resulting in the need to pay a very large reorganization energy. Simply calculating interaction energies with relatively fixed species cannot capture this crucial point, and considering it may help in advancing rational enzyme design.

The energetics of the primary proton transfer in bacteriorhodopsin revisited: it is a sequential light-induced charge separation after all.

The light-induced proton transport in bacteriorhodopsin has been considered as a model for other light-induced proton pumps. However, the exact nature of this process is still unclear. For example, it is not entirely clear what the driving force of the initial proton transfer is and, in particular, whether it reflects electrostatic forces or other effects. The present work simulates the primary proton transfer (PT) by a specialized combination of the EVB and the QCFF/PI methods. This combination allows us to obtain sufficient sampling and a quantitative free energy profile for the PT at different protein configurations. The calculated profiles provide new insight about energetics of the primary PT and its coupling to the protein conformational changes. Our finding confirms the tentative analysis of an earlier work (A. Warshel, Conversion of light energy to electrostatic energy in the proton pump of Halobacterium halobium, Photochem. Photobiol. 30 (1979) 285-290) and determines that the overall PT process is driven by the energetics of the charge separation between the Schiff base and its counterion Asp85. Apparently, the light-induced relaxation of the steric energy of the chromophore leads to an increase in the ion-pair distance, and this drives the PT process. Our use of the linear response approximation allows us to estimate the change in the protein conformational energy and provides the first computational description of the coupling between the protein structural changes and the PT process. It is also found that the PT is not driven by twist-modulated changes of the Schiff base's pKa, changes in the hydrogen bond directionality, or other non-electrostatic effects. Overall, based on a consistent use of structural information as the starting point for converging free energy calculations, we conclude that the primary event should be described as a light-induced formation of an unstable ground state, whose relaxation leads to charge separation and to the destabilization of the ion-pair state. This provides the driving force for the subsequent PT steps.

Accelerating QM/MM free energy calculations: representing the surroundings by an updated mean charge distribution.

Reliable studies of enzymatic reactions by combined quantum mechanical/molecular mechanics (QM(ai)/MM) approaches with an ab initio description of the quantum region presents a major challenge to computational chemists. The main problem is the need for very large computer time to evaluate the QM energy, which in turn makes it extremely challenging to perform proper configurational sampling. One of the most obvious options for accelerating QM/MM simulations is the use of an average solvent potential. In fact, the idea of using an average solvent potential is rather obvious and has implicitly been used in Langevin dipole/QM calculations. However, in the case of explicit solvent models the practical implementations are more challenging, and the accuracy of the averaging approach has not been validated. The present study introduces the average effect of the fluctuating solvent charges by using equivalent charge distributions, which are updated every m steps. Several models are evaluated in terms of the resulting accuracy and efficiency. The most effective model divides the system into an inner region with N explicit solvent atoms and an external region with two effective charges. Different models are considered in terms of the division of the solvent system and the update frequency. Another key element of our approach is the use of the free energy perturbation (FEP) and/or linear response approximation treatments that guarantees the evaluation of the rigorous solvation free energy. Special attention is paid to the convergence of the calculated solvation free energies and the corresponding solute polarization. The performance of the method is examined by evaluating the solvation of a water molecule and a formate ion in water and also the dipole moment of water in water solution. Remarkably, it is found that different averaging procedures eventually converge to the same value but some protocols provide optimal ways of obtaining the final QM(ai)/MM converged results. The current method can provide computational time saving of 1000 for properly converging simulations relative to calculations that evaluate the QM(ai)/MM energy every time step. A specialized version of our approach that starts with a classical FEP charging and then evaluates the free energy of moving from the classical potential to the QM/MM potential appears to be particularly effective. This approach should provide a very powerful tool for QM(ai)/MM evaluation of solvation free energies in aqueous solutions and proteins.

Path Integral Simulation of the H/D Kinetic Isotope Effect in Monoamine Oxidase B Catalyzed Decomposition of Dopamine.

Brain monoamines regulate many centrally mediated body functions, and can cause adverse symptoms when they are out of balance. A starting point to address challenges raised by the increasing burden of brain diseases is to understand, at atomistic level, the catalytic mechanism of an essential amine metabolic enzyme-monoamine oxidase B (MAO B). Recently, we demonstrated that the rate-limiting step of MAO B catalyzed conversion of amines into imines represents the hydride anion transfer from the substrate α-CH2 group to the N5 atom of the flavin cofactor moiety. In this article we simulated for MAO B catalyzed dopamine decomposition the effects of nuclear tunneling by the calculation of the H/D kinetic isotope effect. We applied path integral quantization of the nuclear motion for the methylene group and the N5 atom of the flavin moiety in conjunction with the QM/MM treatment on the empirical valence bond (EVB) level for the rest of the enzyme. The calculated H/D kinetic isotope effect of 12.8 ± 0.3 is in a reasonable agreement with the available experimental data for closely related biogenic amines, which gives strong support for the proposed hydride mechanism. The results are discussed in the context of tunneling in enzyme centers and advent of deuterated drugs into clinical practice.

Simulating the catalytic effect of a designed mononuclear zinc metalloenzyme that catalyzes the hydrolysis of phosphate triesters.

One of the greatest challenges in biotechnology and in biochemistry is the ability to design efficient enzymes. In fact, such an ability would be one of the most convincing manifestations of a full understanding of the origin of enzyme catalysis. Despite some progress on this front, most of the advances have been made by placing the reacting fragments in the proper places rather than by optimizing the preorganization of the environment, which is the key factor in enzyme catalysis. A rational improvement of the preorganization and a consistent assessment of the effectiveness of different design options require approaches capable of evaluating reliably the actual catalytic effect. In this work we examine the ability of the empirical valence bond (EVB) to reproduce the results of directed evolution improvements of the catalysis of diethyl 7-hydroxycoumarinyl by a designed mononuclear zinc metalloenzyme. Encouragingly, our study reproduced the catalytic effect obtained by directed evolution and offers a good start for further studies of this system.

Quantifying the mechanism of phosphate monoester hydrolysis in aqueous solution by evaluating the relevant ab initio QM/MM free-energy surfaces.

Understanding the nature of the free-energy surfaces for phosphate hydrolysis is a prerequisite for understanding the corresponding key chemical reactions in biology. Here, the challenge has been to move to careful ab initio QM/MM (QM(ai)/MM) free-energy calculations, where obtaining converging results is very demanding and computationally expensive. This work describes such calculations, focusing on the free-energy surface for the hydrolysis of phosphate monoesters, paying special attention to the comparison between the one water (1W) and two water (2W) paths for the proton-transfer (PT) step. This issue has been explored before by energy minimization with implicit solvent models and by nonsystematic QM/MM energy minimization, as well as by nonsystematic free-energy mapping. However, no study has provided the needed reliable 2D (3D) surfaces that are necessary for reaching concrete conclusions. Here we report a systematic evaluation of the 2D (3D) free-energy maps for several relevant systems, comparing the results of QM(ai)/MM and QM(ai)/implicit solvent surfaces, and provide an advanced description of the relevant energetics. It is found that the 1W path for the hydrolysis of the methyl diphosphate (MDP) trianion is 6-9 kcal/mol higher than that the 2W path. This difference becomes slightly larger in the presence of the Mg(2+) ion because this ion reduces the pKa of the conjugated acid form of the phosphate oxygen that accepts the proton. Interestingly, the BLYP approach (which has been used extensively in some studies) gives a much smaller difference between the 1W and 2W activation barriers. At any rate, it is worth pointing out that the 2W transition state for the PT is not much higher that the common plateau that serves as the starting point of both the 1W and 2W PT paths. Thus, the calculated catalytic effects of proteins based on the 2W PT mechanistic model are not expected to be different from the catalytic effects predicted using the 1W PT mechanistic model, which was calibrated on the observed barrier in solution and in which the TS charge distribution was similar to the that of the plateau (as was done in all of our previous EVB studies).

Electrostatic contributions to binding of transition state analogues can be very different from the corresponding contributions to catalysis: phenolates binding to the oxyanion hole of ketosteroid isomerase.

The relationship between binding of transition state analogues (TSAs) and catalysis is an open problem. A recent study of the binding of phenolate TSAs to ketosteroid isomerase (KSI) found a small change in the binding energy with a change in charge delocalization of the TSAs. This has been taken as proof that electrostatic effects do not contribute in a major way to catalysis. Here we reanalyze the relationship between the binding of the TSAs and the chemical catalysis by KSI as well as the binding of the transition state (TS), by computer simulation approaches. Since the simulations reproduce the relevant experimental results, they can be used to quantify the different contributions to the observed effects. It is found that the binding of the TSA and the chemical catalysis represent different thermodynamic cycles with very different electrostatic contributions. While the binding of the TSA involves a small electrostatic contribution, the chemical catalysis involves a charge transfer process and a major electrostatic contribution due to the preorganization of the active site. Furthermore, it is found that the electrostatic preorganization contributions to the binding of the enolate intermediate of KSI and the TS are much larger than the corresponding effect for the binding of the TSAs. This reflects the dependence of the preorganization on the orientation of the nonpolar form of the TSAs relative to the oxyanion hole. It seems to us that this work provides an excellent example of the need for computational studies in analyzing key experimental findings about enzyme catalysis.

A new paradigm for electrostatic catalysis of radical reactions in vitamin B12 enzymes.

The catalytic power of enzymes containing coenzyme B(12) cofactor has been, in some respects, the "last bastion" for the strain hypothesis. The present work explores the origin of this effect by using simulation methods that overcome the sampling difficulties of previous energy minimization studies. It is found that the major part of the catalytic effect is due to the electrostatic interaction between the ribose and the protein, and that the strain contribution is very small. Remarkably, enzymes can use electrostatic effects even in a radical process, when the charge distribution of the reacting fragments does not change significantly during the reaction. Electrostatic catalysis can, in such cases, be obtained by attaching a polar group to the leaving fragment and designing an active site that interacts more strongly with this group in the product state than in the reactant state. The finding that evolution had to use this trick provides further evidence to the observation that it is extremely hard to catalyze enzymatic reactions by nonelectrostatic factors. The trick used by B(12) enzymes may, in fact, be a very powerful new strategy in enzyme design.