Absolute Free Energy Research Articles

Free Energies of Solvation with Surface, Volume, and Local Electrostatic Effects and Atomic Surface Tensions to Represent the First Solvation Shell.

Building on the SVPE (surface and volume polarization for electrostatics) model for electrostatic contributions to the free energy of solvation with explicit consideration of both surface and volume polarization effects, on the SMx approach to including first-solvation-shell contributions, and on the linear relationship between the electric field and short-range electrostatic contributions found by Chipman, we have developed a new method for computing absolute aqueous solvation free energies by combining the SVPE method with semiempirical terms that account for effects beyond bulk electrostatics. The new method is called SMVLE, and the elements it contains are denoted by SVPE-CDSL where SVPE denotes accounting for bulk electrostatic interactions between solute and solvent with both surface and volume contributions, CDS denotes the inclusion of solvent cavitation, changes in dispersion energy, and possible changes in local solvent structure by a semiempirical term utilizing geometry-dependent atomic surface tensions as implemented in SMx models, and L represents the local electrostatic effect derived from the outward-directed normal electric field on the cavity surface. The semiempirical CDS and L terms together represent the deviation of short-range contributions to the free energy of solvation from those accounted for by the SVPE term based on the bulk solvent dielectric constant. A solute training set containing a broad range of molecules used previously in the development of SM6 is used here for SMVLE model calibration. The aqueous solvation free energies predicted by the parameterized SMVLE model correlate exceedingly well with experimental values. The square of the correlation coefficient is 0.9949 and the slope is 1.0079. Comparison of the final SMVLE model against the earlier SMx solvation model shows that the parameterized SMVLE model not only yields good accuracy for neutrals but also significantly increases the accuracy for ions, making it the best implicit solvation model to date for aqueous solvation free energies of ions. The semiempirical terms associated with the outward-directed electric field account in a physical way for the improvement in the predictive accuracy for ions. The SMVLE method greatly decreases the need to include explicit water molecules for accurate modeling of solvation free energies of ions.

Computing absolute free energies of disordered structures by molecular simulation

We present a Monte Carlo simulation technique by which the free energy of disordered systems can be computed directly. It is based on thermodynamic integration. The central idea is to construct an analytically solvable reference system from a configuration which is representative for the state of interest. The method can be applied to lattice models (e.g., the Ising model) as well as off-lattice molecular models. We focus mainly on the more challenging off-lattice case. We propose a Monte Carlo algorithm, by which the thermodynamic integration path can be sampled efficiently. At the examples of the hard sphere liquid and a hard disk solid with a defect, we discuss several properties of the approach.

Absolute free energies and equilibrium ensembles of dense fluids computed from a nondynamic growth method

We demonstrate a nondynamical Monte Carlo method to compute free energies and generate equilibrium ensembles of dense fluids. In this method, based on step-by-step polymer growth algorithms, an ensemble of n+1 particles is obtained from an ensemble of n particles by generating configurations of the n+1st particle. A statistically rigorous resampling scheme is utilized to remove configurations with low weights and to avoid a combinatorial explosion; the free energy is obtained from the sum of the weights. In addition to the free energy, the method generates an equilibrium ensemble of the full system. We consider two different system sizes for a Lennard-Jones fluid and compare the results with conventional Monte Carlo methods.

Spatial updating in the great grand canonical ensemble

In spatial updating grand canonical Monte Carlo, particle transfers are implemented by examining the local environment around a point in space. In the present work, these algorithms are extended to very high densities by allowing the volume to fluctuate, thus forming a great grand canonical ensemble. Since fluctuations are unbounded, a constraint must be imposed. The constrained ensemble may be viewed as a superposition of either constant-pressure or grand canonical ensembles. Each simulation of the constrained ensemble requires a set of weights that must be determined iteratively. The outcome of a single simulation is the density of states in terms of all its independent variables. Since all extensive variables fluctuate, it is also possible to estimate absolute free energies and entropies from a single simulation. The method is tested on a system of hard spheres and the transition from the fluid to a face-centered cubic crystal is located with high precision.

Polarizable Simulations with Second-Order Interaction Model—Force Field and Software for Fast Polarizable Calculations: Parameters for Small Model Systems and Free Energy Calculations

We are presenting POSSIM (POlarizable Simulations with Second order Interaction Model) - a software package and a set of parameters designed for molecular simulations. The key feature of POSSIM is that the electrostatic polarization is taken into account using a previously introduced fast formalism. This permits cutting computational cost of using the explicit polarization by about an order of magnitude. In this article, parameters for water, methane, ethane, propane, butane, methanol and NMA are introduced. These molecules are viewed as model systems for protein simulations. We have achieved our goal of ca. 0.5 kcal/mol accuracy for gas-phase dimerization energies and no more than 2% deviations in liquid state heats of vaporization and densities. Moreover, free energies of hydration of the polarizable methane, ethane and methanol have been calculated using the statistical perturbation theory. These calculations are a model for calculating protein pKa shifts and ligand binding affinities. The free energies of hydration were found to be 2.12 kcal/mol, 1.80 kcal/mol and -4.95 kcal/mol for methane, ethane and methanol, respectively. The experimentally determined literature values are 1.91 kcal/mol, 1.83 kcal/mol and -5.11 kcal/mol. The POSSIM average error in these absolute free energies of hydration is only about 0.13 kcal/mol. Using the statistical perturbation theory with polarizable force fields is not widespread, and we believe that this work opens road to further development of the POSSIM force field and its applications for obtaining accurate energies in protein-related computer modeling.

Computation of Absolute Hydration and Binding Free Energy with Free Energy Perturbation Distributed Replica-Exchange Molecular Dynamics

Distributed Replica (REPDSTR) is a powerful parallelization technique enabling simulations of a group of replicas in a parallel/parallel fashion, where each replica is distributed to different nodes of a large cluster [Theor. Chem. Acc. 109: 140 (2003)]. Here, we use the framework provided by REPDSTR to combine a staged free energy perturbation protocol with different values of the thermodynamic coupling parameters with replica-exchange molecular dynamics (FEP/REMD. The structure of REPDSTR, which allows multiple parallel input/output (I/O), facilitates the treatment of replica-exchange to couple the N window simulations. As a result, each of the N synchronous window simulations benefit from the sampling carried out by the N-1 others. As illustrative examples of the FEP/REMD strategy, calculations of the absolute hydration and binding free energy of small molecules were performed using the biomolecular simulation program CHARMM adapted for the IBM Blue Gene/P platform. The computations show that a FEP/REMD strategy significantly improves the sampling and accelerate the convergence of absolute free energy computations.

1-Octanol/Water Partition Coefficients of n-Alkanes from Molecular Simulations of Absolute Solvation Free Energies.

The 1-octanol/water partition coefficient is an important thermodynamic variable usually employed to understand and quantify the partitioning of solutes between aqueous and organic phases. It finds widespread use in many empirical correlations to evaluate the environmental fate of pollutants as well as in the design of pharmaceuticals. The experimental evaluation of 1-octanol/water partition coefficients is an expensive and time-consuming procedure, and thus, theoretical estimation methods are needed, particularly when a physical sample of the solute may not yet be available, such as in pharmaceutical screening. 1-Octanol/water partition coefficients can be obtained from Gibbs free energies of solvation of the solute in both the aqueous and the octanol phases. The accurate evaluation of free energy differences remains today a challenging problem in computational chemistry. In order to study the absolute solvation Gibbs free energies in 1-octanol, a solvent that can mimic many properties of important biological systems, free energy calculations for n-alkanes in the range C1-C8 were performed using molecular simulation techniques, following the thermodynamic integration approach. In the first part of this paper, we test different force fields by evaluating their performance in reproducing pure 1-octanol properties. It is concluded that all-atom force fields can provide good accuracy but at the cost of a higher computational time compared to that of the united-atom force fields. Recent versions of united-atom force fields, such as Gromos and TraPPE, provide satisfactory results and are, thus, useful alternatives to the more expensive all-atom models. In the second part of the paper, the Gibbs free energy of solvation in 1-octanol is calculated for several n-alkanes using three force fields to describe the solutes, namely Gromos, TraPPE, and OPLS-AA. Generally, the results obtained are in excellent agreement with the available experimental data and are of similar accuracy to commonly used QSPR models. Moreover, we have estimated the Gibbs free energy of hydration for the different compounds with the three force fields, reaching average deviations from experimental data of less than 0.2 kcal/mol for the case of the Gromos force field. Finally, we systematically compare different strategies to obtain the 1-octanol/water partition coefficient from the simulations. It is shown that a fully predictive method combining the Gromos force field in the aqueous phase and the OPLS-AA/TraPPE force field for the organic phase can give excellent predictions for n-alkanes up to C8 with an absolute average deviation of 0.1 log P units to the experimental data.

Ab initio derivation of a dataset of real temperature thermodynamic properties: Case study with SiC

A wholly ab initio method for deriving chemical thermodynamic data for inorganic crystals is a useful tool for theoretical research in material science. In this work, all electron density functional computations together with phonon vibrational analysis were used to derive thermodynamic profiles of two polytypes of SiC (3C and 2H) at real temperatures. This included standard Gibbs free energy of formation () and enthalpy of formation () values from 0 to 1000 K. Ground state properties were computed by the full potential linearized augmented plane wave method. Temperature effects were obtained by direct method computation of harmonic vibrations. To obtain and similar computations were also carried out for the reference state elements, Si (diamond) and C (graphite). Derived values include entropy, specific heat, enthalpy increment with temperature, absolute total enthalpy, internal energy and free energy, standard enthalpy of formation and standard free energy of formation for temperatures from 0 to 1000 K. Computed values are in very good agreement with the experimental data. This study demonstrates the utility of ab initio methods for deriving a thermodynamic data set.

Methods for calculating the absolute entropy and free energy of biological systems based on ideas from polymer physics.

The commonly used simulation techniques, Metropolis Monte Carlo (MC) and molecular dynamics (MD) are of a dynamical type which enables one to sample system configurations i correctly with the Boltzmann probability, P(i)(B), while the value of P(i)(B) is not provided directly; therefore, it is difficult to obtain the absolute entropy, S approximately -ln P(i)(B), and the Helmholtz free energy, F. With a different simulation approach developed in polymer physics, a chain is grown step-by-step with transition probabilities (TPs), and thus their product is the value of the construction probability; therefore, the entropy is known. Because all exact simulation methods are equivalent, i.e. they lead to the same averages and fluctuations of physical properties, one can treat an MC or MD sample as if its members have rather been generated step-by-step. Thus, each configuration i of the sample can be reconstructed (from nothing) by calculating the TPs with which it could have been constructed. This idea applies also to bulk systems such as fluids or magnets. This approach has led earlier to the "local states" (LS) and the "hypothetical scanning" (HS) methods, which are approximate in nature. A recent development is the hypothetical scanning Monte Carlo (HSMC) (or molecular dynamics, HSMD) method which is based on stochastic TPs where all interactions are taken into account. In this respect, HSMC(D) can be viewed as exact and the only approximation involved is due to insufficient MC(MD) sampling for calculating the TPs. The validity of HSMC has been established by applying it first to liquid argon, TIP3P water, self-avoiding walks (SAW), and polyglycine models, where the results for F were found to agree with those obtained by other methods. Subsequently, HSMD was applied to mobile loops of the enzymes porcine pancreatic alpha-amylase and acetylcholinesterase in explicit water, where the difference in F between the bound and free states of the loop was calculated. Currently, HSMD is being extended for calculating the absolute and relative free energies of ligand-enzyme binding. We describe the whole approach and discuss future directions.

Absolute Hydration Entropies of Alkali Metal Ions from Molecular Dynamics Simulations

Molecular dynamics simulations in combination with the free energy perturbation technique are used in this work to calculate absolute ion hydration entropies. The hydration entropies for five alkali metal ions are estimated from van't Hoff plots using hydration free energies calculated at eight different temperatures. Considering that the ion-water potentials were parametrized only on absolute hydration free energies and ionic radii, the absolute hydration entropies agree very well with experimental data. Simulation lengths of about 3 ns at each temperature were required to achieve an uncertainty below 1 kcal/mol for the entropic contribution to the hydration free energy (-TDeltaShyd). The uncertainties for the calculated entropies are typically four times larger than for the free energies. The possibility to use approximate approaches to calculate hydration entropies is also investigated. The entropy of creating the uncharged van der Waals spheres in water correlates well with the solvent accessible surface area of the ions. The Born continuum model and the linear response approximation cannot be used to predict the entropy of charging the van der Waals spheres in water without introducing temperature dependent empirical parameters.

Calculation of free-energy differences by confinement simulations. Application to peptide conformers.

Conformational free-energy differences are key quantities for understanding important phenomena in molecular biology that involve large structural changes of macromolecules. In this paper, an improved version of the confinement approach, which is based on earlier developments, is used to determine the free energy of individual molecular states by progressively restraining the corresponding molecular structures to pure harmonic basins, whose absolute free energy can be computed by normal-mode analysis. The method is used to calculate the free-energy difference between two conformational states of the alanine dipeptide in vacuo, and of the beta-hairpin from protein G with an implicit solvation model. In all cases, the confinement results are in excellent agreement with the ones obtained from converged equilibrium molecular dynamics simulations, which have a much larger computational cost. The systematic and statistical errors of the results are evaluated and the origin of the errors is identified. The sensitivity of the calculated free-energy differences to structure-based definitions of the molecular states is discussed. A variant of the method, which closes the thermodynamic cycle by a quasi-harmonic rather than harmonic analysis, is introduced. The latter is proposed for possible use with explicit solvent simulations.

Absolute free energies estimated by combining precalculated molecular fragment libraries

The absolute free energy--or partition function, equivalently--of a molecule can be estimated computationally using a suitable reference system. Here, we demonstrate a practical method for staging such calculations by growing a molecule based on a series of fragments. Significant computer time is saved by precalculating fragment configurations and interactions for reuse in a variety of molecules. We use such fragment libraries and interaction tables for amino acids and capping groups to estimate free energies for small peptides. Equilibrium ensembles for the molecules are generated at no additional computational cost and are used to check our results by comparison to standard dynamics simulation. We explain how our work can be extended to estimate relative binding affinities.

Ab initio molecular dynamics calculations of ion hydration free energies.

We apply ab initio molecular dynamics (AIMD) methods in conjunction with the thermodynamic integration or "lambda-path" technique to compute the intrinsic hydration free energies of Li(+), Cl(-), and Ag(+) ions. Using the Perdew-Burke-Ernzerhof functional, adapting methods developed for classical force field applications, and with consistent assumptions about surface potential (phi) contributions, we obtain absolute AIMD hydration free energies (DeltaG(hyd)) within a few kcal/mol, or better than 4%, of Tissandier et al.'s [J. Phys. Chem. A 102, 7787 (1998)] experimental values augmented with the SPC/E water model phi predictions. The sums of Li(+)/Cl(-) and Ag(+)/Cl(-) AIMD DeltaG(hyd), which are not affected by surface potentials, are within 2.6% and 1.2 % of experimental values, respectively. We also report the free energy changes associated with the transition metal ion redox reaction Ag(+)+Ni(+)-->Ag+Ni(2+) in water. The predictions for this reaction suggest that existing estimates of DeltaG(hyd) for unstable radiolysis intermediates such as Ni(+) may need to be extensively revised.

Absolute free energy and entropy of a mobile loop of the enzyme acetylcholinesterase.

The loop 287-290 (Ile, Phe, Arg, and Phe) of the protein acetylcholinesterase (AChE) changes its structure upon interaction of AChE with diisopropylphosphorofluoridate (DFP). Reversible dissociation measurements suggest that the free-energy (F) penalty for the loop displacement is DeltaF=Ffree-Fbound approximately -4 kcal/mol. Therefore, this loop has been the target of two studies by Olson's group for testing the efficiency of procedures for calculating F. In this paper, we test for the first time the performance of our "hypothetical scanning molecular dynamics" (HSMD) method and the validity of the related modeling for a loop with bulky side chains in explicit water. Thus, we consider only atoms of the protein that are the closest to the loop (they constitute the "template"), where the rest of the atoms are ignored. The template's atoms are fixed in the X-ray coordinates of the free protein, and the loop is capped with a sphere of TIP3P water molecules; also, the X-ray structure of the bound loop is attached to the free template. We carry out two separate MD simulations starting from the free and bound X-ray structures, where only the atoms of the loop and water are allowed to move while the template-water and template-loop (AMBER) interactions are considered. The absolute Ffree and Fbound (of the loop and water) are calculated from the corresponding trajectories. A main objective of this paper is to assess the reliability of this model, and for this several template sizes are studied capped with 80-220 water molecules. We find that consistent results for the free energy (which also agree with the experimental data above) require a template larger than a minimal size and a number of water molecules approximately equal to the experimental density of bulk water. For example, we obtain DeltaFtotal=DeltaFwater+DeltaFloop=-3.1+/-2.5 and -3.6+/-4 kcal/mol for a template consisting of 944 atoms and a sphere containing 160 and 180 waters, respectively. Our calculations demonstrate the important contribution of water to the total free energy. Namely, for water densities close to the experimental value, DeltaFwater is always negative leading thereby to a negative DeltaFtotal (while DeltaFloop is always positive). Also, the contribution of the water entropy TDeltaSwater to DeltaFtotal is significant. Various aspects related to the efficiency of HSMD are tested and improved, and plans for future studies are discussed.

Computations of Absolute Solvation Free Energies of Small Molecules Using Explicit and Implicit Solvent Model.

Accurate determination of absolute solvation free energy plays a critical role in numerous areas of biomolecular modeling and drug discovery. A quantitative representation of ligand and receptor desolvation, in particular, is an essential component of current docking and scoring methods. Furthermore, the partitioning of a drug between aqueous and nonpolar solvents is one of the important factors considered in pharmacokinetics. In this study, the absolute hydration free energy for a set of 239 neutral ligands spanning diverse chemical functional groups commonly found in drugs and drug-like candidates is calculated using the molecular dynamics free energy perturbation method (FEP/MD) with explicit water molecules, and compared to experimental data as well as its counterparts obtained using implicit solvent models. The hydration free energies are calculated from explicit solvent simulations using a staged FEP procedure permitting a separation of the total free energy into polar and nonpolar contributions. The nonpolar component is further decomposed into attractive (dispersive) and repulsive (cavity) components using the Weeks-Chandler-Anderson (WCA) separation scheme. To increase the computational efficiency, all of the FEP/MD simulations are generated using a mixed explicit/implicit solvent scheme with a relatively small number of explicit TIP3P water molecules, in which the influence of the remaining bulk is incorporated via the spherical solvent boundary potential (SSBP). The performances of two fixed-charge force fields designed for small organic molecules, the General Amber force field (GAFF), and the all-atom CHARMm-MSI, are compared. Because of the crucial role of electrostatics in solvation free energy, the results from various commonly used charge generation models based on the semiempirical (AM1-BCC) and QM calculations [charge fitting using ChelpG and RESP] are compared. In addition, the solvation free energies of the test set are also calculated using Poisson-Boltzmann (PB) and Generalized Born model of solvation (GB), which are two widely used continuum electrostatic implicit solvent models. The protocol for running the absolute solvation free energy calculations used throughout is automated as much as possible, with minimum user intervention, so that it can be used in large-scale analysis and force field optimization.

A theoretical study on the mechanism of the base-promoted decomposition of N-chloro,N-methylethanolamine

The first step of the base-promoted decomposition of N-chloro,N-methylethanolamine in aqueous solution (CH3N(Cl)CH2CH2OH + HO- --> imine + Cl- + H2O (+ CH2O) --> amine + aldehyde) is investigated at the MP2/6-31++G(d,p) computing level. Solvation is included by using both a microsolvated model, in which two explicit water molecules simulate the specific solvent effects, and a hybrid cluster-continuum model, by applying a polarized continuum on the previous results, to account for the bulk effect of the solvent. Four alternative pathways (bimolecular fragmentation, Hofmann, Zaitsev and intramolecular eliminations) are possible for the rate-limiting step of this base-promoted decomposition. These reactive processes are bimolecular asynchronous concerted reactions. The common feature of the four pathways is the proton transfer to HO- being more advanced than all other molecular events, whereas imine formation is delayed. Non-reactive cyclic arrangements involving one of the explicit water molecules are found at transition structures of Hofmann and Zaitsev eliminations, such water molecule acting both as H+ donor and acceptor. Although MP2 calculations misjudge the absolute activation Gibbs free energy values, this computational level adequately predicts the enhancement in the decomposition rate due to the presence of the -OH group.

A unified hydrophobicity scale for multispan membrane proteins

The concept of hydrophobicity is critical to our understanding of the principles of membrane protein (MP) folding, structure, and function. In the last decades, several groups have derived hydrophobicity scales using both experimental and statistical methods that are optimized to mimic certain natural phenomena as closely as possible. The present work adds to this toolset the first knowledge-based scale that unifies the characteristics of both alpha-helical and beta-barrel multispan MPs. This unified hydrophobicity scale (UHS) distinguishes between amino acid preference for solution, transition, and trans-membrane states. The scale represents average hydrophobicity values of amino acids in folded proteins, irrespective of their secondary structure type. We furthermore present the first knowledge-based hydrophobicity scale for mammalian alpha-helical MPs (mammalian hydrophobicity scale--MHS). Both scales are particularly useful for computational protein structure elucidation, for example as input for machine learning techniques, such as secondary structure or trans-membrane span prediction, or as reference energies for protein structure prediction or protein design. The knowledge-based UHS shows a striking similarity to a recent experimental hydrophobicity scale introduced by Hessa and coworkers (Hessa T et al., Nature 2007;450:U1026-U1032). Convergence of two very different approaches onto similar hydrophobicity values consolidates the major differences between experimental and knowledge-based scales observed in earlier studies. Moreover, the UHS scale represents an accurate absolute free energy measure for folded, multispan MPs--a feature that is absent from many existing scales. The utility of the UHS was demonstrated by analyzing a series of diverse MPs. It is further shown that the UHS outperforms nine established hydrophobicity scales in predicting trans-membrane spans along the protein sequence. The accuracy of the present hydrophobicity scale profits from the doubling of the number of integral MPs in the PDB over the past four years. The UHS paves the way for an increased accuracy in the prediction of trans-membrane spans.

Quantum chemical identification of blue and red forms of protonated pheophytin- a dianion

Quantum mechanics/molecular mechanics (QM/MM) calculations were employed so as to determine the most stable structures of protonated forms of Pheophytin- a (Pheo) dianion in DMF. The protonated forms include PheoH − and PheoH 2. Electrostatic potential-derived charges (ESP charges) employing the CHelpG scheme were obtained for Pheo 2− and for possible blue and red forms of PheoH − by ROB3LYP methodology using 6-31G(d) basis. The ESP charges and the most stable structures obtained by QM/MM calculations provide a clue to the probable site of protonation in Pheo 2− and PheoH −. A total of 18 different possible structures were investigated. From these calculations, we identify the structures of the blue and red forms of PheoH − and PheoH 2. Furthermore, density functional treatment (DFT) along with dielectric polarizable continuum model (DPCM) calculations using 6-311+G(2d,2p) basis were done to determine the absolute free energy of reduction of Pheo to PheoH − and PheoH 2 in DMF. The calculated mid-point reduction potentials relative to the standard hydrogen electrode are −0.63 V (blue), −0.71 V (red) for reduction of Pheo to PheoH −, and 0.15 V (blue), 0.11 V (red) for reduction to PheoH 2. The calculated numbers are in good agreement with observed values, −0.68, −0.72, 0.15 and 0.05 V, respectively.

Modeling molecular and ionic absolute solvation free energies with quasichemical theory bounds

A recently developed statistical mechanical quasichemical theory (QCT) has led to significant insights into solvation phenomena for both hydrophilic and hydrophobic solutes. The QCT exactly partitions solvation free energies into three components: (1) Inner-shell chemical, (2) outer-shell packing, and (3) outer-shell long-ranged contributions. In this paper, we discuss efficient methods for computing each of the three parts of the free energy. A Bayesian estimation approach is developed to compute the inner-shell chemical and outer-shell packing contributions. We derive upper and lower bounds on the outer-shell long-ranged portion of the free energy by expressing this component in two equivalent ways. Local, high-energy contacts between the solute and solvent are eliminated by spatial conditioning in this free energy piece, leading to near-Gaussian distributions of solute-solvent interaction energies. Thus, the average of the two mean-field bounds yields an accurate and efficient free energy estimate. Aqueous solvation free energy results are presented for several solutes, including methane, perfluoromethane, water, and sodium and chloride ions. The results demonstrate the accuracy and efficiency of the methods. The approach should prove useful in computing solvation free energies in inhomogeneous, restricted environments.

Theoretical determination of absolute free energy of reduction of plastoquinone-9 in photosystem II and of plastoquinone- n in DMF

We employed static continuum electrostatics and multi-conformation continuum electrostatics (MCCE) methods to determine the reduction potential ( E red 0 ) of PQ-9 in a section of Photosystem II (PSII). Both methods relied on the finite difference Poisson–Boltzmann (FDPB) solution. The static method brings out a E red 0 value (0.01 V) that is close to the experimental one (0.05 V), thereby demonstrating that the surrounding environment critically decides the net free energy change. The E red 0 value obtained from MCCE (0.04 V) is even closer to the observed value, thereby indicating the importance of protein side-chain and proton motions in the electron transfer process. Furthermore, density functional theory-dielectric polarisable continuum model (DFT-DPCM) was employed to calculate the absolute free energy of reduction of plastoquinone- n (PQ- n, where n is the number of isoprenoid units) in N, N dimethyl formamide (DMF) solvent. The DFT-DPCM method produced reduction potential values of −0.59 and −0.65 V for PQ-1 and PQ-9, respectively. These are more or less in agreement with the experimentally reported values of −0.64 and −0.62 V, respectively.