A molecular level picture of the stabilization of A-DNA in mixed ethanol–water solutions

Advances in computer power, methodology, and empirical force fields now allow routine “stable” nanosecond-length molecular dynamics simulations of DNA in water. The accurate representation of environmental influences on structure remains a major, unresolved issue. In contrast to simulations of A-DNA in water (where an A-DNA to B-DNA transition is observed) and in pure ethanol (where disruption of the structure is observed), A-DNA in ≈85% ethanol solution remains in a canonical A-DNA geometry as expected. The stabilization of A-DNA by ethanol is likely due to disruption of the spine of hydration in the minor groove and the presence of ion-mediated interhelical bonds and extensive hydration across the major groove.

Molecular switch in signal transduction: Reaction paths of the conformational changes in ras p21

Conformational changes in ras p21 triggered by the hydrolysis of GTP play an essential role in the signal transduction pathway. The path for the conformational change is determined by molecular dynamics simulation with a holonomic constraint directing the system from the known GTP-bound structure (with the γ-phosphate removed) to the GDP-bound structure. The simulation is done with a shell of water molecules surrounding the protein. In the switch I region, the side chain of Tyr-32, which undergoes a large displacement, moves through the space between loop 2 and the rest of the protein, rather than on the outside of the protein. As a result, the charged residues Glu-31 and Asp-33, which interact with Raf in the homologous RafRBD–Raps complex, remain exposed during the transition. In the switch II region, the conformational changes of α2 and loop 4 are strongly coupled. A transient hydrogen bonding complex between Arg-68 and Tyr-71 in the switch II region and Glu-37 in switch I region stabilizes the intermediate conformation of α2 and facilitates the unwinding of a helical turn of α2 (residues 66–69), which in turn permits the larger scale motion of loop 4. Hydrogen bond exchange between the protein and solvent molecules is found to be important in the transition. Possible functional implications of the results are discussed.

A transgenic mouse model of metastatic prostate cancer originating from neuroendocrine cells

A transgenic mouse model of metastatic prostate cancer has been developed that is 100% penetrant in multiple pedigrees. Nucleotides −6500 to +34 of the mouse cryptdin-2 gene were used to direct expression of simian virus 40 T antigen to a subset of neuroendocrine cells in all lobes of the FVB/N mouse prostate. Transgene expression is initiated between 7 and 8 weeks of age and leads to development of prostatic intraepithelial neoplasia within a week. Prostatic intraepithelial neoplasia progresses rapidly to local invasion. Metastases to lymph nodes, liver, lung, and bone are common by 6 months. Tumorigenesis is not dependent on androgens. This model indicates that the neuroendocrine cell lineage of the prostate is exquisitely sensitive to transformation and provides insights about the significance of neuroendocrine differentiation in human prostate cancer.

PCR candidate region mismatch scanning: adaptation to quantitative, high-throughput genotyping

Linkage and association analyses were performed to identify loci affecting disease susceptibility by scoring previously characterized sequence variations such as microsatellites and single nucleotide polymorphisms. Lack of markers in regions of interest, as well as difficulty in adapting various methods to high-throughput settings, often limits the effectiveness of the analyses. We have adapted the Escherichia coli mismatch detection system, employing the factors MutS, MutL and MutH, for use in PCR-based, automated, high-throughput genotyping and mutation detection of genomic DNA. Optimal sensitivity and signal-to-noise ratios were obtained in a straightforward fashion because the detection reaction proved to be principally dependent upon monovalent cation concentration and MutL concentration. Quantitative relationships of the optimal values of these parameters with length of the DNA test fragment were demonstrated, in support of the translocation model for the mechanism of action of these enzymes, rather than the molecular switch model. Thus, rapid, sequence-independent optimization was possible for each new genomic target region. Other factors potentially limiting the flexibility of mismatch scanning, such as positioning of dam recognition sites within the target fragment, have also been investigated. We developed several strategies, which can be easily adapted to automation, for limiting the analysis to intersample heteroduplexes. Thus, the principal barriers to the use of this methodology, which we have designated PCR candidate region mismatch scanning, in cost-effective, high-throughput settings have been removed.

Cardiac βARK1 inhibition prolongs survival and augments β blocker therapy in a mouse model of severe heart failure

Chronic human heart failure is characterized by abnormalities in β-adrenergic receptor (βAR) signaling, including increased levels of βAR kinase 1 (βARK1), which seems critical to the pathogenesis of the disease. To determine whether inhibition of βARK1 is sufficient to rescue a model of severe heart failure, we mated transgenic mice overexpressing a peptide inhibitor of βARK1 (βARKct) with transgenic mice overexpressing the sarcoplasmic reticulum Ca2+-binding protein, calsequestrin (CSQ). CSQ mice have a severe cardiomyopathy and markedly shortened survival (9 ± 1 weeks). In contrast, CSQ/βARKct mice exhibited a significant increase in mean survival age (15 ± 1 weeks; P < 0.0001) and showed less cardiac dilation, and cardiac function was significantly improved (CSQ vs. CSQ/βARKct, left ventricular end diastolic dimension 5.60 ± 0.17 mm vs. 4.19 ± 0.09 mm, P < 0.005; % fractional shortening, 15 ± 2 vs. 36 ± 2, P < 0.005). The enhancement of the survival rate in CSQ/βARKct mice was substantially potentiated by chronic treatment with the βAR antagonist metoprolol (CSQ/βARKct nontreated vs. CSQ/βARKct metoprolol treated, 15 ± 1 weeks vs. 25 ± 2 weeks, P < 0.0001). Thus, overexpression of the βARKct resulted in a marked prolongation in survival and improved cardiac function in a mouse model of severe cardiomyopathy that can be potentiated with β-blocker therapy. These data demonstrate a significant synergy between an established heart-failure treatment and the strategy of βARK1 inhibition.

Progress toward chemical accuracy in the computer simulation of condensed phase reactions.

We describe a procedure for the generation of chemically accurate computer-simulation models to study chemical reactions in the condensed phase. The process involves (i) the use of a coupled semiempirical quantum and classical molecular mechanics method to represent solutes and solvent, respectively; (ii) the optimization of semiempirical quantum mechanics (QM) parameters to produce a computationally efficient and chemically accurate QM model; (iii) the calibration of a quantum/classical microsolvation model using ab initio quantum theory; and (iv) the use of statistical mechanical principles and methods to simulate, on massively parallel computers, the thermodynamic properties of chemical reactions in aqueous solution. The utility of this process is demonstrated by the calculation of the enthalpy of reaction in vacuum and free energy change in aqueous solution for a proton transfer involving methanol, methoxide, imidazole, and imidazolium, which are functional groups involved with proton transfers in many biochemical systems. An optimized semiempirical QM model is produced, which results in the calculation of heats of formation of the above chemical species to within 1.0 kcal/mol (1 kcal = 4.18 kJ) of experimental values. The use of the calibrated QM and microsolvation QM/MM (molecular mechanics) models for the simulation of a proton transfer in aqueous solution gives a calculated free energy that is within 1.0 kcal/mol (12.2 calculated vs. 12.8 experimental) of a value estimated from experimental pKa values of the reacting species.

A structural basis for a phosphoramide mustard-induced DNA interstrand cross-link at 5'-d(GAC).

Phosphoramide mustard-induced DNA interstrand cross-links were studied both in vitro and by computer simulation. The local determinants for the formation of phosphoramide mustard-induced DNA interstrand cross-links were defined by using different pairs of synthetic oligonucleotide duplexes, each of which contained a single potentially cross-linkable site. Phosphoramide mustard was found to cross-link dG to dG at a 5'-d(GAC)-3'. The structural basis for the formation of this 1,3 cross-link was studied by molecular dynamics and quantum chemistry. Molecular dynamics indicated that the geometrical proximity of the binding sites also favored a 1,3 dG-to-dG linkage over a 1,2 dG-to-dG linkage in a 5'-d(GCC)-3' sequence. While the enthalpies of 1,2 and 1,3 mustard cross-linked DNA were found to be very close, a 1,3 structure was more flexible and may therefore be in a considerably higher entropic state.