PubMed Health⌕ Search

Biomedical subjects

Alexander D MacKerell

Publications and source records attributed to Alexander D MacKerell.

35 records · Page 2Linked to original sources

Identification of non-phosphate-containing small molecular weight inhibitors of the tyrosine kinase p56 Lck SH2 domain via in silico screening against the pY + 3 binding site.

The protein p56 lymphoid T cell tyrosine kinase (Lck) is predominantly expressed in T lymphocytes where it plays a critical role in T-cell-mediated immune response. Lck participates in phosphotyrosine-dependent protein-protein interactions through its modular binding unit, the Src homology-2 (SH2) domain. Accordingly, virtual screening methods combined with experimental assays were used to identify small molecular weight nonpeptidic compounds that block Lck SH2 domain-dependent interactions. Virtual screening included scoring normalization procedures and postdocking structural clustering that is shown to facilitate the selection of active compounds. By targeting the well-defined hydrophobic binding pocket known to impart specificity on Lck-protein interactions (i.e., pY + 3 site), inhibitors of the Lck SH2 domain were discovered that omit the phosphotyrosine (pY) or related moieties. The 34 out of 196 computationally selected compounds were shown to inhibit Lck SH2 domain association with phosphorylated immunoreceptor tyrosine based activation motifs peptide. Twenty-four of the active compounds were further tested for their ability to modulate biological function. Thirteen of these compounds showed inhibitory activity in mixed lymphocyte culture assay. Fluorescence titration experiments on four of these active compounds further verified their binding to the SH2 domain. Because of their simple chemical structures, these small organic compounds have the potential to act as lead compounds for the development of novel immunosuppressant drugs.

Animals↗

3-chloropropanoic acid (UMB66): a ligand for the gamma-hydroxybutyric acid receptor lacking a 4-hydroxyl group.

Gamma-Hydroxybutyric acid (GHB) has gained in notoriety in recent years due to its association with sexual assaults. GHB is an endogenous ligand for GHB receptors, but its complete pharmacological mechanism of action in vivo remains unclear due to apparent GABAergic components. It has been proposed that the hydroxyl group in the 4-position acts as a hydrogen bond donor to the GHB receptor. Herein we show that 3-chloropropanoic acid possesses significant affinity for the GHB receptor, has no affinity for GABA receptors, and cannot undergo metabolism to GABAergic compounds. UMB66 is thus a selective agent for the study of GHB in vivo. These results, in combination with data from quantum mechanical calculations, suggest that the hydroxyl group of GHB actually acts as a hydrogen bond acceptor in contrast to the currently accepted model. This finding is anticipated to facilitate the rational design of novel agents with selectivity for GHB receptors that may be used to elucidate the mechanism of action of this common drug of abuse.

Animals↗

Influence of conformation on the EPR spectrum of 5,5-dimethyl-1-hydroperoxy-1-pyrrolidinyloxyl: a spin trapped adduct of superoxide.

Spin trapping, a technique used to characterize short-lived free radicals, consists of using a nitrone or nitroso compound to "trap" an unstable free radical as a long-lived aminoxyl that can be characterized by EPR spectroscopy. The resultant aminoxyl exhibits hyperfine splitting constants that are dependent on the spin trap and the free radical. Such is the case with 2,2-dimethyl-5-hydroxy-1-pyrrolidinyloxyl (DMPO-OH) and 2,2-dimethyl-5-hydroperoxy-1-pyrrodinyloxyl (DMPO-OOH) whose hyperfine splitting constants, A(N) = A(H) = 14.9 G and A(N) = 14.3 G, A(H)(beta) = 11.7 G, and A(H)(gamma) = 1.25 G, respectively, have been used to demonstrate the generation of HO(*) and O(2)(*)(-). However, to date, the source of the apparent A(H)(gamma) hyperfine splitting in DMPO-OOH is not known. We consider three possible explanations to account for the unique EPR spectrum of DMPO-OOH. The first is that the gamma-splitting arises from one of the hydrogen atoms at either carbon 3 or carbon 4 of DMPO-OOH. The second is that the gamma-splitting originates from the hydrogen atom of DMPO-OOH. The third is that the conformational properties of DMPO-R change upon going from DMPO-OH to DMPO-OOH. Experimental and theoretical chemical approaches as well as EPR spectral modeling were used to investigate which of these hypotheses may explain the asymmetric EPR spectrum of DMPO-OOH. From these studies it is shown that the 12-line EPR spectrum of DMPO-OOH results not from any proximal hydrogen, but from additional conformers of DMPO-OOH. Thus, the 1.25 G hyperfine splitting, which has been assigned as a gamma-splitting, is actually from two individual EPR spectra associated with different conformers of DMPO-OOH.

Electron Spin Resonance Spectroscopy↗

Improved treatment of the protein backbone in empirical force fields.

Empirical force field-based calculations of proteins, including protein-folding studies, have improved our understanding of the relationship of their structure to their biological function. However, limitations in the accuracy of empirical force fields in the treatment of the peptide backbone exist. Presented is a grid correction approach to improve the treatment of the peptide backbone phi/psi conformational energies. Inclusion of this correction with the CHARMM22 all-atom protein force field is shown to lead to significant improvement in the treatment of the conformational energies of both the peptide model compound, the alanine dipeptide, and of proteins in their crystal environment. The developed approach is suggested to lead to significant improvements in the accuracy of empirical force fields to treat peptides and proteins.

Dipeptides↗

Altered structural fluctuations in duplex RNA versus DNA: a conformational switch involving base pair opening.

DNA and RNA are known to have different structural properties. In the present study, molecular dynamics (MD) simulations on a series of RNA and DNA duplexes indicate differential structural flexibility for the two classes of oligonucleotides. In duplex RNA, multiple base pairs experienced local opening events into the major groove on the nanosecond time scale, while such events were not observed in the DNA simulations. Three factors are indicated to be responsible for the base opening events in RNA: solvent-base interactions, 2'OH(n)-O4'(n+1) intra-strand hydrogen bonding, and enhanced rigid body motion of RNA at the nucleoside level. Water molecules in the major groove of RNA contribute to initiation of base pair opening. Stabilization of the base pair open state is due to a 'conformational switch' comprised of 2'OH(n)-O4'(n+1) hydrogen bonding and a rigid body motion of the nucleoside moiety in RNA. This rigid body motion is associated with decreased flexibility of the glycosyl linkage and sugar moieties in A-form structures. The observed opening rates in RNA are consistent with the imino proton exchange experiments for AU base pairs, although not for GC base pairs, while structural and flexibility changes associated with the proposed conformational switch are consistent with survey data of RNA and DNA crystal structures. The possible relevance of base pair opening events in RNA to its many biological functions is discussed.

Base Pairing↗

2D conformationally sampled pharmacophore: a ligand-based pharmacophore to differentiate delta opioid agonists from antagonists.

Pharmacophores are widely used for rational drug design and include those based on receptor binding sites or on known ligands. To date, ligand-based pharmacophores have typically used one or a small number of conformers of known receptor ligands. However, this method does not take into account the inherent dynamic nature of molecules, which sample a wide range of conformations, any of which could be the bound form. In the present study, molecular dynamics (MD) simulations were used as a means to sample the conformational space of ligands to include all accessible conformers at room temperature in the development of a pharmacophore. On the basis of these conformers, probability distributions of selected distances and angles in a series of delta specific opioid ligands were obtained and correlated with agonist versus antagonist activities. Individually, the distributions did not allow for unique agonist and antagonist pharmacophores to be identified. However, by extending the conformational analysis to two dimensions, a 2D conformationally sampled pharmacophore (CSP) for distinguishing delta receptor agonists and antagonists was developed. Application of this model to the compound DPI2505 suggests that it may have agonist activity. It is anticipated that the CSP method, which does not require alignment of compounds during pharmacophore development, will be a useful tool for obtaining structure-function relationships of ligands particularly in systems where the receptor 3D structure is not known.

Computer Simulation↗

Rearrangement of 5-trimethylsilylthebaine on treatment with L-selectride: an efficient synthesis of (+)-bractazonine.

Treatment of 5-trimethylsilylthebaine with L-Selectride gave rise to a rearrangement to 10-trimethylsilylbractazonine through migration of the phenyl group, whereas treatment of thebaine with strong Lewis acids is known to lead to a similar rearrangement through migration of the alkyl bridge to give, after reduction, (+)-neodihydrothebaine. It is suggested that the rearrangement of the alkyl group of thebaine is favored due to the formation of a tertiary benzylic cation. However, for 5-trimethylsilylthebaine, the lithium ion of L-Selectride acts as the Lewis acid and the beta-silyl effect dominates in the stabilization of any positive charge. This rearrangement provides a clear example of the greater relative migratory aptitude of phenyl groups over alkyl groups, and provides an efficient synthesis of (+)-bractazonine from thebaine.

Alkaloids↗

Protein-facilitated base flipping in DNA by cytosine-5-methyltransferase.

DNA methylation, various DNA repair mechanisms, and possibly early events in the opening of DNA as required for transcription and replication are initiated by flipping of a DNA base out of the DNA double helix. The energetics and structural mechanism of base flipping in the presence of the DNA-processing enzyme, cytosine 5-methyltransferase from HhaI (M.HhaI), were obtained through molecular dynamics based upon free-energy calculations. Free-energy profiles for base flipping show that, when in the closed conformation, M.HhaI lowers the free-energy barrier to flipping by 17 kcalmol and stabilizes the fully flipped state. Flipping is shown to occur via the major groove of the DNA. Structural analysis indicates that flipping is facilitated by destabilization of the DNA double-helical structure and substitution of DNA base-pairing and base-stacking interactions with DNA-protein interactions. The fully flipped state is stabilized by DNA-protein interactions that are enhanced upon binding of coenzyme. This study represents an atomic detail description of the mechanism by which a protein facilitates specific structural distortion in DNA.

Base Composition↗

Free energy and structural pathways of base flipping in a DNA GCGC containing sequence.

Structural distortions of DNA are essential for its biological function due to the genetic information of DNA not being physically accessible in the duplex state. Base flipping is one of the simplest structural distortions of DNA and may represent an initial event in strand separation required to access the genetic code. Flipping is also utilized by DNA-modifying and repair enzymes to access specific bases. It is typically thought that base flipping (or base-pair opening) occurs via the major groove whereas minor groove flipping is only possible when mediated by DNA-binding proteins. Here, umbrella sampling with a novel center-of-mass pseudodihedral reaction coordinate was used to calculate the individual potentials of mean force (PMF) for flipping of the Watson-Crick (WC) paired C and G bases in the CCATGCGCTGAC DNA dodecamer. The novel reaction coordinate allowed explicit investigation of the complete flipping process via both the minor and major groove pathways. The minor and major groove barriers to flipping are similar for C base flipping while the major groove barrier is slightly lower for G base flipping. Minor groove flipping requires distortion of the WC partner while the flipping base pulls away from its partner during major groove flipping. The flipped states are represented by relatively flat free energy surfaces, with a small, local minimum observed for the flipped G base. Conserved patterns of phosphodiester backbone dihedral distortions during flipping indicate their essential role in the flipping process. During flipping, the target base tracks along the respective grooves, leading to hydrogen-bonding interactions with neighboring base-pairs. Such hydrogen-bonding interactions with the neighboring sequence suggest a novel mechanism of sequence dependence in DNA dynamics.

Base Pairing↗

Combined ab initio/empirical approach for optimization of Lennard-Jones parameters for polar-neutral compounds.

The study of small functionalized organic molecules in aqueous solution is a useful step toward gaining a basic understanding of the behavior of biomolecular systems in their native aqueous environment. Interest in studying amines and fluorine-substituted compounds has risen from their intrinsic physicochemical properties and their prevalence in biological and pharmaceutical compounds. In the present study, a previously developed approach which optimizes Lennard-Jones (LJ) parameters via the use of rare gas atoms combined with the reproduction of experimental condensed phase properties was extended to polar-neutral compounds. Compounds studied included four amines (ammonia, methylamine, dimethylamine, and trimethylamine) and three fluoroethanes (1-fluoroethane, 1,1-difluoroethane, and 1,1,1-trifluoroethane). The resulting force field yielded heats of vaporization and molecular volumes in excellent agreement with the experiment, with average differences less than 1%. The current amine CHARMM parameters successfully reproduced experimental aqueous solvation data where methylamine is more hydrophilic than ammonia, with hydrophobicity increasing with additional methylation on the nitrogen. For both the amines and fluoroethanes the parabolic relationship of the extent of methylation or fluorination, respectively, to the heats of vaporization were reproduced by the new parameters. The present results are also discussed with respect to the impact of parameterization approach to molecular details obtained from computer simulations and to the unique biological properties of fluorine in pharmaceutical compounds.

Amines↗

Polyunsaturated fatty acids in lipid bilayers: intrinsic and environmental contributions to their unique physical properties.

Polyunsaturated lipids are an essential component of biological membranes, influencing order and dynamics of lipids, protein-lipid interaction, and membrane transport properties. To gain an atomic level picture of the impact of polyunsaturation on membrane properties, quantum mechanical (QM) and empirical force field based calculations have been undertaken. The QM calculations of the torsional energy surface for rotation about vinyl-methylene bonds reveal low barriers to rotation, indicating an intrinsic propensity toward flexibility. Based on QM and experimental data, empirical force field parameters were developed for polyunsaturated lipids and applied in a 16 ns molecular dynamics (MD) simulation of a 1-stearoyl-2-docosahexaenoyl-sn-glyerco-3-phosphocholine (SDPC) lipid bilayer. The simulation results are in good agreement with experimental data, suggesting an unusually high degree of conformational flexibility of polyunsaturated hydrocarbon chains in membranes. The detailed analysis of chain conformation and dynamics by simulations is aiding the interpretation of experimental data and is useful for understanding the unique role of polyunsaturated lipids in biological membranes. The complete force field is included as Supporting Information and is available from http://www.pharmacy.umaryland.edu/faculty/amackere/research.html.

Computer Simulation↗

Intrinsic conformational energetics associated with the glycosyl torsion in DNA: a quantum mechanical study.

The glycosyl torsion (chi) in nucleic acids has long been recognized to be a major determinant of their conformational properties. chi torsional energetics were systematically mapped in deoxyribonucleosides using high-level quantum mechanical methods, for north and south sugar puckers and with gamma in the g(+) and trans conformations. In all cases, the syn conformation is found higher in energy than the anti. When gamma is changed from g(+) to trans, the anti orientation of the base is strongly destabilized, and the energy difference and barrier between anti and syn are significantly decreased. The barrier between anti and syn in deoxyribonucleosides is found to be less than 10 kcal/mol and tends to be lower with purines than with pyrimidines. With gamma = g(+)/chi = anti, a south sugar yields a significantly broader energy well than a north sugar with no energy barrier between chi values typical of A or B DNA. Contrary to the prevailing view, the syn orientation is not more stable with south puckers than with north puckers. The syn conformation is significantly more energetically accessible with guanine than with adenine in 5-nucleotides but not in nucleosides. Analysis of nucleic acid crystal structures shows that gamma = trans/chi = anti is a minor but not negligible conformation. Overall, chi appears to be a very malleable structural parameter with the experimental chi distributions reflecting, to a large extent, the associated intrinsic torsional energetics.

Adenine↗

Structure-based inhibitor design targeting HIV-1 integrase.

HIV integrase (IN) is a viral-encoded protein that catalyzes the breaking and joining reactions that mediate integration of viral DNA into the host genome. Therefore, IN offers a unique target for the development of novel anti-HIV and anti-AIDS therapeutics. To take advantage of this potential, drug discovery efforts via structure-based design approaches have been undertaken. Presented is a review of computer-aided drug design efforts targeting HIV IN. Included is an overview of the life-cycle of HIV, with emphasis on the mechanism of action of IN, biological assays for measuring IN activity and identifying IN inhibitors, and the appropriate cell-based assays required for determining the antiviral activity of IN inhibitors. This is followed by a review of the available three-dimensional structures of HIV IN. Structure-based drug design efforts are then critiqued, including both ligand-based (e.g. pharmacophore) and target-based (e.g. docking) methods. Results from recent computational chemistry studies of IN are also discussed.

Catalysis↗

DNA bending induced by carbocyclic sugar analogs constrained to the north conformation.

DNA bending caused by introduction of carbocyclic sugars constrained to the north conformation was studied, using explicit solvent molecular dynamic (MD) simulations. The native Drew-Dickerson (DD) dodecamer and its three modifications containing north carbocyclic sugars in the 7th (T7*), 8th (T8*) or both 7th and 8th (T7T8*) nucleotide positions were examined. Introduction of the carbocyclic sugar results in A-form conformations for the alpha, beta, chi, zeta, and sugar pucker backbone parameters in the modified nucleotides. Increased steric repulsion between the sugar and its parent base in the modified oligonucleotides impacts the roll and cup dinucleotide step parameters, increasing the bending of the oligomer axis. Increased buckling of the substituted nucleotides disrupts the usual stabilizing base stacking interactions. The level of overall bending depends on the number and position of carbocyclic sugars introduced in the DNA sequence. Single sugar substitutions are unable to induce substantial bending due to the neighboring unmodified nucleotides counterbalancing the distortion. Significant bending can, however, be induced by two consecutive north sugars (T7T8*), which is in agreement with experimental results. The modified oligomers populate a wide range of bend angles, indicating that they maintain flexibility in the bent state. The present results suggest that insertion of carbocyclic sugars into DNA or RNA duplexes can be used to engineer bending of the duplexes without impacting the electrostatic or chemical properties of the phosphodiester backbone, thereby serving as excellent tools for experimental elucidation of nucleic acid structure-function relationships.

Base Pairing↗

Consideration of molecular weight during compound selection in virtual target-based database screening.

Virtual database screening allows for millions of chemical compounds to be computationally selected based on structural complimentary to known inhibitors or to a target binding site on a biological macromolecule. Compound selection in virtual database screening when targeting a biological macromolecule is typically based on the interaction energy between the chemical compound and the target macromolecule. In the present study it is shown that this approach is biased toward the selection of high molecular weight compounds due to the contribution of the compound size to the energy score. To account for molecular weight during energy based screening, we propose normalization strategies based on the total number of heavy atoms in the chemical compounds being screened. This approach is computationally efficient and produces molecular weight distributions of selected compounds that can be selected to be (1) lower than that of the original database used in the virtual screening, which may be desirable for selection of leadlike compounds or (2) similar to that of the original database, which may be desirable for the selection of drug-like compounds. By eliminating the bias in target-based database screening toward higher molecular weight compounds it is anticipated that the proposed procedure will enhance the success rate of computer-aided drug design.

Computer-Aided Design↗

Lead validation and SAR development via chemical similarity searching; application to compounds targeting the pY+3 site of the SH2 domain of p56lck.

Compound selection based on chemical similarity has been used to validate active "parent" compounds identified via database searching as viable lead compounds and to obtain initial structure-activity relationships for those leads. Twelve parent compounds that have inhibitory activity against the SH2 domain of the p56 T-cell tyrosine kinase (Lck) are the focus of this study. Lck is involved in the T-cell mediated immune response, and inhibitors of Lck protein-protein interactions could potentially be used to develop novel immunosuppressants. Similarity searches for each parent compound were performed using 2D structural fingerprints on a database containing 1,300,000 commercially available compounds. The inhibitory activity of the selected compounds was assessed using enzyme immunoassay (EIA). In general, the most active parent compounds yield the most high activity similar compounds; however, in two cases low activity parent compounds (i.e. inhibitory activity < 25% at 100 microM) yielded multiple similar compounds with activities > 60%. Such compounds may, therefore, be considered as viable lead compounds for optimization. Structure-activity relationships were explored by examining both ligand structures and their computed bound conformations to the protein. Functional groups common to the active compounds as well as key amino acid residues that form hydrogen bonds with the active compounds were identified. This information will act as the basis for the rational optimization of the lead compounds.

Computer Simulation↗

Chemical substituent effect on pyridine permeability and mechanistic insight from computational molecular descriptors.

The objective was (1) to evaluate the chemical substituent effect on Caco-2 permeability, using a congeneric series of pyridines, and (2) compare molecular descriptors from a computational chemistry approach against molecular descriptors from the Hansch approach for their abilities to explain the chemical substituent effect on pyridine permeability. The passive permeability of parent pyridine and 14 monosubstituted pyridines were measured across Caco-2 monolayers. Computational chemistry analysis was used to obtain the following molecular descriptions: solvation free energies, solvent accessible surface area, polar surface area, and cavitation energy. Results indicate that the parent pyridine was highly permeable and that chemical substitution was able to reduce pyridine permeability almost 20-fold. The substituent effect on permeability provided the following rank order: 3-COO- < 4-NH2 < 3-CONH2 < 3-Cl < 3-CHO < 3-OH < 3-CH2OH < 3-C6H5 < 3-NH2 < 3-CH2C6H5 < 3-C2H5 < 3-H < 3-CH3 < 3-F < 4-C6H5. This substituent effect was better explained via molecule descriptors from the computational chemistry approach than explained by classic descriptors from Hansch. Computational descriptors indicate that aqueous desolvation, but not membrane partitioning per se, dictated substituent effect on permeability.

Caco-2 Cells↗