We present the newest version of the GROningen MOlecular Simulation program package, GROMOS96. GROMOS96 has been developed for the dynamic modelling of (bio)molecules using the methods of molecular dynamics, stochastic dynamics, and energy minimization as well as the path-integral formalism. An overview of its functionality is given, highlighting methodology not present in the last major release, GROMOS87. The organization of the code is outlined, and reliability, testing, and efficiency issues involved in the design of this large (73 000 lines of FORTRAN77 code) and complex package are discussed. Finally, we present two applications illustrating new functionality: local elevation simulation and molecular dynamics in four spatial dimensions.

This curled-in conformational change has never been previously described. Previous models of this movement, with the flaps as rigid levers, are not consistent with the experimental data. The residues that participate in this hydrophobic cluster as a result of the conformational change are highly sensitive to mutation and often contribute to drug resistance when they do change. However, several of these residues are not part of the active site cavity, and their essential role in causing drug resistance could possibly be rationalized if this conformational change actually occurs. Trapping HIV-1 protease in this inactive conformation would provide a unique opportunity for future drug design.

The structure and function of the synthetic innate defense regulator peptide 1018 was investigated. This 12 residue synthetic peptide derived by substantial modification of the bovine cathelicidin bactenecin has enhanced innate immune regulatory and moderate direct antibacterial activities. The solution state NMR structure of 1018 in zwitterionic dodecyl phosphocholine (DPC) micelles indicated an α-helical conformation, while secondary structures, based on circular dichroism measurements, in anionic sodium dodecyl sulfate (SDS) and phospholipid vesicles (POPC/PG in a 1:1 molar ratio) and simulations revealed that 1018 can adopt a variety of folds, tailored to its different functions. The structural data are discussed in light of the ability of 1018 to potently induce chemokine responses, suppress the LPS-induced TNF-α response, and directly kill both Gram-positive and Gram-negative bacteria.

No abstract

No abstract

Hydrophobic residues outside the active site of HIV-1 protease frequently mutate in patients undergoing protease inhibitor therapy; however, the mechanism by which these mutations confer drug resistance is not understood. From analysis of molecular dynamics simulations, 19 core hydrophobic residues appear to facilitate the conformational changes that occur in HIV-1 protease. The hydrophobic core residues slide by each other, exchanging one hydrophobic van der Waal contact for another, with little energy penalty, while maintaining many structurally important hydrogen bonds. Such hydrophobic sliding may represent a general mechanism by which proteins undergo conformational changes. Mutation of these residues in HIV-1 protease would alter the packing of the hydrophobic core, affecting the conformational flexibility of the protease. Therefore these residues impact the dynamic balance between processing substrates and binding inhibitors, and thus contribute to drug resistance.

Daptomycin, a cyclic anionic lipopeptide antibiotic, whose three-dimensional structure was recently solved using solution state NMR (Ball et al. 2004; Jung et al. 2004; Rotondi and Gierasch 2005), requires calcium for function. To date, the exact nature of the interaction between divalent cations, such as Ca(2+) or Mg(2+), has not been fully characterized. It has, however, been suggested that addition of Ca(2+) to daptomycin in a 1:1 molar ratio induces aggregation. Moreover, it has been suggested that certain residues, e.g. Asp3 and Asp7, which are essential for activity (Grunewald et al. 2004; Kopp et al. 2006), may also be important for Ca(2+) binding (Jung et al. 2004). In this work, we have tried: (1) to further pinpoint how Ca(2+) affects daptomycin structure/oligomerization using analytical ultracentrifugation; and (2) to determine whether a specific calcium binding site exists, based on one-dimensional (13)C NMR spectra and molecular dynamics (MD) simulations. The centrifugation results indicated that daptomycin formed micelles of between 14 and 16 monomers in the presence of a 1:1 molar ratio of Ca(2+) and daptomycin. The (13)C NMR data indicated that addition of calcium had a significant effect on the Trp1 and Kyn13 residues, indicating that either calcium binds in this region or that these residues may be important for oligomerization. Finally, the molecular dynamics simulation results indicated that the conformational change of daptomycin upon calcium binding might not be as significant as originally proposed. Similar studies on the divalent cation Mg(2+) are also presented. The implication of these results for the biological function of daptomycin is discussed.

Daptomycin is a cyclic anionic lipopeptide that exerts its rapid bactericidal effect by perturbing the bacterial cell membrane, a mode of action different from most other currently commercially available antibiotics (except e.g. polymyxin and gramicidin). Recent work has shown that daptomycin requires calcium in the form of Ca2+ to form a micellar structure in solution and to bind to bacterial model membranes. This evidence sheds light on the initial steps in the mechanism of action of this novel antibiotic. To understand how daptomycin goes on to perturb bacterial membranes, its three-dimensional structure has been determined in the presence of 1,2-dihexanoyl-sn-glycero-3-phosphocholine (DHPC) micelles. NMR spectra of daptomycin in DHPC were obtained under two conditions, namely in the presence of Ca2+ as used by Jung et al. [D. Jung, A. Rozek, M. Okon, R.E.W. Hancock, Structural transitions as determinants of the action of the calcium-dependent antibiotic daptomycin, Chem. Biol. 11 (2004) 949-57] to solve the calcium-conjugated structure of daptomycin in solution and in a phosphate buffer as used by Rotondi and Gierasch [K.S. Rotondi, L.M. Gierasch, A well-defined amphipathic conformation for the calcium-free cyclic lipopeptide antibiotic, daptomycin, in aqueous solution, Biopolymers 80 (2005) 374-85] to solve the structure of apo-daptomycin. The structures were calculated using molecular dynamics time-averaged refinement. The different sample conditions used to obtain the NMR spectra are discussed in light of fluorescence data, lipid flip-flop and calcein release assays in PC liposomes, in the presence and absence of Ca2+ [D. Jung, A. Rozek, M. Okon, R.E.W. Hancock, Structural transitions as determinants of the action of the calcium-dependent antibiotic daptomycin, Chem. Biol. 11 (2004) 949-57]. The implications of these results for the membrane perturbation mechanism of daptomycin are discussed.