PubMed Health⌕ Search

Biomedical subjects

Alexander D MacKerell

Publications and source records attributed to Alexander D MacKerell.

At least 19 recordsLinked to original sources
PubMed2007

CHARMM force field parameters for simulation of reactive intermediates in native and thio-substituted ribozymes.

Force field parameters specifically optimized for residues important in the study of RNA catalysis are derived from density-functional calculations, in a fashion consistent with the CHARMM27 all-atom empirical force field. Parameters are presented for residues that model reactive RNA intermediates and transition state analogs, thio-substituted phosphates and phosphoranes, and bound Mg(2+) and di-metal bridge complexes. Target data was generated via density-functional calculations at the B3LYP/6-311++G(3df,2p)// B3LYP/6-31++G(d,p) level. Partial atomic charges were initially derived from CHelpG electrostatic potential fitting and subsequently adjusted to be consistent with the CHARMM27 charges. Lennard-Jones parameters were determined to reproduce interaction energies with water molecules. Bond, angle, and torsion parameters were derived from the density-functional calculations and renormalized to maintain compatibility with the existing CHARMM27 parameters for standard residues. The extension of the CHARMM27 force field parameters for the nonstandard biological residues presented here will have considerable use in simulations of ribozymes, including the study of freeze-trapped catalytic intermediates, metal ion binding and occupation, and thio effects.

Catalysis↗
PubMed2006

Characterization of ATP-independent ERK inhibitors identified through in silico analysis of the active ERK2 structure.

The extracellular signal-regulated kinases (ERK1 and ERK2) are important mediators of cell proliferation. Constitutive activation of the ERK proteins plays a critical role in the proliferation of many human cancers. Taking advantage of recently identified substrate docking domains on ERK2, we have used computer-aided drug design (CADD) to identify novel low molecular weight compounds that interact with ERK2 in an ATP-independent manner and disrupt substrate-specific interactions. In the current study, a CADD screen of the 3D structure of active phosphorylated ERK2 protein was used to identify inhibitory compounds. We tested 13 compounds identified by the CADD screen in ERK-specific phosphorylation, cell proliferation, and binding assays. Of the 13 compounds tested, 4 compounds strongly inhibited ERK-mediated phosphorylation of ribosomal S6 kinase-1 (Rsk-1) and/or the transcription factor Elk-1 and inhibited the proliferation of HeLa cervical carcinoma cells with IC(50) values in the 2-10 microM range. These studies demonstrate that CADD can be used to identify lead compounds for development of novel non-ATP-dependent inhibitors selective for active ERK and its interactions with substrates involved in cancer cell proliferation.

Adenosine Triphosphate↗
PubMed2006

Conserved patterns in backbone torsional changes allow for single base flipping from duplex DNA with minimal distortion of the double helix.

Base flipping is a structural mechanism common to many DNA processing and repair enzymes. Changes in the local backbone torsions that occur during base flipping and the effect of environment on their behavior are of particular interest in understanding different base flipping mechanisms. In the present study, structures sampled during umbrella sampling molecular dynamics (MD) simulations of base flipping in aqueous and protein-bound environments, carried out with two different MD simulation strategies, are analyzed to find the most significant phosphodiester backbone distortions in the vicinity of the flipping base. Torsional sampling on the 5' side of the flipping base during flipping through the major groove shows similarities to the torsional sampling on the 3' side during flipping through the minor groove and vice versa. In differing environments, this behavior varies only marginally. These compensating torsional changes in the DNA backbone on 5' and 3' sides of the flipping base limit overall distortion of the DNA double helix during single base flipping. Rotameric intermediate states observed during base flipping are identified and postulated to be metastable states implicated in both large-scale structural changes and functional effects of chemical modifications in DNA.

DNA↗
PubMed2005

Importance of the CMAP correction to the CHARMM22 protein force field: dynamics of hen lysozyme.

The recently developed CMAP correction to the CHARMM22 force field (C22) is evaluated from 25 ns molecular dynamics simulations on hen lysozyme. Substantial deviations from experimental backbone root mean-square fluctuations and N-H NMR order parameters obtained in the C22 trajectories (especially in the loops) are eliminated by the CMAP correction. Thus, the C22/CMAP force field yields improved dynamical and structural properties of proteins in molecular dynamics simulations.

Animals↗
PubMed2005

Conformationally sampled pharmacophore for peptidic delta opioid ligands.

Opioids represent the frontline treatment for acute pain, despite their side effects, motivating efforts toward developing novel opioid analgesics. To facilitate these efforts, a novel modeling approach, the conformationally sampled pharmacophore (CSP), has been developed that increases the probability of including the receptor bound form in the model. This method, originally used for developing a nonpeptidic delta opioid efficacy pharmacophore, is extended to peptidic ligands using replica exchange molecular dynamics simulation for conformational sampling. The developed 2D CSP indicates that the spatial relationship of the basic nitrogen and the hydrophobic moiety in the delta opioid ligands differentiates activity. In addition, results indicate that both peptidic and nonpeptidic ligands have the same binding mode with the receptor. Thus, the CSP approach distinguishes both peptidic and nonpeptidic delta opioid agonists and antagonists and is anticipated to be of general utility for the development of pharmacophores for species with multiple rotatable bonds.

Analgesics, Opioid↗
PubMed2005

Polarizable empirical force field for alkanes based on the classical Drude oscillator model.

Recent extensions of potential energy functions used in empirical force field calculations have involved the inclusion of electronic polarizability. To properly include this extension into a potential energy function it is necessary to systematically and rigorously optimize the associated parameters based on model compounds for which extensive experimental data are available. In the present work, optimization of parameters for alkanes in a polarizable empirical force field based on a classical Drude oscillator is presented. Emphasis is placed on the development of parameters for CH3, CH2, and CH moieties that are directly transferable to long chain alkanes, as required for lipids and other biomolecules. It is shown that a variety of quantum mechanical and experimental target data are reproduced by the polarizable model. Notable is the proper treatment of the dielectric constant of pure alkanes by the polarizable force field, a property essential for the accurate treatment of, for example, hydrophobic solvation in lipid bilayers. The present alkane force field will act as the basis for the aliphatic moieties in an extensive empirical force field for biomolecules that includes the explicit treatment of electronic polarizability.

Alkanes↗
PubMed2005

Cooperative binding of DNA and CBFbeta to the Runt domain of the CBFalpha studied via MD simulations.

The Runt domain (RD) is the DNA-binding region of the Runx genes. A related protein, known as core binding factor beta (CBFbeta) also binds to the RD to enhance RD-DNA interaction by 6- to 10-fold. Here, we report results from molecular dynamics (MD) simulations of RD alone, as a dimer in complexes with DNA and CBFbeta and in a ternary complex with DNA and CBFbeta. Consistent with the experimental findings, in the presence of CBFbeta the estimated free energy of binding of RD to the DNA is more favorable, which is shown to be due to more favorable intermolecular interactions and desolvation contributions. Also contributing to the enhanced binding are favorable intramolecular interactions between the 'wing' residues (RD residues 139-145) and the 'wing1' residues (RD residues 104-116). The simulation studies also indicate that the RD-CBFbeta binding is more favorable in the presence of DNA due to a more favorable RD-CBFbeta interaction energy. In addition, it is predicted that long-range interactions involving ionic residues contribute to binding cooperativity. Results from the MD calculations are used to interpret a variety of experimental mutagenesis data. A novel role for RD Glu116 to the RD-CBFbeta interaction is predicted.

Amino Acids↗
PubMed2005

Position of coordination of the lithium ion determines the regioselectivity of demethylations of 3,4-dimethoxymorphinans with L-selectride.

[reaction: see text] L-Selectride is an efficient agent for the 3-O-demethylation of opioids and is known to cleave the least hindered methoxyl group in a molecule. The treatment of a 3,4-dimethoxymorphinan containing a 6-ketal with L-Selectride gave selective 4-O-demethylation, rather than cleavage of the less hindered 3-methoxyl. In contrast, a 3,4-dimethoxymorphinan lacking a 6-ketal gave selective 3-O-demethylation, suggesting that the regiochemistry of L-Selectride-mediated O-demethylation can be manipulated through altering the position of coordination of the lithium ion.

Catalysis↗
PubMed2005

Unexpected relative aqueous solubilities of a phosphotyrosine analogue and two phosphonate derivatives.

Phosphotyrosine (pTyr) is an essential component of biological signaling, often being a determinant of protein-protein interactions. Accordingly, a number of drug discovery efforts targeting signal transduction pathways have included phosphotyrosine and analogues as essential components of the lead compounds. Toward the goal of improved biological efficacy, the phosphonate and difluoro phosphonate analogues of pTyr have been employed in inhibitor design because of their stability to hydrolysis and enhanced binding affinity in certain cases. To quantitate the contribution of aqueous solubility of pTyr, phosphonomethyl phenylalanine (Pmp), and difluorophosphonomethyl phenylalanine (F(2)Pmp) to their relative binding affinities, free energy perturbation calculations were undertaken on the mimetics phenol phosphate (PP), benzyl phosphonate (BP), and difluorobenzyl phosphonate (F(2)BP), including development of empirical force field parameters compatible with the CHARMM all-atom force fields. Notably, it is shown that the most favorably solvated compound of the series is BP, followed by PP, with F(2)BP the least favorably solvated for both the mono- and dianionic forms of the compounds. The molecular origin of this ordering is shown to be due to changes in charge distribution, in the comparatively larger size of the fluorine atoms, as well as in differences of local solvation between PP and BP. The implications of the differences in aqueous solubility toward the relative binding potencies of pTyr-, Pmp-, and F(2)Pmp-containing peptide ligands are discussed. Our results indicate that one general principle explaining the efficacy of selective fluorination to enhance binding affinities may lie in the ability of fluorine atoms to increase the hydrophobicity of a ligand while maintaining its capability to form hydrogen bonds.

Computer Simulation↗
PubMed2005

An ab initio study on the torsional surface of alkanes and its effect on molecular simulations of alkanes and a DPPC bilayer.

Energies of 119 conformations of normal alkanes from butane to heptane were calculated at approximately the CCSD(T)/cc-pVQZ level. Energies of gauche (g) conformers relative to trans (t) decrease as chain length increases. In what is termed the "positive pentane effect", adjacent gauche conformers of the same sign are stabilized compared to nonadjacent conformers; e.g., for hexane the energies of tgt, tgg, and gtg are 0.600, 0.930, and 1.18 kcal/mol, respectively. Torsional terms in the CHARMM27 (C27) force field were fit to the calculated QM energies to yield a revised potential, C27r. Molecular dynamics simulations of normal alkanes (heptane, decane, tridecane, and pentadecane) with C27r yield higher populations of gauche states, increased transition rates, and improved agreement with experiment as compared to C27. In addition, C27r simulations of a hydrated DPPC lipid bilayer yield improved agreement with the experimental NMR deuterium order parameters for the aliphatic chain ends.

1,2-Dipalmitoylphosphatidylcholine↗
PubMed2005

Functionalization of the 6,14-bridge of the orvinols. 2. Preparation of 18- and 19-hydroxyl-substituted thevinols and their treatment with benzyl bromide.

The etheno bridge of a thevinone was treated with BH3 and H2O2 to give both the 18- and 19-hydroxyl- substituted thevinols. Selective benzylation of the primary 20-hydroxyl over the 19-hydroxyl was successful; however, benzylation of the 18-hydroxylated product led to a reaction at the more hindered alcohol. Thus, the 6,14-bridge of the Diels-Alder products of thebaine can be hydroxylated, which opens up these positions for further chemical manipulation.

Benzyl Compounds↗
PubMed2005

Conformational determinants of tandem GU mismatches in RNA: insights from molecular dynamics simulations and quantum mechanical calculations.

Structure and energetic properties of base pair mismatches in duplex RNA have been the focus of numerous investigations due to their role in many important biological functions. Such efforts have contributed to the development of models for secondary structure prediction of RNA, including the nearest-neighbor model. In RNA duplexes containing GU mismatches, 5'-GU-3' tandem mismatches have a different thermodynamic stability than 5'-UG-3' mismatches. In addition, 5'-GU-3' mismatches in some sequence contexts do not follow the nearest-neighbor model for stability. To characterize the underlying atomic forces that determine the structural and thermodynamic properties of GU tandem mismatches, molecular dynamics (MD) simulations were performed on a series of 5'-GU-3' and 5'-UG-3' duplexes in different sequence contexts. Overall, the MD-derived structural models agree well with experimental data, including local deviations in base step helicoidal parameters in the region of the GU mismatches and the model where duplex stability is associated with the pattern of GU hydrogen bonding. Further analysis of the simulations, validated by data from quantum mechanical calculations, suggests that the experimentally observed differences in thermodynamic stability are dominated by GG interstrand followed by GU intrastrand base stacking interactions that dictate the one versus two hydrogen bonding scenarios for the GU pairs. In addition, the inability of 5'-GU-3' mismatches in different sequence contexts to all fit into the nearest-neighbor model is indicated to be associated with interactions of the central four base pairs with the surrounding base pairs. The results emphasize the role of GG and GU stacking interactions on the structure and thermodynamics of GU mismatches in RNA.

Base Pair Mismatch↗
PubMed2005

Specificity in protein-DNA interactions: energetic recognition by the (cytosine-C5)-methyltransferase from HhaI.

Sequence-specific interactions between proteins and DNA are essential for a variety of biological functions. The (cytosine-C5)-methyltransferase from HhaI (M.HhaI) specifically modifies the second base in GCGC sequences, employing a base flipping mechanism to access the target base being chemically modified. The mechanism of sequence-specific recognition of M.HhaI is not evident based on crystallographic structures, leading to the suggestion that recognition is linked to the flipping event itself, a process that may be referred to as energetic recognition. Using computational methods, it is shown that the free energy barriers to flipping are significantly higher in non-cognate versus the cognate sequence, supporting the energetic recognition mechanism. Energetic recognition is imparted by two protein "selectivity filters" that function via a "web" of protein-DNA interactions in short-lived, high energy states present along the base flipping pathway. Other sequence-specific DNA binding proteins whose function involves significant distortion of DNA's conformation may use a similar recognition mechanism.

Crystallography, X-Ray↗
PubMed2004

Identification and characterization of small molecule inhibitors of the calcium-dependent S100B-p53 tumor suppressor interaction.

The binding of S100B to p53 down-regulates wild-type p53 tumor suppressor activity in cancer cells such as malignant melanoma, so a search for small molecules that bind S100B and prevent S100B-p53 complex formation was undertaken. Chemical databases were computationally searched for potential inhibitors of S100B, and 60 compounds were selected for testing on the basis of energy scoring, commercial availability, and chemical similarity clustering. Seven of these compounds bound to S100B as determined by steady state fluorescence spectroscopy (1.0 microM < or = K(D) < or = 120 microM) and five inhibited the growth of primary malignant melanoma cells (C8146A) at comparable concentrations (1.0 microM < or = IC(50) < or = 50 microM). Additionally, saturation transfer difference (STD) NMR experiments confirmed binding and qualitatively identified protons from the small molecule at the small molecule-S100B interface. Heteronuclear single quantum coherence (HSQC) NMR titrations indicate that these compounds interact with the p53 binding site on S100B. An NMR-docked model of one such inhibitor, pentamidine, bound to Ca(2+)-loaded S100B was calculated using intermolecular NOE data between S100B and the drug, and indicates that pentamidine binds into the p53 binding site on S100B defined by helices 3 and 4 and loop 2 (termed the hinge region).

Antineoplastic Agents↗
PubMed2004

The structure of aqueous guanidinium chloride solutions.

The combination of neutron diffraction with isotopic substitution (NDIS) experiments and molecular dynamics (MD) simulations to characterize the structuring in an aqueous solution of the denaturant guanidinium chloride is described. The simulations and experiments were carried out at a concentration of 3 m at room temperature, allowing for an examination of any propensity for ion association in a realistic solution environment. The simulations satisfactorily reproduced the principal features of the neutron scattering and indicate a bimodal hydration of the guanidinium ions, with the N-H groups making well-ordered hydrogen bonds in the molecular plane, but with the planar faces relatively deficient in interactions with water. The most striking feature of these solutions is the rich ion-ion ordering observed around the guanidinium ion in the simulations. The marked tendency of the guanidinium ions to stack parallel to their water-deficient surfaces indicates that the efficiency of this ion as a denaturant is due to its ability to simultaneously interact favorably with both water and hydrophobic side chains of proteins.

Cations↗
PubMed2004

Caught in the act: visualization of an intermediate in the DNA base-flipping pathway induced by HhaI methyltransferase.

Rotation of a DNA or RNA nucleotide out of the double helix and into a protein pocket ('base flipping') is a mechanistic feature common to some DNA/RNA-binding proteins. Here, we report the structure of HhaI methyltransferase in complex with DNA containing a south-constrained abasic carbocyclic sugar at the target site in the presence of the methyl donor byproduct AdoHcy. Unexpectedly, the locked south pseudosugar appears to be trapped in the middle of the flipping pathway via the DNA major groove, held in place primarily through Van der Waals contacts with a set of invariant amino acids. Molecular dynamics simulations indicate that the structural stabilization observed with the south-constrained pseudosugar will not occur with a north-constrained pseudosugar, which explains its lowered binding affinity. Moreover, comparison of structural transitions of the sugar and phosphodiester backbone observed during computational studies of base flipping in the M.HhaI-DNA-AdoHcy ternary complex indicate that the south-constrained pseudosugar induces a conformation on the phosphodiester backbone that corresponds to that of a discrete intermediate of the base-flipping pathway. As previous crystal structures of M.HhaI ternary complex with DNA displayed the flipped sugar moiety in the antipodal north conformation, we suggest that conversion of the sugar pucker from south to north beyond the middle of the pathway is an essential part of the mechanism through which flipping must proceed to reach its final destination. We also discuss the possibility of the south-constrained pseudosugar mimicking a transition state in the phosphodiester and sugar moieties that occurs during DNA base flipping in the presence of M.HhaI.

Binding Sites↗
PubMed2004

Atomistic view of base flipping in DNA.

Base flipping is essential for the enzyme-catalysed methylation of DNA. In our previous studies, the flipping of bases out of duplex DNA in DNA alone and when bound to the (cytosine-C5)-methyltransferase from HhaI (M.HhaI) were investigated via potential of mean force calculations. Insights into various experimental observations were obtained. In the present paper we present an overview of previous computational studies of base flipping along with new detailed structural and energetic analysis on atomic events that contribute to the free energy surfaces. The contributions from different intrinsic and environmental effects to the base-flipping process are explored, and experimental data derived from a variety of methods are reconciled. A detailed protein-facilitated base-flipping mechanism is proposed. Ground-state destabilization is achieved via disruption of the target base Watson-Crick interactions by substitution with favourable DNA-protein interactions. In addition, specific DNA-protein interactions and favourable solvation effects further promote target base flipping along the major groove through the protein matrix, and maximal interactions occur between the DNA and the protein upon reaching the fully flipped state. Other DNA binding proteins that involve distortion of DNA's conformation may use a similar mechanism to that by which M.HhaI facilitates base flipping.

Base Pairing↗