Biomedical subjects

Michael Feig

Publications and source records attributed to Michael Feig.

At least 19 recordsLinked to original sources

PubMed2006

Implicit solvent simulations of DNA and DNA-protein complexes: agreement with explicit solvent vs experiment.

Molecular dynamics simulations of biomolecules with implicit solvent reduce the computational cost and complexity of such simulations so that longer time scales and larger system sizes can be reached. While implicit solvent simulations of proteins have become well established, the success of implicit solvent in the simulation of nucleic acids has not been fully established to date. Results obtained in this study demonstrate that stable and efficient simulations of DNA and a protein-DNA complex can be achieved with an implicit solvent model based on continuum dielectric electrostatics. Differences in conformational sampling of DNA with two sets of atomic radii that are used to define the dielectric interface between the solute and the continuum dielectric model of the solvent are investigated. Results suggest that depending on the choice of atomic radii agreement is either closer to experimental data or to explicit solvent simulations. Furthermore, partial conformational transitions toward A-DNA conformations when salt is added within the implicit solvent framework are observed.

Computer Simulation↗

PubMed2006

The unorthodox SNAP50 zinc finger domain contributes to cooperative promoter recognition by human SNAPC.

Human small nuclear RNA gene transcription by RNA polymerases II and III depends upon promoter recognition by the SNAPC general transcription factor. DNA binding by SNAPC involves direct DNA contacts by the SNAP190 subunit in cooperation with SNAP50 and SNAP43. The data presented herein shows that SNAP50 plays an important role in DNA binding by SNAPC through its zinc finger domain. The SNAP50 zinc finger domain contains 15 cysteine and histidine residues configured in two potential zinc coordination arrangements. Individual alanine substitution of each cysteine and histidine residue demonstrated that eight sites are important for DNA binding by SNAPC. However, metal binding studies revealed that SNAPC contains a single zinc atom indicating that only one coordination site functions as a zinc finger. Of the eight residues critical for DNA binding, four cysteine residues were also essential for both U1 and U6 transcription by RNA polymerase II and III, respectively. Surprisingly, the remaining four residues, although critical for U1 transcription could support partial U6 transcription. DNA binding studies showed that defects in DNA binding by SNAPC alone could be suppressed through cooperative DNA binding with another member of the RNA polymerase III general transcription machinery, TFIIIB. These results suggest that these eight cysteine and histidine residues perform different functions during DNA binding with those residues involved in zinc coordination likely performing a dominant role in domain stabilization and the others involved in DNA binding. These data further define the unorthodox SNAP50 zinc finger region as an evolutionarily conserved DNA binding domain.

Amino Acid Sequence↗

PubMed2006

Balancing an accurate representation of the molecular surface in generalized born formalisms with integrator stability in molecular dynamics simulations.

Different integrator time steps in NVT and NVE simulations of protein and nucleic acid systems are tested with the GBMV (Generalized Born using Molecular Volume) and GBSW (Generalized Born with simple SWitching) methods. The simulation stability and energy conservation is investigated in relation to the agreement with the Poisson theory. It is found that very close agreement between generalized Born methods and the Poisson theory based on the commonly used sharp molecular surface definition results in energy drift and simulation artifacts in molecular dynamics simulation protocols with standard 2-fs time steps. New parameters are proposed for the GBMV method, which maintains very good agreement with the Poisson theory while providing energy conservation and stable simulations at time steps of 1 to 1.5 fs.

Computer Simulation↗

PubMed2006

A correlation-based method for the enhancement of scoring functions on funnel-shaped energy landscapes.

A correlation-based approach is introduced for enhancing the ability of structure-scoring methods to identify and distinguish native-like conformations. The proposed method relies on a funnel-shaped scoring function that decreases steadily toward the native state. It takes advantage of the idea that the structure from a given ensemble that is closest to the native basin leads to the highest correlation coefficient between a given score and distance to that structure as an approximation of the native state for the entire ensemble. The method is applied successfully to a number of different test cases that demonstrate substantial improvements in the correlation of the score with the distance from the true native state but also result in the selection of more native-like structures compared to the original score.

Algorithms↗

PubMed2006

Molecular dynamics simulations of large integral membrane proteins with an implicit membrane model.

The heterogeneous dielectric generalized Born (HDGB) methodology is an the extension of the GBMV model for the simulation of integral membrane proteins with an implicit membrane environment. Three large integral membrane proteins, the bacteriorhodopsin monomer and trimer and the BtuCD protein, were simulated with the HDGB model in order to evaluate how well thermodynamic and dynamic properties are reproduced. Effects of the truncation of electrostatic interactions were examined. For all proteins, the HDGB model was able to generate stable trajectories that remained close to the starting experimental structures, in excellent agreement with explicit membrane simulations. Dynamic properties evaluated through a comparison of B-factors are also in good agreement with experiment and explicit membrane simulations. However, overall flexibility was slightly underestimated with the HDGB model unless a very large electrostatic cutoff is employed. Results with the HDGB model are further compared with equivalent simulations in implicit aqueous solvent, demonstrating that the membrane environment leads to more realistic simulations.

Bacteriorhodopsins↗

PubMed2006

Structural basis for the enantiospecificities of R- and S-specific phenoxypropionate/alpha-ketoglutarate dioxygenases.

(R)- and (S)-dichlorprop/alpha-ketoglutarate dioxygenases (RdpA and SdpA) catalyze the oxidative cleavage of 2-(2,4-dichlorophenoxy)propanoic acid (dichlorprop) and 2-(4-chloro-2-methyl-phenoxy)propanoic acid (mecoprop) to form pyruvate plus the corresponding phenol concurrent with the conversion of alpha-ketoglutarate (alphaKG) to succinate plus CO2. RdpA and SdpA are strictly enantiospecific, converting only the (R) or the (S) enantiomer, respectively. Homology models were generated for both enzymes on the basis of the structure of the related enzyme TauD (PDB code 1OS7). Docking was used to predict the orientation of the appropriate mecoprop enantiomer in each protein, and the predictions were tested by characterizing the activities of site-directed variants of the enzymes. Mutant proteins that changed at residues predicted to interact with (R)- or (S)-mecoprop exhibited significantly reduced activity, often accompanied by increased Km values, consistent with roles for these residues in substrate binding. Four of the designed SdpA variants were (slightly) active with (R)-mecoprop. The results of the kinetic investigations are consistent with the identification of key interactions in the structural models and demonstrate that enantiospecificity is coordinated by the interactions of a number of residues in RdpA and SdpA. Most significantly, residues Phe171 in RdpA and Glu69 in SdpA apparently act by hindering the binding of the wrong enantiomer more than the correct one, as judged by the observed decreases in Km when these side chains are replaced by Ala.

2-Methyl-4-chlorophenoxyacetic Acid↗

PubMed2006

Synthetic curcuminoids modulate the arachidonic acid metabolism of human platelet 12-lipoxygenase and reduce sprout formation of human endothelial cells.

Platelet 12-lipoxygenase (P-12-LOX) is overexpressed in different types of cancers, including prostate cancer, and the level of expression is correlated with the grade of this cancer. Arachidonic acid is metabolized by 12-LOX to 12(S)-hydroxyeicosatetraenoic acid [12(S)-HETE], and this biologically active metabolite is involved in prostate cancer progression by modulating cell proliferation in multiple cancer-related pathways inducing angiogenesis and metastasis. Thus, inhibition of P-12-LOX can reduce these two processes. Several lipoxygenase inhibitors are known, including plant and mammalian lipoxygenases, but only a few of them are known inhibitors of P-12-LOX. Curcumin is one of these lipoxygenase inhibitors. Using a homology model of the three-dimensional structure of human P-12-LOX, we did computational docking of synthetic curcuminoids (curcumin derivatives) to identify inhibitors superior to curcumin. Docking of the known inhibitors curcumin and NDGA to P-12-LOX was used to optimize the docking protocol for the system in study. Over 75% of the compounds of interest were successfully docked into the active site of P-12-LOX, many of them sharing similar binding modes. Curcuminoids that did not dock into the active site did not inhibit P-12-LOX. From a set of the curcuminoids that were successfully docked and selected for testing, two were found to inhibit human lipoxygenase better than curcumin. False-positive curcuminoids showed high LogP (theoretical) values, indicating poor water solubility, a possible reason for lack of inhibitory activity or/and nonrealistic binding. Additionally, the curcuminoids inhibiting P-12-LOX were tested for their ability to reduce sprout formation of endothelial cells (in vitro model of angiogenesis). We found that only curcuminoids inhibiting human P-12-LOX and the known inhibitor NDGA reduced sprout formation. Only limited inhibition of sprout formation at approximately IC(50) concentrations has been seen. At IC(50), a substantial amount of 12-HETE can be produced by lipoxygenase, providing a stimulus for angiogenic sprouting of endothelial cells. Increasing the concentration of lipoxygenase inhibitors above IC(50), thus decreasing the concentration of 12(S)-HETE produced, greatly reduced sprout formation for all inhibitors tested. This universal event for all tested lipoxygenase inhibitors suggests that the inhibition of sprout formation was most likely due to the inhibition of human P-12-LOX but not other cancer-related pathways.

Amino Acid Sequence↗

PubMed2005

Dynamic error correction and regulation of downstream bubble opening by human RNA polymerase II.

The nucleotide triphosphate (NTP)-driven translocation hypothesis posits that NTP substrates bind to templated DNA sites prior to translocation into the active site. By using millisecond phase kinetics, we demonstrate this prediction in three different ways. First, we show that, in the presence of the translocation blocker alpha-amanitin, NTPs (but not deoxynucleotide triphosphate [dNTPs]) templated at downstream sites (i + 2 and i + 3) dislodge an active site (i + 1) NTP, which was otherwise fated to complete bond synthesis. Second, we show that NTPs templated at i + 2 and/or i + 3 downstream sites suppress misincorporation errors. Third, we show that NTPs templated at downstream sites stabilize the posttranslocated elongation complex at a stall position. Therefore, at least two NTP substrates pair to DNA templated sites downstream of the active site. These results demonstrate the mechanisms of NTP loading and transcriptional efficiency and fidelity for human RNA polymerase II and indicate regulation of downstream bubble opening by NTPs.

Catalytic Domain↗

PubMed2005

A generalized Born formalism for heterogeneous dielectric environments: application to the implicit modeling of biological membranes.

Reliable computer simulations of complex biological environments such as integral membrane proteins with explicit water and lipid molecules remain a challenging task. We propose a modification of the standard generalized Born theory of homogeneous solvent for modeling the heterogeneous dielectric environments such as lipid/water interfaces. Our model allows the representation of biological membranes in the form of multiple layered dielectric regions with dielectric constants that are different from the solute cavity. The proposed new formalism is shown to predict the electrostatic component of solvation free energy with a relative error of 0.17% compared to exact finite-difference solutions of the Poisson equation for a transmembrane helix test system. Molecular dynamics simulations of melittin and bacteriorhodopsin are carried out and performed over 10 ns and 7 ns of simulation time, respectively. The center of melittin along the membrane normal in these stable simulations is in excellent agreement with the relevant experimental data. Simulations of bacteriorhodopsin started from the experimental structure remained stable and in close agreement with experiment. We also examined the free energy profiles of water and amino acid side chain analogs upon membrane insertion. The results with our implicit membrane model agree well with the experimental transfer free energy data from cyclohexane to water as well as explicit solvent simulations of water and selected side chain analogs.

Bacteriorhodopsins↗

PubMed2005

The requirement for mechanical coupling between head and S2 domains in smooth muscle myosin ATPase regulation and its implications for dimeric motor function.

A combination of experimental structural data, homology modelling and elastic network normal mode analysis is used to explore how coupled motions between the two myosin heads and the dimerization domain (S2) in smooth muscle myosin II determine the domain movements required to achieve the inhibited state of this ATP-dependent molecular motor. These physical models rationalize the empirical requirement for at least two heptads of non-coiled alpha-helix at the junction between the myosin heads and S2, and the dependence of regulation on S2 length. The results correlate well with biochemical data regarding altered conformational-dependent solubility and stability. Structural models of the conformational transition between putative active states and the inhibited state show that torsional flexibility of the S2 alpha-helices is a key mechanical requirement for myosin II regulation. These torsional motions of the myosin heads about their coiled coil alpha-helices affect the S2 domain structure, which reciprocally affects the motions of the myosin heads. This inter-relationship may explain a large body of data on function of molecular motors that form dimers through a coiled-coil domain.

Animals↗

PubMed2005

NTP-driven translocation and regulation of downstream template opening by multi-subunit RNA polymerases.

Multi-subunit RNA polymerases bind nucleotide triphosphate (NTP) substrates in the pretranslocated state and carry the dNMP-NTP base pair into the active site for phosphoryl transfer. NTP-driven translocation requires that NTP substrates enter the main-enzyme channel before loading into the active site. Based on this model, a new view of fidelity and efficiency of RNA synthesis is proposed. The model predicts that, during processive elongation, NTP-driven translocation is coupled to a protein conformational change that allows pyrophosphate release: coupling the end of one bond-addition cycle to substrate loading and translocation for the next. We present a detailed model of the RNA polymerase II elongation complex based on 2 low-affinity NTP binding sites located in the main-enzyme channel. This model posits that NTP substrates, elongation factors, and the conserved Rpb2 subunit fork loop 2 cooperate to regulate opening of the downstream transcription bubble.

Humans↗

PubMed2004

Performance comparison of generalized born and Poisson methods in the calculation of electrostatic solvation energies for protein structures.

This study compares generalized Born (GB) and Poisson (PB) methods for calculating electrostatic solvation energies of proteins. A large set of GB and PB implementations from our own laboratories as well as others is applied to a series of protein structure test sets for evaluating the performance of these methods. The test sets cover a significant range of native protein structures of varying size, fold topology, and amino acid composition as well as nonnative extended and misfolded structures that may be found during structure prediction and folding/unfolding studies. We find that the methods tested here span a wide range from highly accurate and computationally demanding PB-based methods to somewhat less accurate but more affordable GB-based approaches and a few fast, approximate PB solvers. Compared with PB solvation energies, the latest, most accurate GB implementations were found to achieve errors of 1% for relative solvation energies between different proteins and 0.4% between different conformations of the same protein. This compares to accurate PB solvers that produce results with deviations of less than 0.25% between each other for both native and nonnative structures. The performance of the best GB methods is discussed in more detail for the application for force field-based minimizations or molecular dynamics simulations.

Hydrogen Bonding↗

PubMed2004

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↗

PubMed2004

Implicit solvation based on generalized Born theory in different dielectric environments.

In this paper we are investigating the effect of the dielectric environment on atomic Born radii used in generalized Born (GB) methods. Motivated by the Kirkwood expression for the reaction field of a single off-center charge in a spherical cavity, we are proposing extended formalisms for the calculation of Born radii as a function of external and internal dielectric constants. We demonstrate that reaction field energies calculated from environmentally dependent Born radii lead to much improved agreement with Poisson-Boltzmann solutions for low dielectric external environments, such as biological membranes or organic solvent, compared to previous methods where the calculation of Born radii does not depend on the environment. We also examine how this new approach can be applied for the calculation of transfer free energies from vacuum to a given external dielectric for a system with an internal dielectric larger than one. This has not been possible with standard GB theory but is relevant when scoring minimized or average structures with implicit solvent.

Biopolymers↗

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

Aberrant activity of the DNA repair enzyme AlkB.

Escherichia coli AlkB is a DNA/RNA repair enzyme containing a mononuclear Fe(II) site that couples the oxidative decomposition of alpha-ketoglutarate (alphaKG) to the hydroxylation of 1-methyladenine or 3-methylcytosine lesions in DNA or RNA, resulting in release of formaldehyde and restoration of the normal bases. In the presence of Fe(II), alphaKG, and oxygen, but the absence of methylated DNA, AlkB was found to catalyze an aberrant reaction that generates a blue chromophore. The color is proposed to derive from Fe(III) coordinated by a hydroxytryptophan at position 178 as revealed by mass spectrometric analysis. Protein structural modeling confirms that Trp 178 is reasonably positioned to react with the Fe(IV)-oxo intermediate proposed to form at the active site.

5-Hydroxytryptophan↗

PubMed2004

MMTSB Tool Set: enhanced sampling and multiscale modeling methods for applications in structural biology.

We describe the Multiscale Modeling Tools for Structural Biology (MMTSB) Tool Set (https://mmtsb.scripps.edu/software/mmtsbToolSet.html), which is a novel set of utilities and programming libraries that provide new enhanced sampling and multiscale modeling techniques for the simulation of proteins and nucleic acids. The tool set interfaces with the existing molecular modeling packages CHARMM and Amber for classical all-atom simulations, and with MONSSTER for lattice-based low-resolution conformational sampling. In addition, it adds new functionality for the integration and translation between both levels of detail. The replica exchange method is implemented to allow enhanced sampling of both the all-atom and low-resolution models. The tool set aims at applications in structural biology that involve protein or nucleic acid structure prediction, refinement, and/or extended conformational sampling. With structure prediction applications in mind, the tool set also implements a facility that allows the control and application of modeling tasks on a large set of conformations in what we have termed ensemble computing. Ensemble computing encompasses loosely coupled, parallel computation on high-end parallel computers, clustered computational grids and desktop grid environments. This paper describes the design and implementation of the MMTSB Tool Set and illustrates its utility with three typical examples--scoring of a set of predicted protein conformations in order to identify the most native-like structures, ab initio folding of peptides in implicit solvent with the replica exchange method, and the prediction of a missing fragment in a larger protein structure.

Computer Simulation↗

PubMed2004

Recent advances in the development and application of implicit solvent models in biomolecule simulations.

Advances have recently been made in the development of implicit solvent methodologies and their application to the modeling of biomolecules, particularly with regard to generalized Born approaches, dielectric screening function formulations and models based on solvent-accessible surface areas. Interesting new developments include more refined non-polar solvation energy estimators, and implicit methods for modeling low-dielectric and heterogeneous environments such as membrane systems. These have been successfully applied to molecular dynamics simulations, the scoring of protein conformations, and the calculation of binding affinities and folding free energy landscapes.

Algorithms↗