Investigations of neurotrophic inhibitors of FK506 binding protein via Monte Carlo simulations

Interdomain communication suppressing high intrinsic ATPase activity of Sse1 is essential for its co-disaggregase activity with Ssa1

In eukaryotes, Hsp110s are unambiguous cognates of the Hsp70 chaperones, in primary sequence, domain organization, and structure. Hsp110s function as nucleotide exchange factors (NEFs) for the Hsp70s although their apparent loss of Hsp70-like chaperone activity, nature of interdomain communication, and breadth of domain functions are still puzzling. Here, by combining single-molecule FRET, small angle X-ray scattering measurements (SAXS), and MD simulation, we show that yeast Hsp110, Sse1 lacks canonical Hsp70-like interdomain allostery. However, the protein exhibits unique noncanonical conformational changes within its domains. Sse1 maintains an open-lid substrate-binding domain (SBD) in close contact with its nucleotide-binding domain (NBD), irrespective of its ATP hydrolysis status. To further appreciate such ATP-hydrolysis-independent exhaustive interaction between two domains of Hsp110s, NBD-SBD chimera was constructed between Hsp110 (Sse1) and Hsp70 (Ssa1). In Sse1/Ssa1 chimera, we observed undocking of two domains leading to complete loss of NEF activity of Sse1. Interestingly, chimeric proteins exhibited significantly enhanced ATPase rate of Sse1-NBD compared to wild-type protein, implying that intrinsic ATPase activity of the protein remains mostly repressed. Apart from repressing the high ATPase activity of its NBD, interactions between two domains confer thermal stability to Sse1 and play critical role in the (co)chaperoning function of Sse1 in Ssa1-mediated disaggregation activity. Altogether, Sse1 exhibits a unique interdomain interaction, which is essential for its NEF activity, suppression of high intrinsic ATPase activity, co-chaperoning activity in disaggregase machinery, and stability of the protein.

Structure-affinity properties of a high-affinity ligand of FKBP12 studied by molecular simulations of a binding intermediate

With a view to explaining the structure-affinity properties of the ligands of the protein FKBP12, we characterized a binding intermediate state between this protein and a high-affinity ligand. Indeed, the nature and extent of the intermolecular contacts developed in such a species may play a role on its stability and, hence, on the overall association rate. To find the binding intermediate, a molecular simulation protocol was used to unbind the ligand by gradually decreasing the biasing forces introduced. The intermediate was subsequently refined with 17 independent stochastic boundary molecular dynamics simulations that provide a consistent picture of the intermediate state. In this state, the core region of the ligand remains stable, notably because of the two anchoring oxygen atoms that correspond to recurrent motifs found in all FKBP12 ligand core structures. Besides, the non-core regions participate in numerous transient intermolecular and intramolecular contacts. The dynamic aspect of most of the contacts seems important both for the ligand to retain at least a part of its configurational entropy and for avoiding a trapped state along the binding pathway. Since the transient and anchoring contacts contribute to increasing the stability of the intermediate, as a corollary, the dissociation rate constant of this intermediate should be decreased, resulting in an increase of the affinity constant . The present results support our previous conclusions and provide a coherent rationale for explaining the prevalence in high-affinity ligands of (i) the two oxygen atoms found in carbonyl or sulfonyl groups of dissimilar core structures and of (ii) symmetric or pseudo-symmetric mobile groups of atoms found as non-core moieties. Another interesting aspect of the intermediate is the distortion of the flexible 80 s loop of the protein, mainly in its tip region, that promotes the accessibility to the bound state.

Refined OPLS all-atom force field for saturated phosphatidylcholine bilayers at full hydration

We report parametrization of dipalmitoyl-phosphatidylcholine (DPPC) in the framework of the Optimized Parameters for Liquid Simulations all-atom (OPLS-AA) force field. We chose DPPC as it is one of the most studied phospholipid species and thus has plenty of experimental data necessary for model validation, and it is also one of the highly important and abundant lipid types, e.g., in lung surfactant. Overall, PCs have not been previously parametrized in the OPLS-AA force field; thus, there is a need to derive its bonding and nonbonding parameters for both the polar and nonpolar parts of the molecule. In the present study, we determined the parameters for torsion angles in the phosphatidylcholine and glycerol moieties and in the acyl chains, as well the partial atomic charges. In these calculations, we used three methods: (1) Hartree-Fock (HF), (2) second order Møller-Plesset perturbation theory (MP2), and (3) density functional theory (DFT). We also tested the effect of the polar environment by using the polarizable continuum model (PCM), and for acyl chains the van der Waals parameters were also adjusted. In effect, six parameter sets were generated and tested on a DPPC bilayer. Out of these six sets, only one was found to be able to satisfactorily reproduce experimental data for the lipid bilayer. The successful DPPC model was obtained from MP2 calculations in an implicit polar environment (PCM).

Absolute Free Energy of Binding and Entropy of the FKBP12-FK506 Complex: Effects of the Force Field

The hypothetical scanning molecular dynamics (HSMD) method combined with thermodynamic integration (HSMD-TI) has been extended recently for calculating ΔA(0)-the absolute free energy of binding of a ligand to a protein. With HSMD-TI, ΔA(0) is obtained in a new way as a sum of several components, among them is ΔSligand-the change in the conformational entropy as the ligand is transferred from the bulk solvent to the active site-this entropy is obtained by a specific reconstruction procedure. This unique aspect of HSMD (which is useful in rational drug design) is in particular important for treating large ligands, where ΔSligand might be significant. Technically, one should verify that the results for ΔSligand converge-a property that might become more difficult for large ligands; therefore, studying ligands of increasing size would define the range of applicability of HSMD-TI for binding. In this paper, we check the performance of HSMD-TI by applying it to the relatively large ligand FK506 (126 atoms) complexed with the protein FKBP12, where ΔA(0) = -12.8 kcal/mol is known experimentally as well as the crystal structure of the complex. This structure was initially equilibrated by carrying out a 100 ns molecular dynamics trajectory, where the system is modeled by the AMBER force field, TIP3P water, and Particle Mesh Ewald. HSMD-TI calculations were carried out in three conformational regions defined by the intervals [0.2,2], [2,5], and [5,100] ns along the trajectory, where local equilibration of the total energy has been observed; we obtained ΔA(0) = -13.6 ± 1.1, -16.6 ± 1.4, and -16.7 ± 1.4 kcal/mol, respectively indicating the following: (1) The second and third regions belong to the same conformational subspace of the complex, which is different from the [0.2,2] ns subspace. (2) The unsatisfactory result for ΔA(0) obtained in the well equilibrated (hence theoretically preferred) latter regions reflects the nonperfect modeling used, which however (3) has led to the experimental ΔA(0) in the [0.2,2] ns region close to the crystal structure. Keeping the complex near its crystal structure has been a successful approach in the literature. To check this avenue further, we applied harmonic restraints on backbone atoms and obtained unsatisfactory results for ΔA(0), suggesting that implementation of this approach is not straightforward. Converging results for ΔSligand were obtained in all regions, where the result ΔSligand([0.2,2]) = 7.1 ± 1.2 kcal/mol is less region dependent than ΔA(0) and is relatively large probably due to the large ligand.

Parallelized-over-parts computation of absolute binding free energy with docking and molecular dynamics

We present a technique for biomolecular free energy calculations that exploits highly parallelized sampling to significantly reduce the time to results. The technique combines free energies for multiple, nonoverlapping configurational macrostates and is naturally suited to distributed computing. We describe a methodology that uses this technique with docking, molecular dynamics, and free energy perturbation to compute absolute free energies of binding quickly compared to previous methods. The method does not require a priori knowledge of the binding pose as long as the docking technique used can generate reasonable binding modes. We demonstrate the method on the protein FKBP12 and eight of its inhibitors.

A computational analysis of the binding affinities of FKBP12 inhibitors using the MM-PB/SA method

The FK506-binding proteins have been targets of pharmaceutical interests over years. We have studied the binding of a set of 12 nonimmunosuppressive small-molecule inhibitors to FKBP12 through molecular dynamics simulations. Each complex was subjected to 1-ns MD simulation conducted in an explicit solvent environment under constant temperature and pressure. The binding free energy of each complex was then computed by the MM-PB/SA method in the AMBER program. Our MM-PB/SA computation produced a good correlation between the experimentally determined and the computed binding free energies with a correlation coefficient (R(2)) of 0.93 and a standard deviation as low as 0.30 kcal/mol. The vibrational entropy term given by the normal mode analysis was found to be helpful for achieving this correlation. Moreover, an adjustment to one weight factor in the PB/SA model was essential to correct the absolute values of the final binding free energies to a reasonable range. A head-to-head comparison of our MM-PB/SA model with a previously reported Linear Response Approximation (LRA) model suggested that the MM-PB/SA method is more robust in binding affinity prediction for this class of compounds.

Alchemical free energy calculations and multiple conformational substates

Thermodynamic integration (TI) was combined with (adaptive) umbrella sampling to improve the convergence of alchemical free energy simulations in which multiple conformational substates are present. The approach, which we refer to as non-Boltzmann TI (NBTI), was tested by computing the free energy differences between three five-atomic model systems, as well as the free energy difference of solvation between leucine and asparagine. In both cases regular TI failed to give converged results, whereas the NBTI results were free from hysteresis and had standard deviations well below +/-0.7 kcal/mole. We also present theoretical considerations that make it possible to compute free energy differences between simple molecules, such as the five-atomic model systems, by numerical integration of the partition functions at the respective end points.

Molecular modeling of organic and biomolecular systems using BOSS and MCPRO

An overview is provided of the capabilities for the current versions of the BOSS and MCPRO programs for molecular modeling of organic and biomolecular systems. Recent applications are noted, particularly for QM/MM studies of organic and enzymatic reactions and for protein-ligand binding.

HIV-1 reverse transcriptase variants: molecular modeling of Y181C, V106A, L100I, and K103N mutations with nonnucleoside inhibitors using Monte Carlo simulations in combination with a linear response method

The energies and physical descriptors for the binding of 21 novel 1-(2,6-difluorobenzyl)-2-(2,6-difluorophenyl)-benzimidazole (BPBI) analogs to HIV-1 reverse transcriptase (RT) variants Y181C, L100I, V106A, and K103N have been determined using Monte Carlo (MC) simulations. The crystallographic structure of the lead compound, 4-methyl BPBI, was used as a starting point to model the inhibitors in both the mutant bound and the unbound states. The energy terms and physical descriptors obtained from the calculations were reasonably correlated with the respective experimental EC50 values for the inhibitors against the various mutant RTs. Using the linear response correlations from the calculations, 2 novel BPBI inhibitors have been designed and simulations have been carried out. The results show the computed deltaG(binding) values match the experimental data for the analogs. Given the ongoing problem with drug resistance, the ability to predict the activity of novel analogs against variants prior to synthesis is highly advantageous.

Syntheses and EGFR kinase inhibitory activity of 6-substituted-4-anilino [1,7] and [1,8] naphthyridine-3-carbonitriles

The syntheses and EGFR kinase inhibitory activity of a series of 6-substituted-4-anilino [1,7] and [1,8] naphthyridine-3-carbonitriles are described. Both reversible and irreversible binding inhibitors were prepared. These series were compared with each other and with the corresponding 4-anilinoquinoline-3-carbonitriles. Compounds having a 1,7-naphthyridine core structure can retain high potency while those with a 1,8-naphthyridine core are significantly less active. These results are consistent with molecular modeling observations.

Solution structure of a neurotrophic ligand bound to FKBP12 and its effects on protein dynamics

The structure of a recently reported neurotrophic ligand, 3-(3-pyridyl)-1-propyl(2S)-1-(3,3-dimethyl-1, 2-dioxopentyl)-2-pyrrolidinecarboxylate, in complex with FKBP12 was determined using heteronuclear NMR spectroscopy. The inhibitor exhibits a binding mode analogous to that observed for the macrocycle FK506, used widely as an immunosuppressant, with the prolyl ring replacing the pipecolyl moiety and the amide bond in a trans conformation. However, fewer favourable protein-ligand interactions are detected in the structure of the complex, suggesting weaker binding compared with the immunosuppressant drug. Indeed, a micromolar dissociation constant was estimated from the NMR ligand titration profile, in contrast to the previously published nanomolar inhibition activity. Although the inhibitor possesses a remarkable structural simplicity with respect to FK506, 15N relaxation studies show that it induces similar effects on the protein dynamics, stabilizing the conformation of solvent-exposed residues which are important for mediating the interaction of immunophilin/ligand complexes with molecular targets and potentially for the transmission of the neurotrophic action of FKBP12 inhibitors.

Alpha-ketoamides and alpha-ketocarbonyls: conformational analysis and development of all-atom OPLS force field

The molecular structures and barriers for the internal rotation around the OC-CO single bond in four alpha-ketoamides and eight alpha-ketocarbonyls have been determined from the MP3/aug-cc-pVDZ and MP2/aug-cc-pVDZ calculations. Alpha-ketocarbonyls with non-bulky substituents adopt planar conformations with two carbonyl oxygens in s-trans arrangement. The s-cis conformation is significantly less stable due to the electrostatic repulsion between the two carbonyl groups. Primary and secondary alpha-ketoamides are planar when the substituent at the carbonyl carbon is hydrogen or methyl group but tertiary alpha-ketoamides adopt a conformation where the OC-CO unit is significantly bent. Based on current ab initio structural data, a set of OPLS-AA force field parameters has been derived. These parameters can be used for the modeling of a variety of alpha-ketoamide or alpha-ketocarbonyl containing drugs such as novel protease inhibitors or neuroregenerative polyketides.

Monte Carlo calculations on HIV-1 reverse transcriptase complexed with the non-nucleoside inhibitor 8-Cl TIBO: contribution of the L100I and Y181C variants to protein stability and biological activity

A computational model of the non-nucleoside inhibitor 8-Cl TIBO complexed with HIV-1 reverse transcriptase (RT) was constructed in order to determine the binding free energies. Using Monte Carlo simulations, both free energy perturbation and linear response calculations were carried out for the transformation of wild-type RT to two key mutants, Y181C and L100I. The newer linear response method estimates binding free energies based on changes in electrostatic and van der Waals energies and solvent-accessible surface areas. In addition, the change in stability of the protein between the folded and unfolded states was estimated for each of these mutations, which are known to emerge upon treatment with the inhibitor. Results from the calculations revealed that there is a large hydrophobic contribution to protein stability in the native, folded state. The calculated absolute free energies of binding from both the linear response, and also the more rigorous free energy perturbation method, gave excellent agreement with the experimental differences in activity. The success of the relatively rapid linear response method in predicting experimental activities holds promise for estimating the activity of the inhibitors not only against the wild-type RT, but also against key protein variants whose emergence undermines the efficacy of the drugs.

Enzymes that catalyse the restructuring of proteins

The large enzyme families of protein disulfide isomerases and peptidyl prolyl cis/trans isomerases have been shown to assist polypeptide restructuring. Various folding states of polypeptides may serve as substrates of the catalysed reaction. Our understanding of the cellular function of these enzymes is increasing as a result of the availability of more specific inhibitors, the discovery of natural substrates and the use of genetically modified organisms. Further highlights of these studies include insights into the three-dimensional structures of enzyme-ligand complexes, as well as into the mechanism of slow folding phases on the atomic level.

Probing the hydrophobic pocket of farnesyltransferase: aromatic substitution of CAAX peptidomimetics leads to highly potent inhibitors

Cysteine farnesylation at the carboxylate terminal tetrapeptide CAAX of Ras protein is catalyzed by farnesyltransferase. This lipid modification is necessary for regulatory function of both normal and oncogenic Ras. The high frequency of Ras mutation in human cancers has prompted an intensive study on finding ways of controlling oncogenic Ras function. Inhibition of farnesyltransferase is among the most sought after targets for cancer chemotherapy. We report here the design, synthesis and biological characterization of a series of peptidomimetics as farnesyltransferase inhibitors. These compounds are extremely potent towards farnesyltransferase with IC50 values ranging from subnanomolar to low nanomolar concentrations. They have a high selectivity for farnesyltransferase over the closely related geranylgeranyltransferase-I. Structure-activity relationship studies demonstrated that a properly positioned hydrophobic group significantly enhanced inhibition potency, reflecting an improved complementarity to the large hydrophobic pocket in the CAAX binding site.

Estimation of the binding affinities of FKBP12 inhibitors using a linear response method

A series of non-immunosuppressive inhibitors of FK506 binding protein (FKBP12) are investigated using Monte Carlo statistical mechanics simulations. These small molecules may serve as scaffolds for chemical inducers of protein dimerization, and have recently been found to have FKBP12-dependent neurotrophic activity. A linear response model was developed for estimation of absolute binding free energies based on changes in electrostatic and van der Waals energies and solvent-accessible surface areas, which are accumulated during simulations of bound and unbound ligands. With average errors of 0.5 kcal/mol, this method provides a relatively rapid way to screen the binding of ligands while retaining the structural information content of more rigorous free energy calculations.