Minimalist explicit solvation models for surface loops in proteins

Open and Closed Form of Maltose Binding Protein in Its Native and Molten Globule State As Studied by Electron Paramagnetic Resonance Spectroscopy

An intensively investigated intermediate state of protein folding is the molten globule (MG) state, which contains secondary but hardly any tertiary structure. In previous work, we have determined the distances between interacting spins within maltose binding protein (MBP) in its native state using continuous wave and double electron-electron resonance (DEER) electron paramagnetic resonance (EPR) spectroscopy. Seven double mutants had been employed to investigate the structure within the two domains of MBP. DEER data nicely corroborated the previously available X-ray data. Even in its MG state, MBP is known to still bind its ligand maltose. We therefore hypothesized that there must be a defined structure around the binding pocket of MBP already in the absence of tertiary structure. Here we have investigated the functional and structural difference between native and MG state in the open and closed form with a new set of MBP mutants. In these, the spin-label positions were placed near the active site. Binding of its ligands leads to a conformational change from open to closed state, where the two domains are more closely together. The complete set of MBP mutants was analyzed at pH 3.2 (MG) and pH 7.4 (native state) using double-quantum coherence EPR. The values were compared with theoretical predictions of distances between the labels in biradicals constructed by molecular modeling from the crystal structures of MBP in open and closed form and were found to be in excellent agreement. Measurements show a defined structure around the binding pocket of MBP in MG, which explains maltose binding. A new and important finding is that in both states ligand-free MBP can be found in open and closed form, while ligand-bound MBP appears only in closed form because of maltose binding.

Conformational energy range of ligands in protein crystal structures: The difficult quest for accurate understanding

In this review, we address a fundamental question: What is the range of conformational energies seen in ligands in protein-ligand crystal structures? This value is important biophysically, for better understanding the protein-ligand binding process; and practically, for providing a parameter to be used in many computational drug design methods such as docking and pharmacophore searches. We synthesize a selection of previously reported conflicting results from computational studies of this issue and conclude that high ligand conformational energies really are present in some crystal structures. The main source of disagreement between different analyses appears to be due to divergent treatments of electrostatics and solvation. At the same time, however, for many ligands, a high conformational energy is in error, due to either crystal structure inaccuracies or incorrect determination of the reference state. Aside from simple chemistry mistakes, we argue that crystal structure error may mainly be because of the heuristic weighting of ligand stereochemical restraints relative to the fit of the structure to the electron density. This problem cannot be fixed with improvements to electron density fitting or with simple ligand geometry checks, though better metrics are needed for evaluating ligand and binding site chemistry in addition to geometry during structure refinement. The ultimate solution for accurately determining ligand conformational energies lies in ultrahigh-resolution crystal structures that can be refined without restraints.

Relative stability of the open and closed conformations of the active site loop of streptavidin

The eight-residue surface loop, 45-52 (Ser, Ala, Val, Gly, Asn, Ala, Glu, Ser), of the homotetrameric protein streptavidin has a "closed" conformation in the streptavidin-biotin complex, where the corresponding binding affinity is one of the strongest found in nature (ΔG ∼ -18 kcal∕mol). However, in most of the crystal structures of apo (unbound) streptavidin, the loop conformation is "open" and typically exhibits partial disorder and high B-factors. Thus, it is plausible to assume that the loop structure is changed from open to closed upon binding of biotin, and the corresponding difference in free energy, ΔF = F(open) - F(closed) in the unbound protein, should therefore be considered in the total absolute free energy of binding. ΔF (which has generally been neglected) is calculated here using our "hypothetical scanning molecular-dynamics" (HSMD) method. We use a protein model in which only the atoms closest to the loop are considered (the "template") and they are fixed in the x-ray coordinates of the free protein; the x-ray conformation of the closed loop is attached to the same (unbound) template and both systems are capped with the same sphere of TIP3P water. Using the force field of the assisted model building with energy refinement (AMBER), we carry out two separate MD simulations (at temperature T = 300 K), starting from the open and closed conformations, where only the atoms of the loop and water are allowed to move (the template-water and template-loop interactions are considered). The absolute F(open) and F(closed) (of loop + water) are calculated from these trajectories, where the loop and water contributions are obtained by HSMD and a thermodynamic integration (TI) process, respectively. The combined HSMD-TI procedure leads to total (loop + water) ΔF = -27.1 ± 2.0 kcal∕mol, where the entropy TΔS constitutes 34% of ΔF, meaning that the effect of S is significant and should not be ignored. Also, ΔS is positive, in accord with the high flexibility of the open loop observed in crystal structures, while the energy ΔE is unexpectedly negative, thus also adding to the stability of the open loop. The loop and the 250 capped water molecules are the largest system studied thus far, which constitutes a test for the efficiency of HSMD-TI; this efficiency and technical issues related to the implementation of the method are also discussed. Finally, the result for ΔF is a prediction that will be considered in the calculation of the absolute free energy of binding of biotin to streptavidin, which constitutes our next project.

Entropy and Free Energy of a Mobile Loop Based on the Crystal Structures of the Free and Bound Proteins

A mobile loop changes its conformation from "open" (free enzyme) to "closed" upon ligand binding. The difference in the Helmholtz free energy, ΔF(loop) between these states sheds light on the mechanism of binding. With our "hypothetical scanning molecular dynamics" (HSMD-TI) method ΔF(loop) = F(free) - F(bound) where F(free) and F(bound) are calculated from two MD samples of the free and bound loop states; the contribution of water is obtained by a thermodynamic integration (TI) procedure. In previous work the free and bound loop structures were both attached to the same "template" which was "cut" from the crystal structure of the free protein. Our results for loop 287-290 of AcetylCholineEsterase agree with the experiment, ΔF(loop)~ -4 kcal/mol if the density of the TIP3P water molecules capping the loop is close to that of bulk water, i.e., N(water) = 140 - 180 waters in a sphere of a 18 Å radius. Here we calculate ΔF(loop) for the more realistic case, where two templates are "cut" from the crystal structures, 2dfp.pdb (bound) and 2ace.pdb (free), where N(water) = 40 - 160; this requires adding a computationally more demanding (second) TI procedure. While the results for N(water) ≤ 140 are computationally sound, ΔF(loop) is always positive (18 ± 2 kcal/mol for N(water) = 140). These (disagreeing) results are attributed to the large average B-factor, 41.6 of 2dfp (23.4 Å(2) for 2ace). While this conformational uncertainty is an inherent difficulty, the (unstable) results for N(water) = 160 suggest that it might be alleviated by applying different (initial) structural optimizations to each template.

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.

Entropy and free energy of a mobile protein loop in explicit water

Estimation of the energy from a given Boltzmann sample is straightforward since one just has to average the contribution of the individual configurations. On the other hand, calculation of the absolute entropy, S (hence the absolute free energy F) is difficult because it depends on the entire (unknown) ensemble. We have developed a new method called "the hypothetical scanning molecular dynamics" (HSMD) for calculating the absolute S from a given sample (generated by any simulation technique). In other words, S (like the energy) is "written" on the sample configurations, where HSMD provides a prescription of how to "read" it. In practice, each sample conformation, i, is reconstructed with transition probabilities, and their product leads to the probability of i, hence to the entropy. HSMD is an exact method where all interactions are considered, and the only approximation is due to insufficient sampling. In previous studies HSMD (and HS Monte CarloHSMC) has been extended systematically to systems of increasing complexity, where the most recent is the seven-residue mobile loop, 304-310 (Gly-His-Gly-Ala-Gly-Gly-Ser) of the enzyme porcine pancreatic alpha-amylase modeled by the AMBER force field and AMBER with the implicit solvation GB/SA (paper I, Cheluvaraja, S.; Meirovitch, H. J. Chem. Theory Comput. 2008, 4, 192). In the present paper we make a step further and extend HSMD to the same loop capped with TIP3P explicit water at 300 K. As in paper I, we are mainly interested in entropy and free energy differences between the free and bound microstates of the loop, which are obtained from two separate MD samples of these microstates. The contribution of the loop to S and F is calculated by HSMD and that of water by a particular thermodynamic integration procedure. As expected, the free microstate is more stable than the bound microstate by a total free energy difference, Ffree-Fbound=-4.8+/-1, as compared to -25.5 kcal/mol obtained with GB/SA. We find that relatively large systematic errors in the loop entropies, Sfree(loop) and Sbound(loop) are cancelled in their difference which is thus obtained efficiently and with high accuracy, i.e., with a statistical error of 0.1 kcal/mol. This cancellation, which has been observed in previous HSMD studies, is in accord with theoretical arguments given in paper I.