OPLS-AA force field parameters have been developed and validated for use in the simulation of 68 unique combinations of room temperature ionic liquids featuring 1-alkyl-3-methylimidazolium [RMIM] (R = Me, Et, Bu, Hex, Oct), N-alkylpyridinium [RPyr], and choline cations, along with Cl(-), PF6(-), BF4(-), NO3(-), AlCl4(-), Al2Cl7(-), TfO(-), saccharinate, and acesulfamate anions. The new parameters were fit to conformational profiles from gas-phase ab initio calculations at the LMP2/cc-pVTZ(-f)//HF/6-31G(d) theory level and compared to experimental condensed-phase structural and thermodynamic data. Monte Carlo simulations of the ionic liquids gave relative deviations from experimental densities of ca. 1-3% at 25 °C for most combinations and also yielded close agreement over a temperature range of 5 to 90 °C. Predicted heats of vaporization compared well with available experimental data and estimates. Transferability of the new parameters to multiple alkyl side-chain lengths for [RMIM] and [RPyr] was determined to give excellent agreement with charges and torsion potentials developed specific to desired alkyl lengths in 35 separate ionic liquid simulations. As further validation of the newly developed parameters, the Kemp elimination reaction of benzisoxazole via piperidine was computed in 1-butyl-3-methylimidazolium hexafluorophosphate [BMIM][PF6] using mixed quantum and molecular mechanics (QM/MM) simulations and was found to give close agreement with the experimental free energy of activation.

The 3M 825EW perfluorosulfonic acid (PFSA) ionomer was doped with the heteropoly acids (HPAs), H 3 PW 12 O 40 (HPW) and H 4 SiW 12 O 40 (HSiW). Dynamic vapor sorption measurements at 95% RH showed a decrease in water content as a result of HPA doping from λ = 8.05 for the undoped 825EW 3M ionmer to λ = 4.40 for a 5% HSiW doped film. FTIR measurement revealed strong interactions between the HPAs, ionomer, and H 2 O. Irrespective of hydration level, it was found that the PFSA films showed tortuous proton diffusion behavior. At maximum hydration and 25 °C, the self-diffusion coefficient of water was found to decrease upon addition of HPA from 4.97 to 2.19 (×10 −6 cm 2 / s) for the undoped 825EW 3M ionomer and 1% HPW, respectively, in excellent agreement with computation. The model at low HPA loadings revealed the decreased diffusion coefficient was due to the water preferentially residing near the HPA as opposed to the SO 3 − groups. The addition of HPA generally improved overall conductivity due to additional formation of hydrogen-bonding networks between the HPA particles. At high RH, it was observed that the proton conductivity increased while the water diffusion coefficient decreased as HPA was incorporated. In addition, a mismatch between the conductivity values and those calculated from the Nernst−Einstein equation using the self-diffusion coefficient of water in the system indicated an enhancement in Grotthuss hopping mechanism upon addition of HPA.

Ionic liquids have been proposed to induce a mechanistic change in the reaction pathway for the fundamentally important base-induced β-elimination class compared to conventional solvents. The role of the reaction medium in the elimination of 1,1,1-tribromo-2,2-bis(3,4-dimethoxyphenyl)ethane via two bases, piperidine and pyrrolidine, has been computationally investigated using methanol and the ionic liquids 1-butyl-3-methylimidazolium tetrafluoroborate and hexafluorophosphate [BMIM][BF(4)] and [BMIM][PF(6)], respectively. QM/MM Monte Carlo simulations utilizing free-energy perturbation theory found the ionic liquids did produce a reaction pathway change from an E1cB-like mechanism in methanol to a pure E2 route that is consistent with experimental observations. The origin of the ionic liquid effect has been found as: (1) a combination of favorable electrostatic interactions, for example, bromine-imidazolium ion, and (2) π-π interactions that enhance the coplanarity between aromatic rings maximizing the electronic effects exerted on the reaction route. Solute-solvent interaction energies have been analyzed and show that liquid clathrate solvation of the transition state is primarily responsible for the observed mechanistic changes. This work provides the first theoretical evidence of an ionic liquid dependent mechanism and elucidates the interplay between sterics and electrostatics crucial to understanding the effect of these unique solvents upon chemical reactions.

Solution pH and the pKa values of ionizable residues are critical factors known to influence enzyme catalysis, structural stability, and dynamical fluctuations. Presented here is an exhaustive computational study utilizing long time constant pH molecular dynamics, pH replica exchange simulations, and kinetic modeling to evaluate pH-dependent conformations, charge dynamics, residue pKa values, and the catalytic activity-pH profile for cellobiohydrolase Cel7B from Melanocarpus albomyces . The predicted pKa values support the role of Glu212 as the catalytic nucleophile and Glu217 as the acid-base residue. The presence of a charge-correlated active site and an extensive hydrogen bonding network is found to be critical in enabling favorable residue orientations for catalysis and shuttling excess protons around the active site. Clusters of amino acids are identified that act in concert to effectively modulate the optimal pH for catalysis while elevating the overall catalytic rate with respect to a noncoupled system. The work presented here demonstrates the complex and critical role of coupled ionizable residues to the proper functioning of cellobiohydrolase Cel7B, functionally related glycosyl hydrolases, and enzymes in general. The simulations also support the use of the CpHMD for the accurate prediction of residue pKa values and to evaluate the impact of pH on protein structure and charge dynamics.

Scaling factors for atomic charges derived from the RM1 semiempirical quantum mechanical wavefunction in conjunction with CM1 and CM3 charge models have been optimized by minimizing errors in absolute free energies of hydration, ΔG(hyd) , for a set of 40 molecules. Monte Carlo statistical mechanics simulations and free energy perturbation theory were used to annihilate the solutes in gas and in a box of TIP4P water molecules. Lennard-Jones parameters from the optimized potentials for liquid simulations-all atom (OPLS-AA) force field were utilized for the organic compounds. Optimal charge scaling factors have been determined as 1.11 and 1.14 for the CM1R and CM3R methods, respectively, and the corresponding unsigned average errors in ΔG(hyd) relative to experiment were 2.05 and 1.89 kcal/mol. Computed errors in aniline and two derivatives were particularly large for RM1 and their removal from the data set lowered the overall errors to 1.61 and 1.75 kcal/mol for CM1R and CM3R. Comparisons are made to the AM1 method which yielded total errors in ΔG(hyd) of 1.50 and 1.64 kcal/mol for CM1A*1.14 and CM3A*1.15, respectively. This work is motivated by the need for a highly efficient yet accurate quantum mechanical (QM) method to study condensed-phase and enzymatic chemical reactions via mixed QM and molecular mechanical (QM/MM) simulations. As an initial test, the Menshutkin reaction between NH(3) and CH(3) Cl in water was computed using a RM1/TIP4P-Ew/CM3R procedure and the resultant ΔG(‡) , ΔG(rxn) , and geometries were in reasonable accord with other computational methods; however, some potentially serious shortcomings in RM1 are discussed.

α-Conotoxin MII (α-CTxMII) is a 16 amino acid peptide with the sequence GCCSNPVCHLEHSNLC containing disulfide bonds between Cys2-Cys8 and Cys3-Cys16. This peptide, isolated from the venom of the marine cone snail Conus magus, is a potent and selective antagonist of neuronal nicotinic acetylcholine receptors (nAChRs). To evaluate the impact of channel-ligand interactions on ligand binding affinity, homology models of the heteropentameric α3β2-nAChR were constructed. The models were created in MODELLER using crystal structures of the Torpedo marmorata-nAChR (Tm-nAChR, PDB ID: 2BG9) and the Aplysia californica-acetylcholine binding protein (Ac-AChBP, PDB ID: 2BR8) as templates for the α3 and β2 subunit isoforms derived from rat neuronal nAChR primary amino acid sequences. Molecular docking calculations were performed with AutoDock to evaluate interactions of the heteropentameric nAChR homology models with the ligands acetylcholine (ACh) and α-CTxMII. The nAChR homology models described here bind ACh with commensurate binding energies to previously reported systems, and identify critical interactions that facilitate both ACh and α-CTxMII ligand binding. The docking calculations revealed an increased binding affinity of the α3β2-nAChR for α-CTxMII with ACh bound to the receptor, which was confirmed through two-electrode voltage clamp experiments on oocytes from Xenopus laevis. These findings provide insights into the inhibition and mechanism of electrostatically driven antagonist properties of the α-CTxMIIs on nAChRs.

The impact of loading the heteropoly acid, 12-phosphotungstic acid (HPW), on a perfluorosulfonic acid (PFSA) proton exchange membrane's morphology was evaluated by means of molecular dynamics (MD) simulations and small-angle X-ray scattering (SAXS) experiments. It is found that the addition of HPW significantly modifies the solvent structure and dynamics in the PFSA membrane, which favors the formation of interconnected proton conducting networks. It is hypothesized that these HPW induced solvent modifications account for the enhanced proton conducting characteristics of these doped membranes. Radial distribution functions and water cluster analysis indicate that the HPW organizes the local solvent water and attracts the nearby excess protons thereby creating localized "nodes" of ordered water and hydronium ions. The "nodes" are found to connect surrounding water wires/channels resulting in a more efficient proton conducting network. This redistribution of solvent and hydronium ions upon addition of HPW creates a shift in the hydrophilic cluster size distribution and the overall membrane morphology. Hydrophilic cluster size analysis indicates that a high percentage of small clusters (d < 15 Å) exist in low HPW doped systems (i.e., 1%), while larger clusters (d > 15 Å) exist for the high HPW doped systems (i.e., 5%). At low hydration levels, the water domains are found to be spheroidal inverted micelles embedded in an ionomer matrix, while at high hydration levels the solvent morphology shifts to a parallel spheroidal elongated cylinder. It is also observed that for the high HPW doping levels the SAXS pattern changes intensity at the low q region and Bragg peaks become present, which indicates the presence of crystalline HPW. These morphological changes create a more interconnected pathway through which the hydrated excess protons may transverse thereby enhancing the PFSA membrane's conductivity

Cyclophilins (Cyp) are a family of cellular enzymes possessing peptidyl-prolyl isomerase activity, which catalyze the cis-trans interconversion of proline-containing peptide bonds. The two most abundant family members, CypA and CypB, have been identified as valid drug targets for a wide range of diseases, including HCV, HIV, and multiple cancers. However, the development of small molecule inhibitors that possess nM potency and high specificity for a particular Cyp is difficult given the complete conservation of all active site residues between the enzymes. Monte Carlo statistical sampling coupled to free energy perturbation theory (MC/FEP) calculations have been carried out to elucidate the origin of the experimentally observed nM inhibition of CypA by acylurea-based derivatives and the >200-fold in vitro selectivity between CypA and CypB from aryl 1-indanylketone-based μM inhibitors. The computed free-energies of binding were in close accord with those derived from experiments. Binding affinity values for the inhibitors were determined to be dependent upon the stabilization strength of the nonbonded interactions provided toward two catalytic residues: Arg55 and Asn102 in CypA and the analogous Arg63 and Asn110 residues in CypB. Fine-tuning of the hydrophobic interactions allowed for enhanced potency among derivatives. The aryl 1-indanylketones are predicted to differentiate between the cyclophilins by using distinct binding motifs that exploit subtle differences in the active site arrangements. Ideas for the development of new selective compounds with the potential for advancement to low-nanomolar inhibition are presented.