Investigations on human immunodeficiency virus type 1 integrase/DNA binding interactions via molecular dynamics and electrostatics calculations.

The complete three-dimensional structure of the active site region of the human immunodeficiency virus type 1 (HIV-1) integrase (IN) is not unambiguously known. This region includes a flexible loop comprising residues 141-148 and the N-terminal portion of the helix alpha-4, which contains E152, the third catalytic residue, and Y143, which plays a secondary role in catalysis. Relatively high B-factors exist for most of the residues in the aforementioned region. The HIV-1 IN belongs to the polynucleotidyl transferase superfamily, whose members have been proposed to use two divalent metal ions for catalysis. Although only the position of the first metal ion has been determined crystallographically for the HIV-1 IN, we recently have proposed a binding site for the second metal ion. Based on this information, we have performed two 500-psec molecular dynamics simulations of the catalytic domain of the HIV-1 IN containing two Mg(2)+ ions. In one of the simulations, we included a dianionic phosphate group (HPO(4)(2)-) in the active site to mimic a portion of the DNA backbone of a substrate for the integration reaction. Electrostatics calculations and ionization state predictions were carried out on representative structures taken from the molecular dynamics simulations. Different conformational behaviors of the enzyme were observed, depending upon whether two Mg(2)+ ions were bound or two Mg(2)+ ions plus phosphate. The electrostatic calculations performed on the dynamical structures provide a further refinement about which regions of the catalytic domain of the HIV-1 IN may be involved in the DNA binding.

Developing a dynamic pharmacophore model for HIV-1 integrase.

We present the first receptor-based pharmacophore model for HIV-1 integrase. The development of "dynamic" pharmacophore models is a new method that accounts for the inherent flexibility of the active site and aims to reduce the entropic penalties associated with binding a ligand. Furthermore, this new drug discovery method overcomes the limitation of an incomplete crystal structure of the target protein. A molecular dynamics (MD) simulation describes the flexibility of the uncomplexed protein. Many conformational models of the protein are saved from the MD simulations and used in a series of multi-unit search for interacting conformers (MUSIC) simulations. MUSIC is a multiple-copy minimization method, available in the BOSS program; it is used to determine binding regions for probe molecules containing functional groups that complement the active site. All protein conformations from the MD are overlaid, and conserved binding regions for the probe molecules are identified. Those conserved binding regions define the dynamic pharmacophore model. Here, the dynamic model is compared to known inhibitors of the integrase as well as a three-point, ligand-based pharmacophore model from the literature. Also, a "static" pharmacophore model was determined in the standard fashion, using a single crystal structure. Inhibitors thought to bind in the active site of HIV-1 integrase fit the dynamic model but not the static model. Finally, we have identified a set of compounds from the Available Chemicals Directory that fit the dynamic pharmacophore model, and experimental testing of the compounds has confirmed several new inhibitors.

Brownian and essential dynamics studies of the HIV-1 integrase catalytic domain.

The three-dimensional structure of the active site region of the enzyme HIV-1 integrase is not unambiguously known. This region includes a flexible peptide loop that cannot be well resolved in crystallographic determinations. Here we present two different computational approaches with different levels of resolution and on different time-scales to understand this flexibility and to analyze the dynamics of this part of the protein. We have used molecular dynamics simulations with an atomic model to simulate the region in a realistic and reliable way for 1 ns. It is found that parts of the loop wind up after 300 ps to extend an existing helix. This indicates that the helix is longer than in the earlier crystal structures that were used as basis for this study. Very recent crystal data confirms this finding, underlining the predictive value of accurate MD simulations. Essential dynamics analysis of the MD trajectory yields an anharmonic motion of this loop. We have supplemented the MD data with a much lower resolution Brownian dynamics simulation of 600 ns length. It provides ideas about the slow-motion dynamics of the loop. It is found that the loop explores a conformational space much larger than in the MD trajectory, leading to a "gating"-like motion with respect to the active site.

Solvent accessibility to aspartyl and succinimidyl residues at positions 7 and 23 in the amyloid beta 1-28 peptide.

The water accessibilities to aspartyl residues at positions 7 and 23 in the amyloid beta 1-28 peptide associated with Alzheimer's Disease have been calculated using different techniques. These accessibilities of water were compared to those of the succinimidyl residues (SUC) replacing the aspartyl ones (ASP). It has been possible to ascertain that these modifications (ASP--->SUC) lead to a significant increase in the water accessibility to the backbone and alpha-carbon atom of the SUC7 and SUC23 residues. It is suggested that the spontaneous transformation of the ASP--->SUC might lead to an increase of the racemization rates due to the higher accessibility of water at these sites. It is also proposed that the behavior of the adjacent residues in the selectivity of the racemization is to control the water accessibility at the reactive residue.

Plural origins of molecular homochirality in our biota Part II. The relative stabilities of homochiral and mixed oligoribotides and peptides.

By computer simulations--molecular mechanics and molecular dynamics with the amber force field (Weiner et al, (1986), J. Comp. Chem. 7, 230-252)--we have determined the stabilities of oligoribotide strands built with D- and L-riboses, and of peptide chains with D- and L-amino acid residues. In particular, complementary double-chains of oligoribotides were studied, since they are an important feature of the growing mechanism of modern nucleic acids. Peptide chains on the other hand, grow without need of a template. We found that mixed oligoribotides are less stable than homochiral ones, and that this chiral effect is less noticeable in peptide chains. The results support the interpretation that L-riboses act as terminators to the template-assisted growth of oligo-r-GD (enantiometric cross-inhibition; Joyce et al., (1987), Proc. Natl. Acad. Sci. USA 84, 4398-4402). Based on this effect, a chemical pathway is proposed which could, under assumed prebiotic conditions, bypass the hindrance of homochiral growth.

Computer simulation of the rough lipopolysaccharide membrane of Pseudomonas aeruginosa.

Lipopolysaccharides (LPSs) form the major constituent of the outer membrane of Gram-negative bacteria, and are believed to play a key role in processes that govern microbial metal binding, microbial adsorption to mineral surfaces, and microbe-mediated oxidation/reduction reactions at the bacterial exterior surface. A computational modeling capability is being developed for the study of geochemical reactions at the outer bacterial envelope of Gram-negative bacteria. A molecular model for the rough LPS of Pseudomonas aeruginosa has been designed based on experimentally determined structural information. An electrostatic model was developed based on Hartree-Fock SCF calculations of the complete LPS molecule to obtain partial atomic charges. The exterior of the bacterial membrane was assembled by replication of a single LPS molecule and a single phospholipid molecule. Molecular dynamics simulations of the rough LPS membrane of P. aeruginosa were carried out and trajectories were analyzed for the energetic and structural factors that determine the role of LPS in processes at the cell surface.

Similarities in the HIV-1 and ASV integrase active sites upon metal cofactor binding.

The HIV-1 integrase, which is essential for viral replication, catalyzes the insertion of viral DNA into the host chromosome thereby recruiting host cell machinery into making viral proteins. It represents the third main HIV enzyme target for inhibitor design, the first two being the reverse transcriptase and the protease. We report here a fully hydrated 2 ns molecular dynamics simulation performed using parallel NWChem3.2.1 with the AMBER95 force field. The HIV-1 integrase catalytic domain previously determined by crystallography (1B9D) and modeling including two Mg(2+) ions placed into the active site based on an alignment against an ASV integrase structure containing two divalent metals (1VSH), was used as the starting structure. The simulation reveals a high degree of flexibility in the region of residues 140-149 even in the presence of a second divalent metal ion and a dramatic conformational change of the side chain of E152 when the second metal ion is present. This study shows similarities in the behavior of the catalytic residues in the HIV-1 and ASV integrases upon metal binding. The present simulation also provides support to the hypothesis that the second metal ion is likely to be carried into the HIV-1 integrase active site by the substrate, a strand of DNA.

Internal dynamics and ionization states of the macrophage migration inhibitory factor: comparison between wild-type and mutant forms.

The macrophage migration inhibitory factor (MIF) is a cytokine that shares a common structural architecture and catalytic strategy with three isomerases: 4-oxalocrotonate tautomerase, 5-carboxymethyl-2-hydroxymuconate isomerase, and D-dopachrome tautomerase. A highly conserved N-terminal proline acts as a base-acid during the proton transfer reaction catalyzed by these enzymes. Such unusual catalytic strategy appears to be possible only due to the N-terminal proline pK(a) shifted to 5.0-6.0 units. Mutations of this residue result in a significant decrease of the catalytic activity of MIF. Two hypotheses have been proposed to explain the catalytic inefficiency of MIF: the lower basicity of primary amines with regard to secondary ones and the increased flexibility resulting from the replacement of a proline by residues like glycine. To investigate that, we have performed molecular dynamics simulations of MIF wild-type and its mutant P1G, as well as calculated the protonation properties of several mutant forms. It was found that the N-terminal glycine does not show larger fluctuations compared to proline, but the former residue is more exposed to the solvent throughout the simulations. The apparent pK(a) of these residues displays very little change (as expected from the structural rigidity of MIF) and is not significantly affected by the surrounding ionizable residues. Instead, the hydrophobic character of the active site seems to be the main factor in determining the pKa of the N-terminal residue and the catalytic efficiency of MIF.

Molecular dynamics studies on the HIV-1 integrase catalytic domain.

The HIV-1 integrase, which is essential for viral replication, catalyzes the insertion of viral DNA into the host chromosome, thereby recruiting host cell machinery into making viral proteins. It represents the third main HIV enzyme target for inhibitor design, the first two being the reverse transcriptase and the protease. Two 1-ns molecular dynamics simulations have been carried out on completely hydrated models of the HIV-1 integrase catalytic domain, one with no metal ions and another with one magnesium ion in the catalytic site. The simulations predict that the region of the active site that is missing in the published crystal structures has (at the time of this work) more secondary structure than previously thought. The flexibility of this region has been discussed with respect to the mechanistic function of the enzyme. The results of these simulations will be used as part of inhibitor design projects directed against the catalytic domain of the enzyme.