Monte Carlo Free Energy Perturbation Research Articles

Effect of Fluorination on the Competition of Halogen Bonding and Hydrogen Bonding: Complexes of Fluoroiodomethane with Dimethyl Ether and Trimethylamine.

To further rationalize the competition between halogen and hydrogen bonding, a combined experimental and theoretical study on the weakly bound molecular complexes formed between the combined halogen bond/hydrogen bond donor fluoroiodomethane and the Lewis bases dimethyl ether and trimethylamine (in standard and fully deuterated form) is presented. The experimental data are obtained by recording infrared and Raman spectra of mixtures of the compounds in liquid krypton, at temperatures between 120 and 156 K. The experiments are supported by ab initio calculations at the MP2/aug-cc-pVDZ-PP level, statistical thermodynamics and Monte Carlo free energy perturbation calculations. For the mixtures containing fluoroiodomethane and dimethyl ether a hydrogen-bonded complex with an experimental complexation enthalpy of -7.0(2) kJ mol-1 is identified. Only a single weak spectral feature is observed which can be tentatively assigned to the halogen-bonded complex. For the mixtures involving trimethylamine, both halogen- and hydrogen-bonded complexes are observed, the experimental complexation enthalpies being -12.5(1) and -9.6(2) kJ mol-1 respectively. To evaluate the influence of fluorination on the competition between halogen and hydrogen bonding, the results obtained for fluoroiodomethane are compared with those of a previous study involving difluoroiodomethane.

Hydration Properties and Solvent Effects for All-Atom Solutes in Polarizable Coarse-Grained Water.

Due to the importance of water in chemical and biological systems, a coarse-grained representation of the solvent can greatly simplify the description of the system while retaining key thermodynamic properties of the medium. A multiscale solvation model that couples all-atom solutes and polarizable Martini coarse-grained water (AAX/CGS) is developed to reproduce free energies of hydration of organic solutes. Using Monte Carlo/free energy perturbation (MC/FEP) calculations, results from multiscale and all-atom simulations are compared. Improved accuracy is obtained with the AAX/CGS approach for hydrophobic and sulfur- or halogen-containing solutes, but larger deviations are found for polar solute molecules where hydrogen bonding is featured. Furthermore, solvent effects on conformational and tautomeric equilibria of AA solutes were investigated using AA, CG, and GB/SA solvent models. It is found that the CG solvent model can reproduce well the medium effects from experiment and AA simulations; however, the GB/SA solvent model fails in some cases. A 7-30-fold reduction in computational cost is found for the present AAX/CGS multiscale simulations compared to the AA alternative.

Combining X-ray crystallography and molecular modeling toward the optimization of pyrazolo[3,4-d]pyrimidines as potent c-Src inhibitors active in vivo against neuroblastoma.

c-Src is a tyrosine kinase belonging to the Src-family kinases. It is overexpressed and/or hyperactivated in a variety of cancer cells, thus its inhibition has been predicted to have therapeutic effects in solid tumors. Recently, the pyrazolo[3,4-d]pyrimidine 3 was reported as a dual c-Src/Abl inhibitor. Herein we describe a multidisciplinary drug discovery approach for the optimization of the lead 3 against c-Src. Starting from the X-ray crystal structure of c-Src in complex with 3, Monte Carlo free energy perturbation calculations were applied to guide the design of c-Src inhibitors with improved activities. As a result, the introduction of a meta hydroxyl group on the C4 anilino ring was computed to be particularly favorable. The potency of the synthesized inhibitors was increased with respect to the starting lead 3. The best identified compounds were also found active in the inhibition of neuroblastoma cell proliferation. Furthermore, compound 29 also showed in vivo activity in xenograft model using SH-SY5Y cells.

A cryospectroscopic infrared and Raman study of the C[sbnd]X⋯π halogen bonding motif: Complexes of the CF3Cl, CF3Br, and CF3I with ethyne, propyne and 2-butyne

Experimental information on the CX⋯π halogen bonding motif was obtained by studying the formation of molecular complexes of CF3Cl, CF3Br and CF3I with ethyne, propyne and 2-butyne in liquid krypton, using FTIR and Raman spectroscopy. For CF3Br, experimental evidence was found for the formation of 1:1 complexes with propyne and 2-butyne only, while for CF3I spectroscopic features confirming the existence of the halogen bonded complexes were observed for ethyne, propyne and 2-butyne. In addition, at higher concentrations of CF3I and 2-butyne, weak absorptions due to a 2:1 complex were also observed. The experimental complexation enthalpies, obtained by using spectra recorded at temperatures between 120K and 140K, are −5.9(3) kJmol−1 for CF3I.ethyne, −5.6(3)kJmol−1 for CF3Br.propyne, −8.1(2)kJmol−1 for CF3I.propyne, −7.3(2)kJmol−1 for CF3Br.2-butyne, −10.9(2)kJmol−1 for CF3I.2-butyne and −20.9(7)kJmol−1 for (CF3I)2.2-butyne. The experimental study is supported by theoretical data obtained from ab initio calculations at the MP2/aug-cc-pVDZ(-PP) and MP2/aug-cc-pVTZ(-PP) levels, and Monte Carlo Free Energy Perturbation (MC-FEP) simulations. The experimental and theoretical values on the CX⋯π halogen-bonding motifs studied are compared with previously reported data for the complexes with ethene and propene and with preliminary results obtained for benzene and toluene.

Phenol–Quinone Tautomerism in (Arylazo)naphthols and the Analogous Schiff Bases: Benchmark Calculations

Tautomerization energies of a series of isomeric [(4-R-phenyl)azo]naphthols and the analogous Schiff bases (R = N(CH3)2, OCH3, H, CN, NO2) are calculated by LPNO-CEPA/1-CBS using the def2-TZVPP and def2-QZVPP basis sets for extrapolation. The performance of various density functionals (B3LYP, M06-2X, PW6B95, B2PLYP, mPW2PLYP, PWPB95) as well as MP2 and SCS-MP2 is evaluated against these results. M06-2X and SCS-MP2 yield results close to the LPNO-CEPA/1-CBS values. Solvent effects (CCl4, CHCl3, CH3CN, and CH3OH) are treated by a variety of bulk solvation models (SM8, IEFPCM, COSMO, PBF, and SMD) as well as explicit solvation (Monte Carlo free energy perturbation using the OPLSAA force field).

Exploring the C-X…π Halogen Bonding Motif: An Infrared and Raman Study of the Complexes of CF3X (X = Cl, Br and I) with the Aromatic Model Compounds Benzene and Toluene

The formation of halogen bonded complexes formed between the trifluorohalomethanes CF3Cl, CF3Br and CF3I and the Lewis bases benzene and toluene at temperatures below 150K was investigated using FTIR and Raman spectroscopy. Experiments using liquid krypton as solvent show that for both CF3Br and CF3I substantial fractions of the monomers can be involved in 1:1 complexes. In addition, weak absorptions illustrating the formation of 2:1 complexes between CF3I and benzene are observed. Using spectra recorded at temperatures between 120 and 140 K, observed information on the relative stability was obtained for all complexes by determining the complexation enthalpies in solution. The resulting values for CF3Br·benzene, CF3I·benzene and (CF3I)2·benzene are −6.5(3), −7.6(2) and −14.5(9) kJ mol−1. The values for CF3Br·toluene and CF3I·toluene are −6.2(5) and −7.4(5) kJ mol−1. The experimental complexation enthalpies are compared with theoretical data obtained by combining results from MP2/aug-cc-pVDZ(-PP) and MP2/aug-cc-pVTZ(-PP) ab initio calculations, from statistical thermodynamical calculations and from Monte Carlo Free Energy Perturbation simulations. The data are also compared with results derived for other C-X···π halogen bonded complexes involving unsaturated Lewis bases such as ethene and ethyne.

On the Interaction of Aliphatic Amines and Ammonium Ions with Carboxylic Acids in Solution and in Receptor Pockets

Association energies of the acetate ion with cationic amines bearing one to three methyl groups were calculated in the range of -14 to -17 kcal/mol in aqueous solution by means of the IEF-PCM method at the CCSD(T)/CBS//MP2/aug-cc-pvdz and DFT/B97D/CBS//B97D/aug-cc-pvtz levels. The main stabilization factor for the association is the possibility for the formation of an ionic intermolecular hydrogen bond between the elements of the complex. For a quaternary ammonium ion, the favorable electrostatic interaction energy is the only driving force, and the stabilization energy for the complex is reduced to -4 kcal/mol. The internal free energies of the ion-pair tautomers of the studied species are higher by 10-15 kcal/mol in water than those for the neutral, hydrogen-bonded forms. Monte Carlo free energy perturbation calculations at T = 298 K and p = 1 atm predict -11 to -16 kcal/mol relative solvation free energy in favor of the corresponding ionic form. As a result, the ion-pair tautomer is the prevailing form in aqueous solution and on the extracellular surface of a receptor. Modeling the complex of a protonated ligand interacting with an Asp/Glu carboxylate side-chain in the binding cavity of a receptor, two strongly bound water molecules were considered so as to form hydrogen-bonded water bridges between the elements of the ion-pair. Nonetheless, the low polarity environment mimicked by a chloroform solvent cannot stabilize the ionic tautomer. A proton jump was predicted, which suggests that acetylcholine, an inherent cation by structure, might have evolved as the natural agonist for muscarinic receptors because a quaternary ammonium system assures the maintenance of the ion-pair form with a carboxylate side-chain in a protein cavity, the latter perhaps then being needed for receptor activation.

Methyl Effects on Protein–Ligand Binding

The effects of addition of a methyl group to a lead compound on biological activity are examined. A literature analysis of >2000 cases reveals that an activity boost of a factor of 10 or more is found with an 8% frequency, and a 100-fold boost is a 1 in 200 event. Four cases in the latter category are analyzed in depth to elucidate any unusual aspects of the protein-ligand binding, distribution of water molecules, and changes in conformational energetics. The analyses include Monte Carlo/free-energy perturbation (MC/FEP) calculations for methyl replacements in inhibitor series for p38α MAP kinase, ACK1, PTP1B, and thrombin. Methyl substitutions ortho to an aryl ring can be particularly effective at improving activity by inducing a propitious conformational change. The greatest improvements in activity arise from coupling the conformational gain with the burial of the methyl group in a hydrophobic region of the protein.

Solute–solvent interactions in cryosolutions: a study of halothane–ammonia complexes

The formation of C-H···N bonded complexes of halothane with ammonia has been studied using infrared and Raman spectroscopy of solutions in the liquid rare gases argon, krypton and xenon, of supersonic jet expansions and of room temperature vapor phase mixtures. For the solutions and for the vapor phase experiments, the formation of complexes with 1:1 and 1:2 stoichiometry was observed. The complexation enthalpy for the 1:1 complex was determined to be -20 (1) kJ mol(-1) in the vapor phase, -17.0 (5) kJ mol(-1) in liquid xenon and -17.3 (6) kJ mol(-1) in liquid krypton. For the 1:2 complex in liquid xenon, the complexation enthalpy was determined to be -31.5 (12) kJ mol(-1). Using the complexation enthalpies for the vapor phase and for the solutions in liquid xenon and krypton, a critical assessment is made of the Monte Carlo Free Energy Perturbation approach to model solvent influences on the thermodynamical properties of the cryosolutions. The influences of temperature and solvent on the complexation shifts of the halothane C-H stretching mode are discussed.

On the weakly C–H⋯π hydrogen bonded complexes of sevoflurane and benzene

A vibrational assignment of the anaesthetic sevoflurane, (CF(3))(2)CHOCH(2)F, is proposed and its interaction with the aromatic model compound benzene is studied using vibrational spectroscopy of supersonic jet expansions and of cryosolutions in liquid xenon. Ab initio calculations, at the MP2/cc-pVDZ and MP2/aug-cc-pVDZ levels, predict two isomers for the 1 : 1 complex, one in which the near-cis, gauche conformer of sevoflurane is hydrogen bonded through its isopropyl-hydrogen atom, the other in which the same conformer is bonded through a bifurcated hydrogen bond with the fluoromethyl hydrogen atoms. From the experiments it is shown that the two isomers are formed, however with a strong population dominance of the isopropyl-bonded species, both in the jet and liquid phase spectra. The experimental complexation enthalpy in liquid xenon, ΔH(o)(LXe), of this species equals -10.9(2) kJ mol(-1), as derived from the temperature dependent behaviour of the cryosolution spectra. Theoretical complexation enthalpies in liquid xenon were obtained by combining the complete basis set extrapolated complexation energies at the MP2/aug-cc-pVXZ (X = D,T) level with corrections derived from statistical thermodynamics and Monte Carlo Free Energy Perturbation calculations, resulting in a complexation enthalpy of -11.2(3) kJ mol(-1) for the isopropyl-bonded complex, in very good agreement with the experimental value, and of -11.4(4) kJ mol(-1), for the fluoromethyl-bonded complex. The Monte Carlo calculations show that the solvation entropy of the isopropyl-bonded species is considerably higher than that of the fluoromethyl-bonded complex, which assists in explaining its dominance in the liquid phase spectra.

Effects of Water Placement on Predictions of Binding Affinities for p38α MAP Kinase Inhibitors

Monte Carlo free energy perturbation (MC/FEP) calculations have been applied to compute the relative binding affinities of 17 congeneric pyridazo-pyrimidinone inhibitors of the protein p38α MAP kinase. Overall correlation with experiment was found to be modest when the complexes were hydrated using a traditional procedure with a stored solvent box. Significant improvements in accuracy were obtained when the MC/FEP calculations were repeated using initial solvent distributions optimized by the water placement algorithm JAWS. The results underscore the importance of accurate placement of water molecules in a ligand binding site for the reliable prediction of relative free energies of binding.

CX⋅⋅⋅O Halogen Bonding: Interactions of Trifluoromethyl Halides with Dimethyl Ether

The formation of weakly bound molecular complexes between dimethyl ether (DME) and the trifluoromethyl halides CF(3)Cl, CF(3)Br and CF(3)I dissolved in liquid argon and in liquid krypton is investigated, using Raman and FTIR spectroscopy. For all halides evidence is found for the formation of C-X...O halogen-bonded 1:1 complexes. At higher concentrations of CF(3)Br, a weak absorption due to a 1:2 complex is also observed. Using spectra recorded at temperatures between 87 and 125 K, the complexation enthalpies for the complexes are determined to be -6.8(3) kJ mol(-1) (DME x CF(3)Cl), -10.2(1) kJ mol(-1) (DME x CF(3)Br), -15.5(1) kJ mol(-1) (DME x CF(3)I), and -17.8(5) kJ mol(-1) [DME(x CF(3)Br)(2)]. Structural and spectral information on the complexes is obtained from ab initio calculations at the MP2/6-311++G(d,p) and MP2/6-311++G(d,p)+LanL2DZ* levels. By applying Monte Carlo free energy perturbation calculations to account for the solvent influences, and statistical thermodynamics to estimate the zero-point vibrational and thermal influences, the ab initio complexation energies are converted into complexation enthalpies for the solutions in liquid argon. The resulting values are compared with the experimental data deduced from the cryosolutions.

Vancomycin analogs: Seeking improved binding of d-Ala- d-Ala and d-Ala- d-Lac peptides by side-chain and backbone modifications

In order to seek vancomycin analogs with improved performance against VanA and VanB resistant bacterial strains, extensive computational investigations have been performed to examine the effects of side-chain and backbone modifications. Changes in binding affinities for tripeptide cell-wall precursor mimics, Ac 2- l-Lys- d-Ala- d-Ala ( 3) and Ac 2- l-Lys- d-Ala- d-Lac ( 4), with vancomycin analogs were computed with Monte Carlo/free energy perturbation (MC/FEP) calculations. Replacements of the 3-hydroxyl group in residue 7 with small alkyl or alkoxy groups, which improve contacts with the methyl side chain of the ligands’ d-Ala residue, are predicted to be the most promising to enhance binding for both ligands. The previously reported amine backbone modification as in 5 is shown to complement the hydrophobic modifications for binding monoacetylated tripeptides. In addition, replacement of the hydroxyl groups in residues 5 and 7 by fluorine is computed to have negligible impact on binding the tripeptides, though it may be pharmacologically advantageous.

Vancomycin resistance: Modeling backbone variants with d-Ala- d-Ala and d-Ala- d-Lac peptides

To seek vancomycin analogs with broader antibacterial activity, effects of backbone modifications for the agylcon 2 on binding with d-Ala- d-Ala- and d-Ala- d-Lac-containing peptides were investigated by Monte Carlo/free energy perturbation (MC/FEP) calculations. The experimental trend in binding affinities for 2 with three tripeptides was well reproduced. Possible modifications of the peptide bond between residues 4 and 5 were then considered, specifically for conversion of the O C NH linkage to CH 2NH 2 + ( 6), FC CH ( 7), HC CH ( 8), and HN C O ( 9). The MC/FEP results did not yield binding improvements for 7, 8, and 9, though the fluorovinyl replacement is relatively benign. The previously reported analog 6 remains as the only variant that exhibits improved affinity for the d-Ala- d-Lac sequence and acceptable affinity for the d-Ala- d-Ala sequence.

Improved Interaction Potentials for Charged Residues in Proteins

Electrostatic interactions dominate the structure and free energy of biomolecules. To obtain accurate free energies involving charged groups from molecular simulations, OPLS-AA parameters have been reoptimized using Monte Carlo free energy perturbation. New parameters fit a self-consistent, experimental set of hydration free energies for acetate (Asp), propionate (Glu), 4-methylimidazolium (Hip), n-butylammonium (Lys), and n-propylguanidinium (Arg), all resembling charged residue side chains, including beta-carbons. It is shown that OPLS-AA free energies depend critically on the type of water model, TIP4P or TIP3P; i.e., each water model requires specific water-charged molecule interaction potentials. New models (models 1 and 3) are thus described for both water models. Uncertainties in relative free energies of charged residues are approximately 2 kcal/mol with the new parameters, due to variations in system setup (MAEs of ca. 1 kcal/mol) and noise from simulations (ca. 1 kcal/mol). The latter error of approximately 1 kcal/mol contrasts MAEs from standard OPLS-AA of up to 13 kcal/mol for the entire series of charged residues or up to 5 kcal/mol for the cationic series Lys, Arg, and Hip. The new parameters can be used directly in molecular simulations with no modification of neutral residues needed and are envisioned to be particular important in simulations where charged residues change environment.

Infrared Spectra of the Complexes of Trifluoroethene with Dimethyl Ether, Acetone, and Oxirane: A Cryosolution Study

Infrared spectra of solutions of trifluoroethene and dimethyl ether, acetone, or oxirane in liquid krypton and liquid argon have been studied. For each Lewis base the formation of a 1:1 complex with the Lewis acid was observed. The C-H stretching of trifluoroethene being perturbed by a strong Fermi resonance, the complexes with trifuloroethene-d were also investigated and showed that in each case the hydrogen bond between the acid and base is of the traditional, red-shifting type. The structures of the complexes were investigated using ab initio calculations. These indicate that with dimethyl ether and acetone two different isomeres can be formed, but with a single one detected in the solution in each case. The Fermi resonance in the complex with unlabeled trifluoroethene is discussed using data derived form ab initio potential and dipole hypersurface calculations. The complexation enthalpies of the complexes were obtained from temperature dependent studies of the solutions and are discussed in relation to the ab initio complexation energies and Monte Carlo free energy perturbation calculations of solvent effects.

Benzoderivatives of Nucleic Acid Bases as Modified DNA Building Blocks

The tautomeric properties of benzoderivatives of the canonical nucleic acid bases have been studied by using different computational approaches. Attention has been paid to the impact of the benzene group in altering the tautomeric preferences of the canonical bases both in the gas phase and in aqueous solution. To this end, relative solvation free energies of the tautomers determined from Self-Consistent Reaction Field continuum calculations and Monte Carlo-Free Energy Perturbation are combined with gas-phase tautomerization free energies determined from quantum mechanical calculations. The results provide a detailed picture of the tautomeric preferences of the benzoderivatives of nucleic acid bases. This information is used to examine the recognition properties of the preferred tautomers of the benzo-fused derivatives, paying particular attention to the ability to form Watson-Crick hydrogen-bonding and stacking interactions as well as to the hydrophobic nature of the modified bases. The implications of present results on the potential use of benzo-fused bases as potential building blocks in modified DNA duplexes are examined.

Theoretical Conformational Analysis for Neurotransmitters in the Gas Phase and in Aqueous Solution. Serotonin.

Conformational analyses have been performed for protonated serotonin in the gas phase, aqueous solution, and in the binding cavity of a 5-HT2A receptor model. DFT geometry optimizations have been performed in the gas phase at the B3LYP/6-31G* levels. Optimized calculations up to the B3LYP/6-311++G** level find two low-energy gauche conformations separated by 8-10 kcal/mol barriers from a trans conformation with relative energy of about 6 kcal/mol. In aqueous solution as concluded from IEF-PCM/B3LYP/6-31G* and IEF-PCM/MP2/6-31G*//IEF-PCM/B3LYP/6-31G* continuum solvent calculations as well as Monte Carlo free energy perturbation simulations with explicit solvent molecules, those barriers decrease to 2-7 kcal/mol, while the two gauche and one trans conformers are within a 3 kcal/mol relative free energy range. The solute is strongly hydrated by about three water molecules around the -NH3(+) group and by one water molecule for each of the pyrrole and phenolic hydrogen atoms. Docking studies of the protonated ligand predicted both gauche and trans ligand conformers to favorably interact with the 5-HT2A receptor in its hypothesized binding cavity. The theoretical studies confirm the experimental results regarding strong interactions with the Asp155 and Ser159 residues (TM helix III) and the interactions of the indole ring with Phe, Trp, and Tyr side chains in TM V, VI, and VII helices within a 24 kcal/mol range for the relative interaction energies.

Free energy and entropy calculations in the conformational equilibrium of 1,2-dichloroethane in water

Using Monte Carlo free energy perturbation simulations, the entropy variation of the trans– gauche conformational equilibrium of 1,2-dichloroethane in water has been computed. Thermodynamic cycles have been used to evaluate the contribution of the van der Waals, electrostatic and polarization forces. Although the experimental free energy is reproduced, the simultaneous reproduction of the experimental free energy and entropy thermodynamic quantities cannot be obtained. The scaling of the solute atomic charge and polarizability values is not able to improve the results. To get proper entropy values for this conformational equilibrium, a new reparametrization involving both solute and solvent molecules is required.

Solvent Effect on the Conformational Equilibrium of 1,2-Dichloroethane in Water. The Role of Solute Polarization

The conformational equilibrium of the 1,2-dichloroethane molecule in water has been studied with Monte Carlo free energy perturbation simulations. A polarizable model for the solute molecule combined with the nonpolarizable TIP4P model for the water molecules has been employed. The approach of only taking into account the polarization of solute produces a rather small increase of computational time with respect to the nonpolarizable simulations. Relative solvation free energies of the trans and gauche 1,2-dichloroethane conformations have been calculated through three perturbation paths which correspond to the progressive incorporation of three different intermolecular solute−solvent energy contributions (van der Waals, Coulombic, and polarization). This stepped procedure allows the separate analysis of the different interaction forces. The solute polarization is studied by introducing atom polarizabilities in the 1,2-dichloroethane molecule. Three polarization procedures which differ in the treatment of the intramolecular interactions among the polarization sites of the solute molecule have been analyzed. It is found that the inclusion of solute polarization produces a different reorganization of solvent molecules around the 1,2-dichloroethane conformations.