Abundant Perithecial Protein (APP) from Neurospora is a primitive functional analog of ocular crystallins.

The eye arose during the Cambrian explosion from pre-existing proteins that would have been recruited for the formation of the specialized components of this organ, such as the transparent lens. Proteins suitable for the role of lens crystallins would need to possess unusual physical properties and the study of such earliest analogs of ocular crystallins would add to our understanding of the nature of recruitment of proteins as lens/corneal crystallins. We show that the Abundant Perithecial Protein (APP) of the fungi Neurospora and Sordaria fulfils the criteria for an early crystallin analog. The perithecia in these fungal species are phototropic, and APP accumulates at a high concentration in the neck of the pitcher-shaped perithecium. Spores are formed at the base of the perithecium, and light contributes to their maturation. The hydrodynamic properties of APP appear to exclude dimer formation or aggregation at high protein concentrations. APP is also deficient in Ca2+ binding, a property seen in its close homolog, the calcium-binding cell adhesion molecule (DdCAD-1) from Dictyostelium discoidum. Comparable to crystallins, APP occurs in high concentrations and seems to have dispensed with Ca2+ binding in exchange for greater stability. These crystallin-like attributes of APP lead us to demonstrate that it is a primitive form of ocular crystallins.

Structural study of the G57W mutant of human gamma-S-crystallin, associated with congenital cataract.

PURPOSE: Human γS-crystallin (CrygS) is an important component of the human eye lens nucleus and cortex. The mutation G57W in the molecule is reported to be associated with congenital cataract in children. We compare the conformational features and aggregation properties of the mutant protein G57W with the wild-type CrygS to understand how the structural changes in the mutant are related to the mechanism of opacification. METHODS: Wild-type and mutant proteins were cloned, expressed, and purified, and their structural properties were studied in solution. Conformational features and the structural stability of the proteins were compared in solution, using circular dichroism (CD) and fluorescence spectroscopic analysis, and the proteins' tendencies to aggregate were compared using extrinsic spectral probes. In addition, we analyzed the proteins' structural differences with extensive molecular modeling in silico. RESULTS: CD and intrinsic fluorescence analysis suggested the secondary and tertiary structures of the mutant are slightly altered. Experiments using extrinsic spectral probes revealed that the compact close-packed structure is loosened somewhat, and the mutant tends to self-aggregate. Denaturation (both thermal and chemical) studies indicate that the replacement of glycine (G) in position 57 by tryptophan (W) lowered the structural stability of the molecule. Further, the mutant had a tendency to precipitate and scatters light more easily than the wild-type. CONCLUSIONS: The replacement of glycine at position 57 by the tryptophan residue in human γS-crystallin weakens the stability of the mutant molecule and causes the molecule to self-aggregate, thus generating light-scattering particles. This set of changes in the mutant offers a molecular insight into the mechanism of opacification.

Gamma crystallins of the human eye lens.

BACKGROUND: Protein crystallins co me in three types (α, ß and γ) and are found predominantly in the eye, and particularly in the lens, where they are packed into a compact, plastic, elastic, and transparent globule of proper refractive power range that aids in focusing incoming light on to the retina. Of these, the γ-crystallins are found largely in the nuclear region of the lens at very high concentrations (>400 mg/ml). The connection between their structure and inter-molecular interactions and lens transparency is an issue of particular interest. SCOPE OF REVIEW: We review the origin and phylogeny of the gamma crystallins, their special structure involving the use of Greek key supersecondary structural motif, and how they aid in offering the appropriate refractive index gradient, intermolecular short range attractive interactions (aiding in packing them into a transparent ball), the role that several of the constituent amino acid residues play in this process, the thermodynamic and kinetic stability and how even single point mutations can upset this delicate balance and lead to intermolecular aggregation, forming light-scattering particles which compromise transparency. We cite several examples of this, and illustrate this by cloning, expressing, isolating and comparing the properties of the mutant protein S39C of human γS-crystallin (associated with congenital cataract-microcornea), with those of the wild type molecule. In addition, we note that human γ-crystallins are also present in other parts of the eye (e.g., retina), where their functions are yet to be understood. MAJOR CONCLUSIONS: There are several 'crucial' residues in and around the Greek key motifs which are essential to maintain the compact architecture of the crystallin molecules. We find that a mutation that replaces even one of these residues can lead to reduction in solubility, formation of light-scattering particles and loss of transparency in the molecular assembly. GENERAL SIGNIFICANCE: Such a molecular understanding of the process helps us construct the continuum of genotype-molecular structural phenotype-clinical (pathological) phenotype. This article is part of a Special Issue entitled Crystallin Biochemistry in Health and Disease.

Structural integrity of the Greek key motif in ßγ-crystallins is vital for central eye lens transparency.

BACKGROUND: We highlight an unrecognized physiological role for the Greek key motif, an evolutionarily conserved super-secondary structural topology of the ßγ-crystallins. These proteins constitute the bulk of the human eye lens, packed at very high concentrations in a compact, globular, short-range order, generating transparency. Congenital cataract (affecting 400,000 newborns yearly worldwide), associated with 54 mutations in ßγ-crystallins, occurs in two major phenotypes nuclear cataract, which blocks the central visual axis, hampering the development of the growing eye and demanding earliest intervention, and the milder peripheral progressive cataract where surgery can wait. In order to understand this phenotypic dichotomy at the molecular level, we have studied the structural and aggregation features of representative mutations. METHODS: Wild type and several representative mutant proteins were cloned, expressed and purified and their secondary and tertiary structural details, as well as structural stability, were compared in solution, using spectroscopy. Their tendencies to aggregate in vitro and in cellulo were also compared. In addition, we analyzed their structural differences by molecular modeling in silico. RESULTS: Based on their properties, mutants are seen to fall into two classes. Mutants A36P, L45PL54P, R140X, and G165fs display lowered solubility and structural stability, expose several buried residues to the surface, aggregate in vitro and in cellulo, and disturb/distort the Greek key motif. And they are associated with nuclear cataract. In contrast, mutants P24T and R77S, associated with peripheral cataract, behave quite similar to the wild type molecule, and do not affect the Greek key topology. CONCLUSION: When a mutation distorts even one of the four Greek key motifs, the protein readily self-aggregates and precipitates, consistent with the phenotype of nuclear cataract, while mutations not affecting the motif display 'native state aggregation', leading to peripheral cataract, thus offering a protein structural rationale for the cataract phenotypic dichotomy "distort motif, lose central vision".

Structural analysis of the mutant protein D26G of human γS-crystallin, associated with Coppock cataract.

PURPOSE: To analyze the protein structural features responsible for the aggregation properties of the mutant protein D26G human γS-crystallin (HGSC) associated with congenital Coppock-type cataract. METHODS: cDNAs of wild-type (WT) and D26G mutant HGSC were cloned and expressed in BL21 (DE3) pLysS cells and the proteins isolated and purified. Their secondary and tertiary structural features, aggregation tendencies, and structural stabilities were compared using spectroscopic (circular dichroism, intrinsic and extrinsic fluorescence), molecular modeling, and dynamics methods. RESULTS: No difference was observed between the conformational (secondary and tertiary structural) features and aggregation properties between the WT and D26G proteins. The mutant, however, was structurally less stable; it denatured at a slightly lower concentration of the added chemical denaturant (at 2.05 M guanidinium chloride, cf. 2.20 M for the WT) and at a slightly lower temperature (at 70.8 °C, cf. 72.0 °C for the WT). The mutant also self-aggregated more readily (it turned turbid upon standing; at 65 °C, it started precipitating beyond 200 s, while the WT did not, even after 900 s). Molecular modeling showed that the Asp26-Arg84 contact (and the related Arg84-Asn54 interaction) was disturbed in the mutant, making the latter less compact around the mutation site. CONCLUSIONS: The cataract-associated mutant D26G of HGSC is remarkably close to the WT molecule in structural features, with only a microenvironmental change in the packing around the mutation site. This alteration appears sufficient to promote self-aggregation, resulting in peripheral cataract.

Structural model of the Y-Family DNA polymerase V/RecA mutasome.

To synthesize past DNA damaged by chemicals or radiation, cells have lesion bypass DNA polymerases (DNAPs), most of which are in the Y-Family. One class of Y-Family DNAPs includes DNAP η in eukaryotes and DNAP V in bacteria, which have low fidelity when replicating undamaged DNA. In Escherchia coli, DNAP V is carefully regulated to insure it is active for lesion bypass only, and one mode of regulation involves interaction of the polymerase subunit (UmuC) and two regulatory subunits (UmuD') with a RecA-filament bound to ss-DNA. Taking a docking approach, ∼150,000 unique orientations involving UmuC, UmuD' and RecA were evaluated to generate models, one of which was judged best able to rationalize the following published findings. (1) In the UmuD'(2)C/RecA-filament model, R64-UmuC interacts with S117-RecA, which is known to be at the UmuC/RecA interface. (2) At the model's UmuC/RecA interface, UmuC has three basic amino acids (K59/R63/R64) that anchor it to RecA. No other Y-Family DNAP has three basic amino acids clustered in this region, making it a plausible site for UmuC to form its unique interaction with RecA. (3) In the model, residues N32/N33/D34 of UmuC form a second interface with RecA, which is consistent with published findings. (4) Active UmuD' is generated when 24 amino acids in the N-terminal tail of UmuD are proteolyzed, which occurs when UmuD(2)C binds the RecA-filament. When UmuD is included in an UmuD(2)C/RecA-filament model, plausible UmuD/RecA contacts guide the UmuD cleavage site (C24/G25) into the UmuD proteolysis active site (S60/K97). One contact involves E11-UmuD interacting with R243-RecA, where the latter is known to be important for UmuD cleavage. (5) The UmuD(2)C/RecA-filament model rationalizes published findings that at least some UmuD-to-UmuD' cleavage occurs intermolecularly. (6) Active DNAP V is known to be the heterotetramer UmuD'(2)C/RecA, a model of which can be generated by a simple rearrangement of the RecA monomer at the 3'-end of the RecA-filament. The rearranged UmuD'(2)C/RecA model rationalizes published findings about UmuD' residues in proximity to RecA. In summary, docking and molecular simulations are used to develop an UmuD'(2)C/RecA model, whose structure rationalizes much of the known properties of the active form of DNA polymerase V.

The mutation V42M distorts the compact packing of the human gamma-S-crystallin molecule, resulting in congenital cataract.

BACKGROUND: Human γS-crystallin is an important component of the human eye lens nucleus and cortex. The mutation V42M in the molecule causes severe congenital cataract in children. We compare the structure of the mutant protein with that of the wild type in order to understand how structural changes in the mutant relate to the mechanism of opacification. METHODS: Both proteins were made using conventional cloning and expression procedures. Secondary and tertiary structural features of the proteins were analyzed using spectral methods. Structural stabilities of the proteins were analyzed using chemical and thermal denaturation methods. Self-aggregation was monitored using extrinsic spectral probes. Molecular modeling was used to compare the structural features of the two proteins. RESULTS: While the wild type and mutant have the same secondary structure, molecular modeling and fluorescence analysis suggest the mutant to have a more open tertiary structure, with a larger hydrophobic surface. Experiments using extrinsic probes reveal that the mutant readily self-aggregates, with the suggestion that the aggregates might be similar to amyloidogenic fibrils. Chemical denaturation indicates that while the wild type exhibits the classic two-state transition, V42M goes through an intermediate state, and has a distinctly lower stability than the wild type. The temperature of thermal unfolding of the mutant is also distinctly lower. Further, the mutant readily precipitates and scatters light more easily than the wild type. CONCLUSION: The replacement of valine in position 42 by the longer and bulkier methionine in human γS-crystallin perturbs the compact ß-sheet core packing topology in the N-terminal domain of the molecule, exposes nonpolar residues thereby increasing the surface hydrophobicity and weakens the stability of the protein, thus promoting self-aggregation leading to light scattering particles. This set of changes in the properties of the mutant offers a molecular insight into the mechanism of opacification.

Two polymorphic residues account for the differences in DNA binding and transcriptional activation by NF-κB proteins encoded by naturally occurring alleles in Nematostella vectensis.

The NF-κB family of transcription factors is activated in response to many environmental and biological stresses, and plays a key role in innate immunity across a broad evolutionary expanse of animals. A simple NF-κB pathway is present in the sea anemone Nematostella vectensis, an important model organism in the phylum Cnidaria. Nematostella has previously been shown to have two naturally occurring NF-κB alleles (Nv-NF-κB-C and Nv-NF-κB-S) that encode proteins with different DNA-binding and transactivation abilities. We show here that polymorphic residues 67 (Cys vs. Ser) and 269 (Ala vs. Glu) play complementary roles in determining the DNA-binding activity of the NF-κB proteins encoded by these two alleles and that residue 67 is primarily responsible for the difference in their transactivation ability. Phylogenetic analysis indicates that Nv-NF-κB-S is the derived allele, consistent with its restricted geographic distribution. These results define polymorphic residues that are important for the DNA-binding and transactivating activities of two naturally occurring variants of Nv-NF-κB. The implications for the appearance of the two Nv-NF-κB alleles in natural populations of sea anemones are discussed.

Stabilization of the predominant disease-causing aldolase variant (A149P) with zwitterionic osmolytes.

Hereditary fructose intolerance (HFI) is a disease of carbohydrate metabolism that can result in hyperuricemia, hypoglycemia, liver and kidney failure, coma, and death. Currently, the only treatment for HFI is a strict fructose-free diet. HFI arises from aldolase B deficiency, and the most predominant HFI mutation is an alanine to proline substitution at position 149 (A149P). The resulting aldolase B with the A149P substitution (AP-aldolase) has activity that is <100-fold that of the wild type. The X-ray crystal structure of AP-aldolase at both 4 and 18 °C reveals disordered adjacent loops of the (α/ß)(8) fold centered around the substitution, which leads to a dimeric structure as opposed to the wild-type tetramer. The effects of osmolytes were tested for restoration of structure and function. An initial screen of osmolytes (glycerol, sucrose, polyethylene glycol, 2,4-methylpentanediol, glutamic acid, arginine, glycine, proline, betaine, sarcosine, and trimethylamine N-oxide) reveals that glycine, along with similarly structured compounds, betaine and sarcosine, protects AP-aldolase structure and activity from thermal inactivation. The concentration and functional moieties required for thermal protection show a zwitterion requirement. The effects of osmolytes in restoring structure and function of AP-aldolase are described. Testing of zwitterionic osmolytes of increasing size and decreasing fractional polar surface area suggests that osmolyte-mediated AP-aldolase stabilization occurs neither primarily through excluded volume effects nor through transfer free energy effects. These data suggest that AP-aldolase is stabilized by binding to the native structure, and they provide a foundation for developing stabilizing compounds for potential therapeutics for HFI.

Architecture of y-family DNA polymerases relevant to translesion DNA synthesis as revealed in structural and molecular modeling studies.

DNA adducts, which block replicative DNA polymerases (DNAPs), are often bypassed by lesion-bypass DNAPs, which are mostly in the Y-Family. Y-Family DNAPs can do non-mutagenic or mutagenic dNTP insertion, and understanding this difference is important, because mutations transform normal into tumorigenic cells. Y-Family DNAP architecture that dictates mechanism, as revealed in structural and modeling studies, is considered. Steps from adduct blockage of replicative DNAPs, to bypass by a lesion-bypass DNAP, to resumption of synthesis by a replicative DNAP are described. Catalytic steps and protein conformational changes are considered. One adduct is analyzed in greater detail: the major benzo[a]pyrene adduct (B[a]P-N(2)-dG), which is bypassed non-mutagenically (dCTP insertion) by Y-family DNAPs in the IV/κ-class and mutagenically (dATP insertion) by V/η-class Y-Family DNAPs. Important architectural differences between IV/κ-class versus V/η-class DNAPs are discussed, including insights gained by analyzing ~400 sequences each for bacterial DNAPs IV and V, along with sequences from eukaryotic DNAPs kappa, eta and iota. The little finger domains of Y-Family DNAPs do not show sequence conservation; however, their structures are remarkably similar due to the presence of a core of hydrophobic amino acids, whose exact identity is less important than the hydrophobic amino acid spacing.

Amino acid architecture that influences dNTP insertion efficiency in Y-family DNA polymerase V of E. coli.

Y-family DNA polymerases (DNAPs) are often required in cells to synthesize past DNA-containing lesions, such as [+ta]-B[a]P-N(2)-dG, which is the major adduct of the potent mutagen/carcinogen benzo[a]pyrene. The current model for the non-mutagenic pathway in Escherichia coli involves DNAP IV inserting deoxycytidine triphosphate opposite [+ta]-B[a]P-N(2)-dG and DNAP V doing the next step(s), extension. We are investigating what structural differences in these related Y-family DNAPs dictate their functional differences. X-ray structures of Y-family DNAPs reveal a number of interesting features in the vicinity of the active site, including (1) the "roof-amino acid" (roof-aa), which is the amino acid that lies above the nucleobase of the deoxynucleotide triphosphate (dNTP) and is expected to play a role in dNTP insertion efficiency, and (2) a cluster of three amino acids, including the roof-aa, which anchors the base of a loop, whose detailed structure dictates several important mechanistic functions. Since no X-ray structures existed for UmuC (the polymerase subunit of DNAP V) or DNAP IV, we previously built molecular models. Herein, we test the accuracy of our UmuC(V) model by investigating how amino acid replacement mutants affect lesion bypass efficiency. A ssM13 vector containing a single [+ta]-B[a]P-N(2)-dG is transformed into E. coli carrying mutations at I38, which is the roof-aa in our UmuC(V) model, and output progeny vector yield is monitored as a measure of the relative efficiency of the non-mutagenic pathway. Findings show that (1) the roof-aa is almost certainly I38, whose beta-carbon branching R-group is key for optimal activity, and (2) I38/A39/V29 form a hydrophobic cluster that anchors an important mechanistic loop, aa29-39. In addition, bypass efficiency is significantly lower both for the I38A mutation of the roof-aa and for the adjacent A39T mutation; however, the I38A/A39T double mutant is almost as active as wild-type UmuC(V), which probably reflects the following. Y-family DNAPs fall into several classes with respect to the [roof-aa/next amino acid]: one class has [isoleucine/alanine] and includes UmuC(V) and DNAP eta (from many species), while the second class has [alanine (or serine)/threonine] and includes DNAP IV, DNAP kappa (from many species), and Dpo4. Thus, the high activity of the I38A/A39T double mutant probably arises because UmuC(V) was converted from the V/eta class to the IV/kappa class with respect to the [roof-aa/next amino acid]. Structural and mechanistic aspects of these two classes of Y-family DNAPs are discussed.

Y-Family DNA polymerases may use two different dNTP shapes for insertion: a hypothesis and its implications.

Chemicals and radiation can damage DNA leading to the formation of adducts/lesions, which - if not removed by DNA repair pathways - usually block replicative DNA polymerases (DNAPs). To overcome such potentially lethal blockage, cells have lesion bypass DNAPs, which are often in the Y-Family and include several classes. One class includes human DNAP kappa and E. coli DNAP IV, and they insert dCTP in the non-mutagenic pathway opposite [+ta]-B[a]P-N(2)-dG, which is the major adduct formed by the environmental carcinogen benzo[a]pyrene. Another class includes hDNAP eta and ecDNAP V, and they insert dATP opposite [+ta]-B[a]P-N(2)-dG in the dominant G-->T mutagenic pathway. Herein we develop a hypothesis for why the IV/kappa-class preferentially does cellular dCTP insertion. On the minor groove side of the active site, Y-Family DNAPs have a cleft/hole that can be analyzed based on an analogy to a "chimney." Our models of DNAP IV show a large chimney opening from which the pyrene of [+ta]-B[a]P-N(2)-dG can protrude, which allows canonical adduct-dG:dCTP pairing. In contrast, our models of DNAP V have small chimney openings that forces adduct-dG downward in the active site such that canonical adduct-dG:dCTP pairing is not possible. Based on X-ray structures, sequence alignment and our modeled structures of Y-Family DNAPs, chimney opening size seems primarily controlled by one amino acid ("flue-handle"), which dictates whether nearby amino acids ("flue") plug the chimney or not. Based on this analysis, a correlation is apparent: the flue is closed in V/eta-class DNAPs giving small chimney openings, while the flue is open for the IV/kappa-class giving large chimney openings. Secondarily, a hypothesis is developed for why the V/eta-class might preferentially do cellular dATP insertion opposite [+ta]-B[a]P-N(2)-dG: the small chimney forces adduct-dG lower in the active site, possibly leading to catalysis using a non-canonical dNTP shape that permits syn-adenine:adduct-dG base pairing. In summary, a hypothesize is developed that the pyrene moiety of [+ta]-B[a]P-N(2)-dG protrudes from the large chimney opening of DNAP IV, thus permitting canonical dCTP:adduct-dG pairing, while the small chimney opening of DNAP V forces [+ta]-B[a]P-N(2)-dG lower down in the active site, in which syn-adenine can pair with adduct-dG via a non-canonical dNTP shape.

Intermolecular disulfide bond formation in the NEMO dimer requires Cys54 and Cys347.

NEMO is an essential regulatory component of the IkappaB kinase (IKK) complex, which controls activation of the NF-kappaB signaling pathway. Herein, we show that NEMO exists as a disulfide-bonded dimer when isolated from several cell types and analyzed by SDS-polyacrylamide gel electrophoresis under non-reducing conditions. Treatment of cells with hydrogen peroxide (H(2)O(2)) induces further formation of NEMO dimers. Disulfide bond-mediated formation of NEMO dimers requires Cys54 and Cys347. The ability of these residues to form disulfide bonds is consistent with their location in a NEMO dimer structure that we generated by molecular modeling. We also show that pretreatment with H(2)O(2) decreases TNFalpha-induced IKK activity in NEMO-reconstituted cells, and that TNFalpha has a diminished ability to activate NF-kappaB DNA binding in cells reconstituted with NEMO mutant C54/347A. This study implicates NEMO as a target of redox regulation and presents the first structural model for the NEMO protein.

Molecular modeling benzo[a]pyrene N2-dG adducts in the two overlapping active sites of the Y-family DNA polymerase Dpo4.

The potent, ubiquitous environmental mutagen/carcinogen benzo[a]pyrene (B[a]P) induces a single major adduct [+ta]-B[a]P-N2-dG, whose bypass in most cases results in either no mutation (dCTP insertion) or a G-->T mutation (dATP insertion). Translesion synthesis (TLS) of [+ta]-B[a]P-N2-dG generally requires DNA polymerases (DNAPs) in the Y-family, which exist in cells to bypass DNA damage caused by chemicals and radiation. A molecular dynamics (MD) study is described with dCTP opposite [+ta]-B[a]P-N2-dG in Dpo4, which is the best studied Y-family DNAP from a structural point of view. Two orientations of B[a]P-N2-dG (BPmi5 and BPmi3) are considered, along with two orientations of the dCTP (AS1 and AS2), as outlined next. Based on NMR studies, the pyrene moiety of B[a]P-N2-dG is in the minor groove, when paired with dC, and can point toward either the base on the 5'-side (BPmi5) or the 3'-side (BPmi3). Based on published X-ray structures, Dpo4 appears to have two partially overlapping active sites. The architecture of active site 1 (AS1) is similar to all other families of DNAPs (e.g., the shape of the dNTP). Active site 2 (AS2), however, is non-canonical (e.g., the beta- and gamma-phosphates in AS2 are approximately where the alpha- and beta-phosphates are in AS1). In the Dpo4 models generated herein, using the BPmi3 orientation the pyrene moiety of [+ta]-B[a]P-N2-dG points toward the duplex region of the DNA, and is accommodated without distortions in AS1, but with distortions in AS2. Considering the BPmi5 orientation, the pyrene moiety points toward the ss-region of DNA in Dpo4, and sits in a hole defined by the fingers and little fingers domain ("chimney"); BPmi5 is accommodated in AS2 without significant distortions, but poorly in AS1. In summary, when dCTP is paired with [+ta]-B[a]P-N2-dG in the two overlapping active sites in Dpo4, the pyrene in the BPmi3 orientation is accommodated better in active site 1 (AS1), while the pyrene in the BPmi5 orientation is accommodated better in AS2. Finally, we discuss why Y-family DNAPs might have two catalytic active sites.

Genetic effects of oxidative DNA damages: comparative mutagenesis of the imidazole ring-opened formamidopyrimidines (Fapy lesions) and 8-oxo-purines in simian kidney cells.

Fapy.dG and 8-oxo-7,8-dihydro-2'-deoxyguanosine (8-oxo-dG) are formed in DNA by hydroxyl radical damage. In order to study replication past these lesions in cells, we constructed a single-stranded shuttle vector containing the lesion in 5'-TGT and 5'-TGA sequence contexts. Replication of the modified vector in simian kidney (COS-7) cells showed that Fapy.dG is mutagenic inducing primarily targeted Fapy.G-->T transversions. In the 5'-TGT sequence mutational frequency of Fapy.dG was approximately 30%, whereas in the 5'-TGA sequence it was approximately 8%. In parallel studies 8-oxo-dG was found to be slightly less mutagenic than Fapy.dG, though it also exhibited a similar context effect: 4-fold G-->T transversions (24% versus 6%) occurred in the 5'-TGT sequence relative to 5'-TGA. To investigate a possible structural basis for the higher G-->T mutations induced by both lesions when their 3' neighbor was T, we carried out a molecular modeling investigation in the active site of DNA polymerase beta, which is known to incorporate both dCTP (no mutation) and dATP (G-->T substitution) opposite 8-oxo-G. In pol beta, the syn-8-oxo-G:dATP pair showed greater stacking with the 3'-T:A base pair in the 5'-TGT sequence compared with the 3'-A:T in the 5'-TGA sequence, whereas stacking for the anti-8-oxo-G:dCTP pair was similar in both 5'-TGT and 5'-TGA sequences. Similarly, syn-Fapy.G:dATP pairing showed greater stacking in the 5'-TGT sequence compared with the 5'-TGA sequence, while stacking for anti-Fapy.G:dCTP pairs was similar in the two sequences. Thus, for both lesions less efficient base stacking between the lesion:dATP pair and the 3'-A:T base pair in the 5'-TGA sequence might cause lower G-->T mutational frequencies in the 5'-TGA sequence compared to 5'-TGT. The corresponding lesions derived from 2'-deoxyadenosine, Fapy.dA and 8-oxo-dA, were not detectably mutagenic in the 5'-TAT sequence, and were only weakly mutagenic (<1%) in the 5'-TAA sequence context, where both lesions induced targeted A-->C transversions. To our knowledge this is the first investigation using extrachromosomal probes containing a Fapy.dG or Fapy.dA site-specifically incorporated, which showed unequivocally that in simian kidney cells Fapy.G-->T substitutions occur at a higher frequency than 8-oxo-G-->T and that Fapy.dA is very weakly mutagenic, as is 8-oxo-dA.

Mirror image stereoisomers of the major benzo[a]pyrene N2-dG adduct are bypassed by different lesion-bypass DNA polymerases in E. coli.

The potent mutagen/carcinogen benzo[a]pyrene (B[a]P) is metabolically activated to (+)-anti-B[a]PDE, which induces a full spectrum of mutations (e.g., G-to-T, G-to-A, -1 frameshifts, etc.) via its major adduct [+ta]-B[a]P-N2-dG. We recently showed that the dominant G-to-T mutation depends on DNA polymerase V (DNAP V), but not DNAPs IV or II, when studied in a 5'-TG sequence in E. coli. Herein we investigate what DNAPs are responsible for non-mutagenic bypass with [+ta]-B[a]P-N2-dG, along with its mirror image adduct [-ta]-B[a]P-N2-dG. Each adduct is built into a 5'-TG sequence in a single stranded M13 phage vector, which is then transformed into eight different E. coli strains containing all combinations of proficiency and deficiency in the three lesion-bypass DNAPs II, IV and V. Based on M13 progeny output, non-mutagenic bypass with [-ta]-B[a]P-N2-dG depends on DNAP IV. In contrast, non-mutagenic bypass with [+ta]-B[a]P-N2-dG depends on both DNAPs IV and V, where arguments suggest that DNAP IV is involved in dCTP insertion, while DNAP V is involved in extension of the adduct-G:C base pair. Numerous findings indicate that DNAP II has a slight inhibitory effect on the bypass of [+ta]- and [-ta]-B[a]P-N2-dG in the case of both DNAPs IV and V. In conclusion, for efficient non-mutagenic bypass (dCTP insertion) in E. coli, [+ta]-B[a]P-N2-dG requires DNAPs IV and V, [-ta]-B[a]P-N2-dG requires only DNAP IV, while DNAP II is inhibitory to both, and experiments to investigate these differences should provide insights into the mechanism and purpose of these lesion-bypass DNAPs.

Homology modeling of four Y-family, lesion-bypass DNA polymerases: the case that E. coli Pol IV and human Pol kappa are orthologs, and E. coli Pol V and human Pol eta are orthologs.

Y-family DNA polymerases (DNAPs) are a superfamily of evolutionarily related proteins that exist in cells to bypass DNA damage caused by both radiation and chemicals. Cells have multiple Y-family DNAPs, presumably to conduct translesion synthesis (TLS) on DNA lesions of varying structure and conformation. The potent, ubiquitous environmental mutagen/carcinogen benzo[a]pyrene (B[a]P) induces all classes of mutations with G-->T base substitutions predominating. We recently showed that a G-->T mutagenesis pathway for the major adduct of B[a]P ([+ta]-B[a]P-N2-dG) in Escherichia coli depends on Y-family member DNAP V. Since no X-ray crystal study for DNAP V has been reported, no structure is available to help in understanding the structural basis for dATP insertion associated with G-->T mutations from [+ta]-B[a]P-N2-dG. Herein, we do homology modeling to construct a model for UmuC, which is the polymerase subunit of DNAP V. The sequences of eight Y-family DNAPs were aligned based on the positioning of conserved amino acids and an analysis of conserved predicted secondary structure, as well as insights gained from published X-ray structures of five Y-family members. Starting coordinates for UmuC were generated from the backbone coordinates for the Y-family polymerase Dpo4 for reasons discussed, and were refined using molecular dynamics with CHARMM 27. A survey of the literature revealed that E. coli DNAP V and human DNAP eta show a similar pattern of dNTP insertion opposite a variety of DNA lesions. Furthermore, E. coli DNAP IV and human DNAP kappa show a similar dNTP insertional pattern with these same DNA lesions, although the insertional pattern for DNAP IV/kappa differs from the pattern for DNAPs V/eta. These comparisons prompted us to construct and refine models for E. coli DNAP IV and human DNAPs eta and kappa as well. The dNTP/template binding pocket of all four DNAPs was inspected, focusing on the array of seven amino acids that contact the base of the incoming dNTP, as well as the template base. DNAPs V and eta show similarities in this array, and DNAPs IV and kappa also show similarities, although the arrays are different for the two pairs of DNAPs. Thus, there is a correlation between structural similarities and insertional similarities for the pairs DNAPs V/eta and DNAPs IV/kappa. Although the significance of this correlation remains to be elucidated, these observations point the way for future experimental studies.

Free-energy perturbation methods to study structure and energetics of DNA adducts: results for the major N2-dG adduct of benzo[a]pyrene in two conformations and different sequence contexts.

The potent mutagen/carcinogen benzo[a]pyrene (B[a]P) is activated to (+)-anti-B[a]PDE, which induces a variety of mutations (e.g., G --> T, G --> A, etc.) via its major adduct [+ta]-B[a]P-N2-dG. One hypothesis is that adducts (such as [+ta]-B[a]P-N2-dG) induce different mutations via different conformations, probably when replicated by different lesion-bypass DNA polymerases (DNAPs). We showed that Escherichia coli DNAP V was responsible for G --> T mutations with [+ta]-B[a]P-N2-dG in a 5'-TGT sequence (Yin et al., (2004) DNA Repair 3, 323), so we wish to study conformations of this adduct/sequence context by molecular modeling. The development of a CHARMM-based molecular dynamics (MD) simulations protocol with free-energy calculations in the presence of solvent and counterions is described. A representative base-pairing and base-displaced conformation of [+ta]-B[a]P-N2-dG in the 5'-TGT sequence are used: (1) BPmi5, which has the B[a]P moiety in the minor groove pointing toward the base on the 5'-side of the adduct, and (2) Gma5, which has the B[a]P moiety stacked with the surrounding base pairs and the dG moiety displaced into the major groove. The MD output structures are reasonable when compared to known NMR structures. Changes in DNA sequence context dramatically affect the biological consequences (e.g., mutagenesis) of [+ta]-B[a]P-N2-dG. Consequently, we also developed a MD-based free-energy perturbation (FEP) protocol to study DNA sequence changes. FEP involves the gradual "fading-out" of atoms in a starting structure (A) and "fading-in" of atoms in a final structure (B), which allows a realistic assessment of the energetic and structural changes when two structures A and B are closely related. Two DNA sequence changes are described: (1) 5'-TGT --> 5'-TGG, which involves two steps [T:A --> T:C --> G:C], and (2) 5'-TGT --> 5'-TGC, which involves three steps [T:A --> T:2AP --> C:2AP --> C:G], where 2AP (2-aminopurine) is included, because T:2AP and C:2AP retain more-or-less normal pairing orientations between complementary bases. FEP is also used to evaluate the impact that a 5'-TGT to 5'-UGT sequence change might have on mutagenesis with [+ta]-B[a]P-N2-dG. In summary, we developed (1) a CHARMM-based molecular dynamics (MD) simulations protocol with free-energy calculations in the presence of solvent and counterions to study B[a]P-N2-dG adducts in DNA duplexes, and (2) a MD-based free-energy perturbation (FEP) protocol to study DNA sequence context changes around B[a]P-N2-dG adducts.

Molecular modeling of four stereoisomers of the major B[a]PDE adduct (at N(2)-dG) in five cases where the structure is known from NMR studies: molecular modeling is consistent with NMR results.

The potent mutagen/carcinogen benzo[a]pyrene (B[a]P) is metabolically activated to (+)-anti-B[a]PDE, which is known to induce a variety of mutations (e.g., GC --> TA, GC --> AT, etc.). One hypothesis for this complexity is that different mutations are induced by different conformations of its major adduct [+ta]-B[a]P-N(2)-dG when bypassed during DNA replication (perhaps by different DNA polymerases). Our previous molecular modeling studies have suggested that conformational complexity might be extensive in that B[a]P-N(2)-dG adducts appeared capable of adopting at least sixteen potential conformational classes in ds-DNA [e.g., Kozack and Loechler (1999) Carcinogenesis 21, 1953], although only eight seemed likely to be relevant to base substitution mutagenesis. Such molecular modeling studies are only likely to be valuable for the interpretation of mutagenesis results if global minimum energy conformations for adducts are found and if the differences in the energies of these different conformations can be computed reasonably accurately. One approach to assessing the reliability of our molecular modeling techniques is considered herein. Using a five-step molecular modeling protocol, which importantly included a molecular dynamics version of simulated annealing, eight conformations are studied in each of five cases. (The five cases are listed below, and were chosen because in each case the preferred solution conformation is known from a NMR study.) Of the eight conformations studied, the one computed to be lowest in energy is the same conformation as the one observed by NMR in four of the five cases: 5'-CGC sequence with [+ta]-, [-ta]-, and [+ca]-B[a]P-N(2)-dG, and 5'-TGC sequence with [+ta]-B[a]P-N(2)-dG. In the fifth case (5'-CGC sequence with [-ca]-B[a]P-N(2)-dG), the known NMR conformation is computed to be second lowest in energy, but it is within approximately 1.7 kcal of the computed lowest energy conformation. These results suggest that molecular modeling is surprisingly accurate in computing lowest energy conformations and that it should be useful in assessing the relative energies of different conformations. This is especially important given that currently molecular modeling is the only means available to study the energetics of minor conformations of DNA adducts.

Molecular and cellular assessment of ginkgo biloba extract as a possible ophthalmic drug.

We have investigated the biochemical and cell biological basis of the reported beneficiary effects of the leaf extracts of the plant Ginkgo biloba, which has been used as a possible ophthalmic drug. The antioxidant, antimicrobial, anti-apoptotic and cytoprotective properties of the standardized extract called EGb761 were assayed. Chemical stresses were induced in cells using alloxan or dexamethasone, and the effect of EGb761 on them was studied using the MTT and TUNEL assays. Its ability to modulate the activities of some antioxidant enzymes was tested in vitro. In addition, cataract was induced in rats through selenite injection, and the effect of EGb761 administration on the progression of cataract was studied using slit lamp examination. Ginkgo biloba was found to be an excellent antioxidant. It readily scavenges reactive oxygen and nitrogen radicals and inhibits oxidative modifications that occur to proteins in vitro. It enters intact cells and protects them from alloxan-mediated and light-mediated stress, and the nuclear DNA from single strand breaks. It also effectively inhibits chemically induced apoptosis. It does not modulate the activities of endogenous antioxidant enzymes, nor does it have any significant antimicrobial activity. Unlike some other plant extracts, it is not phototoxic. In experiments wherein selenite cataract was induced in laboratory rats, treatment with the extract significantly retards the progression of lens opacification in vivo. Ginkgo biloba's inherent antioxidant, antiapoptotic and cytoprotective action and potential anticataract ability appear to be some of the factors responsible for its beneficial effects.