Molecular dynamics simulations with a combined quantum mechanical and molecular mechanical (QM/MM) potential have been carried out to investigate the squalene-to-hopene carbocation cyclization mechanism in squalene-hopene cyclase (SHC). The present study is based on free energy simulations by constructing the free energy surface for the cyclization steps along the reaction pathway. The picture that emerges for the carbocation cyclization cascade is a delicate balance of thermodynamic and kinetic control that ultimately favors the formation of the final hopanoids carbon skeleton. A key finding is that the five- to six-membered ring expansion process is not a viable reaction pathway for either C- or D-ring formation in the cyclization reaction. The only significant intermediate is the A/B-bicyclic cyclohexyl cation (III), from which two asynchronous concerted reaction pathways lead to, respectively, the 6,6,6,5-tetracyclic carbon skeleton and the 6,6,6,6,5-pentacyclic hopanoids. Experimentally, these two products are observed to have 1% and 99% yields, respectively, in the wild-type enzyme. We conclude that the product distribution in the wild-type enzyme is dictated by kinetic control of these two reaction pathways.

The discovery of asunaprevir (BMS-650032, 24) is described. This tripeptidic acylsulfonamide inhibitor of the NS3/4A enzyme is currently in phase III clinical trials for the treatment of hepatitis C virus infection. The discovery of 24 was enabled by employing an isolated rabbit heart model to screen for the cardiovascular (CV) liabilities (changes to HR and SNRT) that were responsible for the discontinuation of an earlier lead from this chemical series, BMS-605339 (1), from clinical trials. The structure-activity relationships (SARs) developed with respect to CV effects established that small structural changes to the P2* subsite of the molecule had a significant impact on the CV profile of a given compound. The antiviral activity, preclincial PK profile, and toxicology studies in rat and dog supported clinical development of BMS-650032 (24).

Beta-secretase (BACE) is a critical enzyme in the production of beta-amyloid, a protein that has been implicated as a potential cause of Alzheimer's disease (AD). There are two aspartic acid residues (Asp 32 and Asp 228) present in the catalytic region of BACE that can adopt multiple protonation states. The protonation state and precise location of the protons for these two residues, particularly in the presence of an inhibitor, are subjects of great interest since they have a direct bearing on the mechanism of aspartyl proteases and efforts to model beta-secretase. We have carried out full liner-scaling quantum mechanical (QM) calculations that include Poisson-Boltzmann solvation in order to identify the preferred protonation state and proton location in the presence and absence of an inhibitor. These calculations favor the monoprotonated state in the presence of ligand, and di-deprotonated state in the absence of ligand. Further the proton in the monoprotonated state is located on the inner oxygen of Asp 228. These results have implications for the catalytic mechanism of BACE and related aspartyl proteases. They also provide a reference state for the protein in structure-based modeling studies of this therapeutically important target.

The discovery of BMS-605339 (35), a tripeptidic inhibitor of the NS3/4A enzyme, is described. This compound incorporates a cyclopropylacylsulfonamide moiety that was designed to improve the potency of carboxylic acid prototypes through the introduction of favorable nonbonding interactions within the S1' site of the protease. The identification of 35 was enabled through the optimization and balance of critical properties including potency and pharmacokinetics (PK). This was achieved through modulation of the P2* subsite of the inhibitor which identified the isoquinoline ring system as a key template for improving PK properties with further optimization achieved through functionalization. A methoxy moiety at the C6 position of this isoquinoline ring system proved to be optimal with respect to potency and PK, thus providing the clinical compound 35 which demonstrated antiviral activity in HCV-infected patients.

We describe the implementation of an adaptive umbrella sampling method, making use of the weighted histogram analysis method, for computing multidimensional potential of mean force for chemical reaction in solution. The approach is illustrated by investigating the effect of aqueous solution on the free energy surface for the proton transfer reaction of [H(3)N-H-NH(3)](+) using a combined quantum mechanical and molecular mechanical AM1/TIP3P potential.

Combined quantum mechanical and molecular mechanical (QM/MM) calculations and molecular dynamics simulations of bacteriorhodopsin (bR) in the membrane matrix have been carried out to determine the factors that make significant contributions to the opsin shift. We found that both solvation and interactions with the protein significantly shifts the absorption maximum of the retinal protonated Schiff base, but the effects are much more pronounced in polar solvents such as methanol, acetonitrile, and water than in the protein environment. The differential solvatochromic shifts of PSB in methanol and in bR leads to a bathochromic shift of about 1800 cm(-1). Because the combined QM/MM configuration interaction calculation is essentially a point charge model, this contribution is attributed to the extended point-charge model of Honig and Nakanishi. The incorporation of retinal in bR is accompanied by a change in retinal conformation from the 6-s-cis form in solution to the 6-s-trans configuration in bR. The extension of the pi-conjugated system further increases the red-shift by 2400 cm(-1). The remaining factors are due to the change in dispersion interactions. Using an estimate of about 1000 cm(-1) in the dispersion contribution by Houjou et al., we obtained a theoretical opsin shift of 5200 cm(-1) in bR, which is in excellent agreement with the experimental value of 5100 cm(-1). Structural analysis of the PSB binding site revealed the specific interactions that make contributions to the observed opsin shift. The combined QM/MM method used in the present study provides an opportunity to accurately model the photoisomerization and proton transfer reactions in bR.

No abstract

Hyperpolarization-activated, cyclic nucleotide-modulated (HCN) channels mediate rhythmic electrical activity of neural and cardiac pacemaker cells. Drugs that block these channels slow the beating rate of the heart and are used to treat angina. Here, we characterized the effect of the HCN channel blocker, ZD7288 [4-(N-ethyl-N-phenylamino)-1,2-dimethyl-6-(methylamino) pyrimidinium chloride] on HCN2 channels that were heterologously expressed in Xenopus oocytes. A site-directed mutagenesis approach was used to identify specific residues of the mouse HCN2 channel pore that interact with ZD7288. Two residues (Ala425 and Ile432) located in the S6 transmembrane domain were found to be the primary determinants for block of HCN2 channels by ZD7288. I432A mutant HCN2 channels were ϳ100-fold less sensitive to block by ZD7288. Substitution of Ile432 with more hydrophobic residues (Phe, Leu, or Val) caused only modest shifts in the IC 50 for the drug. HCN1 channels have a Val (Val390) in the equivalent position of Ile432 and are less sensitive to block by ZD7288. Accordingly, mutation of this Val390 to Ile in HCN1 increased the sensitivity of these channels to drug block. Mutation of Ala425 and Ile432 also attenuated the block of HCN2 by the more potent blocker cilobradine. An HCN2 homology model based on the bacterial KcsA K ϩ channel predicts that the phenyl ring of ZD7288 occupies a hydrophobic cavity formed by Ala425 and Ile432 and that the charged ring aligns with the axis of the inner pore closely corresponding to the localization of K ϩ ions observed in the KcsA crystal structure.Pacemaker channels are activated by membrane hyperpolarization and conduct an inward cation current that contributes to spontaneous activity of specialized pacemaker cells in the heart (Baruscotti et al., 2005) and synaptic integration and network rhythmicity of central neurons (Biel et al