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.

A set of nearly 100 crystallographic structures was analyzed using ab initio methods in order to verify the effect of the conformational variability of Watson-Crick guanine-cytosine and adenine-thymine base pairs on the intermolecular interaction energy and its components. Furthermore, for the representative structures, a potential energy scan of the structural parameters describing mutual orientation of the base pairs was carried out. The results were obtained using the hybrid variational-perturbational interaction energy decomposition scheme. The electron correlation effects were estimated by means of the second-order Møller-Plesset perturbation theory and coupled clusters with singles and doubles method adopting AUG-cc-pVDZ basis set. Moreover, the characteristics of hydrogen bonds in complexes, mimicking those appearing in B-DNA, were evaluated using topological analysis of the electron density. Although the first-order electrostatic energy is usually the largest stabilizing component, it is canceled out by the associated exchange repulsion in majority of the studied crystallographic structures. Therefore, the analyzed complexes of the nucleic acid bases appeared to be stabilized mainly by the delocalization component of the intermolecular interaction energy which, in terms of symmetry adapted perturbation theory, encompasses the second- and higher-order induction and exchange-induction terms. Furthermore, it was found that the dispersion contribution, albeit much smaller in terms of magnitude, is also a vital stabilizing factor. It was also revealed that the intermolecular interaction energy and its components are strongly influenced by four (out of six) structural parameters describing mutual orientation of bases in Watson-Crick pairs, namely shear, stagger, stretch, and opening. Finally, as a part of a model study, much of the effort was devoted to an extensive testing of the UBDB databank. It was shown that the databank quite successfully reproduces the electrostatic energy determined with the aid of ab initio methods.

No abstract

The interaction of natural lipopeptide pseudofactin II with a series of doubly charged metal cations was examined by matrix-assisted laser-desorption ionization-time of flight (MALDI-TOF) mass spectrometry and molecular modelling. The molecular modelling for metal-pseudofactin II provides information on the metal-peptide binding sites. Overall, Mg(2+), Ca(2+) and Zn(2+) favor the association with oxygen atoms spanning the peptide backbone, whereas Cu(2+) is coordinated by three nitrogens. Circular dichroism (CD) results confirmed that Zn(2+) and Cu(2+) can disrupt the secondary structure of pseudofactin II at high concentrations, while Ca(2+) and Mg(2+) did not essentially affect the structure of the lipopeptide. Interestingly, our results showed that the addition of Zn(2+) and Cu(2+) helped smaller micelles to form larger micellar aggregates. Since pseudofactin II binds metals, we tested whether this phenomena was somehow related to its antimicrobial activity against Staphylococcus epidermidis and Proteus mirabilis. We found that the antimicrobial effect of pseudofactin II was increased by supplementation of culture media with all tested divalent metal ions. Finally, by using Gram-positive and Gram-negative bacteria we showed that the higher antimicrobial activity of metal complexes of pseudofactin II is attributed to the disruption of the cytoplasmic membrane.

No abstract

No abstract

The stacking interactions of two guanine molecules were analyzed detail at the DF-MP2/aug-cc-pVDZ level of theory for conformations appearing B-DNA. The dependence of intermolecular interaction energies on the pairs of step parameters (shift, slide, rise, tilt, roll and twist) was determined. The values of these parameters were chosen to cover the whole range of variability appearing crystallographic data. The scanning procedure was performed by subsequent changes of two variables with fixed values of the remaining base-pair and base-step BDNA parameters. Additionally, the hybrid variational-perturbational scheme was applied for the decomposition of the interaction energy into physically meaningful contributions at the MP2 level of theory. The significant impact of the mutual orientations of guanine bases was observed not only on the total intermolecular energy but also on its components. The second-order dispersion interaction is the most significant contribution to stabilization energy and is about eight times larger compared to the first-order electrostatic term with relaxation effects, which is also of stabilizing character. The dispersion interactions may vary up to 9.6 kcal mol(-1) between different guanine-guanine conformations. The parameters defining the mutual orientation of stacked guanine molecules have a different impact on the stabilization of the investigated complex. The following base-step parameters have only a minor impact on the stabilization energies: shift-slide, shift-roll, shift-twist, slide-twist and roll-twist. On the other hand, parameters such as rise and tilt significantly influence intermolecular interactions, i.e. strong attraction occurs only for a limited range of their values.

Twenty hydrogen-bonded complexes composed of nucleic acid base and amino acid side chain have been analyzed using ab initio quantum chemistry methods with the aim of gaining insights into the nature of molecular interactions in these systems. The intermolecular interaction energies were estimated using the second-order Møller-Plesset perturbation theory and coupled clusters approach with single and double excitations, while their components have been determined by means of a hybrid variational-perturbational decomposition scheme. Additionally, the topological analysis of an electron density distribution of the studied complexes has been performed. In the case of all of the studied neutral complexes, the main source of stabilization is the delocalizaction energy associated with the electron density deformation upon the interaction which contributes almost half of the total interaction energy. Furthermore, analysis of the interaction induced difference density maps of complexes containing neutral amino acid side chains reveals that the delocalization component involves the electron density changes localized in the double-hydrogen-bonded ring structures. A relatively good correlation between the sum of densities at hydrogen-bond critical points and the Hartree-Fock intermolecular interaction energy components (electrostatic, delocalization, and exchange) has been observed for the two independently considered sets of complexes, containing positively charged and neutral amino acid side chains.