Molecular dynamics studies on the HIV-1 integrase catalytic domain

Analogue and structure based approaches for modelling HIV-1 integrase inhibitors

A set of 220 inhibitors belonging to different structure classes and having HIV-1 integrase activity were collected along with their experimental pIC50 values. Geometries of all the inhibitors were fully optimized using B3LYP/6-31 + G(d) level of theory. These ligands were docked against 4 different HIV-1 integrase receptors (PDB IDs: 4LH5, 5KRS, 3ZSQ and 3ZSV). 30 docked poses were generated for all 220 inhibitors and ligand interaction of the first docked pose and the docked pose with the highest score were analysed. Residue GLU170 of 4LH5 receptor shows the highest number of interactions followed by ALA169, GLN168, HIS171 and ASP167 residues. Hydrogen bonding and stacking are mainly responsible for the interactions of these inhibitors with the receptor. We performed Molecular Dynamics (MD) simulation to observe the root-mean-square deviation (RMSD), for measure the average change of displacement between the atoms for a particular frame with respect to a reference and The Root Mean Square Fluctuation (RMSF) for characterization of local changes along the protein chain of the docked complexes. Analogue based models were generated to predict the pIC50 values for integrase inhibitors using various types of descriptors such as constitutional, geometrical, topological, quantum chemical and docking based descriptors. The best models were selected on the basis of statistical parameters and were validated by training and test set division. A few new inhibitors were designed on the basis of structure activity relationship and their pIC50 values were predicted using the generated models. All the designed new inhibitors a very high potential and may be used as potent inhibitors of HIV integrase. These models may be useful for further design and development of new and potent HIV integrase inhibitors.Communicated by Ramaswamy H. Sarma.

Rigidity and flexibility characteristics of DD[E/D]-transposases Mos1 and Sleeping Beauty

DD[E/D]-transposases catalyze the multistep reaction of cut-and-paste DNA transposition. Structurally, several DD[E/D]-transposases have been characterized, revealing a multi-domain structure with the catalytic domain possessing the RNase H-like structural motif that brings three catalytic residues (D, D, and E or D) into close proximity for the catalysis. However, the dynamic behavior of DD[E/D]-transposases during transposition remains poorly understood. Here, we analyze the rigidity and flexibility characteristics of two representative DD[E/D]-transposases Mos1 and Sleeping Beauty (SB) using the minimal distance constraint model (mDCM). We find that the catalytic domain of both transposases is globally rigid, with the notable exception of the clamp loop being flexible in the DNA-unbound form. Within this globally rigid structure, the central β-sheet of the RNase H-like motif is much less rigid in comparison to its surrounding α-helices, forming a cage-like structure. The comparison of the original SB transposase to its hyperactive version SB100X reveals the region where the change in flexibility/rigidity correlates with increased activity. This region is found to be within the RNase H-like structural motif and comprise the loop leading from beta-strand B3 to helix H1, helices H1 and H2, which are located on the same side of the central beta-sheet, and the loop between helix H3 and beta-strand B5. We further identify the RKEN214-217DAVQ mutations of the set of hyperactive mutations within the catalytic domain of SB transposase to be the driving factor that induces change in residue-pair rigidity correlations within SB transposase. Given that a signature RNase H-like structural motif is found in DD[E/D]-transposases and, more broadly, in a large superfamily of polynucleotidyl transferases, our results are relevant to these proteins as well.

The conformational feasibility for the formation of reaching dimer in ASV and HIV integrase: a molecular dynamics study

Retroviral integrases are reported to form alternate dimer assemblies like the core-core dimer and reaching dimer. The core-core dimer is stabilized predominantly by an extensive interface between two catalytic core domains. The reaching dimer is stabilized by N-terminal domains that reach to form intermolecular interfaces with the other subunit's core and C-terminal domains (CTD), as well as CTD-CTD interactions. In this study, molecular dynamics (MD), Brownian dynamics (BD) simulations, and free energy analyses, were performed to elucidate determinants for the stability of the reaching dimer forms of full-length Avian Sarcoma Virus (ASV) and Human Immunodeficiency Virus (HIV) IN, and to examine the role of the C-tails (the last ~16-18 residues at the C-termini) in their structural dynamics. The dynamics of an HIV reaching dimer derived from small angle X-ray scattering and protein crosslinking data, was compared with the dynamics of a core-core dimer model derived from combining the crystal structures of two-domain fragments. The results showed that the core domains in the ASV reaching dimer express free dynamics, whereas those in the HIV reaching dimer are highly stable. BD simulations suggest a higher rate of association for the HIV core-core dimer than the reaching dimer. The predicted stability of these dimers was therefore ranked in the following order: ASV reaching dimer < HIV reaching dimer < composite core-core dimer. Analyses of MD trajectories have suggested residues that are critical for intermolecular contacts in each reaching dimer. Tests of these predictions and insights gained from these analyses could reveal a potential pathway for the association and dissociation of full-length IN multimers.

A targeted DNA substrate mechanism for the inhibition of HIV-1 integrase by inhibitors with antiretroviral activity

We recently reported that viral DNA could be the primary target of raltegravir (RAL), an efficient anti‐HIV‐1 drug, which acts by inhibiting integrase. To elucidate this mechanism, we conducted a comparative analysis of RAL and TB11, a diketoacid abandoned as an anti‐HIV‐1 drug for its weak efficiency and marked toxicity, and tested the effects of the catalytic cofactor Mg²⁺ (5 mm) on drug‐binding properties. We used circular dichroism and fluorescence to determine drug affinities for viral DNA long terminal repeats (LTRs) and peptides derived from the integrase active site and DNA retardation assays to assess drug intercalation into DNA base pairs. We found that RAL bound more tightly to LTR ends than did TB11 (a diketo acid bearing an azido group) and that Mg²⁺ significantly increased the affinity of both RAL and TB11. We also observed a good relationship between drug binding with processed LTR and strand transfer inhibition. This unusual type of inhibition was caused by Mg²⁺‐assisted binding of drugs to DNA substrate, rather than to enzyme. Notably, while RAL bound exclusively to the cleavable/cleaved site, TB11 further intercalated into DNA base pairs and interacted with the integrase‐derived peptides. These unwanted binding sites explain the weaker bioavailability and higher toxicity of TB11 compared with the more effective RAL.

Assessing protein conformational sampling and structural stability via de novo design and molecular dynamics simulations

Molecular dynamics and de novo techniques, associated to quality parameter sets, have excelled at determining the structure of small proteins with high accuracy. To achieve a detailed description of protein conformations, these methods must critically assess the thermodynamic features of the molecular ensembles. Here, a comparison of the conformational ensemble generated by molecular dynamics and de novo techniques were carried out for six Top7-based proteins carrying gp41 HIV-1 epitopes. The native Top7, a highly stable computationally designed protein, was used as benchmark. Structural stability, flexibility, and secondary structure content were assessed. The consistency of the latter was compared to experimental circular dichroism spectra for all proteins. While both methods are capable to identify the stable from unstable chimeric proteins, the sampled conformational space and flexibility differ significantly in both methods. Molecular dynamics simulations seem to better describe secondary structure content and identify regions responsible for conformational instability. The de novo method, as implemented in Rosetta-a prime tool for protein design, overestimates secondary structure content. On the other hand, its empirical energy function is capable to predict the threshold for protein stability.

Study of the Conformational Dynamics of the Catalytic Loop of WT and G140A/G149A HIV-1 Integrase Core Domain Using Reversible Digitally Filtered Molecular Dynamics

The HIV-1 IN enzyme is one of three crucial virally encoded enzymes (HIV-1 IN, HIV-1 PR, and HIV-1 RT) involved in the life-cycle of the HIV-1 virus, making it an attractive target in the development of drugs against the AIDS virus. The structure and mechanism of the HIV-1 IN enzyme is the least understood of the three enzymes due to the lack of three-dimensional structural information. X-ray cystallographic studies have not yet been able to resolve the full-length structure, and studies have been mainly focused on the catalytic domain. This central domain possesses an important catalytic loop observed to overhang the active site, and experimental studies have shown that its dynamics affects the catalytic activity of mutant HIV-1 IN enzymes. In this study, the enhanced sampling technique, Reversible Digitally Filtered Molecular Dynamics (RDFMD), has been applied to the catalytic domain of the WT and G140A/G149A HIV-1 IN enzymes and has highlighted significant differences between the behavior of the catalytic loop which may explain the decrease of activity observed in experimental studies for this mutant.

Structural dynamics of native and V260E mutant C-terminal domain of HIV-1 integrase

The C-terminal domain (CTD) of HIV-1 integrase is a five stranded β-barrel resembling an SH3 fold. Mutational studies on isolated CTD and full-length IN have reported V260E mutant as either homo-dimerization defective or affecting the stability and folding of CTD. In this study, molecular dynamics simulation techniques were used to unveil the effect of V260E mutation on isolated CTD monomer and dimer. Both monomeric and dimeric forms of wild type and V260E mutant are highly stable during the simulated period. However, the stabilizing π-stacking interaction between Trp243 and Trp243' at the dimer interface is highly disturbed in CTD-V260E (>6 Å apart). The loss in entropy for dimerization is -30 and -25 kcal/mol for CTD-wt and CTD-V260E respectively signifying a weak hydrophobic interaction and its perturbation in CTD-V260E. The mutant Glu260 exhibits strong attraction/repulsion with all the basic/acidic residues of CTD. In addition to this, the dynamics of CTD-wild type and V260E monomers at 498 K was analyzed to elucidate the effect of V260E mutation on CTD folding. Increase in SASA and reduction in the number of contacts in CTD-V260E during simulation highlights the instability caused by the mutation. In general, V260E mutation affects both multimerization and protein folding with a pronounced effect on protein folding rather than multimerization. This study emphasizes the importance of the hydrophobic nature and SH3 fold of CTD in proper functioning of HIV integrase and perturbing this nature would be a rational approach toward designing more selective and potent allosteric anti-HIV inhibitors.

Predicting protein-DNA interactions by full search computational docking

Protein-DNA interactions are essential for many biological processes. X-ray crystallography can provide high-resolution structures, but protein-DNA complexes are difficult to crystallize and typically contain only small DNA fragments. Thus, there is a need for computational methods that can provide useful predictions to give insights into mechanisms and guide the design of new experiments. We used the program DOT, which performs an exhaustive, rigid-body search between two macromolecules, to investigate four diverse protein-DNA interactions. Here, we compare our computational results with subsequent experimental data on related systems. In all cases, the experimental data strongly supported our structural hypotheses from the docking calculations: a mechanism for weak, nonsequence-specific DNA binding by a transcription factor, a large DNA-binding footprint on the surface of the DNA-repair enzyme uracil-DNA glycosylase (UNG), viral and host DNA-binding sites on the catalytic domain of HIV integrase, and a three-DNA-contact model of the linker histone bound to the nucleosome. In the case of UNG, the experimental design was based on the DNA-binding surface found by docking, rather than the much smaller surface observed in the crystallographic structure. These comparisons demonstrate that the DOT electrostatic energy gives a good representation of the distinctive electrostatic properties of DNA and DNA-binding proteins. The large, favourably ranked clusters resulting from the dockings identify active sites, map out large DNA-binding sites, and reveal multiple DNA contacts with a protein. Thus, computational docking can not only help to identify protein-DNA interactions in the absence of a crystal structure, but also expand structural understanding beyond known crystallographic structures.

Computational design of a full-length model of HIV-1 integrase: modeling of new inhibitors and comparison of their calculated binding energies with those previously studied

A full-length model of integrase (IN) of the human immunodeficiency virus type 1 (HIV-1) was constructed based on the distinctly resolved X-ray crystal structures of its three domains, named N-terminal, catalytic core and C-terminal. Thirty-one already known inhibitors with varieties of structural differences as well as nine newly tested ones were docked into the catalytic core. The molecular dynamic (MD) and binding properties of these complexes were obtained by MD calculations. The binding energies calculated by molecular mechanic/Poisson Boltzmann solvation area were significantly correlationed with available IC50. Four inhibitors including two newly designed were also docked into the full-length model and their MD behaviors and binding properties were calculated. It was found that one of the newly designed compounds forms a better complex with HIV-1 IN compared to the rest including raltegravir. MD calculations were performed with AMBER suite of programs using ff99SB force field for the proteins and the general Amber force field for the ligands. In conclusion, the results have produced a promising standpoint not only in the construction of the full-length model but also in development of new drugs against it. However, the role of multimer formation and the involvement of DNAs, and their subsequent effect on the complexation and inhibition, are required to arrive at a conclusive decision.

Emerging methods for ensemble-based virtual screening

Ensemble based virtual screening refers to the use of conformational ensembles from crystal structures, NMR studies or molecular dynamics simulations. It has gained greater acceptance as advances in the theoretical framework, computational algorithms, and software packages enable simulations at longer time scales. Here we focus on the use of computationally generated conformational ensembles and emerging methods that use these ensembles for discovery, such as the Relaxed Complex Scheme or Dynamic Pharmacophore Model. We also discuss the more rigorous physics-based computational techniques such as accelerated molecular dynamics and thermodynamic integration and their applications in improving conformational sampling or the ranking of virtual screening hits. Finally, technological advances that will help make virtual screening tools more accessible to a wider audience in computer aided drug design are discussed.

Communications: Electron polarization critically stabilizes the Mg2+ complex in the catalytic core domain of HIV-1 integrase

In this paper, we present a detailed dynamics study of the catalytic core domain (CCD) of HIV-1 integrase using both polarized and nonpolarized force fields. The numerical results reveal the critical role of protein polarization in stabilizing Mg(2+) coordination complex in CCD. Specifically, when nonpolarized force field is used, a remarkable drift of the Mg(2+) complex away from its equilibrium position is observed, which causes the binding site blocked by the Mg(2+) complex. In contrast, when polarized force field is employed in MD simulation, HIV-1 integrase CCD structure is stabilized and both the position of the Mg(2+) complex and the binding site are well preserved. The detailed analysis shows the transition of alpha-helix to 3(10)-helix adjacent to the catalytic loop (residues 139-147), which correlates with the dislocation of the Mg(2+) complex. The current study demonstrates the importance of electronic polarization of protein in stabilizing the metal complex in the catalytic core domain of HIV-1 integrase.

Computer tools in the discovery of HIV-1 integrase inhibitors

Computer-aided drug design (CADD) methodologies have made great advances and contributed significantly to the discovery and/or optimization of many clinically used drugs in recent years. CADD tools have likewise been applied to the discovery of inhibitors of HIV-1 integrase, a difficult and worthwhile target for the development of efficient anti-HIV drugs. This article reviews the application of CADD tools, including pharmacophore search, quantitative structure-activity relationships, model building of integrase complexed with viral DNA and quantum-chemical studies in the discovery of HIV-1 integrase inhibitors. Different structurally diverse integrase inhibitors have been identified by, or with significant help from, various CADD tools.

Raltegravir: molecular basis of its mechanism of action

Integration of the HIV-1 viral DNA generated by reverse transcription of the RNA genome into the host cell chromosomes is a key step of viral replication, catalyzed by the viral integrase. In October 2007, the first integrase inhibitor, raltegravir, was approved for clinical use under the name of Isentress™. The results of the various clinical trials that have evaluated raltegravir have been very encouraging with regard to the immunological and virological efficacy and tolerance. However, as observed for other anti-retrovirals, specific resistance mutations have been identified in patients failing to respond to treatment with raltegravir. Although knowledge of the integrase structural biology remains fragmentary, the structures and modeling data available might provide relevant clues on the origin of the emergence of these resistance mutations. In this review, we describe the mechanism of action of this drug and the main data relating to its use in vivo, together with recent structural data important to our understanding of the origin of viral resistance.

Solution conformation and dynamics of the HIV-1 integrase core domain

The human immunodeficiency virus type 1 (HIV-1) integrase (IN) is a critical enzyme involved in infection. It catalyzes two reactions to integrate the viral cDNA into the host genome, 3' processing and strand transfer, but the dynamic behavior of the active site during catalysis of these two processes remains poorly characterized. NMR spectroscopy can reveal important structural details about enzyme mechanisms, but to date the IN catalytic core domain has proven resistant to such an analysis. Here, we present the first NMR studies of a soluble variant of the catalytic core domain. The NMR chemical shifts are found to corroborate structures observed in crystals, and confirm prior studies suggesting that the alpha4 helix extends toward the active site. We also observe a dramatic improvement in NMR spectra with increasing MgCl(2) concentration. This improvement suggests a structural transition not only near the active site residues but also throughout the entire molecule as IN binds Mg(2+). In particular, the stability of the core domain is linked to the conformation of its C-terminal helix, which has implications for relative domain orientation in the full-length enzyme. (15)N relaxation experiments further show that, although conformationally flexible, the catalytic loop of IN is not fully disordered in the absence of DNA. Indeed, automated chemical shift-based modeling of the active site loop reveals several stable clusters that show striking similarity to a recent crystal structure of prototype foamy virus IN bound to DNA.

HIV-1 IN alternative molecular recognition of DNA induced by raltegravir resistance mutations

Virologic failure during treatment with raltegravir, the first effective drug targeting HIV integrase, is associated with two exclusive pathways involving either Q148H/R/K, G140S/A or N155H mutations. We carried out a detailed analysis of the molecular and structural effects of these mutations. We observed no topological change in the integrase core domain, with conservation of a newly identified Omega-shaped hairpin containing the Q148 residue, in particular. In contrast, the mutations greatly altered the specificity of DNA recognition by integrase. The native residues displayed a clear preference for adenine, whereas the mutant residues strongly favored pyrimidines. Raltegravir may bind to N155 and/or Q148 residues as an adenine bioisoster. This may account for the selected mutations impairing raltegravir binding while allowing alternative DNA recognition by integrase. This study opens up new opportunities for the design of integrase inhibitors active against raltegravir-resistant viruses.

In pursuit of virtual lead optimization: pruning ensembles of receptor structures for increased efficiency and accuracy during docking

Representing receptors as ensembles of protein conformations during docking is a powerful method to approximate protein flexibility and increase the accuracy of the resulting ranked list of compounds. Unfortunately, docking compounds against a large number of ensemble members can increase computational cost and time investment. In this article, we present an efficient method to evaluate and select the most contributive ensemble members prior to docking for targets with a conserved core of residues that bind a ligand moiety. We observed that ensemble members that preserve the geometry of the active site core are most likely to place ligands in the active site with a conserved orientation, generally rank ligands correctly and increase interactions with the receptor. A relative distance approach is used to quantify the preservation of the three-dimensional interatomic distances of the conserved ligand-binding atoms and prune large ensembles quickly. In this study, we investigate dihydrofolate reductase as an example of a protein with a conserved core; however, this method for accurately selecting relevant ensemble members a priori can be applied to any system with a conserved ligand-binding core, including HIV-1 protease, kinases, and acetylcholinesterase. Representing a drug target as a pruned ensemble during in silico screening should increase the accuracy and efficiency of high-throughput analyses of lead analogs.

HIV-1 integrase inhibitors: 2005-2006 update

HIV-1 integrase (IN) catalyzes the integration of proviral DNA into the host genome, an essential step for viral replication. Inhibition of IN catalytic activity provides an attractive strategy for antiretroviral drug design. Currently two IN inhibitors, MK-0518 and GS-9137, are in advanced stages of human clinical trials. The IN inhibitors in clinical evaluation demonstrate excellent antiretroviral efficacy alone or in combination regimens as compared to previously used clinical antiretroviral agents in naive and treatment-experienced HIV-1 infected patients. However, the emergence of viral strains resistant to clinically studied IN inhibitors and the dynamic nature of the HIV-1 genome demand a continued effort toward the discovery of novel inhibitors to keep a therapeutic advantage over the virus. Continued efforts in the field have resulted in the discovery of compounds from diverse chemical classes. In this review, we provide a comprehensive report of all IN inhibitors discovered in the years 2005 and 2006.

Changes to the HIV long terminal repeat and to HIV integrase differentially impact HIV integrase assembly, activity, and the binding of strand transfer inhibitors

Human immunodeficiency virus (HIV) integrase enzyme is required for the integration of viral DNA into the host cell chromosome. Integrase complex assembly and subsequent strand transfer catalysis are mediated by specific interactions between integrase and bases at the end of the viral long terminal repeat (LTR). The strand transfer reaction can be blocked by the action of small molecule inhibitors, thought to bind in the vicinity of the viral LTR termini. This study examines the contributions of the terminal four bases of the nonprocessed strand (G(2)T(1)C(-1)A(-2)) of the HIV LTR on complex assembly, specific strand transfer activity, and inhibitor binding. Base substitutions and abasic replacements at the LTR terminus provided a means to probe the importance of each nucleotide on the different functions. An approach is described wherein the specific strand transfer activity for each integrase/LTR variant is derived by normalizing strand transfer activity to the concentration of active sites. The key findings of this study are as follows. 1) The G(2):C(2) base pair is necessary for efficient assembly of the complex and for maintenance of an active site architecture, which has high affinity for strand transfer inhibitors. 2) Inhibitor-resistant enzymes exhibit greatly increased sensitivity to LTR changes. 3) The strand transfer and inhibitor binding defects of a Q148R mutant are due to a decreased affinity of the complex for magnesium. 4) Gln(148) interacts with G(2), T(1), and C(-1) at the 5' end of the viral LTR, with these four determinants playing important and overlapping roles in assembly, strand transfer catalysis and high affinity inhibitor binding.

Hybrid quantum mechanical/molecular mechanical molecular dynamics simulations of HIV-1 integrase/inhibitor complexes

Human immunodeficiency virus (HIV)-1 integrase (IN) is an attractive target for development of acquired immunodeficiency syndrome chemotherapy. In this study, conventional and coupled quantum mechanical and molecular mechanical (QM/MM) molecular dynamics (MD) simulations of HIV-1 IN complexed with 5CITEP (IN-5CITEP) were carried out. In addition to differences in the bound position of 5CITEP, significant differences at the two levels of theory were observed in the metal coordination geometry and the areas involving residues 116-119 and 140-166. In the conventional MD simulation, the coordination of Mg(2+) was found to be a near-perfect octahedral geometry whereas a distorted octahedral complex was observed in QM/MM. All of the above reasons lead to a different pattern of protein-ligand salt link formation that was not observed in the classical MD simulation. Furthermore to provide a theoretical understanding of inhibition mechanisms of 5CITEP and its derivative (DKA), hybrid QM/MM MD simulations of the two complexes (IN-5CITEP and IN-DKA) have been performed. The results reveal that areas involving residues 60-68, 116-119, and 140-149 were substantially different among the two systems. The two systems show similar pattern of metal coordination geometry, i.e., a distorted octahedron. In IN-DKA, both OD1 and OD2 of Asp-64 coordinate the Mg(2+) in a monodentate fashion whereas only OD1 is chelated to the metal as observed in IN-5CITEP. The high potency of DKA as compared to 5CITEP is supported by a strong salt link formed between its carboxylate moiety and the ammonium group of Lys-159. Detailed comparisons between HIV-1 IN complexed with DKA and with 5CITEP provide information about ligand structure effects on protein-ligand interactions in particular with the Lys-159. This is useful for the design of new selective HIV-1 IN inhibitors.

A comparative molecular dynamics study on thermostability of human and chicken prion proteins

To compare the thermostabilities of human and chicken normal cellular prion proteins (HuPrP(C) and CkPrP(C)), molecular dynamics (MD) simulations were performed for both proteins at an ensemble level (10 parallel simulations at 400 K and 5 parallel simulations at 300 K as a control). It is found that the thermostability of HuPrP(C) is comparable with that of CkPrP(C), which implicates that the non-occurrence of prion diseases in non-mammals cannot be completely attributed to the thermodynamic properties of non-mammalian PrP(C).

Structural and Dynamical Properties of a Full-length HIV-1 Integrase: Molecular Dynamics Simulations

The structural and dynamical properties of the complete full-length structure of HIV-1 integrase were investigated using Molecular Dynamics approach. Simulations were carried out for the three systems, core domain only (CORE), full-length structure without (FULL) and with a Mg2+ (FULL+ION) in its active site, aimed to investigate the difference in the molecular properties of the full-length models due to their different construction procedures as well as the effects of the two ends, C- and N-terminal, on those properties in the core domain. The full-length structure was prepared from the two experimental structures of two-domain fragment. The following properties were observed to differ significantly from the previous reports: (i) relative topology formed by an angle between the three domains; (ii) the cavity size defined by the catalytic triad, Asp64, Asp116, and Glu152; (iii) distances and solvation of the Mg2+; and (iv) conformation of the catalytic residues. In addition, the presence of the two terminal domains decreases the mobility of the central core domain significantly.

Novel bifunctional quinolonyl diketo acid derivatives as HIV-1 integrase inhibitors: design, synthesis, biological activities, and mechanism of action

The virally encoded integrase protein is an essential enzyme in the life cycle of the HIV-1 virus and represents an attractive and validated target in the development of therapeutics against HIV infection. Drugs that selectively inhibit this enzyme, when used in combination with inhibitors of reverse transcriptase and protease, are believed to be highly effective in suppressing the viral replication. Among the HIV-1 integrase inhibitors, the beta-diketo acids (DKAs) represent a major lead for anti-HIV-1 drug development. In this study, novel bifunctional quinolonyl diketo acid derivatives were designed, synthesized, and tested for their inhibitory ability against HIV-1 integrase. The compounds are potent inhibitors of integrase activity. Particularly, derivative 8 is a potent IN inhibitor for both steps of the reaction (3'-processing and strand transfer) and exhibits both high antiviral activity against HIV-1 infected cells and low cytotoxicity. Molecular modeling studies provide a plausible mechanism of action, which is consistent with ligand SARs and enzyme photo-cross-linking experiments.

Molecular basis for the Cu2+ binding-induced destabilization of beta2-microglobulin revealed by molecular dynamics simulation

According to experimental data, binding of the Cu(2+) ions destabilizes the native state of beta2-microglobulin (beta2m). The partial unfolding of the protein was generally considered an early step toward fibril formation in dialysis-related amyloidosis. Recent NMR studies have suggested that the destabilization of the protein might be achieved through increased flexibility upon Cu(2+) binding. However, the molecular mechanism of destabilization due to Cu(2+), its role in amyloid formation, and the relative contributions of different potential copper-binding sites remain unclear. To elucidate the effect of ion ligation at atomic detail, a series of molecular dynamics simulations were carried out on apo- and Cu(2+)-beta2m systems in explicit aqueous solutions, with varying numbers of bound ions. Simulations at elevated temperatures (360 K) provide detailed pictures for the process of Cu(2+)-binding-induced destabilization of the native structure at the nanosecond timescale, which are in agreement with experiments. Conformational transitions toward partially unfolded states were observed in protein solutions containing bound copper ions at His-31 and His-51, which is marked by an increase in the protein vibrational entropy, with TDeltaS(vibr) ranging from 30 to 69 kcal/mol. The binding of Cu(2+) perturbs the secondary structure and the hydrogen bonding pattern disrupts the native hydrophobic contacts in the neighboring segments, which include the beta-strand D2 and part of the beta-strand E, B, and C and results in greater exposure of the D-E loop and the B-C loop to the water environment. Analysis of the MD trajectories suggests that the changes in the hydrophobic environment near the copper-binding sites lower the barrier of conformational transition and stabilize the more disordered conformation. The results also indicate that the binding of Cu(2+) at His-13 has little effect on the conformational stability, whereas the copper-binding site His-31, and to a lesser extent His-51, are primarily responsible for the observed changes in the protein conformation and dynamics.

Optimization of the GB/SA solvation model for predicting the structure of surface loops in proteins

Implicit solvation models are commonly optimized with respect to experimental data or Poisson-Boltzmann (PB) results obtained for small molecules, where the force field is sometimes not considered. In previous studies, we have developed an optimization procedure for cyclic peptides and surface loops in proteins based on the entire system studied and the specific force field used. Thus, the loop has been modeled by the simplified solvation function E(tot) = E(FF) (epsilon = 2r) + Sigma(i) sigma(i)A(i), where E(FF) (epsilon = nr) is the AMBER force field energy with a distance-dependent dielectric function, epsilon = nr, A(i) is the solvent accessible surface area of atom i, and sigma(i) is its atomic solvation parameter. During the optimization process, the loop is free to move while the protein template is held fixed in its X-ray structure. To improve on the results of this model, in the present work we apply our optimization procedure to the physically more rigorous solvation model, the generalized Born with surface area (GB/SA) (together with the all-atom AMBER force field) as suggested by Still and co-workers (J. Phys. Chem. A 1997, 101, 3005). The six parameters of the GB/SA model, namely, P(1)-P(5) and the surface area parameter, sigma (programmed in the TINKER package) are reoptimized for a "training" group of nine loops, and a best-fit set is defined from the individual sets of optimized parameters. The best-fit set and Still's original set of parameters (where Lys, Arg, His, Glu, and Asp are charged or neutralized) were applied to the training group as well as to a "test" group of seven loops, and the energy gaps and the corresponding RMSD values were calculated. These GB/SA results based on the three sets of parameters have been found to be comparable; surprisingly, however, they are somewhat inferior (e.g, of larger energy gaps) to those obtained previously from the simplified model described above. We discuss recent results for loops obtained by other solvation models and potential directions for future studies.

A three-dimensional model of the human immunodeficiency virus type 1 integration complex

While the general features of HIV-1 integrase function are understood, there is still uncertainty about the composition of the integration complex and how integrase interacts with viral and host DNA. We propose an improved model of the integration complex based on current experimental evidence including a comparison with the homologous Tn5 transposase containing bound DNA and an analysis of DNA binding sites using Goodford's GRID. Our model comprises a pair of integrase dimers, two strands of DNA to represent the viral DNA ends and a strand of bent DNA representing the host chromosome. In our model, the terminal four base pairs of each of the viral DNA strands interact with the integrase dimer providing the active site, while bases one turn away interact with a flexible loop (residues 186-194) on the second integrase dimer. We propose that residues E152, Q148 and K156 are involved in the specific recognition of the conserved CA dinucleotide and that the active site mobile loop (residues 140-149) stabilises the integration complex by acting as a barrier to separate the two viral DNA ends. In addition, the residues responsible for DNA binding in our model show a high level of amino acid conservation.

Active site binding modes of the beta-diketoacids: a multi-active site approach in HIV-1 integrase inhibitor design

Predicting a bioactive conformation of a ligand is of paramount importance in rational drug design. The task becomes very difficult when the receptor site possesses a region with unusual conformational flexibility. Significant conformational differences are present in the active site regions in the available crystal structures of the core domains of HIV-1 integrase (IN). Among all reported IN inhibitors, the beta-diketoacid class of compounds has proved to be of most promise and indeed S-1360 was the first IN inhibitor to enter clinical studies. With an aim to predict the bioactive (active site bound) conformation of S-1360, we performed extensive docking studies using three different reported crystal structures where the active site or partial active site region was resolved. For comparison we extended our studies to include 5CITEP (the first compound cocrystallized with IN core domain) and a bis-diketoacid (BDKA). We found that the conformation of S-1360 when bound in one of the active sites matches that of the experimentally observed results of IN escape mutants resistant to S-1360. Therefore, we propose that this active site conformation is the biologically relevant conformation and can be used for the future structure-based drug design studies selectively targeting IN.

Comparative molecular dynamics simulations of HIV-1 integrase and the T66I/M154I mutant: Binding modes and drug resistance to a diketo acid inhibitor

HIV-1 IN is an essential enzyme for viral replication and an interesting target for the design of new pharmaceuticals for use in multidrug therapy of AIDS. L-731,988 is one of the most active molecules of the class of beta-diketo acids. Individual and combined mutations of HIV-1 IN at residues T66, S153, and M154 confer important degrees of resistance to one or more inhibitors belonging to this class. In an effort to understand the molecular mechanism of the resistance of T66I/M154I IN to the inhibitor L-731,988 and its specific binding modes, we have carried out docking studies, explicit solvent MD simulations, and binding free energy calculations. The inhibitor was docked against different protein conformations chosen from prior MD trajectories, resulting in 2 major orientations within the active site. MD simulations have been carried out for the T66I/M154I DM IN, DM IN in complex with L-731,988 in 2 different orientations, and 1QS4 IN in complex with L-731,988. The results of these simulations show a similar dynamical behavior between T66I/M154I IN alone and in complex with L-731,988, while significant differences are observed in the mobility of the IN catalytic loop (residues 138-149). Water molecules bridging the inhibitor to residues from the active site have been identified, and residue Gln62 has been found to play an important role in the interactions between the inhibitor and the protein. This work provides information about the binding modes of L-731,988, as well as insight into the mechanism of inhibitor-resistance in HIV-1 integrase.

The Role of Electrostatic Interaction in Triggering the Unraveling of Stable Helix 1 in Normal Prion Protein. A Molecular Dynamics Simulation Investigation

The conversion of normal prion protein (PrPC) into scrapie isoform (PrPSc) is a key event in the pathogenesis of prion diseases. However, the conversion mechanism has given rise to much controversy. For instance, there is much debate on the behavior of helix 1 (H1) in the conversion. A series of experiments demonstrated that H1 in isolated state was very stable under a variety of conditions. But, other experiments indicated that helices 2 and 3 rather than H1 were retained in PrPSc. In this paper, molecular dynamics (MD) simulation is employed to investigate the dynamic behavior of H1. It is revealed that although the helix 1 of Human PrPC (HuPrPC) is very stable in the isolated state, it becomes unstable when incorporated into native HuPrPC, which likely results from the long-range electrostatic interaction between Asp147 and Arg208 located in the helices 1 and 3, respectively. This explanation is supported by experimental evaluation and MD simulation on D147N mutant of HuPrPC that the mutant becomes a little more stable than the wild type HuPrPC. This finding not only help to reconcile the existing debate on the role of helix 1 in the PrPC-->PrPSc transition, but also reveals a possible mechanism for triggering the PrPC-->PrPSc conversion.

Comparison of multiple molecular dynamics trajectories calculated for the drug-resistant HIV-1 integrase T66I/M154I catalytic domain

HIV-1 integrase (IN) is an essential enzyme for the viral replication and an interesting target for the design of new pharmaceuticals for multidrug therapy of AIDS. Single and multiple mutations of IN at residues T66, S153, or M154 confer degrees of resistance to several inhibitors that prevent the enzyme from performing its normal strand transfer activity. Four different conformations of IN were chosen from a prior molecular dynamics (MD) simulation on the modeled IN T66I/M154I catalytic core domain as starting points for additional MD studies. The aim of this article is to understand the dynamic features that may play roles in the catalytic activity of the double mutant enzyme in the absence of any inhibitor. Moreover, we want to verify the influence of using different starting points on the MD trajectories and associated dynamical properties. By comparison of the trajectories obtained from these MD simulations we have demonstrated that the starting point does not affect the conformational space explored by this protein and that the time of the simulation is long enough to achieve convergence for this system.

Large-scale conformational dynamics of the HIV-1 integrase core domain and its catalytic loop mutants

HIV-1 integrase is one of the three essential enzymes required for viral replication and has great potential as a novel target for anti-HIV drugs. Although tremendous efforts have been devoted to understanding this protein, the conformation of the catalytic core domain around the active site, particularly the catalytic loop overhanging the active site, is still not well characterized by experimental methods due to its high degree of flexibility. Recent studies have suggested that this conformational dynamics is directly correlated with enzymatic activity, but the details of this dynamics is not known. In this study, we conducted a series of extended-time molecular dynamics simulations and locally enhanced sampling simulations of the wild-type and three loop hinge mutants to investigate the conformational dynamics of the core domain. A combined total of >480 ns of simulation data was collected which allowed us to study the conformational changes that were not possible to observe in the previously reported short-time molecular dynamics simulations. Among the main findings are a major conformational change (>20 A) in the catalytic loop, which revealed a gatinglike dynamics, and a transient intraloop structure, which provided a rationale for the mutational effects of several residues on the loop including Q(148), P(145), and Y(143). Further, clustering analyses have identified seven major conformational states of the wild-type catalytic loop. Their implications for catalytic function and ligand interaction are discussed. The findings reported here provide a detailed view of the active site conformational dynamics and should be useful for structure-based inhibitor design for integrase.

Retroviral DNA integration--mechanism and consequences

Integration of retroviral cDNA into the host cell chromosome is an essential step in its replication. This process is catalyzed by the retroviral integrase protein, which is conserved among retroviruses and retrotransposons. Integrase binds viral and host DNA in a complex, called the preintegration complex (PIC), with other viral and cellular proteins. While the PIC is capable of directing integration of the viral DNA into any chromosomal location, different retroviruses have clear preferences for integration in or near particular chromosomal features. The determinants of integration site selection are under investigation but may include retrovirus-specific interactions between integrase and tethering factors bound to the host cell chromosomes. Research into the mechanisms of retroviral integration site selection has shed light on the phenomena of insertional mutagenesis and viral latency.

Molecular dynamics studies of the wild-type and double mutant HIV-1 integrase complexed with the 5CITEP inhibitor: mechanism for inhibition and drug resistance

The human immunodeficiency virus type 1 (HIV-1) integrase (IN) is an essential enzyme in the life cycle of the virus and is an attractive target for the development of new drugs useful in acquired immunodeficiency syndrome multidrug therapy. Starting from the crystal structure of the 5CITEP inhibitor bound to the active site in the catalytic domain of the HIV-1 IN, two different molecular dynamics simulations in water have been carried out. In the first simulation the wild-type IN was used, whereas in the second one the double mutation T66I/M154I, described to lead to drug resistance, was introduced in the protein. Compelling differences have been observed in these two structures during analyses of the molecular dynamics trajectories, particularly in the inhibitor binding modes and in the conformational flexibility of the loop (residues 138-149) located near the three catalytic residues in the active site (Asp(64), Asp(116), Glu(152)). Because the conformational flexibility of this region is important for efficient biological activity and its behavior is quite different in the two models, we suggest a hypothetical mechanism for the inhibition and drug resistance of HIV-1 IN. These results can be useful for the rational design of more potent and selective integrase inhibitors and may allow for the design of inhibitors that will be more robust against known resistance mutations.

Enzymatic conformational fluctuations along the reaction coordinate of cytidine deaminase

Analysis of the crystal structures for cytidine deaminase complexed with substrate analog 3-deazacytidine, transition-state analog zebularine 3,4-hydrate, and product uridine establishes significant changes in the magnitude of atomic-scale fluctuations along the (approximate) reaction coordinate of this enzyme. Differences in fluctuations between the substrate analog complex, transition-state analog complex, and product complex are monitored via changes in corresponding crystallographic temperature factors. Previously, we reported that active-site conformational disorder is substantially reduced in the transition-state complex relative to the two ground-state complexes. Here, this result is statistically corroborated by crystallographic data for fluorinated zebularine 3,4-hydrate, a second transition-state analog, and by multiple regression analysis. Multiple regression explains 70% of the total temperature factor variation through a predictive model for the average B-value of an amino acid as a function of the catalytic state of the enzyme (substrate, transition state, product) and five other physical and structural descriptors. Furthermore, correlations of atomic fluctuation magnitudes throughout the body of each complex are quantified through an auto-correlation function. The transition-state analog complex shows the greatest correlations between temperature factor magnitudes for spatially separated atoms, underscoring the strong ability of this reaction-coordinate species to "organize" enzymatic fluctuations. The catalytic significance for decreased atomic-scale motions in the transition state is discussed. A thermodynamic argument indicates that the significant decreases in local enzymatic conformational entropy at the transition state result in enhanced energetic stabilization there.

Revisiting the structural flexibility of the complex p21(ras)-GTP: the catalytic conformation of the molecular switch II

The hydrolysis of GTP in p21(ras) triggers conformational changes that regulate the ras/ERK signaling pathway. An important active site residue is Gln61, which has been found to be mutated in 30% of human tumors. The dynamics of the active site conformation is studied by using molecular dynamics simulation of two independent structures of the GTP-bound uncomplexed enzyme. Two distinct conformations of the enzyme are observed, in which the side-chain residue Gln61 is in different orientations. Essential dynamics analysis is used to describe the essential motions in the transition between the two conformations. Results are compared with earlier simulations of p21(ras) and its complex with GTPase activating protein p21-GAP.

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.

HIV-1 integrase catalytic core: molecular dynamics and simulated fluorescence decays

Two molecular dynamics simulations have been carried out on the HIV-1 integrase catalytic core starting from fully determined crystal structures. During the first one, performed in the absence of divalent cation (6-ns long), the catalytic core took on two main conformations. The conformational transition occurs at approximately 3.4 ns. In contrast, during the second one, in the presence of Mg(2+) (4-ns long), there were no such changes. The molecular dynamics simulations were used to compute the fluorescence intensity decays emitted by the four tryptophan residues considered as the only chromophores. The decay was computed by following, frame by frame, the amount of chromophores that remained excited at a certain time after light absorption. The simulation took into account the quenching through electron transfer to the peptide bond and the fluorescence resonance energy transfer between the chromophores. The fit to the experimental intensity decays obtained at 5 degrees C and at 30 degrees C is very good. The fluorescence anisotropy decays were also simulated. Interestingly, the fit to the experimental anisotropy decay was excellent at 5 degrees C and rather poor at 30 degrees C. Various hypotheses such as dimerization and abnormal increase of uncorrelated internal motions 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.