Photoswitchable Molecular Units with Tunable Nonlinear Optical Activity: A Theoretical Investigation.

The first-, second-, and third-order molecular nonlinear optical properties, including two-photon absorption of a series of derivatives, involving two dithienylethene (DTE) groups connected by several molecular linkers (bis(ethylene-1,2-dithiolato)Ni- (NiBDT), naphthalene, quasilinear oligothiophene chains), are investigated by employing density functional theory (DFT). These properties can be efficiently controlled by DTE switches, in connection with light of appropriate frequency. NiBDT, as a linker, is associated with a greater contrast, in comparison to naphthalene, between the first and second hyperpolarizabilities of the "open-open" and the "closed-closed" isomers. This is explained by invoking the low-lying excited states of NiBDT. It is shown that the second hyperpolarizability can be used as an index, which follows the structural changes induced by photochromism. Assuming a Förster type transfer mechanism, the intramolecular excited-state energy transfer (EET) mechanism is studied. Two important parameters related to this are computed: the electronic coupling (VDA) between the donor and acceptor fragments as well as the overlap between the absorption and emission spectra of the donor and acceptor groups. NiBDT as a linker is associated with a low electronic coupling, VDA, value. We found that VDA is affected by molecular geometry. Our results predict that the linker strongly influences the communication between the open-closed DTE groups. The sensitivity of the molecular nonlinear optical properties could assist with identification of molecular isomers.

Tailoring Colors by O Annulation of Polycyclic Aromatic Hydrocarbons.

The synthesis of O-doped polyaromatic hydro- carbons in which two polycyclic aromatic hydrocarbon sub units are bridged through one or two O atoms has been achieved. This includes high-yield ring-closure key steps that, depending on the reaction conditions, result in the formation of furanyl or pyranopyranyl linkages through intramolecular C-O bond formation. Comprehensive photophysical measurements in solution showed that these compounds have exceptionally high emission yields and tunable absorption properties throughout the UV/Vis spectral region. Electrochemical investigations showed that in all cases O annulation increases the electron-donor capabilities by raising the HOMO energy level, whereas the LUMO energy level is less affected. Moreover, third-order nonlinear optical (NLO) measurements on solutions or thin films containing the dyes showed very good values of the second hyperpolarizability. Importantly, poly(methyl methacrylate) films containing the pyranopyranyl derivatives exhibited weak linear absorption and NLO absorption compared to the nonlinearity and NLO refraction, respectively, and thus revealed them to be exceptional organic materials for photonic devices.

Reverse Micelles As Antioxidant Carriers: An Experimental and Molecular Dynamics Study.

Water-in-oil microemulsions with biocompatible components were formulated to be used as carriers of natural antioxidants, such as hydroxytyrosol (HT) and gallic acid (GA). The system was composed of a mixture of natural surfactants, lecithin and monoglycerides, medium chain triglycerides, and aqueous phase. A dual approach was undertaken to study the structure and dynamics of these complicated systems. First, experimental data were collected by using adequate techniques, such as dynamic light scattering (DLS) and electron paramagnetic resonance (EPR) spectroscopy. Following this, a coarse-grained molecular dynamics (CGMD) study based on the experimental composition using the MARTINI force field was conducted. The simulations revealed the spontaneous formation of reverse micelles (RMs) starting from completely random initial conformations, underlying their enhanced thermodynamic stability. The location of the bioactive molecules, as well as the structure of the RM, were in accordance with the experimental findings. Furthermore, GA molecules were found to be located inside the water core, in contrast to the HT ones, which seem to lie at the surfactant interfacial layer. The difference in the antioxidants' molecular location was only revealed in detail from the computational analysis and explains the RM's swelling observed by GA in DLS measurements.

Compilation of Data and Modelling of Nanoparticle Interactions and Toxicity in the NanoPUZZLES Project.

The particular properties of nanomaterials have led to their rapidly increasing use in diverse fields of application. However, safety assessment is not keeping pace and there are still gaps in the understanding of their hazards. Computational models predicting nanotoxicity, such as (quantitative) structure-activity relationships ((Q)SARs), can contribute to safety evaluation, in line with general efforts to apply alternative methods in chemical risk assessment. Their development is highly dependent on the availability of reliable and high quality experimental data, both regarding the compounds' properties as well as the measured toxic effects. In particular, "nano-QSARs" should take the nano-specific characteristics into account. The information compiled needs to be well organized, quality controlled and standardized. Integrating the data in an overarching, structured data collection aims to (a) organize the data in a way to support modelling, (b) make (meta)data necessary for modelling available, and (c) add value by making a comparison between data from different sources possible.Based on the available data, specific descriptors can be derived to parameterize the nanomaterial-specific structure and physico-chemical properties appropriately. Furthermore, the interactions between nanoparticles and biological systems as well as small molecules, which can lead to modifications of the structure of the active nanoparticles, need to be described and taken into account in the development of models to predict the biological activity and toxicity of nanoparticles. The EU NanoPUZZLES project was part of a global cooperative effort to advance data availability and modelling approaches supporting the characterization and evaluation of nanomaterials.

Development of a potent 2-oxoamide inhibitor of secreted phospholipase A2 guided by molecular docking calculations and molecular dynamics simulations.

Inhibition of group IIA secreted phospholipase A2 (GIIA sPLA2) has been an important objective for medicinal chemists. We have previously shown that inhibitors incorporating the 2-oxoamide functionality may inhibit human and mouse GIIA sPLA2s. Herein, the development of new potent inhibitors by molecular docking calculations using the structure of the known inhibitor 7 as scaffold, are described. Synthesis and biological evaluation of the new compounds revealed that the long chain 2-oxoamide based on (S)-valine GK241 led to improved activity (IC50=143 nM and 68 nM against human and mouse GIIA sPLA2, respectively). In addition, molecular dynamics simulations were employed to shed light on GK241 potent and selective inhibitory activity.

A Computational Study of the Interaction and Polarization Effects of Complexes Involving Molecular Graphene and C60 or a Nucleobases.

A systematic analysis of the molecular structure, energetics, electronic (hyper)polarizabilities and their interaction-induced counterparts of C60 with a series of molecular graphene (MG) models, CmHn, where m = 24, 84, 114, 222, 366, 546 and n = 12, 24, 30, 42, 54, 66, was performed. All the reported data were computed by employing density functional theory and a series of basis sets. The main goal of the study is to investigate how alteration of the size of the MG model affects the strength of the interaction, charge rearrangement, and polarization and interaction-induced polarization of the complex, C60-MG. A Hirshfeld-based scheme has been employed in order to provide information on the intrinsic polarizability density representations of the reported complexes. It was found that the interaction energy increases approaching a limit of -26.98 kcal/mol for m = 366 and 546; the polarizability and second hyperpolarizability increase with increasing the size of MG. An opposite trend was observed for the dipole moment. Interestingly, the variation of the first hyperpolarizability is relatively small with m. Since polarizability is a key factor for the stability of molecular graphene with nucleobases (NB), a study of the magnitude of the interaction-induced polarizability of C84H24-NB complexes is also reported, aiming to reveal changes of its magnitude with the type of NB. The binding strength of C84H24-NB complexes is also computed and found to be in agreement with available theoretical and experimental data. The interaction involved in C60 B12N12H24-NB complexes has also been considered, featuring the effect of contamination on the binding strength between MG and NBs.

Biological and computational evaluation of resveratrol inhibitors against Alzheimer's disease.

It has been reported that beta amyloid induces production of radical oxygen species and oxidative stress in neuronal cells, which in turn upregulates ß-secretase (BACE-1) expression and beta amyloid levels, thereby propagating oxidative stress and increasing neuronal injury. A series of resveratrol derivatives, known to be inhibitors of oxidative stress-induced neuronal cell death (oxytosis) were biologically evaluated against BACE-1 using homogeneous time-resolved fluorescence (TRF) assay. Correlation between oxytosis inhibitory and BACE-1 inhibitory activity of resveratrol derivatives was statistically significant, supporting the notion that BACE-1 may act as pivotal mediator of neuronal cell oxytosis. Four of the biologically evaluated resveratrol analogs demonstrated considerably higher activity than resveratrol in either assay. The discovery of some "hits" led us to initiate detailed docking studies associated with Molecular Dynamics in order to provide a plausible explanation for the experimental results and understand their molecular basis of action.

Investigation of the interactions of silibinin with 2-hydroxypropyl-ß-cyclodextrin through biophysical techniques and computational methods.

Cyclodextrins (CDs) are a well-known class of supermolecules that have been widely used to protect drugs against conjugation and metabolic inactivation as well as to enhance the aqueous solubility and hence to ameliorate the oral bioavailability of sparingly soluble drug molecules. The hepatoprotectant drug silibinin can be incorporated into CDs, and here we elucidate the interaction between the drug and the host at the molecular level. The complexation product of silibinin with 2-hydroxypropyl-ß-cyclodextrin (HP-ß-CD) is characterized by Differential Scanning Calorimetry, mass spectrometry, solid and liquid high-resolution NMR spectroscopy. The chemical shift changes using (13)C CP/MAS on the complexing of the guest with the host provided significant information on the molecular interactions, and they were in agreement with the 2D NOESY results. These results point out that in both solid and liquid forms, the drug is engulfed and interacts with HP-ß-CD in identical manner. Molecular dynamics calculations have been performed to examine the thermodynamic characteristics associated with the silibinin-HP-ß-CD interactions and to study the stability of the complex. To approximate the physiological conditions, the aqueous solubility and dissolution characteristics of the complex at pH states simulating those of the upper gastrointestinal tract have been applied. To evaluate the antiproliferative activity of silibinin-HP-ß-CD complex comparatively to silibinin in MCF-7 human cancer cells, MTT assays have been performed.

Stability and binding effects of silver(I) complexes at lipoxygenase-1.

An anti-inflammatory complex of Ag(I), namely [Ag(tpp)3(asp)](dmf) [tpp = triphenylphosphine, aspH = aspirin, dmf = N,N-dimethylformamide], was synthesized in an attempt to develop novel metallotherapeutic molecules. STD (1)H NMR experiments were used to examine if this complex binds to LOX-1. The (1)H NMR spectra in buffer Tris/D2O betrayed the existence of two complexes: the complex of aspirin and the complex of salicylic acid produced after deacetylation of aspirin. Nevertheless, the STD spectra showed that only the complex of salicylic acid is bound to the enzyme. Molecular docking and dynamics were used to complement our study. The complexes were stabilized inside a large LOX-1 cavity by establishing a network of hydrogen bonds and steric interactions. The complex formation with salicylic acid was more favorable. The in silico results provide a plausible explanation of the experimental results, which showed that only the complex with salicylic acid enters the binding cavity.

Elucidation of conformational states, dynamics, and mechanism of binding in human κ-opioid receptor complexes.

Opioid G protein-coupled receptors (GPCRs) have been implicated in modulating pain, addiction, psychotomimesis, mood and memory, among other functions. We have employed the recently reported crystal structure of the human κ-opioid receptor (κ-OR) and performed molecular dynamics (MD), free energy, and ab initio calculations to elucidate the binding mechanism in complexes with antagonist JDTic and agonist SalA. The two systems were modeled in water and in DPPC lipid bilayers, in order to investigate the effect of the membrane upon conformational dynamics. MD and Atoms in Molecules (AIM) ab initio calculations for the complexes in water showed that each ligand was stabilized inside the binding site of the receptor through hydrogen bond interactions that involved residues Asp138 (with JDTic) and Gln115, His291, Leu212 (with SalA). The static description offered by the crystal structure was overcome to reveal a structural rearrangement of the binding pocket, which facilitated additional interactions between JDTic and Glu209/Tyr139. The role of Glu209 was emphasized, since it belongs to an extracellular loop that covers the binding site of the receptor and is crucial for ligand entrapment. The above interactions were retained in membrane complexes (SalA forms additional hydrogen bonds with Tyr139/312), except the Tyr139 interaction, which is abolished in the JDTic complex. For the first time, we report that JDTic alternates between a "V-shape" (stabilized via a water-mediated intramolecular interaction) and a more extended conformation, a feature that offers enough suppleness for effective binding. Moreover, MM-PBSA calculations showed that the more efficient JDTic binding to κ-OR compared to SalA (ΔGJDTic = -31.6 kcal mol(-1), ΔGSalA = -9.8 kcal mol(-1)) is attributed mostly to differences in electrostatic contributions. Importantly, our results are in qualitative agreement with the experiments (ΔGJDTic,exp = -14.4 kcal mol(-1), ΔGSalA,exp = -10.8 kcal mol(-1)). This study provides previously unattainable information on the dynamics of human κ-OR and insight on the rational design of drugs with improved pharmacological properties.

Comparative study of the AT1 receptor prodrug antagonist candesartan cilexetil with other sartans on the interactions with membrane bilayers.

Drug-membrane interactions of the candesartan cilexetil (TCV-116) have been studied on molecular basis by applying various complementary biophysical techniques namely differential scanning calorimetry (DSC), Raman spectroscopy, small and wide angle X-ray scattering (SAXS and WAXS), solution ¹H and ¹³C nuclear magnetic resonance (NMR) and solid state ¹³C and ³¹P (NMR) spectroscopies. In addition, ³¹P cross polarization (CP) NMR broadline fitting methodology in combination with ab initio computations has been applied. Finally molecular dynamics (MD) was applied to find the low energy conformation and position of candesartan cilexetil in the bilayers. Thus, the experimental results complemented with in silico MD results provided information on the localization, orientation, and dynamic properties of TCV-116 in the lipidic environment. The effects of this prodrug have been compared with other AT1 receptor antagonists hitherto studied. The prodrug TCV-116 as other sartans has been found to be accommodated in the polar/apolar interface of the bilayer. In particular, it anchors in the mesophase region of the lipid bilayers with the tetrazole group oriented toward the polar headgroup spanning from water interface toward the mesophase and upper segment of the hydrophobic region. In spite of their localization identity, their thermal and dynamic effects are distinct pointing out that each sartan has its own fingerprint of action in the membrane bilayer, which is determined by the parameters derived from the above mentioned biophysical techniques.

On the vibrational linear and nonlinear optical properties of compounds involving noble gas atoms: HXeOXeH, HXeOXeF, and FXeOXeF.

The vibrational (hyper)polarizabilities of some selected Xe derivatives are studied in the context of Bishop-Kirtman perturbation theory (BKPT) and numerical finite field methodology. It was found that for this set of rare gas compounds, the static vibrational properties are quite large, in comparison to the corresponding electronic ones, especially those of the second hyperpolarizability. This also holds for the dc-Pockels ß(-ω;ω,0), Kerr γ(-ω;ω,0,0) and electric field second harmonic generation γ (-2ω;ω,ω,0) effects, although the computed nuclear relaxation (nr) vibrational contributions are smaller in magnitude than the static ones. HXeOXeH was used to study the effects of electron correlation, basis set, and geometry. Geometry effects were found to lead to noticeable changes of the vibrational and electronic second hyperpolarizability. A limited study of the effect of Xe insertion to the nr vibrational properties is also reported. Assessment of the results revealed that Xe insertion has a remarkable effect on the nr (hyper)polarizabilities. In terms of the BKPT, this is associated with a remarkable increase of the electrical and mechanical anharmonicity terms. The latter is consistent with the anharmonic character of several vibrational modes reported for rare gas compounds.

Performance of density functional theory in computing nonresonant vibrational (hyper)polarizabilities.

A set of exchange-correlation functionals, including BLYP, PBE0, B3LYP, BHandHLYP, CAM-B3LYP, LC-BLYP, and HSE, has been used to determine static and dynamic nonresonant (nuclear relaxation) vibrational (hyper)polarizabilities for a series of all-trans polymethineimine (PMI) oligomers containing up to eight monomer units. These functionals are assessed against reference values obtained using the Møller-Plesset second-order perturbation theory (MP2) and CCSD methods. For the smallest oligomer, CCSD(T) calculations confirm the choice of MP2 and CCSD as appropriate for assessing the density functionals. By and large, CAM-B3LYP is the most successful, because it is best for the nuclear relaxation contribution to the static linear polarizability, intensity-dependent refractive index second hyperpolarizability, static second hyperpolarizability, and is close to the best for the electro-optical Pockels effect first hyperpolarizability. However, none of the functionals perform satisfactorily for all the vibrational (hyper)polarizabilities studied. In fact, in the case of electric field-induced second harmonic generation all of them, as well as the Hartree-Fock approximation, yield the wrong sign. We have also found that the Pople 6-31+G(d) basis set is unreliable for computing nuclear relaxation (hyper)polarizabilities of PMI oligomers due to the spurious prediction of a nonplanar equilibrium geometry.

Theoretical modelling of photoswitching of hyperpolarisabilities in ruthenium complexes.

Static excited-state polarisabilities and hyperpolarisabilities of three Ru(II) ammine complexes are computed at the density functional theory (DFT) and several correlated ab initio levels. Most accurate modelling of the low energy electronic absorption spectrum is obtained with the hybrid functionals B3LYP, B3P86 or M06 for the complex [Ru(II)(NH3)5(MeQ(+))](3+) (MeQ(+)=N-methyl-4,4'-bipyridinium, 3) in acetonitrile. The match with experimental data is less good for [Ru(II)(NH3)5L](3+) (L=N-methylpyrazinium, 2; N-methyl-4-{E,E-4-(4-pyridyl)buta-1,3-dienyl}pyridinium, 4). These calculations confirm that the first dipole- allowed excited state (FDAES) has metal-to-ligand charge-transfer (MLCT) character. Both the solution and gas-phase results obtained for 3 by using B3LYP, B3P86 or M06 are very similar to those from restricted active-space SCF second-order perturbation theory (RASPT2) with a very large basis set and large active space. However, the time-dependent DFT λ(max) predictions from the long-range corrected functionals CAM-B3LYP, LC-ωPBE and wB97XB and also the fully ab initio resolution of identity approximate coupled-cluster method (gas-phase only) are less accurate for all three complexes. The ground state (GS) two-state approximation first hyperpolarisability ß(2SA) for 3 from RASPT2 is very close to that derived experimentally via hyper-Rayleigh scattering, whereas the corresponding DFT-based values are considerably larger. The ß responses calculated by using B3LYP, B3P86 or M06 increase markedly as the π-conjugation extends on moving along the series 2→4, for both the GS and FDAES species. All three functionals predict substantial FDAES ß enhancements for each complex, increasing with the π-conjugation, up to about sevenfold for 4. Also, the computed second hyperpolarisabilities γ generally increase in the FDAES, but the results vary between the different functionals.

A contribution to the drug resistance mechanism of darunavir, amprenavir, indinavir, and saquinavir complexes with HIV-1 protease due to flap mutation I50V: a systematic MM-PBSA and thermodynamic integration study.

The emergence of HIV-1 drug-resistant mutations is the major problem against AIDS treatment. We employed molecular dynamics (MD) calculations and free energy (MM-PBSA and thermodynamic integration) analyses on wild-type (WT) and mutated HIV-1 protease (HIV-1 PR) complexes with darunavir, amprenavir, indinavir, and saquinavir to clarify the mechanism of resistance due to the I50V flap mutation. Conformational analysis showed that the protease flaps are increasingly flexible in the I50V complexes. In the WT, stabilization of the HIV-1 PR structure is achieved via interflap and water-mediated hydrogen-bonding interactions between the flaps. Furthermore, hydrogen bonds between drugs and binding cavity residues (Asp29/29'/30/30') are crucial for effective inhibition. All these interactions were significantly diminished (or absent) in the mutated forms, thus denoting their importance toward binding. Thermodynamic integration calculations reproduced the experimental data to within ≈1 kcal mol⁻¹ and showed that the I50V mutation results in weaker binding free energies for all analyzed complexes with respect to the WT. It was observed that the loss in binding energy upon mutation was mostly enthalpically driven in all complexes, with the greatest effect coming from the reduction of van der Waals interactions. Our results motivated us to test two novel compounds that have been synthesized to maximize interactions with HIV-1 PR. MM-PBSA and TI calculations showed that compound 3c (Ghosh et al. Bioorg. Med. Chem. Lett. 2012, 22, 2308) is a promising protease inhibitor, which presents very effective binding to the WT PR (ΔG(MM-PBSA) = -17.2 kcal mol⁻¹, ΔG(exp) = -16.1 kcal mol⁻¹). Upon I50V mutation, the complex binding free energy was weakened by a ΔΔG(TI) of 1.8 kcal mol⁻¹, comparable to the marketed inhibitors. This predicts that I50V may confer low resistance to 3c. This computational comparative study contributes toward elucidation of the I50V drug-resistance mechanism in HIV-1 PR.

Structural and static electric response properties of highly symmetric lithiated silicon cages: theoretical predictions.

It is shown by density functional theory calculations that high symmetry silicon cages can be designed by coating with Li atoms. The resulting highly symmetric lithiated silicon cages (up to D(5d) symmetry) are low-lying true minima of the energy hypersurface with binding energies of the order of 4.6 eV per Si atom and moderate highest occupied molecular orbital-lowest unoccupied molecular orbital gaps. Moreover, relying on a systematic study of the electric response properties obtained by ab initio (Hartree-Fock, MP2, and configuration interaction singles (CIS)) and density functional (B3LYP, B2PLYP, and CAM-B3LYP) methods, it is shown that lithium coating has a large impact on the magnitude of their second hyperpolarizabilities resulting to highly hyperpolarizable species. Such hyperpolarizable character is directly connected to the increase in the density of the low-lying excited states triggered by the interaction between the Si cage and the surrounding Li atoms.

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.

Conformational analysis of two novel cytotoxic C2-substituted pyrrolo[2,3-f]quinolines in aqueous media, organic solvents, membrane bilayers and at the putative active site.

We have performed: (i) conformational analysis of two novel cytotoxic C2-substituted pyrrolo[2,3-f]quinolines 5e and 5g in deuterated dimethylsulfoxide (DMSO-d(6)) utilizing NOE results from NMR spectroscopy; (ii) molecular dynamics (MD) calculations in water, DMSO and dimyristoyl phosphatidylcholine bilayers and (iii) molecular docking and MD calculations on DNA nucleotide sequences. The obtained results for the two similar in structure molecules showed differences in: (i) their conformational properties in silico and in media that reasonably simulate the biological environment; (ii) the way they are incorporated into the lipid bilayers and therefore their diffusion ability and (iii) molecular docking capacity as it is depicted from their different binding scores.

Binding of novel fullerene inhibitors to HIV-1 protease: insight through molecular dynamics and molecular mechanics Poisson-Boltzmann surface area calculations.

The objectives of this study include the design of a series of novel fullerene-based inhibitors for HIV-1 protease (HIV-1 PR), by employing two strategies that can also be applied to the design of inhibitors for any other target. Additionally, the interactions which contribute to the observed exceptionally high binding free energies were analyzed. In particular, we investigated: (1) hydrogen bonding (H-bond) interactions between specific fullerene derivatives and the protease, (2) the regions of HIV-1 PR that play a significant role in binding, (3) protease changes upon binding and (4) various contributions to the binding free energy, in order to identify the most significant of them. This study has been performed by employing a docking technique, two 3D-QSAR models, molecular dynamics (MD) simulations and the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method. Our computed binding free energies are in satisfactory agreement with the experimental results. The suitability of specific fullerene derivatives as drug candidates was further enhanced, after ADMET (absorption, distribution, metabolism, excretion and toxicity) properties have been estimated to be promising. The outcomes of this study revealed important protein-ligand interaction patterns that may lead towards the development of novel, potent HIV-1 PR inhibitors.

On the stability, electronic structure, and nonlinear optical properties of HXeOXeF and FXeOXeF.

The electronic ground state, stability, and linear and nonlinear optical properties of HXeOXeF and FXeOXeF have been studied theoretically by employing complete active space valence bond (CASVB), multistate complete active space perturbation theory (MS-CASPT2), and coupled cluster methods. It is shown that the oxygen inserted between the two Xe atoms significantly modifies the ground-state electronic configuration of the formed derivative by increasing the closed-shell contribution (σ(2)) and removing the diradicaloid character observed in HXe(2)F. The electronic charge distribution has been analyzed by employing the atoms-in-molecules (AIM) method. The dissociation channels of HXeOXeF and FXeOXeF have been studied in detail. It was found that these compounds are metastable, protected by substantial energy barriers and, thus, they can be prepared under appropriate conditions. Both two- and three-body dissociation reactions have been considered. The effects of inserting O in HXe(2)F and substituting H (HXeOXeF) by F, leading to FXeOXeF, on the energy barriers are discussed. The significant effects of the inserted oxygen on the polarizability and even more on the first hyperpolarizability have been computed and rationalized.