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.     

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.     

Analysis of the substrate specificity loop of the HAD superfamily cap domain

The haloacid dehalogenase (HAD) superfamily includes a variety of enzymes that catalyze the cleavage of substrate C-Cl, P-C, and P-OP bonds via nucleophilic substitution pathways. All members possess the alpha/beta core domain, and many also possess a small cap domain. The active site of the core domain is formed by four loops (corresponding to sequence motifs 1-4), which position substrate and cofactor-binding residues as well as the catalytic groups that mediate the "core" chemistry. The cap domain is responsible for the diversification of chemistry within the family. A tight beta-turn in the helix-loop-helix motif of the cap domain contains a stringently conserved Gly (within sequence motif 5), flanked by residues whose side chains contribute to the catalytic site formed at the domain-domain interface. To define the role of the conserved Gly in the structure and function of the cap domain loop of the HAD superfamily members phosphonoacetaldehyde hydrolase and beta-phosphoglucomutase, the Gly was mutated to Pro, Val, or Ala. The catalytic activity was severely reduced in each mutant. To examine the impact of Gly substitution on loop 5 conformation, the X-ray crystal structure of the Gly50Pro phosphonoacetaldehyde hydrolase mutant was determined. The altered backbone conformation at position 50 had a dramatic effect on the spatial disposition of the side chains of neighboring residues. Lys53, the Schiff Base forming lysine, had rotated out of the catalytic site and the side chain of Leu52 had moved to fill its place. On the basis of these studies, it was concluded that the flexibility afforded by the conserved Gly is critical to the function of loop 5 and that it is a marker by which the cap domain substrate specificity loop can be identified within the amino acid sequence of HAD family members.     

Distinct Conformations of a Putative Translocation Element in Poliovirus Polymerase

The mechanism whereby RNA is translocated by the single subunit viral RNA-dependent RNA polymerases is not yet understood. These enzymes lack homologs of the “O-helix” structures and associated fingers domain movements thought to be responsible for translocation in many DNA-templated polymerases. The structures of multiple picornavirus polymerase elongation complexes suggest that these enzymes use a different molecular mechanism where translocation is not strongly coupled to the opening of the active site following catalysis. Here we present the 2.0- to 2.6-Å-resolution crystal structures and biochemical data for 12 poliovirus polymerase mutants that together show how proper enzyme functions and translocation activity requires conformational flexibility of a loop sequence in the palm domain B-motif. Within the loop, the Ser288-Gly289-Cys290 sequence is shown to play a major role in the catalytic cycle based on RNA binding, processive elongation activity, and single nucleotide incorporation assays. The structures show that Ser288 forms a key hydrogen bond with Asp238, the backbone flexibility of Gly289 is required for translocation competency, and Cys290 modulates the overall elongation activity of the enzyme. Some conformations of the loop represent likely intermediates on the way to forming the catalytically competent closed active site, while others are consistent with a role in promoting translocation of the nascent base pair out of the active site. The loop structure and key residues surrounding it are highly conserved, suggesting that the structural dynamics we observe in poliovirus 3D^pol are a common feature of viral RNA-dependent RNA polymerases.     

Membrane structure of the human immunodeficiency virus gp41 fusion peptide by molecular dynamics simulation. II. The glycine mutants

In this work, molecular dynamics (MD) simulation of the interaction of three mutants, G3V, G5V and G10V, of the human immunodeficiency virus (HIV) gp41 16-residue fusion peptide (FP) with an explicit palmitoyloleoylphosphatidyl-ethanolamine (POPE) lipid bilayer was performed. The goals of this work are to study the correlation of the fusogenic activity of the FPs with the mode of their interaction with the bilayer and to examine the roles of the many glycine residues in the FP in the fusion process. The results of this work corroborate the main conclusion of our earlier MD work of the WT FP and several mutants with polar substitution. These two studies provide correlation between the mode of insertion and the fusogenic activity of these peptides and support the hypothesis that an oblique insertion of the fusion domain of the viral protein is required for fusogenic activity. Inactive mutants interact with the bilayer by a surface-binding mode. The results of this work, combined with the results of our earlier work, show that, while the secondary structures of the wild-type FP and its mutants do not affect the fusogenic activities, the conformational flexibility appears to be an important factor. The active WT FP and its partially active mutants, G3V and G5V, all have significant conformational transitions at one of the glycine sites. They occur at Gly(5) in FP-wt, at Gly(10) in FP-G5V and at Gly(13) in FP-G3V. Thus, a glycine site in each of these active (or partially active) FPs provides conformational flexibility. On the other hand, the inactive mutants FP-G10V, FP-L9R and FP-V2E do not have any conformational transitions except at either terminus and thus possess no conformational flexibility. Thus, the results of this work support the suggestion that the role of glycine residues in the fusion domain is to provide the necessary conformational flexibility for fusion activity.The glycines also form a "glycine strip" in the FP that locates on one (the less hydrophobic) face of the helix (the "sided helix"). However, whether this "glycine strip" is disrupted or not does not seem to correlate with the retention of fusogenic activities. Finally, although the FLGFL (8-12) motif is absolutely conserved in the HIV fusion domain, a well-structured motif stabilized by hydrogen bonding does not appear to be required for activity. In fact, hydrogen bonding in this motif was found to be missing in FP-G3V and FP-G5V. Both of these mutants are partially active.     

Explaining an unusually fast parasitic enzyme: folate tail-binding residues dictate substrate positioning and catalysis in Cryptosporidium hominis thymidylate synthase

The essential enzyme TS-DHFR from Cryptosporidium hominis undergoes an unusually rapid rate of catalysis at the conserved TS domain, facilitated by two nonconserved residues, Ala287 and Ser290, in the folate tail-binding region. Mutation of these two residues to their conserved counterparts drastically affects multiple steps of the TS catalytic cycle. We have determined the crystal structures of all three mutants (A287F, S290G, and A287F/S290G) in complex with active site ligands dUMP and CB3717. The structural data show two effects of the mutations: an increased distance between the ligands in the active site and increased flexibility of the folate ligand in the partially open enzyme state that precedes conformational change to the active catalytic state. The latter effect is able to be rescued by the mutants containing the A287F mutation. In addition, the conserved water network of TS is altered in each of the mutants. The structural results point to a role of the folate tail-binding residues in closely positioning ChTS ligands and restricting ligand flexibility in the partially open state to allow for a rapid transition to the active closed state and enhanced rate of catalysis. These results provide an explanation on how folate tail-binding residues at one end of the active site affect long-range interactions throughout the TS active site and validate these residues as targets for species-specific drug design.     

Probing conformational plasticity of the activation domain of trypsin: the role of glycine hinges

Trypsin-like serine proteases play essential roles in diverse physiological processes such as hemostasis, apoptosis, signal transduction, reproduction, immune response, matrix remodeling, development, and differentiation. All of these proteases share an intriguing activation mechanism that involves the transition of an unfolded domain (activation domain) of the zymogen to a folded one in the active enzyme. During this conformational change, activation domain segments move around highly conserved glycine hinges. In the present study, hinge glycines were replaced by alanine residues via site directed mutagenesis. The effects of these mutations on the interconversion of the zymogen-like and active conformations as well as on catalytic activity were studied. Mutant trypsins showed zymogen-like structures to varying extents characterized by increased flexibility of some activation domain segments, a more accessible N-terminus and a deformed substrate binding site. Our results suggest that the trypsinogen to trypsin transition is hindered by the mutations, which results in a shift of the equilibrium between the inactive zymogen-like and active enzyme conformations toward the inactive state. Our data also showed, however, that the inactive conformations of the various mutants differ from each other. Binding of substrate analogues shifted the conformational equilibrium toward the active enzyme since inhibited forms of the trypsin mutants showed similar structural features as the wild-type enzyme. The catalytic activity of the mutants correlated with the proper conformation of the active site, which could be supported by varying conformations of the N-terminus and the autolysis loop. Transient kinetic measurements confirmed the existence of an inactive to active conformational transition occurring prior to substrate binding.     

Molecular dynamics studies of the full-length integrase-DNA complex

We have carried out a molecular dynamics (MD) simulation of full-length HIV-1 integrase (IN) dimer complexed with viral DNA with the aim of gaining information about the enzyme motion and investigating the movement of the catalytic flexible loop (residues 140-149) thought to be essential in the catalytic mechanism of IN. During the simulation, we observed quite a different behavior of this region in the presence or absence of the viral DNA. In particular, the MD results underline the crucial role of the residue Tyr143 in the mechanism of integration of viral DNA into the host chromosome. The present findings confirm the experimental data (e.g., site-directed mutagenesis experiments) showing that the loop is involved in the integration reactions and its mobility is correlated with the catalytic activity of HIV-1 integrase.     

Role of the 207-218 peptide region of Moloney murine leukemia virus integrase in enzyme catalysis

X-ray diffraction data on a few retroviral integrases show a flexible loop near the active site. By sequence alignment, the peptide region 207-218 of Mo-MLV IN appears to correspond to this flexible loop. In this study, residues H208, Y211, R212, Q214, S215 and S216 of Mo-MLV IN were mutated to determine their role on enzyme activity. We found that Y211A, R212A, R212K and Q214A decreased integration activity, while disintegration and 3′-processing were not significantly affected. By contrast H208A was completely inactive in all the assays. The core domain of Mo-MLV integrase was modeled and the flexibility of the region 207-216 was analyzed. Substitutions with low integration activity showed a lower flexibility than wild type integrase. We propose that the peptide region 207-216 is a flexible loop and that H208, Y211, R212 and Q214 of this loop are involved in the correct assembly of the DNA-integrase complex during integration.     

The effects of amino acid replacements of glycine 20 on conformational stability and catalysis of staphylococcal nuclease

Staphylococcal nuclease (SNase) is a well-established model for protein folding studies. Its three-dimensional structure has been determined. The enzyme, Ca2+, and DNA or RNA substrate form a ternary complex. Glycine 20 is the second position of the first beta-turn of SNase, which may serve as the folding initiation site for the SNase polypeptide. To study the role of Gly20 in the conformational stability and catalysis of SNase, three mutants, in which Gly20 was replaced by alanine, valine, or isoleucine, were constructed and studied by using circular dichroism spectra, intrinsic and ANS-binding fluorescence spectra, stability and activity assays. The mutations have little effect on the conformational integrity of the mutants. However, the catalytic activity is reduced drastically by the mutations, and the stability of the protein is progressively decreased in the order G20A相似文献   

X-ray structure of simian immunodeficiency virus integrase containing the core and C-terminal domain (residues 50-293)--an initial glance of the viral DNA binding platform

The crystal structure of simian immunodeficiency virus (SIV) integrase that contains in a single polypeptide the core and the C-terminal deoxyoligonucleotide binding domain has been determined at 3 A resolution with an R-value of 0.203 in the space group P2(1)2(1)2(1). Four integrase core domains and one C-terminal domain are found to be well defined in the asymmetric unit. The segment extending from residues 114 to 121 assumes the same position as seen in the integrase core domain of avian sarcoma virus as well as human immunodeficiency virus type-1 (HIV-1) crystallized in the absence of sodium cacodylate. The flexible loop in the active site, composed of residues 141-151, remains incompletely defined, but the location of the essential Glu152 residue is unambiguous. The residues from 210-218 that link the core and C-terminal domains can be traced as an extension from the core with a short gap at residues 214-215. The C(alpha) folding of the C-terminal domain is similar to the solution structure of this domain from HIV-1 integrase. However, the dimeric form seen in the NMR structure cannot exist as related by the non-crystallographic symmetry in the SIV integrase crystal. The two flexible loops of the C-terminal domain, residues 228-236 and residues 244-249, are much better fixed in the crystal structure than in the NMR structure with the former in the immediate vicinity of the flexible loop of the core domain. The interface between the two domains encompasses a solvent-exclusion area of 1500 A(2). Residues from both domains purportedly involved in DNA binding are narrowly distributed on the same face of the molecule. They include Asp64, Asp116, Glu152 and Lys159 from the core and Arg231, Leu234, Arg262, Arg263 and Lys264 from the C-terminal domain. A model for DNA binding is proposed to bridge the two domains by tethering the 228-236 loop of the C-terminal domain and the flexible loop of the core.     

Role of amino acid residues in left-handed helical conformation for the conformational stability of a protein

Our previous study of six non-Gly to Gly/Ala mutant human lysozymes in a left-handed helical region showed that only one non-Gly residue at a rigid site had unfavorable strain energy as compared with Gly at the same position (Takano et al., Proteins 2001; 44:233-243). To further examine the role of left-handed residues in the conformational stability of a protein, we constructed ten Gly to Ala mutant human lysozymes. Most Gly residues in human lysozyme are located in the left-handed helix region. The thermodynamic parameters for denaturation and crystal structures were determined by differential scanning calorimetry and X-ray analysis, respectively. The difference in denaturation Gibbs energy (DeltaDeltaG) for the ten Gly to Ala mutants ranged from + 1.9 to -7.5 kJ/mol, indicating that the effect of the mutation depends on the environment of the residue. We confirm that Gly in a left-handed region is more favorable at rigid sites than non-Gly, but there is little difference in energetic cost between Gly and non-Gly at flexible sites. The present results indicate that dihedral angles in the backbone conformation and also the flexibility at the position should be considered for analyses of protein stability, and protein structural determination, prediction, and design.     

Solution structure of the catalytic domain of gammadelta resolvase. Implications for the mechanism of catalysis

The site-specific DNA recombinase, gammadelta resolvase, from Escherichia coli catalyzes recombination of res site-containing plasmid DNA to two catenated circular DNA products. The catalytic domain (residues 1-105), lacking a C-terminal dimerization interface, has been constructed and the NMR solution structure of the monomer determined. The RMSD of the NMR conformers for residues 2-92 excluding residues 37-45 and 64-73 is 0.41 A for backbone atoms and 0.88 A for all heavy atoms. The NMR solution structure of the monomeric catalytic domain (residues 1-105) was found to be formed by a four-stranded parallel beta-sheet surrounded by three helices. The catalytic domain (residues 1-105), deficient in the C-terminal dimerization domain, was monomeric at high salt concentration, but displayed unexpected dimerization at lower ionic strength. The unique solution dimerization interface at low ionic strength was mapped by NMR. With respect to previous crystal structures of the dimeric catalytic domain (residues 1-140), differences in the average conformation of active-site residues were found at loop 1 containing the catalytic S10 nucleophile, the beta1 strand containing R8, and at loop 3 containing D67, R68 and R71, which are required for catalysis. The active-site loops display high-frequency and conformational backbone dynamics and are less well defined than the secondary structures. In the solution structure, the D67 side-chain is proximal to the S10 side-chain making the D67 carboxylate group a candidate for activation of S10 through general base catalysis. Four conserved Arg residues can function in the activation of the phosphodiester for nucleophilic attack by the S10 hydroxyl group. A mechanism for covalent catalysis by this class of recombinases is proposed that may be related to dimer interface dissociation.     

Molecular modeling of the dimeric structure of human lipoprotein lipase and functional studies of the carboxyl-terminal domain.

Lipoprotein lipase (LPL) plays a key role in lipid metabolism. Molecular modeling of dimeric LPL was carried out using insight ii based upon the crystal structures of human, porcine, and horse pancreatic lipase. The dimeric model reveals a saddle-shaped structure and the key heparin-binding residues in the amino-terminal domain located on the top of this saddle. The models of two dimeric conformations - a closed, inactive form and an open, active form - differ with respect to how surface-loop positions affect substrate access to the catalytic site. In the closed form, the surface loop covers the catalytic site, which becomes inaccessible to solvent. Large conformational changes in the open form, especially in the loop and carboxyl-terminal domain, allow substrate access to the active site. To dissect the structure-function relationships of the LPL carboxyl-terminal domain, several residues predicted by the model structure to be essential for the functions of heparin binding and substrate recognition were mutagenized. Arg405 plays an important role in heparin binding in the active dimer. Lys413/Lys414 or Lys414 regulates heparin affinity in both monomeric and dimeric forms. To evaluate the prediction that LPL forms a homodimer in a 'head-to-tail' orientation, two inactive LPL mutants - a catalytic site mutant (S132T) and a substrate-recognition mutant (W390A/W393A/W394A) - were cotransfected into COS7 cells. Lipase activity could be recovered only when heterodimerization occurred in a head-to-tail orientation. After cotransfection, 50% of the wild-type lipase activity was recovered, indicating that lipase activity is determined by the interaction between the catalytic site on one subunit and the substrate-recognition site on the other.     

Crystal structures of the catalytic domain of human protein kinase associated with apoptosis and tumor suppression

We have determined X-ray crystal structures with up to 1.5 A resolution of the catalytic domain of death-associated protein kinase (DAPK), the first described member of a novel family of pro-apoptotic and tumor-suppressive serine/threonine kinases. The geometry of the active site was studied in the apo form, in a complex with nonhydrolyzable AMPPnP and in a ternary complex consisting of kinase, AMPPnP and either Mg2+ or Mn2+. The structures revealed a previously undescribed water-mediated stabilization of the interaction between the lysine that is conserved in protein kinases and the beta- and gamma-phosphates of ATP, as well as conformational changes at the active site upon ion binding. Comparison between these structures and nucleotide triphosphate complexes of several other kinases disclosed a number of unique features of the DAPK catalytic domain, among which is a highly ordered basic loop in the N-terminal domain that may participate in enzyme regulation.     

Vaccinia topoisomerase mutants illuminate conformational changes during closure of the protein clamp and assembly of a functional active site.

We present a mutational analysis of vaccinia topoisomerase that highlights the contributions of five residues in the catalytic domain (Phe-88 and Phe-101 in helix alpha1, Ser-204 in alpha5, and Lys-220 and Asn-228 in alpha6) to the DNA binding and transesterification steps. When augmented by structural information from exemplary type IB topoisomerases and tyrosine recombinases in different functional states, the results suggest how closure of the protein clamp around duplex DNA and assembly of a functional active site might be orchestrated by internal conformational changes in the catalytic domain. Lys-220 is a constituent of the active site, and a positive charge at this position is required for optimal DNA cleavage. Ser-204 and Asn-228 appear not to be directly involved in reaction chemistry at the scissile phosphodiester. We propose that (i) Asn-228 recruits the Tyr-274 nucleophile to the active site by forming a hydrogen bond to the main chain of the tyrosine-containing alpha8 helix and that (ii) contacts between Ser-204 and the DNA backbone upstream of the cleavage site trigger a separate conformational change required for active site assembly. Mutations of Phe-88 and Phe-101 affect DNA binding, most likely at the clamp closure step, which we posit to entail a distortion of helix alpha1.     

Dynamics and DNA substrate recognition by the catalytic domain of lambda integrase

Bacteriophage lambda integrase (lambda-Int) is the prototypical member of a large family of enzymes that catalyze site-specific DNA recombination via the formation of a Holliday junction intermediate. DNA strand cleavage by lambda-Int is mediated by nucleophilic attack on the scissile phosphate by a conserved tyrosine residue, forming an intermediate with the enzyme covalently attached to the 3'-end of the cleaved strand via a phosphotyrosine linkage. The crystal structure of the catalytic domain of lambda-Int (C170) obtained in the absence of DNA revealed the tyrosine nucleophile at the protein's C terminus to be located on a beta-hairpin far from the other conserved catalytic residues and adjacent to a disordered loop. This observation suggested that a conformational change in the C terminus of the protein was required to generate the active site in cis, or alternatively, that the active site could be completed in trans by donation of the tyrosine nucleophile from a neighboring molecule in the recombining synapse. We used NMR spectroscopy together with limited proteolysis to examine the dynamics of the lambda-Int catalytic domain in the presence and absence of DNA half-site substrates with the goal of characterizing the expected conformational change. Although the C terminus is indeed flexible in the absence of DNA, we find that conformational changes in the tyrosine-containing beta-hairpin are not coupled to DNA binding. To gain structural insights into C170/DNA complexes, we took advantage of mechanistic conservation with Cre and Flp recombinases to model C170 in half-site and tetrameric Holliday junction complexes. Although the models do not reveal the nature of the conformational change required for cis cleavage, they are consistent with much of the available experimental data and provide new insights into the how trans complementation could be accommodated.