Biomedical subjects

Carlos Simmerling

Publications and source records attributed to Carlos Simmerling.

18 recordsLinked to original sources

PubMed2006

Comparison of multiple Amber force fields and development of improved protein backbone parameters.

The ff94 force field that is commonly associated with the Amber simulation package is one of the most widely used parameter sets for biomolecular simulation. After a decade of extensive use and testing, limitations in this force field, such as over-stabilization of alpha-helices, were reported by us and other researchers. This led to a number of attempts to improve these parameters, resulting in a variety of "Amber" force fields and significant difficulty in determining which should be used for a particular application. We show that several of these continue to suffer from inadequate balance between different secondary structure elements. In addition, the approach used in most of these studies neglected to account for the existence in Amber of two sets of backbone phi/psi dihedral terms. This led to parameter sets that provide unreasonable conformational preferences for glycine. We report here an effort to improve the phi/psi dihedral terms in the ff99 energy function. Dihedral term parameters are based on fitting the energies of multiple conformations of glycine and alanine tetrapeptides from high level ab initio quantum mechanical calculations. The new parameters for backbone dihedrals replace those in the existing ff99 force field. This parameter set, which we denote ff99SB, achieves a better balance of secondary structure elements as judged by improved distribution of backbone dihedrals for glycine and alanine with respect to PDB survey data. It also accomplishes improved agreement with published experimental data for conformational preferences of short alanine peptides and better accord with experimental NMR relaxation data of test protein systems.

Alanine↗

PubMed2006

The open structure of a multi-drug-resistant HIV-1 protease is stabilized by crystal packing contacts.

The introduction of HIV-1 protease (HIV-PR) inhibitors has led to a dramatic increase in patient survival; however, these gains are threatened by the emergence of multi-drug-resistant strains. Design of inhibitors that overcome resistance would be greatly facilitated by deeper insight into the mechanistic events associated with binding of substrates and inhibitors, as well as an understanding of the effects of resistance mutations on the structure and dynamic behavior of HIV-PR. We previously reported a series of simulations that provide a model for HIV-PR dynamics, with spontaneous conversions between the bound and unbound crystal forms upon addition or removal of an inhibitor. Importantly, the unbound protease transiently sampled a third fully open state that permits entry to the active site, unlike both crystallographic forms. Recently, a crystal structure of unbound HIV-PR was reported for the MDR 769 isolate (PDB: 1TW7); unlike all previous experimental structures, the binding pocket is open. It is suggested that drug resistance in this strain arises at least in part from the inability of inhibitors to induce closing. We carried out simulations of the MDR 769 HIV-PR mutant and observed that the reported structure is unstable in solution and rapidly adopts the semi-open conformation observed for the unbound wild-type protease in solution. Further analysis suggests that the wide-open structure observed for MDR 769 arises not from sequence variation, but instead is an artifact from crystal packing. Thus, despite being the first experimental structure to reveal flap opening sufficient for substrate access to the active site, this structure may not be directly relevant to studies of inhibitor entry or to the cause of HIV-PR drug resistance.

Chemistry, Pharmaceutical↗

PubMed2006

Structure of acyl carrier protein bound to FabI, the FASII enoyl reductase from Escherichia coli.

Acyl carrier proteins play a central role in metabolism by transporting substrates in a wide variety of pathways including the biosynthesis of fatty acids and polyketides. However, despite their importance, there is a paucity of direct structural information concerning the interaction of ACPs with enzymes in these pathways. Here we report the structure of an acyl-ACP substrate bound to the Escherichia coli fatty acid biosynthesis enoyl reductase enzyme (FabI), based on a combination of x-ray crystallography and molecular dynamics simulation. The structural data are in agreement with kinetic studies on wild-type and mutant FabIs, and reveal that the complex is primarily stabilized by interactions between acidic residues in the ACP helix alpha2 and a patch of basic residues adjacent to the FabI substrate-binding loop. Unexpectedly, the acyl-pantetheine thioester carbonyl is not hydrogen-bonded to Tyr(156), a conserved component of the short chain alcohol dehydrogenase/reductase superfamily active site triad. FabI is a proven target for drug discovery and the present structure provides insight into the molecular determinants that regulate the interaction of ACPs with target proteins.

Acyl Carrier Protein↗

PubMed2006

Computational analysis of the mode of binding of 8-oxoguanine to formamidopyrimidine-DNA glycosylase.

8-Oxoguanine (8OG) is the most prevalent form of oxidative DNA damage. In bacteria, 8OG is excised by formamidopyrimidine glycosylase (Fpg) as the initial step in base excision repair. To efficiently excise this lesion, Fpg must discriminate between 8OG and an excess of guanine in duplex DNA. In this study, we explore the structural basis underlying this high degree of selectivity. Two structures have been reported in which Fpg is bound to DNA, differing with respect to the position of the lesion in the active site, one structure showing 8OG bound in the syn conformation and the other in the anti conformation. Remarkably, the results of our all-atom simulations are consistent with both structures. The syn conformation observed in the crystallographic structure of Fpg obtained from Bacillus stearothermophilus is stabilized through interaction with E77, a nonconserved residue. Replacement of E77 with Ser, creating the Fpg sequence found in Escherichia coli and other bacteria, results in preferred binding of 8OG in the anti conformation. Our calculations provide novel insights into the roles of active site residues in binding and recognition of 8OG by Fpg.

Amino Acid Sequence↗

PubMed2006

Insight through molecular mechanics Poisson-Boltzmann surface area calculations into the binding affinity of triclosan and three analogues for FabI, the E. coli enoyl reductase.

Keeping pace with emerging drug resistance in clinically important pathogens will be greatly aided by inexpensive yet reliable computational methods that predict the binding affinities of ligands for drug targets. We present results using the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method to calculate the affinity of a series of triclosan analogues for the E. coli enoyl reductase FabI, spanning a 450000-fold range of binding affinities. Significantly, a high correlation is observed between the calculated binding energies and those determined experimentally. Further examination indicates that the van der Waals energies are the most correlated component of the total affinity (r2 = 0.74), indicating that the shape of the inhibitor is very important in defining the binding energies for this system. The validation of MM-PBSA for the E coli FabI system serves as a platform for inhibitor design efforts focused on the homologous enzyme in Staphylococcus aureus and Mycobacterium tuberculosis.

Binding Sites↗

PubMed2006

HIV-1 protease flaps spontaneously close to the correct structure in simulations following manual placement of an inhibitor into the open state.

We report unrestrained, all-atom molecular dynamics simulations of HIV-1 protease (HIV-PR) with a continuum solvent model that reproducibly sample closing of the active site flaps following manual placement of a cyclic urea inhibitor into the substrate binding site of the open protease. The open form was obtained from the unbound, semi-open HIV-PR crystal structure, which we recently reported (Hornak, V.; et al. Proc. Natl. Acad. Sci. U.S.A. 2006, 103, 915-920.) to have spontaneously opened during unrestrained dynamics. In those simulations, the transiently open flaps always returned to the semi-open form that is observed in all crystal structures of the free protease. Here, we show that manual docking of the inhibitor reproducibly induces spontaneous conversion to the closed form as seen in all inhibitor-bound HIV-PR crystal structures. These simulations reproduced not only the greater degree of flap closure, but also the striking difference in flap "handedness" between bound and free enzyme. In most of the simulations, the final structures were highly accurate. Root-mean-square deviations (RMSD) from the crystal structure of the complex were approximately 1.5 A (averaged over the last 100 ps) for the inhibitor and each flap despite initial RMSD of 2-5 A for the inhibitors and 6-11 A for the flaps. Key hydrogen bonds were formed between the flap tips and between flaps and inhibitor that match those seen in the crystal structure. The results demonstrate that all-atom simulations have the ability to significantly improve poorly docked ligand conformations and reproduce large-scale receptor conformational changes that occur upon binding.

Azepines↗

PubMed2006

HIV-1 protease flaps spontaneously open and reclose in molecular dynamics simulations.

We report unrestrained, all-atom molecular dynamics simulations of HIV-1 protease that sample large conformational changes of the active site flaps. In particular, the unliganded protease undergoes multiple conversions between the "closed" and "semiopen" forms observed in crystal structures of inhibitor-bound and unliganded protease, respectively, including reversal of flap "handedness." Simulations in the presence of a cyclic urea inhibitor yield stable closed flaps. Furthermore, we observe several events in which the flaps of the unliganded protease open to a much greater degree than observed in crystal structures and subsequently return to the semiopen state. Our data strongly support the hypothesis that the unliganded protease predominantly populates the semiopen conformation, with closed and fully open structures being a minor component of the overall ensemble. The results also provide a model for the flap opening and closing that is considered to be essential to enzyme function.

Acquired Immunodeficiency Syndrome↗

PubMed2005

Dynamic behavior of DNA base pairs containing 8-oxoguanine.

The process by which DNA repair enzymes recognize and selectively excise damaged bases in duplex DNA is fundamental to our mechanistic understanding of these critical biological reactions. 8-Oxoguanine (8-oxoG) is the most common form of oxidative DNA damage; unrepaired, this lesion generates a G:C-->T:A mutation. Central to the recognition and repair of DNA damage is base extrusion, a process in which the damaged base lesion or, in some cases, its partner disengages from the helix and is bound to the enzyme's active site where base excision takes place. The conformation adopted by 8-oxoG in duplex DNA is affected by the base positioned opposite this lesion; conformational changes may also take place when the damaged base binds to its cognate repair enzyme. We performed unrestrained molecular dynamics simulations for several 13-mer DNA duplexes. Oligomers containing G:C and 8oxoG:C pairs adopted Watson-Crick geometries in stable B-form duplexes; 8oxoG showed increased local and global flexibility and a reduced barrier to base extrusion. Duplexes containing the G:A mismatch showed much larger structural fluctuations and failed to adopt a well-defined structure. For the 8oxoG:A mismatch that is recognized by the DNA glycosylase MutY, the damaged nucleoside underwent spontaneous and reproducible anti-->syn transitions. The syn conformation is thermodynamically preferred. Steric hindrance and unfavorable electrostatics associated with the 8oxoG O8 atom in the anti conformation were the major driving forces for this transition. Transition events follow two qualitatively different pathways. The overall anti-->syn transition rate and relative probability of the two transition paths were dependent on local sequence context. These simulations indicate that both the dynamic and equilibrium behavior of the duplex change as a result of oxidation; these differences may provide valuable new insight into the selective action of enzymes on damaged DNA.

Base Pairing↗

PubMed2005

Folding cooperativity in a three-stranded beta-sheet model.

The thermodynamic behavior of a previously designed three-stranded beta-sheet was studied via several microseconds of standard and replica exchange molecular dynamics simulations. The system is shown to populate at least four thermodynamic minima, including two partially folded states in which only a single hairpin is formed. Simulated melting curves show different profiles for the C and N-terminal hairpins, consistent with differences in secondary structure content in published NMR and CD/FTIR measurements, which probed different regions of the chain. Individual beta-hairpins that comprise the three-stranded beta-sheet are observed to form cooperatively. Partial folding cooperativity between the component hairpins is observed, and good agreement between calculated and experimental values quantifying this cooperativity is obtained when similar analysis techniques are used. However, the structural detail in the ensemble of conformations sampled in the simulations permits a more direct analysis of this cooperativity than has been performed on the basis of experimental data. The results indicate the actual folding cooperativity perpendicular to strand direction is significantly larger than the lower bound obtained previously.

Computer Simulation↗

PubMed2005

Modified replica exchange simulation methods for local structure refinement.

Parallel tempering, also known as replica exchange molecular dynamics (REMD), has recently been successfully used to study the structure and thermodynamic properties of biomolecules such as peptides and small proteins. For large systems, however, applying REMD can be costly since the number of replicas needed increases as the square root of the number of degrees of freedom in the system. Often, enhanced sampling is only needed for a subset of atoms, such as a loop region of a large protein or a small ligand binding to a receptor. In such applications, it is often reasonable to assume a weak dependence of the structure of the larger region on the instantaneous conformation of the smaller region of interest. For these cases, we derived two variant replica exchange methods, partial replica exchange molecular dynamics (PREMD) and local replica exchange molecular dynamics (LREMD). The Hamiltonian for the system is separated, with replica exchange carried out only for terms involving the subsystem of interest while the remainder of the system is maintained at a single temperature. The number of replicas required for efficient exchange thus depends on the number of degrees of freedom in the fragment needing refinement rather than on the size of the full system. The method can be applied to much larger systems than was previously practical. This also provides a means to preserve the integrity of the structure outside the refinement region without introduction of restraints. LREMD takes this weak coupling approximation a step further, employing only a single representation of the large fragment that simultaneously interacts with all of the replicas of the subsystem of interest. This is obtained by combining replica exchange with the locally enhanced sampling approximation (LES), reducing the computational expense of replica exchange simulations to near that of a single standard molecular dynamics (MD) simulation. Use of LREMD also permits the use of LES without requiring the specification of a single temperature, a known difficulty for standard LES simulations. We tested these two methods on the loop region of an RNA hairpin model system and find significant advantages over standard MD and REMD simulations.

Algorithms↗

PubMed2005

Structural requirements of the extracellular to transmembrane domain junction for erythropoietin receptor function.

The erythropoietin receptor (EpoR) is crucial for erythrocyte formation. The x-ray crystal structures of the EpoR extracellular domain lack the juxtamembrane (JM) region and the junction to the transmembrane (TM) domain. Yet the JM-TM regions are important for transmitting the conformational change imposed on the receptor dimer by Epo binding. Cysteine-scanning mutagenesis of the JM-TM regions identified three novel constitutively active mutants, demonstrating close disulfide-bonded juxtapositioning of these residues in the JM (L223C) and N-terminal TM domain (L226C, I227C). Chemical cross-linking defined the interface of the active helical TM dimer and revealed that the JM-TM segment encompassing Leu(226)-Leu(230) is non-helical. Molecular dynamics and NMR studies indicated that the TM-JM junction forms an N-terminal helix cap. This structure is important for EpoR function because replacement of this motif by consecutive leucines rendered the receptor constitutively active.

Amino Acid Motifs↗

PubMed2005

The Amber biomolecular simulation programs.

We describe the development, current features, and some directions for future development of the Amber package of computer programs. This package evolved from a program that was constructed in the late 1970s to do Assisted Model Building with Energy Refinement, and now contains a group of programs embodying a number of powerful tools of modern computational chemistry, focused on molecular dynamics and free energy calculations of proteins, nucleic acids, and carbohydrates.

Algorithms↗

PubMed2004

Inhibition of the bacterial enoyl reductase FabI by triclosan: a structure-reactivity analysis of FabI inhibition by triclosan analogues.

To explore the molecular basis for the picomolar affinity of triclosan for FabI, the enoyl reductase enzyme from the type II fatty acid biosynthesis pathway in Escherichia coli, an SAR study has been conducted using a series of triclosan analogues. Triclosan (1) is a slow, tight-binding inhibitor of FabI, interacting specifically with the E.NAD(+) form of the enzyme with a K(1) value of 7 pM. In contrast, 2-phenoxyphenol (2) binds with equal affinity to the E.NAD(+) (K(1) = 0.5 microM) and E.NADH (K(2) = 0.4 microM) forms of the enzyme and lacks the slow-binding step observed for triclosan. Thus, removal of the three triclosan chlorine atoms reduces the affinity of the inhibitor for FabI by 70,000-fold and removes the preference for the E.NAD(+) FabI complex. 5-Chloro-2-phenoxyphenol (3) is a slow, tight-binding inhibitor of FabI and binds to the E.NAD(+) form of the enzyme (K(1) = 1.1 pM) 7-fold more tightly than triclosan. Thus, while the two ring B chlorine atoms are not required for FabI inhibition, replacement of the ring A chlorine increases binding affinity by 450,000-fold. Given this remarkable observation, the SAR study was extended to the 5-fluoro-2-phenoxyphenol (4) and 5-methyl-2-phenoxyphenol (5) analogues to further explore the role of the ring A substituent. While both 4 and 5 are slow, tight-binding inhibitors, they bind substantially less tightly to FabI than triclosan. Compound 4 binds to both E.NAD(+) and E.NADH forms of the enzyme with K(1) and K(2) values of 3.2 and 240 nM, respectively, whereas compound 5 binds exclusively to the E.NADH enzyme complex with a K(2) value of 7.2 nM. Thus, the ring A substituent is absolutely required for slow, tight-binding inhibition. In addition, pK(a) measurements coupled with simple electrostatic calculations suggest that the interaction of the ring A substituent with F203 is a major factor in governing the affinity of analogues 3-5 for the FabI complex containing the oxidized form of the cofactor.

Anti-Bacterial Agents↗

PubMed2004

Development of softcore potential functions for overcoming steric barriers in molecular dynamics simulations.

In this work, we describe the development of softcore potential functions that permit occasional "tunneling" through the regions of conformational space during molecular dynamics (MD) simulations, which would otherwise be sterically prohibited. The modification consists of a truncation of the nonbonded interaction before the steeply repulsive region encountered at short interatomic distances. This modification affects both Lennard-Jones and Coulomb parts of the nonbonded potential. Critical to success is the choice of appropriate pairwise switching distances at which this modification should be made. In the present work, these are calculated based on potential of mean force functions extracted from model system molecular dynamics simulations. We believe that these functions describe the dynamic short-range interactions much better than mean force potentials derived from an ensemble of static structures (e.g. protein data bank (PDB)). Once a set of mean force potentials is obtained, a single empirical parameter, effective barrier height, is employed to determine switching distances for all pairwise atomic interactions. Changing this single parameter allows adjustment of the "softness" of the whole system. We tested the applicability of the new softcore potentials in a loop structure optimization study. The H1 loop in the antibody 17/9 was selected as our test case because substantial repacking of loop residues in the dense protein environment is necessary for successful relaxation of random initial conformations. Softcore simulations converted to correct loop conformations, in contrast to standard simulations which never sampled this structure even after 10 ns. The resulting root mean square deviation (RMSD) values (below 1.3 A for all heavy atoms of the loop) demonstrate the usefulness of the approach based on mean force derived softcore functions.

Amino Acids↗

PubMed2003

Generation of accurate protein loop conformations through low-barrier molecular dynamics.

Prediction and refinement of protein loop structures are important and challenging tasks for which no general solution has been found. In addition to the accuracy of scoring functions, the main problems reside in (1) insufficient statistical sampling and (2) crossing energy barriers that impede conformational rearrangements of the loop. We approach these two issues by using "low-barrier molecular dynamics," a combination of energy smoothing techniques. To address statistical sampling, locally enhanced sampling (LES) is used to produce multiple copies of the loop, thus improving statistics and reducing energy barriers. We introduce a novel extension of LES that can improve local sampling even further through hierarchical subdivision of copies. Even though LES reduces energy barriers, it cannot provide for crossing infinite barriers, which can be problematic when substantial rearrangement of residues is necessary. To permit this kind of loop residue repacking, a "soft-core" potential energy function is introduced, so that atomic overlaps are temporarily allowed. We tested this new combined methodology to a loop in anti-influenza antibody Fab 17/9 (7 residues long) and to another loop in the antiprogesterone antibody DB3 (8 residues). In both cases, starting from random conformations, we were able to locate correct loop structures (including sidechain orientations) with heavy-atom root-mean-square deviation (fit to the nonloop region) of approximately 1.1 A in Fab 17/9 and approximately 1.8 A in DB3. We show that the combination of LES and soft-core potential substantially improves sampling compared to regular molecular dynamics. Moreover, the sampling improvement obtained with this combined approach is significantly better than that provided by either of the two methods alone.

Algorithms↗

PubMed2003

Using PC clusters to evaluate the transferability of molecular mechanics force fields for proteins.

The transferability of molecular mechanics parameters derived for small model systems to larger biopolymers such as proteins can be difficult to assess. Even for small peptides, molecular dynamics simulations are typically too short to sample structures significantly different than initial conformations, making comparison to experimental data questionable. We employed a PC cluster to generate large numbers of native and non-native conformations for peptides with experimentally measured structural data, one predominantly helical and the other forming a beta-hairpin. These atomic-detail sets do not suffer from slow convergence, and can be used to rapidly evaluate important force field properties. In this case a suspected bias toward alpha-helical conformations in the ff94 and ff99 force fields distributed with the AMBER package was verified. The sets provide critical feedback not only on force field transferability, but may also predict modifications for improvement. Such predictions were used to modify the ff99 parameter set, and the resulting force field was used to test stability and folding of model peptides. Structural behavior during molecular dynamics with the modified force field is found to be very similar to expectations, suggesting that these basis sets of conformations may themselves have significant transferability among force fields. We continue to improve and expand this data set and plan to make it publicly accessible. The calculations involved in this process are trivially parallel and can be performed using inexpensive personal computers with commodity components.

Computer Simulation↗

PubMed2002

Conformational heterogeneity observed in simulations of a pyrene-substituted DNA.

NMR studies previously carried out for a DNA system with a pyrene-substituted base did not observe NOEs involving the adenine located 5' to the pyrene, and thus the conformation of the adenine was poorly defined in the resulting family of refined structures. However, chemical shift data suggested that an AT base pair may be present. We have carried out fully unrestrained molecular dynamics simulations starting from several members of the family of structures, and these simulations support the existence of an AT base pair for this region. Simulations in both explicit and implicit solvent were carried out, with each converging to either anti or syn conformation for adenine and base pairing in all cases. During these simulations, large and dramatic conformational changes are observed that suggest pathways for complex conformational changes in the highly packed DNA interior. Our analysis reveals little difference in the energies of these syn and anti conformations, in contrast to control calculations carried out for standard DNA (in the absence of a neighboring pyrene). While no interconversion between the conformations was observed in standard simulations, reversible anti/syn exchange was directly simulated using the locally enhanced sampling approach. No exchange was seen in the non-pyrene control sequence. Together, these results suggest that an increased flexibility is introduced as a consequence of the pyrene substitution, offering an explanation that is consistent with the available experimental data. These results increase our optimism that simulations in atomic detail may provide accurate models for experimental observations in complex systems.

Computer Simulation↗

PubMed2002

All-atom structure prediction and folding simulations of a stable protein.

We present results from all-atom, fully unrestrained ab initio folding simulations for a stable protein with nontrivial secondary structure elements and a hydrophobic core. The construct, "trpcage", is a 20-residue sequence optimized by the Andersen group at University of Washington and is currently the smallest protein that displays two-state folding properties. Compared over the well-defined regions of the experimental structure, our prediction has a remarkably low 0.97 A Calpha root-mean-square-deviation (rmsd) and 1.4 A for all heavy atoms. The simulated structure family displays additional features that are suggested by experimental data, yet are not evident in the family of NMR-derived structures.

Amino Acid Sequence↗