Biomedical subjects

Alexander D Mackerell

Publications and source records attributed to Alexander D Mackerell.

12 recordsLinked to original sources

PubMed2006

Quantum mechanical analysis of 1,2-ethanediol conformational energetics and hydrogen bonding.

A proper understanding of the conformational energetics of 1,2-ethanediol (ethylene glycol) is important to the construction of molecular mechanics force fields for the treatment of carbohydrates since these biologically important molecules have a prevalence of vicinal hydroxyl groups. In the present study, quantum mechanical analysis of the 10 unique minimum-energy conformations of ethylene glycol is performed by using 10 model chemistries ranging from HF/6-311++G(d,p) up to a hybrid method that approximates CCSD(T)/cc-pVQZ. In addition, natural bond orbital (NBO) analysis of these conformations with deletion of pairings of CO bond/antibonding and lone pair/antibonding orbitals is used to investigate contributions from the "gauche" effect to ethylene glycol conformational energetics. MP2 with the "correlation consistent" basis sets and DFT/6-311++G(d,p) do the best job of matching the approximate CCSD(T)/cc-pVQZ energies while MP2/6-31G(d) and Hartree-Fock both fare poorly. NBO analysis shows the conformational energies to be independent of the deletion of matrix elements associated with (i) CO bonding and antibonding orbital interactions and (ii) lone pair and antibonding orbital interactions, whereas the energetic ordering correlates with geometric parameters consistent with internal hydrogen bonds. Thus, the present results suggest that standard molecular mechanics potential energy functional forms, which lack explicit terms to account for stereoelectronic effects, are appropriate for carbohydrates.

Ethylene Glycol↗

PubMed2006

Development of an empirical force field for silica. Application to the quartz-water interface.

Interactions of pulverized crystalline silica with biological systems, including the lungs, cause cell damage, inflammation, and apoptosis. To allow computational atomistic modeling of these pathogenic processes, including interactions between silica surfaces and biological molecules, new parameters for quartz, compatible with the CHARMM empirical force field were developed. Parameters were optimized to reproduce the experimental geometry of alpha-quartz, ab initio vibrational spectra, and interactions between model compounds and water. The newly developed force field was used to study interactions of water with two singular surfaces of alpha-quartz, (011) and (100). Properties monitored and analyzed include the variation of the density of water molecules in the plane perpendicular to the surface, disruption of the water H-bond network upon adsorption, and space-time correlations of water oxygen atoms in terms of Van Hove self-correlation functions. The vibrational density of states spectra of water in confined compartments were also computed and compared with experimental neutron-scattering results. Both the attenuation and shifting to higher frequencies of the hindered translational peaks upon confinement are clearly reproduced by the model. However, an upshift of librational peaks under the conditions of model confinement still remains underrepresented at the current empirical level.

Adsorption↗

PubMed2006

NMR imino proton exchange experiments on duplex DNA primarily monitor the opening of purine bases.

Molecular dynamics simulations were performed to investigate GC and AT base opening events in DNA. Calculated equilibrium constants between the base open (or flipped) and closed states were shown to be in good agreement with experimental data from NMR imino proton exchange experiments. Analysis of the computed results indicates that the equilbrium constants are dominated by the opening of the A and G bases in the AT and GC base pairs, respectively. Thus, the present results predict that NMR imino proton exchange experiments of base opening are primarily monitoring the opening of purine bases.

Adenine↗

PubMed2006

Mitogen activated protein (MAP) kinases: development of ATP and non-ATP dependent inhibitors.

Extracellular signals regulate most of the body's physiological functions through the MAP kinase signaling pathways. These MAP kinase signaling pathways are normally under tight regulation such that activation and inactivation occurs only when needed. However, aberrant regulation observed with naturally occurring mutations in specific signaling proteins often results in constitutive activation of the MAP kinases and is involved in several pathophysiological conditions, such as cancer, neurodegeneration, and inflammation. As such, much effort has been expended to develop inhibitory molecules of the MAP kinase signaling pathways. Several compounds have been identified that inhibit MAP kinase signaling by targeting receptors or other proteins upstream of the MAP kinases. The development of specific inhibitors of the MAP kinases themselves has been less successful and only a few compounds, which interfere with ATP binding, have been identified. A common problem with kinase inhibitors that compete with ATP binding is their lack of specificity. Thus, alternative approaches to inhibit MAP kinase function are being sought. The MAP kinase proteins contain docking domains that direct the interactions with a variety of substrate proteins. Using the 3-dimensional structure of MAP kinases and computer modeling, molecules that target specific docking domains and selectively disrupt substrate interactions are being developed. This non-ATP interfering approach may allow the selective inhibition of MAP kinase substrates involved in disease processes while preserving MAP kinase functions associated with normal cells.

Adenosine Triphosphate↗

PubMed2005

CH/pi interactions involving aromatic amino acids: refinement of the CHARMM tryptophan force field.

High-level ab initio calculations have been carried out to study weak CH/pi interactions and as a check of the CHARMM force field for aromatic amino acids. Comparisons with published data indicate that the MP2/cc-pVTZ level of theory is suitable for calculations of CH/pi interaction, including the T-shape benzene dimer. This level of theory was, therefore, applied to investigate CH/pi interactions between ethene or cis-2-butene and benzene in a variety of orientations. In addition, complexes between ethene and a series of model compounds (toluene, methylindole and p-cresol) representing the aromatic amino acids were studied motivated by the presence of CH/pi interactions in biological systems. Ab initio binding energies were compared to the binding energies obtained with the CHARMM22 force field. In the majority of orientations, CHARMM22 reproduces the preferred binding modes, with excellent agreement for the benzene dimer. Small discrepancies found in the calculations involving methylindole along with a survey of published thermodynamic data for the aromatic amino acids prompted additional optimization of the tryptophan force field. Partial atomic charges, Lennard-Jones parameters, and force constants were improved to obtain better intra- and intermolecular properties, with significant improvements obtained in the reproduction of experimental heats of sublimation for indole and free energies of aqueous solvation for methylindole.

Alkenes↗

PubMed2005

Identification of novel extracellular signal-regulated kinase docking domain inhibitors.

The extracellular signal regulated kinase (ERK1 and ERK2) signal transduction pathways play a critical role in cell proliferation. Hyperactivation of the ERK proteins either through increased expression of membrane-bound growth factor receptors or genetic mutations of upstream proteins is thought to be involved in the pathogenesis of many human cancers. Thus, targeted inhibition of ERK signaling is viewed as a potential approach to prevent cancer cell proliferation. Currently, no specific inhibitors of the ERK proteins exist. Moreover, most kinase inhibitors lack specificity because they target the ATP binding region, which is well conserved among the protein kinase families. Taking advantage of recently identified ERK docking domains, which are reported to facilitate substrate protein interactions, we have used computer-aided drug design (CADD) to identify novel small molecular weight ERK inhibitors. Following a CADD screen of over 800 000 molecules, 80 potential compounds were selected and tested for activity in biological assays. Several compounds inhibited ERK-specific phosphorylation of ribosomal S6 kinase-1 (Rsk-1) or the ternary complex factor Elk-1 (TCF/Elk-1), both of which are involved in promoting cell proliferation. Active compounds showed a dose-dependent reduction in the proliferation of several cancer cell lines as measured by colony survival assays. Direct binding between the active compounds and ERK2 was indicated by fluorescence quenching. These active compounds may serve as lead candidates for development of novel specific inhibitors of ERK-substrate interactions involved in cell proliferation.

Binding Sites↗

PubMed2005

Design of Inhibitors for S100B.

S100B interacts with the p53 protein in a calcium-dependent manner and down-regulates its function as a tumor suppressor. Therefore, inhibiting the S100B-p53 interaction represents a new approach for restoring functional wild-type p53 in cancers with elevated S100B such as found in malignant melanoma. A discussion of the biological rational for targeting S100B and a description of methodologies relevant to the discovery of compounds that inhibit S100B-p53 binding, including computational techniques, structural biology techniques, and cellular assays, is presented.

Biomarkers, Tumor↗

PubMed2004

Molecular-level organization of saturated and polyunsaturated fatty acids in a phosphatidylcholine bilayer containing cholesterol.

Cholesterol's preference for specific fatty acid chains is investigated at the atomic level in a 20 ns molecular dynamics computer simulation of a lipid bilayer membrane consisting of cholesterol and 1-stearoyl-2-docosahexaenoyl-sn-glycero-3-phosphocholine (SDPC) in a 1:3 ratio. These simulations reproduce experimental measurements suggesting that cholesterol prefers to be solvated by saturated acyl chains and has a low affinity for polyunsaturated fatty acids. Analyses of the simulation trajectory provide a detailed picture of both the transverse and lateral structures of the lipid bilayer membrane, along with a description of lipid and cholesterol dynamics at high temporal resolution. Comparison with a previous simulation of a pure phospholipid bilayer allows an atomic-level description of the changes in membrane structure and dynamics resulting from incorporation of cholesterol. The observed differential cholesterol interactions with saturated and polyunsaturated lipids provide a mechanism for the formation of laterally inhomogeneous membranes; thus, the simulation provides molecular-level insight into the formation of lipid rafts.

Cholesterol↗

PubMed2004

Molecular recognition of aldehydes by aldehyde dehydrogenase and mechanism of nucleophile activation.

Experimental structural data on the state of substrates bound to class 3 Aldehyde Dehydrogenases (ALDH3A1) is currently unknown. We have utilized molecular mechanics (MM) simulations, in conjunction with new force field parameters for aldehydes, to study the atomic details of benzaldehyde binding to ALDH3A1. Our results indicate that while the nucleophilic Cys243 must be in the neutral state to form what are commonly called near-attack conformers (NACs), these structures do not correlate with increased complexation energy calculated with the MM-Generalized Born Molecular Volume (GBMV) method. The negatively charged Cys243 (thiolate form) of ALDH3A1 also binds benzaldehyde in a stable conformation but in this complex the sulfur of Cys243 is oriented away from benzaldehyde yet yields the most favorable MM-GBMV complexation energy. The identity of the general base, Glu209 or Glu333, in ALDHs remains uncertain. The MM simulations reveal structural and possible functional roles for both Glu209 and Glu333. Structures from the MM simulations that would support either glutamate residue as the general base were further examined with Hybrid Quantum Mechanical (QM)/MM simulations. These simulations show that, with the PM3/OPLS potential, Glu209 must go through a step-wise mechanism to activate Cys243 through an intervening water molecule while Glu333 can go through a more favorable concerted mechanism for the same activation process.

Aldehyde Dehydrogenase↗

PubMed2004

Extending the treatment of backbone energetics in protein force fields: limitations of gas-phase quantum mechanics in reproducing protein conformational distributions in molecular dynamics simulations.

Computational studies of proteins based on empirical force fields represent a powerful tool to obtain structure-function relationships at an atomic level, and are central in current efforts to solve the protein folding problem. The results from studies applying these tools are, however, dependent on the quality of the force fields used. In particular, accurate treatment of the peptide backbone is crucial to achieve representative conformational distributions in simulation studies. To improve the treatment of the peptide backbone, quantum mechanical (QM) and molecular mechanical (MM) calculations were undertaken on the alanine, glycine, and proline dipeptides, and the results from these calculations were combined with molecular dynamics (MD) simulations of proteins in crystal and aqueous environments. QM potential energy maps of the alanine and glycine dipeptides at the LMP2/cc-pVxZ//MP2/6-31G* levels, where x = D, T, and Q, were determined, and are compared to available QM studies on these molecules. The LMP2/cc-pVQZ//MP2/6-31G* energy surfaces for all three dipeptides were then used to improve the MM treatment of the dipeptides. These improvements included additional parameter optimization via Monte Carlo simulated annealing and extension of the potential energy function to contain peptide backbone phi, psi dihedral crossterms or a phi, psi grid-based energy correction term. Simultaneously, MD simulations of up to seven proteins in their crystalline environments were used to validate the force field enhancements. Comparison with QM and crystallographic data showed that an additional optimization of the phi, psi dihedral parameters along with the grid-based energy correction were required to yield significant improvements over the CHARMM22 force field. However, systematic deviations in the treatment of phi and psi in the helical and sheet regions were evident. Accordingly, empirical adjustments were made to the grid-based energy correction for alanine and glycine to account for these systematic differences. These adjustments lead to greater deviations from QM data for the two dipeptides but also yielded improved agreement with experimental crystallographic data. These improvements enhance the quality of the CHARMM force field in treating proteins. This extension of the potential energy function is anticipated to facilitate improved treatment of biological macromolecules via MM approaches in general.

Algorithms↗

PubMed2004

CHARMM fluctuating charge force field for proteins: II protein/solvent properties from molecular dynamics simulations using a nonadditive electrostatic model.

A fluctuating charge (FQ) force field is applied to molecular dynamics simulations for six small proteins in explicit polarizable solvent represented by the TIP4P-FQ potential. The proteins include 1FSV, 1ENH, 1PGB, 1VII, 1H8K, and 1CRN, representing both helical and beta-sheet secondary structural elements. Constant pressure and temperature (NPT) molecular dynamics simulations are performed on time scales of several nanoseconds, the longest simulations yet reported using explicitly polarizable all-atom empirical potentials (for both solvent and protein) in the condensed phase. In terms of structure, the FQ force field allows deviations from native structure up to 2.5 A (with a range of 1.0 to 2.5 A). This is commensurate to the performance of the CHARMM22 nonpolarizable model and other currently existing polarizable models. Importantly, secondary structural elements maintain native structure in general to within 1 A (both helix and beta-strands), again in good agreement with the nonpolarizable case. In qualitative agreement with QM/MM ab initio dynamics on crambin (Liu et al. Proteins 2001, 44, 484), there is a sequence dependence of average condensed phase atomic charge for all proteins, a dependence one would anticipate considering the differing chemical environments around individual atoms; this is a subtle quantum mechanical feature captured in the FQ model but absent in current state-of-the-art nonpolarizable models. Furthermore, there is a mutual polarization of solvent and protein in the condensed phase. Solvent dipole moment distributions within the first and second solvation shells around the protein display a shift towards higher dipole moments (increases on the order of 0.2-0.3 Debye) relative to the bulk; protein polarization is manifested via the enhanced condensed phase charges of typical polar atoms such as backbone carbonyl oxygens, amide nitrogens, and amide hydrogens. Finally, to enlarge the sample set of proteins, gas-phase minimizations and 1 ps constant temperature simulations are performed on various-sized proteins to compare to earlier work by Kaminsky et al. (J Comp Chem 2002, 23, 1515). The present work establishes the feasibility of applying a fully polarizable force field for protein simulations and demonstrates the approach employed in extending the CHARMM force field to include these effects.

Computational Biology↗

PubMed2004

Empirical force fields for biological macromolecules: overview and issues.

Empirical force field-based studies of biological macromolecules are becoming a common tool for investigating their structure-activity relationships at an atomic level of detail. Such studies facilitate interpretation of experimental data and allow for information not readily accessible to experimental methods to be obtained. A large part of the success of empirical force field-based methods is the quality of the force fields combined with the algorithmic advances that allow for more accurate reproduction of experimental observables. Presented is an overview of the issues associated with the development and application of empirical force fields to biomolecular systems. This is followed by a summary of the force fields commonly applied to the different classes of biomolecules; proteins, nucleic acids, lipids, and carbohydrates. In addition, issues associated with computational studies on "heterogeneous" biomolecular systems and the transferability of force fields to a wide range of organic molecules of pharmacological interest are discussed.

Algorithms↗