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.

Electrostatic and non-electrostatic contributions to the binding free energies of anthracycline antibiotics to DNA.

The knowledge about molecular factors driving simple ligand-DNA interactions is still limited. The aim of the present study was to investigate the electrostatic and non-electrostatic contributions to the binding free energies of anthracycline compounds with DNA. Theoretical calculations based on continuum methods (Poisson-Boltzmann and solvent accessible surface area) were performed to estimate the binding free energies of five selected anthracycline ligands (daunomycin, adriamycin, 9-deoxyadriamycin, hydroxyrubicin, and adriamycinone) to DNA. The free energy calculations also took into account the conformational change that DNA undergoes upon ligand binding. This conformational change appeared to be very important for estimating absolute free energies of binding. Our studies revealed that the absolute values of all computed contributions to the binding free energy were quite large compared to the total free energy of binding. However, the sum of these large positive and negative values produced a small negative value of the free energy around -10 kcal/mol. This value is in good agreement with experimental data. Experimental values for relative binding free energies were also reproduced for charged ligands by our calculations. Together, it was found that the driving force for ligand-DNA complex formation is the non-polar interaction between the ligand and DNA even if the ligand is positively charged.

Electrostatic effects in homeodomain-DNA interactions.

We report here an investigation of the role of electrostatics in homeodomain-DNA interactions using techniques based around the use of the Poisson-Boltzmann equation. In the present case such a study is of particular interest, since in contrast to other proteins previously studied with this method, the homeodomain is a small, highly charged protein that forms extensive ion pairs upon binding DNA. We have investigated the salt dependence of the binding constant for specific association and for a variety of models for non-specific association. The results indicate that, in line with the models proposed by Manning and Record, the entropy of counterion release accounts for a significant fraction of the salt dependence of the binding free energy, though this is perhaps due to fortuitous cancellation of other contributing terms. The thermodynamic effects of a number of specific homeodomain mutants were also investigated, and partly rationalized in terms of favorable electrostatic interactions in the major goove of DNA. Investigation of the temperature-dependence of the free energy of association indicates that the electrostatic contributions become increasingly favorable as the temperature rises. For this particular system, however, there appears to be no significant electrostatic contribution to the delta(delta C(p)) of association. Finally, an analysis of the free energy of interaction when the homeodomain is moved ca one Debye length from the DNA suggests that pure electrostatic forces are able to steer the homeodomain into a partially correct orientation for binding to the DNA.

Discovery of HIV-1 integrase inhibitors through a novel combination of ligand and structure-based drug design.

Over the past 10 years, classical computer-aided molecular design methods have not been frequently applied for the discovery of novel HIV-1 integrase (IN) inhibitors, due to the intrinsic challenges that this enzyme presents. Therefore, a novel approach that combines the chemical information of known integrase inhibitors with the enzyme's detailed 3D structure in a stepwise fashion is proposed: (I) use of a pharmacophore model (PM), which takes into account in a weighted fashion the chemical features of known ligands, in analogous manner to the to search the Maybridge and the NCI 3D databases; (II) drug-likeness optimization; (III) virtual high-throughput screening of the hits matching the PM query against 1QS4 wild-type IN structure using different Docking/Scoring combinations; (IV) visual inspection and selection of the hits in function of: binding free energies; binding mode type within the active site; retrieval among the best 20% hits in more than 6 Docking/Scoring protocols at the same time. This approach aims at a rational selection of new potential HIV-1 integrase inhibitors.

Computational science new horizons and relevance to pharmaceutical design.

Computer methods are used extensively in the design and refinement of drug leads. A short summary is given for several computational methods followed by a description of how some of these methods have been applied to design drugs targeted to the renin-angiotensin system and to cholinergic synapses. These methods include quantitative structure-activity relationship (QSAR) methods, comparative molecular field analyses (CoMFA), 3D database searching, de novo design of ligands, docking, and computational alchemy [free energy perturbation (FEP) and thermodynamic integration (MCTI)]. Most of these methods can be used whether or not detailed structural information about the binding site is available, although without an x-ray structure, the analyses are more qualitative. All of these methods are used extensively in the commercial design of pharmaceuticals. The main problem with most of these methods is in the scoring (ranking) of interactions or matches. Advances in this area and others (methods development and increases in capabilities of computers) will increase the predictive power of these methods and help to speed the time to market of new pharmaceuticals. (Trends Cardiovasc Med 1996;6:198-203).

Conservative and nonconservative mutations in proteins: anomalous mutations in a transport receptor analyzed by free energy and quantum chemical calculations.

Experimental studies on a bacterial sulfate receptor have indicated anomalous relative binding affinities for the mutations Ser130-->Cys,Ser130-->Gly, and Ser130-->Ala. The loss of affinity for sulfate in the former mutation was previously attributed to a greater steric effect on the part of the Cys side chain relative to the Ser side chain, whereas the relatively small loss of binding affinity for the latter two mutations was attributed to the loss of a single hydrogen bond. In this report we present quantum chemical and statistical thermodynamic studies of these mutations. Qualitative results from these studies indicate that for the Ser130-->Cys mutation the large decrease in binding affinity is in part caused by steric effects, but also significantly by the differential work required to polarize the Cys thiol group relative to the Ser hydroxyl group. The Gly mutant cobinds a water molecule in the same location as the Ser side chain resulting in a relatively small decrease in binding affinity. Results for the Ala mutant are in disagreement with experimental results but are likely to be limited by insufficient sampling of configuration space due to physical constraints applied during the simulation.

Calculation of the pKa values for the ligands and side chains of Escherichia coli D-alanine:D-alanine ligase.

Poisson-Boltzmann electrostatics methods have been used to calculate the pKa shifts for the ligands and titratable side chains of D-alanine:D-alanine ligase of the ddlb gene of Escherichia coli (DdlB). The focus of this study is to determine the ionization state of the second D-alanine (D-Ala2) in the active site of DdlB. The pKa of the amine is shifted over 5 pKa units more alkaline in the protein, clearly implying that D-Ala2 is bound to DdlB in its zwitterionic state and not in the free-base form as had been previously suggested. Comparisons are made to the depsipeptide ligase from the vancomycin-resistance cascade, VanA. It is suggested that VanA has different enzymatic properties due to a change in binding specificity rather than altered catalytic behavior and that the specificity of binding D-lactate over D-Ala2 may arise from the difference in ionization characteristics of the ligands.

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.

Theoretical study of inhibition of adenosine deaminase by (8R)-coformycin and (8R)-deoxycoformycin.

Molecular dynamics and free energy simulations were performed to examine the binding of (8R)-deoxycoformycin and (8R)-coformycin to adenosine deaminase. The two inhibitors differ only at the 2' position of the sugar ring; the sugar moiety of conformycin is ribose, while it is deoxyribose for deoxycoformycin. The 100 ps molecular dynamics trajectories reveal that Asp 19 and His 17 interact strongly with the 5' hydroxyl group of the sugar moiety of both inhibitors and appear to play an important role in binding the sugar. The 2' and 3' groups of the sugars are near the protein-water interface and can be stabilized by either protein residues or water. The flexibility of the residues at the opening of the active site helps to explain the modest difference in binding of the two inhibitors and how substrates/inhibitors can enter an otherwise inaccessible binding site.

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.

Controls of nitrogen limitation in tallgrass prairie.

The relationship between fire frequency and N limitation to foliage production in tallgrass prairie was studied with a series of fire and N addition experiments. Results indicated that fire history affected the magnitude of the vegetation response to fire and to N additions. Sites not burned for over 15 years averaged only a 9% increase in foliage biomass in response to N enrichment. In contrast, foliage production increased an average of 68% in response to N additions on annually burned sites, while infrequently burned sites, burned in the year of the study, averaged a 45% increase. These findings are consistent with reports indicating that reduced plant growth on unburned prairie is due to shading and lower soil temperatures, while foliage production on frequently burned areas is constrained by N availability. Infrequent burning of unfertilized prairie therefore results in a maximum production response in the year of burning relative to either annually burned or long-term unburned sites.Foliage biomass of tallgrass prairie is dominated by C4 grasses; however, forb species exhibited stronger production responses to nitrogen additions than did the grasses. After four years of annual N additions, forb biomass exceeded that of grass biomass on unburned plots, and grasses exhibited a negative response to fertilizer, probably due to competition from the forbs. The dominant C4 grasses may out-compete forbs under frequent fire conditions not only because they are better adapted to direct effects of burning, but because they can grow better under low available N regimes created by frequent fire.

Metabolic complications of HIV and AIDs.

Numerous anecdotal reports, clinical observations, and newly published studies provide evidence of the growing interest and concern about body shape changes and metabolic complications seen in HIV infection. At first believed to be a single complication caused by the protease inhibitor class of antiretroviral therapy, the focus of research is shifting to distinct syndromes with multiple causes. In this unfolding story, peripheral fat atrophy, central fat accumulation, dyslipidemia, and glucose disregulation characterize the more commonly recognized syndromes. Osteoporosis and osteopenia have been recently observed. While the etiologies await discovery, the long-term consequences of these metabolic changes demand the expertise of clinicians not formerly considered "front-line" in HIV/AIDS treatment, such as orthopaedic nurses.

Nursing management of symptoms influenced by HIV infection of the endocrine system.

Examining the impact of endocrine disregulation is useful in managing symptoms of HIV infection in its later stages. Current studies indicate insults to the endocrine system from a variety of agents present in HIV infection. Establishing direct pathologic, causal effects is difficult; however, the functional domains of the endocrine system provide a good framework for nursing management of the symptoms of later-stage HIV infection.

Comparative molecular field analysis (coMFA) study of epothilones-tubulin depolymerization inhibitors: pharmacophore development using 3D QSAR methods.

A three-dimensional quantitative structure-activity relationship (3D QSAR) study has been carried out on epothilones based on comparative molecular field analyses (CoMFA) using a large data set of epothilone analogs, which are potent inhibitors of tubulin depolymerization. Microtubules, which are polymers of the a/beta-tubulin heterodimer, need to dissociate in order to form the mitotic spindle, a structure required for cell division. A rational pharmacophore searching method using 3D QSAR procedures was carried out and the results for the epothilones are described herein. One-hundred and sixty-six epothilone analogs and their depolymerization inhibition properties with tubulin were used as a training set. Over a thousand molecular field energies were generated and applied to generate the descriptors of QSAR equations. Using a genetic function algorithm (GFA) method, combined with a least square approach, multiple QSAR models were considered during the search for pharmacophore elements. Each GFA run resulted in 100 QSAR models, which were ranked according to their lack of fit (LOF) scores, with a total of 40 GFA runs having been performed. The 40 best QSAR equations from each run had adequate fitted correlation coefficients (R from 0.813 to 0.863) and were of sufficient statistical significance (F value from 7.2 to 10.9). The pharmacophore elements for epothilones were studied by investigating the hit frequency of descriptors (i.e. the sampling probabilities of grid points from the GFA studies) from the set of the 4,000 top scoring QSAR equations. By comparing the frequency with which each grid point appeared in the QSAR equations, three candidate regions in the epothilones were proposed to be pharmacophore elements. Two of them are completely compatible with the recent model proposed by Ojima et al. [Proc. Natl. Acad. Sci. USA, 96 (1999) 4256], however, one is quite different and is necessary to accurately predict the activities of all 166 epothilone molecules used in our training set. Finally, by visualizing the 35 most probable grid points, it was found that changes related to the C6, C7, C8, C 12, S20, and C21 atoms of the epothilones were highly correlated to their activity.

Molecular modeling of the intrastrand guanine-guanine DNA adducts produced by cisplatin and oxaliplatin.

Intrastrand DNA adducts formed by cisplatin and oxaliplatin were modeled with molecular mechanics minimization and restrained molecular dynamics simulations in a comparative study. A reasonable set of force field parameters for the Pt atom were refined by using the available cisplatinated DNA crystal structure as a guide. This crystal structure was also used as the starting structure for the simulations. Analysis of the resulting structures indicated that the covalent effects of oxaliplatin coordination on DNA structure were very similar to those of cisplatin. The most prominent difference between the two structures resulted from the presence of the 1,2-diaminocyclohexane ring in the oxaliplatin adduct. The modeling indicated that this ring protrudes directly outward into, and fills much of, the narrowed major groove of the bound DNA, forming a markedly altered and less polar major groove in the area of the adduct. The differences in the structure of the adducts produced by cisplatin and oxaliplatin are consistent with the observation that they are differentially recognized by the DNA mismatch repair system.

Structure-based drug design: computational advances.

Structure-based computational methods continue to enhance progress in the discovery and refinement of therapeutic agents. Several such methods and their applications are described. These include molecular visualization and molecular modeling, docking, fragment methods, 3-D database techniques, and free-energy perturbation. Related issues that are discussed include the use of simplified potential energy functions and the determination of the positions of tightly bound waters. Strengths and weaknesses of the various methods are described.

Orientational steering in enzyme-substrate association: ionic strength dependence of hydrodynamic torque effects.

The effect of hydrodynamic torques on the association rate constants for enzyme-ligand complexation is investigated by Brownian dynamics simulations. Our hydrodynamic models of the enzyme and ligand are composed of spherical elements with friction forces acting at their centers. A quantitative measure of hydrodynamic torque orientational effects is introduced by choosing, as a reference system, an enzyme-ligand model with the same average hydrodynamic interactions but without orientational dependence. Our simple models show a 15% increase in the rate constant caused by hydrodynamic torques at physiological ionic strength. For more realistic hydrodynamic models, which are not computationally feasible at present, this effect is probably higher. The most important finding of this work is that hydrodynamic complementarity in shape (i.e. like the fitting together of pieces of a puzzle) is most effective for interactions between molecules at physiological ionic strength.

A model for enzyme-substrate interaction in alanine racemase.

We report on a theoretical model for the complex of the enzyme alanine racemase with its natural substrate (L-alanine) and cofactor (pyridoxal 5'-phosphate). Electrostatic potentials were calculated and ionization states were predicted for all of the ionizable groups in alanine racemase. Some rather unusual charge states were predicted for certain residues. Tyr265' has an unusually low predicted pK(a) of 7.9 and at pH 7.0 has a predicted average charge of -0.37, meaning that 37% of the Tyr265' residues in an ensemble of enzyme molecules are in the phenolate form. At pH 8-9, the majority of Tyr265' side groups will be in the phenolate form. This lends support to the experimental evidence that Tyr265' is the catalytic base involved in the conversion of L-alanine to D-alanine. Residues Lys39 and Lys129 have predicted average charges of +0.91 and +0.14, respectively, at pH 7.0. Lys39 is believed to be the catalytic base for the conversion of D-alanine to L-alanine, and the present results show that, at least some of the time, it is in the unprotonated amine form and thus able to act as a base. Cys311', which is located very close to the active site, has an unusually low predicted pK(a) of 5.8 and at pH 7.0 has a predicted average charge of -0.72. The very low predicted charge for Lys129 is consistent with experimental evidence that it is carbamylated, since an unprotonated amine group is available to act as a Lewis base and form the carbamate with CO(2). Repeating the pK(a) calculations on the enzyme with Lys129 in carbamylated form predicts trends similar to those of the uncarbamylated enzyme. It appears that the enzyme has the ability to stabilize negative charge in the region of the active site. Implications for selective inhibitor design are discussed.

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.