G protein coupled receptors -in silico drug discovery and design.

In silico drug discovery is a complex process requiring flexibility and ingenuity in method selection and a careful validation of work protocols. GPCR in silico drug discovery poses additional challenges due to the paucity of crystallographic data. This paper starts by reviewing selected GPCR in silico screening programs reported in the literature, including both structure-based and ligand-based approaches. Particular emphasis is given to library design, binding mode selection, process validation and compound selection for biological testing. Following literature review, we provide insights into in silico methodologies and process workflows used at EPIX to drive over 20 highly successful screening and lead optimization programs performed since 2001. Applications of the various methodologies discussed are demonstrated by examples from recent programs that have not yet been published.

How much do enzymes really gain by restraining their reacting fragments?

The steric effect, exerted by enzymes on their reacting substrates, has been considered as a major factor in enzyme catalysis. In particular, it has been proposed that enzymes catalyze their reactions by pushing their reacting fragments to a catalytic configuration which is sometimes called near attack configuration (NAC). This work uses computer simulation approaches to determine the relative importance of the steric contribution to enzyme catalysis. The steric proposal is expressed in terms of well defined thermodynamic cycles that compare the reaction in the enzyme to the corresponding reaction in water. The S(N)2 reaction of haloalkane dehalogenase from Xanthobacter autotrophicus GJ10, which was used in previous studies to support the strain concept is chosen as a test case for this proposal. The empirical valence bond (EVB) method provides the reaction potential surfaces in our studies. The reliability and efficiency of this method make it possible to obtain stable results for the steric free energy. Two independent strategies are used to evaluate the actual magnitude of the steric effect. The first applies restraints on the substrate coordinates in water in a way that mimics the steric effect of the protein active site. These restraints are then released and the free energy associated with the release process provides the desired estimate of the steric effect. The second approach eliminates the electrostatic interactions between the substrate and the surrounding in the enzyme and in water, and compares the corresponding reaction profiles. The difference between the resulting profiles provides a direct estimate of the nonelectrostatic contribution to catalysis and the corresponding steric effect. It is found that the nonelectrostatic contribution is about -0.7 kcal/mol while the full "apparent steric contribution" is about -2.2 kcal/mol. The apparent steric effect includes about -1.5 kcal/mol electrostatic contribution. The total electrostatic contribution is found to account for almost all the observed catalytic effect ( approximately -6.1 kcal/mol of the -6.8 calculated total catalytic effect). Thus, it is concluded that the steric effect is not the major source of the catalytic power of haloalkane dehalogenase. Furthermore, it is found that the largest component of the apparent steric effect is associated with the solvent reorganization energy. This solvent-induced effect is quite different from the traditional picture of balance between the repulsive interaction of the reactive fragments and the steric force of the protein.

Remarkable rate enhancement of orotidine 5'-monophosphate decarboxylase is due to transition-state stabilization rather than to ground-state destabilization.

The remarkable rate enhancement of orotidine 5'-phosphate decarboxylase (ODCase) has been attributed to ground-state destabilization (GSD) by desolvation and more recently to GSD by electrostatic stress. Here we reiterate our previous arguments that the GSD mechanisms are not likely to play a major role in enzyme catalysis and analyze quantitatively the origin of the rate enhancement of ODCase. This analysis involves energy considerations and computer simulations. Our energy considerations show that (i) the previously proposed desolvation mechanism is based on an improper reference state; (ii) a nonpolar active site cannot account for the catalytic effect of the enzyme; (iii) the focus on the role of the negatively charged protein residues in the electrostatic stress GSD mechanism overlooks the fact that the positively charged Lys72 strongly stabilizes the substrate; (iv) although the previous calculation of the actual enzymatic reaction correctly reproduced the observed rate enhancement, it could not obtain this rate enhancement from the calculated binding energies (which are the relevant quantities for determining GSD effects); (v) the GSD mechanism is inconsistent with the observed binding energy of the phosphoribosyl part of the substrate; and (vi) the presumably unstable substrate (orotate) can be stabilized, at equilibrium, by accepting a proton from the solvent. Our computer simulation studies involve two set of calculations. First, we study the catalytic reaction by using an empirical valence bond potential surface calibrated by ab initio calculations of the reference solution reaction. This calculation reproduces the observed catalytic effect of the enzyme. Next, we use free-energy perturbation calculations and evaluated the electrostatic contributions to the binding energies of the ground state and transition state (TS). These calculations show that the rate enhancement in ODCase is due to the TS stabilization rather than to GSD. The differences between our own and the previous theoretical analyses stem from both the selection of the reacting system and the treatment of the long-range electrostatic contributions to the binding energy. The reacting system was previously assumed to encompass only the orotate. However, this selection does not allow proper description of the reaction catalyzed by the enzyme (i.e., [Orotate(-) + LysH(+)] if [uracil + Lys + CO(2)]). Therefore, the reacting system should include both orotate and the general acid in the form of the protonated Lys72 protein residue. This selection leads to a simple and consistent interpretation of the catalytic effect where the electrostatic stabilization of the transition state is due to the fact that the two negatively charged aspartic residues are already placed near the reactive lysine so that they do not have to reorganize significantly during the reaction. Interestingly, even calculations with only orotate(-) as the reacting system do not produce sufficient destabilization to account for a GSD mechanism. In summary, we conclude, in agreement with previous workers, that ODCase catalyzes its reaction by electrostatic effects. However, we show that these effects are associated with TS stabilization due to a reduction in the protein-protein reorganization energy and not with protein-substrate destabilization effects.

How important are entropic contributions to enzyme catalysis?

The idea that enzymes accelerate their reactions by entropic effects has played a major role in many prominent proposals about the origin of enzyme catalysis. This idea implies that the binding to an enzyme active site freezes the motion of the reacting fragments and eliminates their entropic contributions, (delta S(cat)(double dagger))', to the activation energy. It is also implied that the binding entropy is equal to the activation entropy, (delta S(w)(double dagger))', of the corresponding solution reaction. It is, however, difficult to examine this idea by experimental approaches. The present paper defines the entropic proposal in a rigorous way and develops a computer simulation approach that determines (delta S(double dagger))'. This approach allows us to evaluate the differences between (delta S(double dagger))' of an enzymatic reaction and of the corresponding reference reaction in solution. Our approach is used in a study of the entropic contribution to the catalytic reaction of subtilisin. It is found that this contribution is much smaller than previously thought. This result is due to the following: (i) Many of the motions that are free in the reactants state of the reference solution reaction are also free at the transition state. (ii) The binding to the enzyme does not completely freeze the motion of the reacting fragments so that (delta S(double dagger))' in the enzymes is not zero. (iii) The binding entropy is not necessarily equal to (delta S(w)(double dagger))'.

Ab initio investigation of the molecular structure of methyl methoxymethyl phosphonate, a promising nuclease-resistant alternative of the phosphodiester linkage.

Conformational flexibility of the methyl methoxymethyl phosphonate anion (CH3-O-PO2-CH2-O-CH3)-, a nuclease resistant alternative to the phosphodiester linkage in DNA, have been investigated by ab initio quantum mechanical calculations. The potential of backbone torsional degrees of freedom of methyl methoxymethyl phosphonate anion (MMP) was determined at the Hartree-Fock (HF) 3-21G* level using the adiabatic mapping technique. Energies, geometries, and effective atomic charges of different conformers were calculated at HF/6-31G* and MP2/6-31G* levels of theory. These were compared to the results obtained for dimethyl phosphate calculated at the same level. The impact on DNA structure from inserting a methylene group between phosphorus and oxygen of the nucleoside sugar moiety was examined via distance and angle-constrained geometry optimizations. Due to its high flexibility, MMP has been shown to be compatible with both A and B forms of DNA.