The molecular structures and nature of interactions in CH3+Arn(n=1–8) complexes

Dative versus electron-sharing bonding in the isoelectronic argon compounds ArR+ (R = CH3, NH2, OH, and F)

For the series of isoelectronic ArR⁺ (R = CH₃, NH₂, OH, and F) complexes, the nature of the bonding between Ar and R shifts from an Ar → R⁺ dative σ bond in ArCH₃⁺ and ArNH₂⁺ to an Ar⁺–R electron-sharing σ bond in ArOH⁺ and ArF⁺.

Preparation and characterization of chemically bonded argon-boroxol ring cation complexes

The cation complexes [ArB3O4]+, [ArB3O5]+, [ArB4O6]+ and [ArB5O7]+ were prepared via a laser vaporization supersonic ion source in the gas phase. Their vibrational spectra were measured via mass-selected infrared photodissociation spectroscopy. Spectroscopy combined with quantum chemical calculations revealed that the [ArB3O5]+, [ArB4O6]+ and [ArB5O7]+ cation complexes have planar structures each involving an aromatic boroxol ring and an argon-boron covalent bond. In contrast, the [ArB3O4]+ cation is characterized to be a weakly bound complex with a B3O4+ chain structure.

Computational studies of darunavir into HIV-1 protease and DMPC bilayer: necessary conditions for effective binding and the role of the flaps

Human immunodeficiency virus type 1 protease (HIV-1 PR) is one of the main targets toward AIDS therapy. We have selected the potent drug darunavir and a weak inhibitor (fullerene analog) as HIV-1 PR substrates to compare protease's conformational features upon binding. Molecular dynamics (MD), molecular mechanics Poisson-Boltzmann surface area (MM-PBSA), and quantum-mechanical (QM) calculations indicated the importance of the stability of HIV-1 PR flaps toward effective binding: a weak inhibitor may induce flexibility to the flaps, which convert between closed and semiopen states. A water molecule in the darunavir-HIV-1 PR complex bridged the two flap tips of the protease through hydrogen bonding (HB) interactions in a stable structure, a feature that was not observed for the fullerene-HIV-1 PR complex. Additionally, despite that van der Waals interactions and nonpolar contribution to solvation favored permanent fullerene entrapment into the cavity, these interactions alone were not sufficient for effective binding; enhanced electrostatic interactions as observed in the darunavir-complex were the crucial component of the binding energy. An alternative pathway to the usual way of a ligand to access the cavity was also observed for both compounds. Each ligand entered the binding cavity through an opening between the one flap of the protease and a neighboring loop. This suggested that access to the cavity is not necessarily regulated by flap opening. Darunavir exerts its biological action inside the cell, after crossing the membrane barrier. Thus, we also initiated a study on the interactions between darunavir and the DMPC bilayer to reveal that the drug was accommodated inside the bilayer in conformations that resembled its structure into HIV-1 PR, being stabilized via HBs with the lipids and water molecules.

IR Spectroscopy of Microsolvated Aromatic Cluster Ions: Ionization-Induced Switch in Aromatic Molecule–Solvent Recognition

IR spectroscopy, mass spectrometry, and quantum chemical calculations are employed to characterize the intermolecular interaction of a variety of aromatic cations (A+) with several types of solvents. For this purpose, isolated ionic complexes of the type A+–L n , in which A+ is microsolvated by a controlled number (n) of ligands (L), are prepared in a supersonic plasma expansion, and their spectra are obtained by IR photodissociation (IRPD) spectroscopy in a tandem mass spectrometer. Two prototypes of aromatic ion–solvent recognition are considered: (i) microsolvation of acidic aromatic cations in a nonpolar hydrophobic solvent and (ii) microsolvation of bare aromatic hydrocarbon cations in a polar hydrophilic solvent. The analysis of the IRPD spectra of A+–L dimers provides detailed information about the intermolecular interaction between the aromatic ion and the neutral solvent, such as ion–ligand binding energies, the competition between different intermolecular binding motifs (H-bonds, π-bonds, charge–dipole bonds), and its dependence on chemical properties of both the A+ cation and the solvent type L. IRPD spectra of larger A+–L n clusters yield detailed insight into the cluster growth process, including the formation of structural isomers, the competition between ion–solvent and solvent–solvent interactions, and the degree of (non)cooperativity of the intermolecular interactions as a function of solvent type and degree of solvation. The systematic A+–L n cluster studies are shown to reveal valuable new information about fundamental chemical properties of the bare A+ cation, such as proton affinity, acidity, and reactivity. Because of the additional attraction arising from the excess charge, the interaction in the A+–L n cation clusters differs largely from that in the corresponding neutral A–L n clusters with respect to both the interaction strength and the most stable structure, implying in most cases an ionization-induced switch in the preferred aromatic molecule–solvent recognition motif. This process causes severe limitations for the spectroscopic characterization of ion–ligand complexes using popular photoionization techniques, due to the restrictions imposed by the Franck–Condon principle. The present study circumvents these limitations by employing an electron impact cluster ion source for A+–L n generation, which generates predominantly the most stable isomer of a given cluster ion independent of its geometry.

Theoretical calculations: Can Gibbs free energy for intermolecular complexes be predicted efficiently and accurately?

The theoretical study has been performed to refine the procedure for calculations of Gibbs free energy with a relative accuracy of less than 1 kcal/mol. Three benchmark intermolecular complexes are examined via several quantum-chemical methods, including the second-order Moller-Plesset perturbation (MP2), coupled cluster (CCSD(T)), and density functional (BLYP, B3LYP) theories augmented by Dunnings correlation-consistent basis sets. The effects of electron correlation, basis set size, and anharmonicity are systematically analyzed, and the results are compared with available experimental data. The results of the calculations suggest that experimental accuracy can be reached only by extrapolation of MP2 and CCSD(T) total energies to the complete basis set. The contribution of anharmonicity to the zero point energy and TDeltaSint values is fairly small. The new, economic way to reach chemical accuracy in the calculations of the thermodynamic parameters of intermolecular interactions is proposed. In addition, interaction energy (De) and free energy change (DeltaA) for considered species have been evaluated by Carr-Parrinello molecular dynamics (CPMD) simulations and static BLYP-plane wave calculations. The free energy change along the reaction paths were determined by the thermodynamic integration/"Blue Moon Ensemble" technique. Comparison between obtained values, and available experimental and conventional ab initio results has been made. We found that the accuracy of CPMD simulations is affected by several factors, including statistical uncertainty and convergence of constrained forces (TD integration), and the nature of DFT (density functional theory) functional. The results show that CPMD technique is capable of reproducing interaction and free energy with an accuracy of 1 kcal/mol and 2-3 kcal/mol respectively.

The Nature of Variations of Ammonia Proton Affinity in an Argon Environment

The presence of solvent molecules, even inert, may significantly influence processes taking place in the gas phase. The reason for the solvent activity may be found in studies of complexes originating from microsolvation. The structure, thermodynamics, and vibrational properties of NH(4)(+)Ar(n) (n = 0-5) and NH(3)Ar(n) (n = 0-4) complexes are presented in this work with the aim to elucidate the effect of microsolvation on protonation/deprotonation processes. The relation between the nature of interactions in cationic and neutral clusters and the proton affinity is studied as a function of the number of ligands in complexes.