Molecular dynamics simulations of HIV-1 protease monomer: Assembly of N-terminus and C-terminus into beta-sheet in water solution

HIV-1 protease (HIV-PR) consists of two identical subunits that are united together through a four-stranded antiparallel beta-sheet formed of the peptide termini of each monomer. Since the active site exists only in the dimer, a strategy that is attracting more and more attention in inhibitor design and which may overcome the serious drug resistance caused by competitive inhibitors is to block the peptide termini of the monomer, thereby interfering with formation of the active dimer. In the present work, we performed several extensive molecular dynamics (MD) simulations of the HIV-PR monomer in water to illustrate its solvated conformation and dynamics behavior. We found that the peptide termini usually assembled into beta-sheet after several nanoseconds' simulation, and became much less flexible. This beta-sheet is stabilized by intramolecular interactions and is not easily disaggregated under the present MD simulation conditions. This transformation may be an important transition during the relaxing and equilibrating of the HIV-PR monomer in aqueous solution, and the terminal beta-sheet may be one of the major conformations of the solvated HIV-PR monomer termini in water. This work may provide new insights into the dynamics behavior and dimerization mechanism of HIV-PR, and more significantly, offer a more rational receptor model for the design and discovery of novel dimerization inhibitors than crystalline structures.     

Molecular dynamics studies of caspase-3

Caspase-3 is a fundamental target for pharmaceutical interventions against a variety of diseases involving disregulated apoptosis. The enzyme is active as a dimer with two symmetry-related active sites, each featuring a Cys-His catalytic dyad and a selectivity loop, which recognizes the characteristic DEVD pattern of the substrate. Here, a molecular dynamics study of the enzyme in complex with two pentapeptide substrates DEVDG is presented, which provides a characterization of the dynamic properties of the active form in aqueous solution. The mobility of the substrate and that of the catalytic residues are rather low indicating a distinct preorganization effect of the Michaelis complex. An essential mode analysis permits us to identify coupled motions between the two monomers. In particular, it is found that the motions of the two active site loops are correlated and tend to steer the substrate toward the reactive center, suggesting that dimerization has a distinct effect on the dynamic properties of the active site regions. The selectivity loop of one monomer turns out to be correlated with the N-terminal region of the p12 subunit of the other monomer, an interaction that is also found to play a fundamental role in the electrostatic stabilization of the quaternary structure. To further characterize the specific influence of dimerization on the enzyme essential motions, a molecular dynamics analysis is also performed on the isolated monomer.     

Solution structure of the mature HIV-1 protease monomer: insight into the tertiary fold and stability of a precursor

We present the first solution structure of the HIV-1 protease monomer spanning the region Phe1-Ala95 (PR1-95). Except for the terminal regions (residues 1-10 and 91-95) that are disordered, the tertiary fold of the remainder of the protease is essentially identical to that of the individual subunit of the dimer. In the monomer, the side chains of buried residues stabilizing the active site interface in the dimer, such as Asp25, Asp29, and Arg87, are now exposed to solvent. The flap dynamics in the monomer are similar to that of the free protease dimer. We also show that the protease domain of an optimized precursor flanked by 56 amino acids of the N-terminal transframe region is predominantly monomeric, exhibiting a tertiary fold that is quite similar to that of PR1-95 structure. This explains the very low catalytic activity observed for the protease prior to its maturation at its N terminus as compared with the mature protease, which is an active stable dimer under identical conditions. Adding as few as 2 amino acids to the N terminus of the mature protease significantly increases its dissociation into monomers. Knowledge of the protease monomer structure and critical features of its dimerization may aid in the screening and design of compounds that target the protease prior to its maturation from the Gag-Pol precursor.     

Substrate and Inhibitor-induced Dimerization and Cooperativity in Caspase-1 but Not Caspase-3

Caspases are intracellular cysteine-class proteases with aspartate specificity that is critical for driving processes as diverse as the innate immune response and apoptosis, exemplified by caspase-1 and caspase-3, respectively. Interestingly, caspase-1 cleaves far fewer cellular substrates than caspase-3 and also shows strong positive cooperativity between the two active sites of the homodimer, unlike caspase-3. Biophysical and kinetic studies here present a molecular basis for this difference. Analytical ultracentrifugation experiments show that mature caspase-1 exists predominantly as a monomer under physiological concentrations that undergoes dimerization in the presence of substrate; specifically, substrate binding shifts the K_D for dimerization by 20-fold. We have created a hemi-active site-labeled dimer of caspase-1, where one site is blocked with the covalent active site inhibitor, benzyloxycarbonyl-Val-Ala-Asp-fluoromethylketone. This hemi-labeled enzyme is about 9-fold more active than the apo-dimer of caspase-1. These studies suggest that substrate not only drives dimerization but also, once bound to one site in the dimer, promotes an active conformation in the other monomer. Steady-state kinetic analysis and modeling independently support this model, where binding of one substrate molecule not only increases substrate binding in preformed dimers but also drives the formation of heterodimers. Thus, the cooperativity in caspase-1 is driven both by substrate-induced dimerization as well as substrate-induced activation. Substrate-induced dimerization and activation seen in caspase-1 and not in caspase-3 may reflect their biological roles. Whereas caspase-1 cleaves a dramatically smaller number of cellular substrates that need to be concentrated near inflammasomes, caspase-3 is a constitutively active dimer that cleaves many more substrates located diffusely throughout the cell.     

Differences in dynamic structure of LC8 monomer, dimer, and dimer-peptide complexes

Dimerization of dynein light chain LC8 creates two symmetric grooves at the dimer interface with diverse binding capabilities. In addition to pH and protein concentration, dimerization is affected by phosphorylation, as illustrated by a phosphomimetic mutation that promotes dissociation of LC8 to a monomer and subsequent dissociation from the dynein complex in vitro. In this work we characterize the dynamic structure and unfolding profiles of an LC8 mutant, H55K, as a model for monomeric LC8 at neutral pH. Backbone (15)N relaxation experiments show that the monomer, while primarily ordered, has more heterogeneous dynamics relative to the LC8 dimer, predominantly in residues that ultimately form the binding groove, particularly those in beta 1 and beta 3 strands. This heterogeneity suggests that conformations that are primed for binding are sampled in the inactive monomer and favored in the active dimer. Further changes of LC8 backbone dynamics upon binding to short peptides from Swallow (Swa) and dynein intermediate chain (IC) were elucidated. The conformational heterogeneity apparent in the LC8 dimer is retained in LC8/IC but is lost in LC8/Swa, suggesting that the degree of ordering upon binding is ligand dependent. The reduced complexity of motion in LC8/Swa correlates with the less favorable entropy of binding of LC8 to Swa relative to IC. We propose that the conformational motility of beta 3 has functional significance in dimerization and in ligand binding. In the latter, beta 3 flexibility apparently accommodates different binding modes for different ligands resulting in ligand-specific conformational dynamics of the binding site that may impact other processes such as accessibility to phosphorylation.     

Structural requirements for catalysis and dimerization of human methionine adenosyltransferase I/III

We have used site-directed mutagenesis to probe the structural requirements for catalysis and dimerization of human hepatic methionine adenosyltransferase (hMAT). We built a homology model of the dimeric hMAT III inferred by the crystal structure of the highly homologous Escherichia coli MAT dimer. The active sites of both enzymes comprise the same amino acids and are located in the inter-subunit interface. All of the amino acids predicted to be in the hMAT III active site were mutated, as well as residues in a conserved ATP binding region. All of the mutations except one severely affected catalytic activity. On the other hand, dimerization was affected only by single mutations of three different residues, all on one monomer. The homology model suggested that the side chains of these residues stabilized the monomer and participated in a bridge between subunits consisting of a network of metal and phosphate ions. In agreement with this observation, we demonstrated that dimerization cannot occur in the absence of phosphate.     

Binding and fluorescence studies of the relationship between neurophysin-peptide interaction and neurophysin self-association: an allosteric system exhibiting minimal cooperativity

The mechanism of peptide-enhanced neurophysin self-association was investigated to address questions raised by the crystal structure of a neurophysin-dipeptide complex. The dependence on protein concentration of the binding of a broad range of peptides to the principal hormone-binding site confirmed that occupancy of this site alone, and not a site that bridges the monomer-monomer interface, is the trigger for enhanced dimerization. For the binding of most peptides to the principal hormone-binding site on bovine neurophysin I, the affinity of each dimer site was at least 10 times that of monomer under the conditions used. No interactions between the two sites of the dimer were evident. Fluorescence polarization studies of pressure-induced dimer dissociation indicated that the volume change for this reaction was almost 4 times greater in the liganded than in the unliganded state, pointing to a significant alteration of the monomer-monomer interface upon peptide binding. Novel conformational changes in the vicinity of the single neurophysin tyrosine, Tyr-49, induced by pressures lower than required for subunit dissociation, were also observed. The bovine neurophysin I dimer therefore appears to represent an allosteric system in which there is thermodynamic and functional communication between each binding site and the monomer-monomer interface, but no communication across the interface to the binding site of the other subunit. A model for the peptide-enhanced dimerization is proposed in which intersubunit contacts between monomers reduce the large unfavorable free energy associated with binding-induced intrasubunit conformational change. Structural origins of the lack of communication across the interface are suggested on the basis of the low volume change associated with dimerization in the unliganded state and monomer-monomer contacts in the crystal structure. Potential roles for the peptide alpha-amino group and position 2 phenyl ring in triggering conformational change are discussed.     

Interaction between D-amino acid oxidase and small molecules.

1. From the standpoint of monomer-dimer equilibrium of hog kidney D-amino acid oxidase [EC 1.4.3.3] and the interaction between the enzyme and small molecules, the effect of pH on the binding of p-aminobenzoate to the monomer and dimer of the enzyme was studied by kinetic methods and spectrophotometric titration. 2. The maximum binding number of p-aminobenzoate to the dimer is two molecules, and there is no interaction between the two active sites of the dimer (i.e., no cooperativity) over the range of pH from 6.5 to 10. 3. The affinity of the dimer for p-aminobenzoate is several times higher than that of the monomer at pH 6.5-10, and consequently p-aminobenzoate induces dimerization in the equilibrium state of D-amino acid oxidase. The interaction energy of two subunits of the dimer is stabilized by the binding of p-aminobenzoate by 1-2 kcal/mole over the pH range studied. 4. The binding sites of the quasi-substrate, p-aminobenzoate, in the dimer and the intersubunit binding site of the dimer are clearly different, because p-aminobenzoate induces dimerization of the enzyme. 5. The pK values of ionizing groups in the free monomer and the free dimer which participate in the binding of the competitive inhibitor, p-aminobenzoate, are approximately the same, 8.7, as determined from the pH dependence of the affinity of the inhibitor for the enzyme. Furthermore, no pK for the enzyme-inhibitor complex in the pH range 6.5-10 was observed. 6. There is no interaction between the two ionizing groups of the dimer during protonation-deprotonation, because a theoretical equation involving no cooperativity between the two ionizing groups in the dimer explains the results well.     

Macromolecular chelation as an improved mechanism of protease inhibition: structure of the ecotin-trypsin complex.

The 2.4 A crystal structure (R = 0.180) of the serine protease inhibitor ecotin was determined in a complex with trypsin. Ecotin's dimer structure provides a second discrete and distal binding site for trypsin and, as shown by modelling experiments, other serine proteases. The second site is approximately 45 A from the reactive/active site of the complex and features 13 hydrogen bonds, including six that involve carbonyl oxygen atoms and four bridged by water molecules. Contacts ecotin makes with trypsin's active site are similar to, though more extensive than, those found between trypsin and basic pancreatic trypsin inhibitor. The side chain of ecotin Met84 is found in the substrate binding pocket of trypsin where it makes few contacts, but also does not disrupt the solvent structure or cause misalignment of the scissile bond. This first case of protein dimerization being used to augment binding energy and allow chelation of a target protein provides a new model for protein-protein interactions and for protease inhibition.     

Theoretical studies of relaxation of a monomeric subunit of HIV-1 protease in water using molecular dynamics

The dynamic behavior of one 99-residue subunit of the dimeric aspartyl protease of HIV-1 was studied in a 160 psec molecular dynamics simulation at 300 K in water. The crystal structure of one of the identical subunits of the dimer was the starting point, with the aqueous phase modeled by 4,331 explicit waters in a restrained spherical droplet Analysis of the simulations showed that the monomer displayed considerable flexibility in the interfacial portions of the flap (the region which folds over the substrate), the N- and C-0termini, and, to a lesser extent, the active site. The flap undergoes significant motion as an independent rigid finger, but without the cantilever previously reported hi a simulation of the dimer. The N-terminus displayed the greatest fluctuational disorder whereas the C-terminus exhibited the greatest root mean square movement from the crystal structure. The central core of the monomer had a heavy-atom root mean square deviation from the initial structure of about 3.0 Å during the latter half of the simulation. Although this is larger than the 1.6 Å found for comparable simulations of typical globular proteins, the general features of the tertiary structure were preserved over the course of the simulation. Overall, these results indicate that the relaxed structure obtained in these simulations may provide a better model for the tertiary structure of the solvated HIV-1 protease monomer than the subunit conformation seen in the X-ray crystallographic structure of the dimer. Except in the flap region, the design of compounds intended to interfere with dimerization should take this relaxation and the flexibility of the solvated monomer, especially at the termini, into account. © 1993 Wiley-Liss, Inc.     

Modulation of dimer stability in yeast pyrophosphatase by mutations at the subunit interface and ligand binding to the active site

Yeast (Saccharomyces cerevisiae) pyrophosphatase (Y-PPase) is a tight homodimer with two active sites separated in space from the subunit interface. The present study addresses the effects of mutation of four amino acid residues at the subunit interface on dimer stability and catalytic activity. The W52S variant of Y-PPase is monomeric up to an enzyme concentration of 300 microm, whereas R51S, H87T, and W279S variants produce monomer only in dilute solutions at pH > or = 8.5, as revealed by sedimentation, gel electrophoresis, and activity measurements. Monomeric Y-PPase is considerably more sensitive to the SH reagents N-ethylmaleimide and p-hydroxymercurobenzosulfonate than the dimeric protein. Additionally, replacement of a single cysteine residue (Cys(83)), which is not part of the subunit interface or active site, with Ser resulted in insensitivity of the monomer to SH reagents and stabilization against spontaneous inactivation during storage. Active site ligands (Mg(2+) cofactor, P(i) product, and the PP(i) analog imidodiphosphate) stabilized the W279S dimer versus monomer predominantly by decreasing the rate of dimer to monomer conversion. The monomeric protein exhibited a markedly increased (5-9-fold) Michaelis constant, whereas k(cat) remained virtually unchanged, compared with dimer. These results indicate that dimerization of Y-PPase improves its substrate binding performance and, conversely, that active site adjustment through cofactor, product, or substrate binding strengthens intersubunit interactions. Both effects appear to be mediated by a conformational change involving the C-terminal segment that generally shields the Cys(83) residue in the dimer.     

Functional consequences of the Kaposi's sarcoma-associated herpesvirus protease structure: regulation of activity and dimerization by conserved structural elements

The structure of Kaposi's sarcoma-associated herpesvirus protease (KSHV Pr), at 2.2 A resolution, reveals the active-site geometry and defines multiple possible target sites for drug design against a human cancer-producing virus. The catalytic triad of KSHV Pr, (Ser114, His46, and His157) and transition-state stabilization site are arranged as in other structurally characterized herpesviral proteases. The distal histidine-histidine hydrogen bond is solvent accessible, unlike the situation in other classes of serine proteases. As in all herpesviral proteases, the enzyme is active only as a weakly associated dimer (K(d) approximately 2 microM), and inactive as a monomer. Therefore, both the active site and dimer interface are potential targets for antiviral drug design. The dimer interface in KSHV Pr is compared with the interface of other herpesviral proteases. Two conserved arginines (Arg209), one from each monomer, are buried within the same region of the dimer interface. We propose that this conserved arginine may provide a destabilizing element contributing to the tuned micromolar dissociation of herpesviral protease dimers.     

DNA target sequence identification mechanism for dimer-active protein complexes

Sequence-specific DNA-binding proteins must quickly and reliably localize specific target sites on DNA. This search process has been well characterized for monomeric proteins, but it remains poorly understood for systems that require assembly into dimers or oligomers at the target site. We present a single-molecule study of the target-search mechanism of protelomerase TelK, a recombinase-like protein that is only active as a dimer. We show that TelK undergoes 1D diffusion on non-target DNA as a monomer, and it immobilizes upon dimerization even in the absence of a DNA target site. We further show that dimeric TelK condenses non-target DNA, forming a tightly bound nucleoprotein complex. Together with theoretical calculations and molecular dynamics simulations, we present a novel target-search model for TelK, which may be generalizable to other dimer and oligomer-active proteins.     

Identification of the Serratia endonuclease dimer: structural basis and implications for catalysis.

The Serratia endonuclease is an extracellularly secreted enzyme capable of cleaving both single- and double-stranded forms of DNA and RNA. It is the first member of a large class of related and usually dimeric endonucleases for which a structure is known. Using X-ray crystallography, the structure of monomer of this enzyme was reported by us previously (Miller MD et al., 1994, Nature Struct Biol 1:461-468). We now confirm the dimeric nature of this enzyme through light-scattering experiments and identify the physiologic dimer interface through crystal packing analysis. This dimerization occurs through an isologous twofold interaction localized to the carboxy-terminal subdomain of the enzyme. The dimer is a prolate ellipsoid with dimensions 30 A x 35 A x 90 A. The dimer interface is flat and contains four salt links, several hydrogen bonds, and nonpolar interactions. Buried water is prominent in this interface and it includes an unusual "cubic" water cluster. The position of the two active sites in the dimer suggests that they can act independently in their cleavage of DNA, but have a geometrical advantage in attacking substrate relative to the monomer.     

Crystal structure of Bacillus subtilis guanine deaminase: the first domain-swapped structure in the cytidine deaminase superfamily

Guanine deaminase, a key enzyme in the nucleotide metabolism, catalyzes the hydrolytic deamination of guanine into xanthine. The crystal structure of the 156-residue guanine deaminase from Bacillus subtilis has been solved at 1.17-A resolution. Unexpectedly, the C-terminal segment is swapped to form an intersubunit active site and an intertwined dimer with an extensive interface of 3900 A(2) per monomer. The essential zinc ion is ligated by a water molecule together with His(53), Cys(83), and Cys(86). A transition state analog was modeled into the active site cavity based on the tightly bound imidazole and water molecules, allowing identification of the conserved deamination mechanism and specific substrate recognition by Asp(114) and Tyr(156'). The closed conformation also reveals that substrate binding seals the active site entrance, which is controlled by the C-terminal tail. Therefore, the domain swapping has not only facilitated the dimerization but has also ensured specific substrate recognition. Finally, a detailed structural comparison of the cytidine deaminase superfamily illustrates the functional versatility of the divergent active sites found in the guanine, cytosine, and cytidine deaminases and suggests putative specific substrate-interacting residues for other members such as dCMP deaminases.     

Verification of the interdomain contact site in the inactive monomer, and the domain-swapped fold in the active dimer of Hsp33 in solution

Upon dimerization by oxidation, Hsp33 functions as a molecular chaperone in prokaryotes. Previously published structures of both the inactive and active species are of doubtful relevance to the solution conformations since the inactive (reduced) crystal structure was dimeric, while the active (oxidized) species was crystallized with a truncation of its regulation domain. The interdomain contact site of the inactive monomer, identified in this work, is consistent with that previously observed in the reduced dimer crystal. In contrast, fluorescence quenching of the active dimer contradicted the results expected from the domain-swapped fold observed in the truncated dimer crystal. The results of this study provide important new information concerning controversial issues in the activation process of Hsp33.     

Molecular mechanism for dimerization to regulate the catalytic activity of human cytomegalovirus protease

Biochemical studies indicate that dimerization is required for the catalytic activity of herpesvirus proteases, whereas structural studies show a complete active site in each monomer, away from the dimer interface. Here we report kinetic, biophysical and crystallographic characterizations of structure-based mutants in the dimer interface of human cytomegalovirus (HCMV) protease. Such mutations can produce a 1,700-fold reduction in the kcat while having minimal effects on the K(m). Dimer stability is not affected by these mutations, suggesting that dimerization itself is insufficient for activity. There are large changes in monomer conformation and dimer organization of the apo S225Y mutant enzyme. However, binding of an activated peptidomimetic inhibitor induced a conformation remarkably similar to the wild type protease. Our studies suggest that appropriate dimer formation may be required to indirectly stabilize the protease oxyanion hole, revealing a novel mechanism for dimerization to regulate enzyme activity.     

Dimerization of Bacterial Diaminopimelate Epimerase Is Essential for Catalysis

Diaminopimelate (DAP) epimerase is involved in the biosynthesis of meso-DAP and lysine, which are important precursors for the synthesis of peptidoglycan, housekeeping proteins, and virulence factors in bacteria. Accordingly, DAP epimerase is a promising antimicrobial target. Previous studies report that DAP epimerase exists as a monomeric enzyme. However, we show using analytical ultracentrifugation, X-ray crystallography, and enzyme kinetic analyses that DAP epimerase from Escherichia coli exists as a functional dimer in solution and the crystal state. Furthermore, the 2.0-Å X-ray crystal structure of the E. coli DAP epimerase dimer shows for the first time that the enzyme exists in an open, active conformation. The importance of dimerization was subsequently probed by using site-directed mutagenesis to generate a monomeric mutant (Y268A). Our studies show that Y268A is catalytically inactive, thus demonstrating that dimerization of DAP epimerase is essential for catalysis. Molecular dynamics simulations indicate that the DAP epimerase monomer is inherently more flexible than the dimer, suggesting that dimerization optimizes protein dynamics to support function. Our findings offer insight into the development of novel antimicrobial agents targeting the dimeric antibiotic target DAP epimerase.     

Herpesvirus protease inhibition by dimer disruption

Kaposi's sarcoma-associated herpesvirus (KSHV), like all herpesviruses, encodes a protease (KSHV Pr), which is necessary for the viral lytic cycle. Herpesvirus proteases function as obligate dimers; however, each monomer has an intact, complete active site which does not interact directly with the other monomer across the dimer interface. Protein grafting of an interfacial KSHV Pr alpha-helix onto a small stable protein, avian pancreatic polypeptide, generated a helical 30-amino-acid peptide designed to disrupt the dimerization of KSHV Pr. The chimeric peptide was optimized through protein modeling of the KSHV Pr-peptide complex. Circular dichroism analysis and gel filtration chromatography revealed that the rationally designed peptide adopts a helical conformation and is capable of disrupting KSHV Pr dimerization, respectively. Additionally, the optimized peptide inhibits KSHV Pr activity by 50% at a approximately 200-fold molar excess of peptide to KSHV Pr, and the dissociation constant was estimated to be 300 microM. Mutagenesis of the interfacial residue M197 to a leucine resulted in an inhibitory concentration which was twofold higher for KSHV Pr M197L than for KSHV Pr, in agreement with the model that the dimer interface is involved in peptide binding. These results indicate that the dimer interface, as well as the active sites, of herpesvirus proteases is a viable target for inhibiting enzyme activity.     

Modulation of glycophorin A transmembrane helix interactions by lipid bilayers: molecular dynamics calculations

Starting from the glycophorin A dimer structure determined by NMR, we performed simulations of both dimer and monomer forms in explicit lipid bilayers with constant normal pressure, lateral area, and temperature using the CHARMM potential. Analysis of the trajectories in four different lipids reveals how lipid chain length and saturation modulate the structural and energetic properties of transmembrane helices. Helix tilt, helix-helix crossing angle, and helix accessible volume depend on lipid type in a manner consistent with hydrophobic matching concepts: the most relevant lipid property appears to be the bilayer thickness. Although the net helix-helix interaction enthalpy is strongly attractive, analysis of residue-residue interactions reveals significant unfavorable electrostatic repulsion between interfacial glycine residues previously shown to be critical for dimerization. Peptide volume is nearly conserved upon dimerization in all lipid types, indicating that the monomeric helices pack equally well with lipid as dimer helices do with one another. Enthalpy calculations indicate that the helix-environment interaction energy is lower in the dimer than in the monomer form, when solvated by unsaturated lipids. In all lipid environments there is a marked preference for lipids to interact with peptide predominantly through one rather than both acyl chains. Although our trajectories are not long enough to allow a full thermodynamic treatment, these results demonstrate that molecular dynamics simulations are a powerful method for investigating the protein-protein, protein-lipid, and lipid-lipid interactions that determine the structure, stability and dynamics of transmembrane alpha-helices in membranes.