Specificity and catalysis of uracil DNA glycosylase. A molecular dynamics study of reactant and product complexes with DNA

The structure of uracil DNA glycosylase (UDG) in complex with a nonamer duplex DNA containing a uracil has been determined only in the product state. The reactant state was constructed by reattaching uracil to the deoxyribose, and both complexes were studied by molecular dynamics simulations. Significant changes in the positions of secondary structural elements in the enzyme are induced by the hydrolysis of the glycosidic bond. The simulations show that the specificity of the uracil pocket in the enzyme is largely retained in both complexes with the exception of Asn-204, which has been identified as a residue that contributes to discrimination between uracil and cytosine. The hydrogen bond between the amide group of Asn-204 and O(4) of uracil is disrupted by fluctuations of the side chain in the reactant state and is replaced by a hydrogen bond to water molecules trapped in the interior of the protein behind the uracil binding pocket. The role of two residues implicated by mutation experiments to be important in catalysis, His-268 and Asp-145, is clarified by the simulations. In the reactant state, His-268 is found 3.45 +/- 0.34 A from the uracil, allowing a water molecule to form a bridge to O(2). The environment in the enzyme raises the pK(a) value of His-268 to 7.1, establishing a protonated residue for assisting in the hydrolysis of the glycosidic bond. In agreement with the crystallographic structure, the DNA backbone retracts after the hydrolysis to allow His-268 to approach the O(2) of uracil with a concomitant release of the bridging water molecule and a reduction in the pK(a) to 5.5, which releases the proton to the product. The side chain of Asp-145 is fully solvated in the reactant state and H-bonded through a water molecule to the 3'-phosphate of uridine. Both the proximity of Asp-145 to the negatively charged phosphate and its pK(a) of 4.4 indicate that it cannot act as a general base catalyst. We propose a mechanism in which the bridging water between Asp-145 and the 3'-phosphate accepts a proton from another water to stabilize the bridge through a hydronium ion as well as to produce the hydroxide anion required for the hydrolytic step. The mechanism is consistent with known experimental data.

Characterization of novel cathepsin K mutations in the pro and mature polypeptide regions causing pycnodysostosis

Cathepsin K, a lysosomal cysteine protease critical for bone remodeling by osteoclasts, was recently identified as the deficient enzyme causing pycnodysostosis, an autosomal recessive osteosclerotic skeletal dysplasia. To investigate the nature of molecular lesions causing this disease, mutations in the cathepsin K gene from eight families were determined, identifying seven novel mutations (K52X, G79E, Q190X, Y212C, A277E, A277V, and R312G). Expression of the first pro region missense mutation in a cysteine protease, G79E, in Pichia pastoris resulted in an unstable precursor protein, consistent with misfolding of the proenzyme. Expression of five mature region missense defects revealed that G146R, A277E, A277V, and R312G precursors were unstable, and no mature proteins or protease activity were detected. The Y212C precursor was activated to its mature form in a manner similar to that of the wild-type cathepsin K. The mature Y212C enzyme retained its dipeptide substrate specificity and gelatinolytic activity, but it had markedly decreased activity toward type I collagen and a cathepsin K-specific tripeptide substrate, indicating that it was unable to bind collagen triple helix. These studies demonstrated the molecular heterogeneity of mutations causing pycnodysostosis, indicated that pro region conformation directs proper folding of the proenzyme, and suggested that the cathepsin K active site contains a critical collagen-binding domain.

Human cathepsin V functional expression, tissue distribution, electrostatic surface potential, enzymatic characterization, and chromosomal localization

Cathepsin V, a thymus and testis-specific human cysteine protease, was expressed in Pichia pastoris, and its physicokinetic properties were determined. Recombinant procathepsin V is autocatalytically activated at acidic pH and is effectively inhibited by various cysteine protease class-specific inhibitors. The S2P2 subsite specificity of cathepsin V was found to be intermediate between those of cathepsins S and L. The substrate binding pocket, S2, accepted both aromatic and nonaromatic hydrophobic residues, whereas cathepsins L and S preferred either an aromatic or nonaromatic hydrophobic residue, respectively. In contrast to cathepsin L, but similar to cathepsin S, cathepsin V exhibited only a very weak collagenolytic activity. Furthermore, cathepsin V was determined to be significantly more stable at mildly acidic and neutral pH than cathepsin L, but distinctly less stable than cathepsin S. A homology structure model of cathepsin V revealed completely different electrostatic potentials on the molecular surface when compared with human cathepsin L. The model-based electrostatic potential of human cathepsin V was neutral to weakly positive at and in the vicinity of the active site cleft, whereas that of cathepsin L was negative over extended regions of the surface. Surprisingly, the electrostatic potential of the human cathepsin V model structure resembled that of the model structure of mouse cathepsin L. These differences in the electrostatic potential at the molecular surfaces provide a reactivity determinant that may be the source of differences in substrate selectivity and pH stability. Cathepsin V was mapped to the chromosomal region 9q22.2, a site adjacent to the cathepsin L locus. The high sequence identity and the overlapping chromosomal gene loci suggest that both proteases evolved from an ancestral cathepsin L-like precursor by gene duplication.

Structure and dynamics of calmodulin in solution

To characterize the dynamic behavior of calmodulin in solution, we have carried out molecular dynamics (MD) simulations of the Ca2+-loaded structure. The crystal structure of calmodulin was placed in a solvent sphere of radius 44 A, and 6 Cl- and 22 Na+ ions were included to neutralize the system and to model a 150 mM salt concentration. The total number of atoms was 32,867. During the 3-ns simulation, the structure exhibits large conformational changes on the nanosecond time scale. The central alpha-helix, which has been shown to unwind locally upon binding of calmodulin to target proteins, bends and unwinds near residue Arg74. We interpret this result as a preparative step in the more extensive structural transition observed in the "flexible linker" region 74-82 of the central helix upon complex formation. The major structural change is a reorientation of the two Ca2+-binding domains with respect to each other and a rearrangement of alpha-helices in the N-terminus domain that makes the hydrophobic target peptide binding site more accessible. This structural rearrangement brings the domains to a more favorable position for target binding, poised to achieve the orientation observed in the complex of calmodulin with myosin light-chain kinase. Analysis of solvent structure reveals an inhomogeneity in the mobility of water in the vicinity of the protein, which is attributable to the hydrophobic effect exerted by calmodulin's binding sites for target peptides.