Can MM-PBSA calculations predict the specificities of protein kinase inhibitors?

An application of the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) protocol to the prediction of protein kinase inhibitor selectivity is presented. Six different inhibitors are placed in equivalent orientations in each of six different receptors. Fully solvated molecular dynamics is then run for 1 ns on each of the 36 complexes, and the resulting trajectories scored, using the implicit solvent model. The results show some correlation with experimentally-determined specificities; anomalies may be attributed to a variety of causes, including difficulties in quantifying induced fit penalties and variabilities in normal modes calculations. Decomposing interaction energies on a per-residue basis yields more useful insights into the natures of the binding modes and suggests that the real value of such calculations lies in understanding interactions rather than outright prediction.

Macroscopic pKa Calculations for Fluorescein and Its Derivatives.

This study describes the calculation of the microscopic dissociation and tautomerization constants of fluorescein and its derivatives, 2',7'-dichlorofluorescein (DCF) and 2',7'-difluorofluorescein (DFF), in an aqueous environment. In vacuo free energies were obtained using complete basis set (CBS) and DFT-based methods, while free energies of solvation were calculated with the CPCM implicit solvation protocol using the UAHF, UAKS, and Pauling radii sets. Our results indicate that the different vacuum protocols give free energy changes upon dissociation within 1 kcal/mol of each other for a given molecule. Therefore, we suggest that the computationally less intensive PBE1PBE/6-311+G(2d,2p)//PBE1PBE/6-31+G(d) model chemistry may reasonably be used in pKa calculations of larger molecules. The calculations also provided a rigorous test of the implicit solvation models. Relative calculations of dissociation constants gave results in good agreement with experiment; absolute values deviated from experimental data by 1-3 pKa units. Consistently better results were obtained with the Pauling radii set. The influence of geometry relaxation on going from vacuum to solvent is negligible for pKa2 and larger for pKa1 but still smaller than the variation due to the radii set. Calculation of tautomerization constants gave more variable results, with none of the solvation methods able to reproduce experimental values consistently, although certain individual constants were correctly calculated.

Incorporation of flexibility into rigid-body docking: applications in rounds 3-5 of CAPRI.

We have submitted models for all 9 targets in Rounds 3-5 of CAPRI and have predicted at least 30% of the correct contacts for 4 of the targets and at least 10% of the correct contacts for another 4 targets. We have employed a variety of techniques but have had the greatest success by combining established rigid-body docking with a variety of initial conformations generated by molecular dynamics.

Recognition of phosphorylated-Smad2-containing complexes by a novel Smad interaction motif.

Transforming growth factor beta (TGF-beta) superfamily members signal via complexes of activated Smads, comprising phosphorylated receptor-regulated Smads, such as Smad2 and Smad3, and Smad4. These complexes are recruited to DNA by specific transcription factors. The forkhead/winged-helix transcription factors, XFast-1/XFoxH1a and XFast-3/XFoxH1b, bind an activated Smad heterotrimer comprising two Smad2s and one Smad4. Here we identify a novel Smad2 interaction motif, the Fast/FoxH1 motif (FM), present in all known Fast/FoxH1 family members, N-terminal to the common Smad interaction motif (SIM). The FM is necessary and sufficient to bind active Smad2/Smad4 complexes. The FM differs from the SIM since it discriminates between Smad2 and Smad3, and moreover only binds phosphorylated Smad2 in the context of activated Smad complexes. It is the first Smad interaction motif with this property. Site-directed mutagenesis indicates that the binding site for the FM on a Smad2/Smad4 heterotrimer is a hydrophobic pocket that incorporates the Smad/Smad interface. We demonstrate that the presence of an FM and SIM in the Fast/FoxH1 proteins allows them to compete efficiently for activated Smad2/Smad4 complexes with transcription factors such as Mixer that only contain a SIM. This establishes a hierarchy of Smad-interacting transcription factors, determined by their affinity for active Smad complexes.

Ground and excited state CASPT2 geometry optimizations of small organic molecules.

A method for computing second-order multiconfigurational perturbation theory (CASPT2) energy gradients numerically has been implemented and applied to a range of elementary organic chromophores, including 1,3 butadiene, acrolein, and two protonated Schiff bases. Geometries of ground and excited states-as well as conical intersections-are compared with the corresponding CASSCF structures, illustrating the effect of including the correction for dynamical electron correlation. It is shown that the differences between the two methods are not readily categorized, but that, while individual changes in bond lengths can be quite large ( approximately 0.01-0.02 A), the natures and CASPT2 energetics of the structures remain similar. Exceptions to this tend to be systems that have a strong ionic character and that are not well described at the CASSCF level. Basis set effects (double- vs. triple-zeta) were examined for a limited number of examples, and found to be quite dramatic at both levels of theory.

Reaction path analysis of the "tunable" photoisomerization selectivity of free and locked retinal chromophores.

Multiconfigurational second-order perturbation theory computations and reaction path mapping for the retinal protonated Schiff base models all-trans-nona-2,4,6,8-tetraeniminium and 2-cis-nona-2,4,6,8-tetraeniminium cation demonstrate that, in isolated conditions, retinal chromophores exhibit at least three competing excited-state double bond isomerization paths. These paths are associated with the photoisomerization of the double bonds in positions 9, 11, and 13, respectively, and are controlled by barriers that favor the position 11. The computations provide a basis for the understanding of the observed excited-state lifetime in both naturally occurring and synthetic chromophores in solution and, tentatively, in the protein environment. In particular, we provide a rationalization of the excited-state lifetimes observed for a group of locked retinal chromophores which suggests that photoisomerization in bacteriorhodopsin is the result of simultaneous specific "catalysis" (all-trans --> 13-cis path) accompanied by specific "inhibition" (all-trans --> 11-cis path). The nature of the S(1) --> S(0) decay channel associated with the three paths has also been investigated at the CASSCF level of theory. It is shown that the energy surfaces in the vicinity of the conical intersection for the photoisomerization about the central double bond of retinal (position 11) and the two corresponding lateral double bonds (positions 9 and 13) are structurally different.