Search

Computer simulations of molecular systems enable structure-energy-function relationships of molecular processes to be described at the sub-atomic, atomic, supra-atomic or supra-molecular level and plays an increasingly important role in chemistry, biology and physics. To interpret the results of such simulations appropriately, the degree of uncertainty and potential errors affecting the calculated properties must be considered. Uncertainty and errors arise from (1) assumptions underlying the molecular model, force field and simulation algorithms, (2) approximations implicit in the interatomic interaction function (force field), or when integrating the equations of motion, (3) the chosen values of the parameters that determine the accuracy of the approximations used, and (4) the nature of the system and the property of interest. In this overview, advantages and shortcomings of assumptions and approximations commonly used when simulating bio-molecular systems are considered. What the developers of bio-molecular force fields and simulation software can do to facilitate and broaden research involving bio-molecular simulations is also discussed.

Biological membranes display a great diversity in lipid composition and lateral structure that is crucial in a variety of cellular functions. Simulations of membranes have contributed significantly to the understanding of the properties, functions and behaviour of membranes and membrane-protein assemblies. This success relies on the ability of the force field used to describe lipid-lipid and lipid-environment interactions accurately, reproducibly and realistically. In this review, we present some recent progress in lipid force-field development and validation strategies. In particular, we highlight how a range of properties obtained from various experimental techniques on lipid bilayers and membranes, can be used to assess the quality of a force field. We discuss the limitations and assumptions that are inherent to both computational and experimental approaches and how these can influence the comparison between simulations and experimental data. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

The free-energy landscape of a small protein, the FBP 28 WW domain, has been explored using molecular dynamics (MD) simulations with alternative descriptions of the molecule. The molecular models used range from coarse-grained to all-atom with either an implicit or explicit treatment of the solvent. Sampling of conformation space was performed using both conventional and temperature-replica exchange MD simulations. Experimental chemical shifts and NOEs were used to validate the simulations, and experimental phi values both for validation and as restraints. This combination of different approaches has provided insight into the free energy landscape and barriers encountered by the protein during folding and enabled the characterization of native, denatured and transition states which are compatible with the available experimental data. All the molecular models used stabilize well defined native and denatured basins; however, the degree of agreement with the available experimental data varies. While the most detailed, explicit solvent model predicts the data reasonably accurately, it does not fold despite a simulation time 10 times that of the experimental folding time. The less detailed models performed poorly relative to the explicit solvent model: an implicit solvent model stabilizes a ground state which differs from the experimental native state, and a structure-based model underestimates the size of the barrier between the two states. The use of experimental phi values both as restraints, and to extract structures from unfolding simulations, result in conformations which, although not necessarily true transition states, appear to share the geometrical characteristics of transition state structures. In addition to characterizing the native, transition and denatured states of this particular system in this work, the advantages and limitations of using varying levels of representation are discussed. (C) 2008 Wiley Periodicals, Inc. J Comput Chem 30: 1059-1068, 2009

A large variety of antimicrobial peptides have been shown to act, at least in vitro, by potation of the lipid membrane. The nanometre size of these pores, however, complicates their structural characterization by experimental techniques. Here we use molecular dynamics simulations, to study the interaction of a specific class of antimicrobial peptides, melittin, with a dipalmitoylphosphatidylcholine bilayer in atomic detail. We show that transmembrane pores spontaneously form above a critical peptide to lipid ratio. The lipid molecules bend inwards to form a toroidally shaped pore but with only one or two peptides lining the pore. This is in strong contrast to the traditional models of toroidal pores in which the peptides are assumed to adopt a transmembrane orientation. We find that peptide aggregation, either prior or after binding to the membrane surface. is a prerequisite to pore formation. The presence of a stable helical secondary structure of the peptide, however is not. Furthermore, results obtained with modified peptides point to the importance of electrostatic interactions in the potation process. Removing the charges of the basic amino-acid residues of melittin prevents pore formation. It was also found that in the absence of counter ions pores not only form more rapidly but lead to membrane rupture. The rupture process occurs via a novel recursive potation pathway, which we coin the Droste mechanism. (c) 2008 Elsevier B.V. All rights reserved.

A method is presented to enhance the efficiency of simulations of lipid vesicles. The method increases computational speed by eliminating water molecules that either surround the vesicle or reside in the interior of the vesicle, without altering the properties of the water at the membrane interface. Specifically, mean field force approximation (MFFA) boundary potentials are used to replace both the internal and external excess bulk solvent. In addition to reducing the cost of simulating preformed vesicles, the molding effect of the boundary potentials also enhances the formation and equilibration of vesicles from random solutions of lipid in water. Vesicles with diameters in the range from 20 to 60 nm were obtained on a nanosecond time scale, without any noticeable effect of the boundary potentials on their structure.

Electrophoresis is widely used to determine the electrostatic potential of colloidal particles. Oil droplets in pure water show negative or positive electrophoretic mobilities depending on the pH. This is commonly attributed to the adsorption of hydroxyl or hydronium ions, resulting in a negative or positive surface charge, respectively. This explanation, however, is not in agreement with the difference in isoelectric point and point of zero charge observed in experiment. Here we present molecular dynamics simulations of oil droplets in water in the presence of an external electric field but in the absence of any ions. The simulations reproduce the negative sign and the order of magnitude of the oil droplet mobilities at the point of zero charge in experiment. The electrostatic potential in the oil with respect to the water phase, induced by anisotropic dipole orientation in the interface, is positive. Our results suggest that electrophoretic mobility does not always reflect the net charge or electrostatic potential of a suspended liquid droplet and, thus, the interpretation of electrophoresis in terms of purely continuum effects may need to be reevaluated. (C) 2007 Elsevier Inc. All rights reserved.

Class I hydrophobins are fungal proteins that self-assemble into robust amphipathic rodlet monolayers on the surface of aerial structures such as spores and fruiting bodies. These layers share many structural characteristics with amyloid fibrils and belong to the growing family of functional amyloid-like materials produced by microorganisms. Although the three-dimensional structure of the soluble monomeric form of a class I hydrophobin has been determined, little is known about the molecular structure of the rodlets or their assembly mechanism. Several models have been proposed, some of which suggest that the Cys3–Cys4 loop has a critical role in the initiation of assembly or in the polymeric structure. In order to provide insight into the relationship between hydrophobin sequence and rodlet assembly, we investigated the role of the Cys3–Cys4 loop in EAS, a class I hydrophobin from Neurospora crassa. Remarkably, deletion of up to 15 residues from this 25-residue loop does not impair rodlet formation or reduce the surface activity of the protein, and the physicochemical properties of rodlets formed by this mutant are indistinguishable from those of its full-length counterpart. In addition, the core structure of the truncation mutant is essentially unchanged. Molecular dynamics simulations carried out on the full-length protein and this truncation mutant binding to an air–water interface show that, although it is hydrophobic, the loop does not play a role in positioning the protein at the surface. These results demonstrate that the Cys3–Cys4 loop does not have an integral role in the formation or structure of the rodlets and that the major determinant of the unique properties of these proteins is the amphipathic core structure, which is likely to be preserved in all hydrophobins despite the high degree of sequence variation across the family.

A series of molecular dynamics simulations in explicit solvent were started from nine structural models of the transition state of the SH3 domain of alpha-spectrin, which were generated by Lindorff Larsen et al. (Nat Struct Mol Biol 2004;11:443-449) using molecular dynamics simulations in which experimental (Phi-values were incorporated as restraints. Two of the nine models were simulated 10 times for 200 ns and the remaining models simulated two times for 200 us. Complete folding was observed in one case, while in the other simulations partial folding or unfolding even ts were observed, which were characterized by a regularization of elements of secondary structure. These results are consistent with recent experimental evidence that the folding of SH3 domains involves low populated intermediate states. (C) 2007 Wiley-Liss, Inc.

The question of why plants release isoprene when heat stressed has been continuously debated for more than half a century. In this work we use molecular dynamics simulation techniques to directly investigate the interaction between isoprene and a model phospholipid membrane in atomic detail. It is found that isoprene partitions preferentially in the center of the membrane and in a dose dependent manner enhances the order within the membrane without significantly changing the dynamical properties of the system. At a concentration of 20 mol% isoprene (16 isoprene molecules per 64 lipid molecules) the effect of the addition of isoprene on the membrane order is equivalent to a reduction in temperature of 10 K, rising to a reduction of 30 K at 43 mol% isoprene. The significance of the work is that it provides for the first time direct evidence that isoprene stabilizes lipid membranes and reduces the likelihood of a phospholipid membrane undergoing a heat induced phase transition. Furthermore it provides a clear mechanistic picture as to why plants specifically utilize isoprene for this purpose.

Hydrophobins are small ( approximately 100 aa) proteins that have an important role in the growth and development of mycelial fungi. They are surface active and, after secretion by the fungi, self-assemble into amphipathic membranes at hydrophobic/hydrophilic interfaces, reversing the hydrophobicity of the surface. In this study, molecular dynamics simulation techniques have been used to model the process by which a specific class I hydrophobin, SC3, binds to a range of hydrophobic/hydrophilic interfaces. The structure of SC3 used in this investigation was modeled based on the crystal structure of the class II hydrophobin HFBII using the assumption that the disulfide pairings of the eight conserved cysteine residues are maintained. The proposed model for SC3 in aqueous solution is compact and globular containing primarily beta-strand and coil structures. The behavior of this model of SC3 was investigated at an air/water, an oil/water, and a hydrophobic solid/water interface. It was found that SC3 preferentially binds to the interfaces via the loop region between the third and fourth cysteine residues and that binding is associated with an increase in alpha-helix formation in qualitative agreement with experiment. Based on a combination of the available experiment data and the current simulation studies, we propose a possible model for SC3 self-assembly on a hydrophobic solid/water interface.

Molecular dynamics simulations of the magainin MG-H2 peptide interacting with a model phospholipid membrane have been used to investigate the mechanism by which antimicrobial peptides act. Multiple copies of the peptide were randomly placed in solution close to the membrane. The peptide readily bound to the membrane, and above a certain concentration, the peptide was observed to cooperatively induce the formation of a nanometer- sized, toroidally shaped pore in the bilayer. In sharp contrast with the commonly accepted model of a toroidal pore, only one peptide was typically found near the center of the pore. The remaining peptides lay close to the edge of the pore, maintaining a predominantly parallel orientation with respect to the membrane.

Molecular dynamics simulations have been used to study the phase behavior of a dipalmitoylphosphatidylcholine (DPPC)/palmitic acid (PA)/water 1:2:20 mixture in atomic detail. Starting from a random solution of DPPC and PA in water, the system adopts either a gel phase at temperatures below approximately 330 K or an inverted hexagonal phase above approximately 330 K in good agreement with experiment. It has also been possible to observe the direct transformation from a gel to an inverted hexagonal phase at elevated temperature (approximately 390 K). During this transformation, a metastable fluid lamellar intermediate is observed. Interlamellar connections or stalks form spontaneously on a nanosecond time scale and subsequently elongate, leading to the formation of an inverted hexagonal phase. This work opens the possibility of studying in detail how the formation of nonlamellar phases is affected by lipid composition and (fusion) peptides and, thus, is an important step toward understanding related biological processes, such as membrane fusion.

Bacterial chaperonin, GroEL, together with its co-chaperonin, GroES, facilitates the folding of a variety of polypeptides. Experiments suggest that GroEL stimulates protein folding by multiple cycles of binding and release. Misfolded proteins first bind to an exposed hydrophobic surface on GroEL. GroES then encapsulates the substrate and triggers its release into the central cavity of the GroEL/ES complex for folding. In this work, we investigate the possibility to facilitate protein folding in molecular dynamics simulations by mimicking the effects of GroEL/ES namely, repeated binding and release, together with spatial confinement. During the binding stage, the (metastable) partially folded proteins are allowed to attach spontaneously to a hydrophobic surface within the simulation box. This destabilizes the structures, which are then transferred into a spatially confined cavity for folding. The approach has been tested by attempting to refine protein structural models generated using the ROSETTA procedure for ab initio structure prediction. Dramatic improvements in regard to the deviation of protein models from the corresponding experimental structures were observed. The results suggest that the primary effects of the GroEL/ES system can be mimicked in a simple coarse-grained manner and be used to facilitate protein folding in molecular dynamics simulations. Furthermore, the results Sur port the assumption that the spatial confinement in GroEL/ES assists the folding of encapsulated proteins.

Extensive molecular dynamic simulations (similar to 240 ns) have been used to investigate the conformational behavior of PrP106-126 prion peptide in four different environments (water, dimethyl sulfoxide, hexane, and trifluoroethanol) and under both neutral and acidic conditions. The conformational polymorphism of PrP106-126 in solution observed in the simulations supports the role of this fragment in the structural transition of the native to the abnormal form of prion protein in response to changes in the local environmental conditions. The peptide in solution is primarily unstructured. The simulations show an increased presence of helical structure in an apolar solvent, in agreement with the results from circular dichroism spectroscopy. In water solution, beta-sheet elements were observed between residues 108-112 and either residues 115-121 or 121-126. An alpha-beta transition was observed under neutral conditions. In DMSO, the peptide adopted an extended conformation, in agreement with nuclear magnetic resonance experiments.

The formation, relative stability, and possible stoichiometries of two (self-)complementary peptide sequences (B and E) designed to form either a parallel homodimeric (B + B) or an antiparallel heterodimeric (B + E) coiled coil have been investigated. Peptide B shows a characteristic coiled coil pattern in circular dichroism spectra at pH 7.4, whereas peptide E is apparently random coiled under these conditions. The peptides are complementary to each other, with peptide E forming a coiled coil when mixed with peptide B. Molecular dynamics simulations show that combinations of B + B and B + E readily form a dimeric coiled coil, whereas E + E does not fall in line with the experimental data. However, the simulations strongly suggest the preferred orientation of the helices in the homodimeric coiled coil is antiparallel, with interactions at the interface quite different to that of the idealized model. In addition, molecular dynamics simulations suggest equilibrium between dimers, trimers, and tetramers of alpha-helices for peptide B.

In this work, we compare the results of molecular dynamics simulations involving the application of three generalized Born (GB) models to 10 different proteins. The three GB models, the Still, HCT, and modified analytical generalized Born models, were implemented in the computationally efficient gromacs package. The performance of each model was assessed from the backbone rms deviation from the native structure, the number of native hydrogen bonds retained in the simulation, and the experimental and calculated radius of gyration. Analysis of variance (ANOVA) was used to analyze the results of the simulations. The rms deviation measure was found to be unable to distinguish the quality of the results obtained with the three different GB models, whereas the number of native hydrogen bonds and radius of gyration yielded a statistically meaningful discrimination among models. Our results suggest that, of the three, modified analytical generalized Born yields the best agreement between calculated and experimentally derived structures. More generally, our study highlights the need both to evaluate the effects of different variables on the results of simulations and to verify that the results of molecular dynamics simulations are statistically meaningful.

Molecular dynamics simulations of lecithin lipid bilayers in water as they are cooled from the liquid crystalline phase show the spontaneous formation of rippled bilayers. The ripple consists of two domains of different length and orientation, connected by a kink. The organization of the lipids in one domain of the ripple is found to be that of a splayed gel; in the other domain the lipids are gel-like and fully interdigitated. In the concave part of the kink region between the domains the lipids are disordered. The results are consistent with the experimental information available and provide an atomic-level model that may be tested by further experiments.

HIV-1 protease is most active under weakly acidic conditions (pH 3.5-6.5), when the catalytic Asp25 and Asp25' residues share 1 proton. At neutral pH, this proton is lost and the stability of the structure is reduced. Here we present an investigation of the effect of pH on the dynamics of HIV-1 protease using MD simulation techniques. MD simulations of the solvated HIV-1 protease with the Asp25/25' residues monoprotonated and deprotonated have been performed. In addition we investigated the effect of the inclusion of Na+ and Cl- ions to mimic physiological salt conditions. The simulations of the monoprotonated form and deprotonated form including Na+ show very similar behavior. In both cases the protein remained stable in the compact, "self-blocked" conformation in which the active site is blocked by the tips of the flaps. In the deprotonated system a Na+ ion binds tightly to the catalytic dyad shielding the repulsion between the COO- groups. Ab initio calculations also suggest the geometry of the active site with the Na+ bound closely resembles that of the monoprotonated case. In the simulations of the deprotonated form (without Na+ ions), a water molecule bound between the Asp25 Asp25' side-chains. This disrupted the dimerization interface and eventually led to a fully open conformation. (C) 2004 Wiley-Liss, Inc.

Azurin from Pseudomonas aeruginosa is a small 128-residue, copper-containing protein. Its redox potential can be modified by mutating the protein. Free-energy calculations based on classical molecular-dynamics simulations of the protein and from mutants in aqueous solution at different pH values were used to compute relative redox potentials. The precision of the free-energy calculations with the coupling-parameter approach is evaluated as function of the number and sequence of lambda values, the sampling time and initial conditions. It is found that the precision is critically dependent on the relaxation of hydrogen-bonding networks when changing the atomic-charge distribution due to a change of redox state or pH value. The errors in the free energies range from 1 to 10 k(B)T, depending on the type of process. Only qualitative estimates of the change in redox potential by protein mutation can be obtained.

This article describes the software suite GROMACS (Groningen MAchine for Chemical Simulation) that was developed at the University of Groningen, The Netherlands, in the early 1990s. The software, written in ANSI C, originates from a parallel hardware project, and is well suited for parallelization on processor clusters. By careful optimization of neighbor searching and of inner loop performance, GROMACS is a very fast program for molecular dynamics simulation. It does not have a force field of its own, but is compatible with GROMOS, OPLS, AMBER, and ENCAD force fields. In addition, it can handle polarizable shell models and flexible constraints. The program is versatile, as force routines can be added by the user, tabulated functions can be specified, and analyses can be easily customized. Nonequilibrium dynamics and free energy determinations are incorporated. Interfaces with popular quantum-chemical packages (MOPAC, GAMES-UK, GAUSSIAN) are provided to perform mixed MM/QM simulations. The package includes about 100 utility and analysis programs. GROMACS is in the public domain and distributed (with source code and documentation) under the GNU General Public License. It is maintained by a group of developers from the Universities of Groningen, Uppsala, and Stockholm, and the Max Planck Institute for Polymer Research in Mainz. Its Web site is http://www.gromacs.org.

A novel method for the refinement of misfolded protein structures is proposed in which the properties of the solvent environment are oscillated in order to mimic some aspects of the role of molecular chaperones play in protein folding in vivo. Specifically, the hydrophobicity of the solvent is cycled by repetitively altering the partial charges on solvent molecules (water) during a molecular dynamics simulation. During periods when the hydrophobicity of the solvent is increased, intramolecular hydrogen bonding and secondary structure formation are promoted. During periods of increased solvent polarity, poorly packed regions of secondary structures are destabilized, promoting structural rearrangement. By cycling between these two extremes, the aim is to minimize the formation of long-lived intermediates. The approach has been applied to the refinement of structural models of three proteins generated by using the ROSETTA procedure for ab initio structure prediction. A significant improvement in the deviation of the model structures from the corresponding experimental structures was observed. Although preliminary, the results indicate computationally mimicking some functions of molecular chaperones in molecular dynamics simulations can promote the correct formation of secondary structure and thus be of general use in protein folding simulations and in the refinement of structural models of small- to medium-size proteins.

The incorporation of explicit ions to mimic the effect of ionic strength or to neutralize the overall charge on a system in free energy calculations using molecular dynamics simulations is investigated. The difference in the free energy of hydration between two triosephosphate isomerase inhibitors calculated at five different ion concentrations is used as an example. We show that the free energy difference can be highly sensitive to the presence of explicit ions even in cases where the mutation itself does not involve a change in the overall charge. The effect is most significant if the molecule carries a net charge close to the site mutated. Furthermore, it is shown that the introduction of a small number of ions can lead to very severe sampling problems suggesting that in practical calculations convergence can best be achieved by incorporating either no counterions or by simulating at high ionic strength to ensure sufficient sampling of the ion distribution. (C) 2004 Wiley Periodicals, Inc.

Successive parameterizations of the GROMOS force field have been used successfully to simulate biomolecular systems over a long period of time. The continuing expansion of computational power with time makes it possible to compute ever more properties for an increasing variety of molecular systems with greater precision. This has led to recurrent parameterizations of the GROMOS force field all aimed at achieving better agreement with experimental data. Here we report the results of the latest, extensive reparameterization of the GROMOS force field. In contrast to the parameterization of other biomolecular force fields, this parameterization of the GROMOS force field is based primarily on reproducing the free enthalpies of hydration and apolar solvation for a range of compounds. This approach was chosen because the relative free enthalpy of solvation between polar and apolar environments is a key property in many biomolecular processes of interest, such as protein folding, biomolecular association, membrane formation, and transport over membranes. The newest parameter sets, 53A5 and 53A6, were optimized by first fitting to reproduce the thermodynamic properties of pure liquids of a range of small polar molecules and the solvation free enthalpies of amino acid analogs in cyclohexane (53A5). The partial charges were then adjusted to reproduce the hydration free enthalpies in water (53A6). Both parameter sets are fully documented, and the differences between these and previous parameter sets are discussed.

Here, we use coarse grained molecular dynamics (MD) simulations to study the spontaneous aggregation of dipalmitoylphosphatidylcholine (DPPC) lipids into small unilamellar vesicles. We show that the aggregation process occurs on a nanosecond time scale, with bicelles and cuplike vesicles formed at intermediate stages. Formation of hemifused vesicles is also observed at higher lipid concentration. With either 25% dipalmitoylphosphatidylethanolamine (DPPE) or lysoPC mixed into the system, the final stages of the aggregation process occur significantly faster. The structure of the spontaneously formed vesicles is analyzed in detail. Microsecond simulations of isolated vesicles reveal significant differences in the packing of the lipids between the inner and outer monolayers, and between PC, PE, and lysoPC. Due to the small size of the vesicles they remain almost perfectly spherical, undergoing very limited shape fluctuations or bilayer undulations. The lipid lateral diffusion rate is found to be faster in the outer than in the inner monolayer. The water permeability coefficient of the pure DPPC vesicles is of the order of 10(-3) cm s(-1), in agreement with experimental measurements.

The initial stages of the gating of the mechanoselective channel of large conductance from Mycobacterium tuberculosis have been studied in atomic detail using molecular dynamics simulation techniques. A truncated form of the protein complex embedded in a palmitoyloleoylphosphatidylcholine lipid bilayer and surrounded by explicit water was simulated under different pressure conditions to mimic the effects of tension and compression within the membrane on the protein. As a direct result of lateral tension being applied to the membrane, an increase in the tilt of a subset of the transmembrane helices was observed. This in turn led to the enlargement of the pore and the disruption of the hydrophobic gate consisting of residues Ile-14 and Val-21. The simulations suggest that opening occurs in a sequential staged process. Such a mechanism could explain the partial opening or staged conductance observed in patch-clamp experiments using related large conductance mechanosensitive channel complexes.

Molecular dynamics simulation techniques together with time-dependent density functional theory calculations have been used to investigate the effect of photon absorption by a 4-hydroxy-cinnamic acid chromophore on the structural properties of the photoactive yellow protein (PYP) from Ectothiorodospira halophila. The calculations suggest that the protein not only modifies the absorption spectrum of the chromophore but also regulates the subsequent isomerization of the chromophore by stabilizing the isomerization transition state. Although signaling from PYP is thought to involve partial unfolding of the protein, the mechanical effects accompanying isomerization do not appear to directly destabilize the protein.

The ability of the GROMOS96 force field to reproduce partition constants between water and two less polar solvents (cyclohexane and chloroform) for analogs of 18 of the 20 naturally occurring amino acids has been investigated. The estimations of the solvation free energies in water, in cyclohexane solution, and chloroform solution are based on thermodynamic integration free energy calculations using molecular dynamics simulations. The calculations show that while the force field reproduces the experimental solvation free energies of nonpolar analogs with reasonable accuracy the solvation free energies of polar analogs in water are systematically overestimated (too positive). The dependence of the calculated free energies on the atomic partial charges was also studied. (C) 2002 Wiley Periodicals, Inc.

The dynamics of the three-stranded beta-sheet peptide Betanova has been studied at four different temperatures (280, 300, 350, and 450 K by molecular dynamics simulation techniques, in explicit water. Two 20-ns simulations at 280 K indicate that the peptide remains very flexible under "folding" conditions sampling a range of conformations that together satisfy the nuclear magnetic resonance (NMR)-derived experimental constraints. Two simulations at 300 K (above the experimental folding temperature) of 20 ns each show partial formation of "native"-like structure, which also satisfies most of the NOE constraints at 280 K. At higher temperature, the presence of compact states, in which a series of hydrophobic contacts remain present, are observed. This is consistent with experimental observations regarding the role of hydrophobic contacts in determining the peptide's stability and in initiating the formation of turns and loops. A set of different structures is shown to satisfy NMR-derived distance restraints and a possible mechanism for the folding of the peptide into the NMR-determined structure is proposed.

Extensive molecular dynamics (MD) simulations of binary systems of phospholipids and bile salts, a model for human bile, have been performed. Recent progress in hardware and software development allows simulation of the spontaneous aggregation of the constituents into small mixed micelles, in agreement with experimental observations. The MD simulations reveal the structure of these micelles at atomic detail. The phospholipids are packed radially with their headgroups at the surface and the hydrophobic tails pointing toward the micellar center. The bile salts act as wedges between the phospholipid headgroups, with their hydrophilic sides exposed to the aqueous environment. The structure of the micelles strongly resembles the previously proposed radial shell model. Simulations including small fractions of cholesterol reveal how cholesterol is solubilized inside these mixed micelles without changing their overall structure.

A principal component analysis has been applied on equilibrium simulations of a b-heptapeptide that shows reversible folding in a methanol solution. The analysis shows that the configurational space contains only three dense sub-states. These states of relatively low free energy correspond to the ‘‘native’’ left-handed helix, a partly helical intermediate, and a hairpin-like structure. The collection of unfolded conformations form a relatively diffuse cloud with little substructure. Internal hydrogen-bonding energies were found to correlate well with the degree of folding. The native helical structure folds from the N terminus; the transition from the major folding intermediate to the native helical structure involves the formation of the two most C-terminal backbone hydrogen bonds. A fourstate Markov model was found to describe transition frequencies between the conformational states within error limits, indicating that memoryeffects are negligible beyond the nanosecond time-scale. The dominant native state fluctuations were found to be very similar to unfolding motions, suggesting that unfolding pathways can be inferred from fluctuations in the native state. The low-dimensional essential subspace, describing 69 % of the collective atomic fluctuations, was found to converge at time-scales of the order of one nanosecond at all temperatures investigated, whereas folding/unfolding takes place at significantly longer time-scales, even above the melting temperature. # 2001 Academic Press

Many viral fusion proteins only become activated under mildly acidic condition (pH 4.5–6.5) close to the pKa of histidine side-chain protonation. Analysis of the sequences and structures of influenza HA (haemagglutinin) and flaviviral envelope glycoproteins has led to the identification of a number of histidine residues that are not only fully conserved themselves but have local environments that are also highly conserved [Kampmann, Mueller, Mark, Young and Kobe (2006) Structure 14, 1481–1487]. Here, we summarize studies aimed at determining the role, if any, that protonation of these potential switch histidine residues plays in the low-pH-dependent conformational changes associated with fusion activation of a flaviviral envelope protein. Specifically, we report on MD (Molecular Dynamics) simulations of the DEN2 (dengue virus type 2) envelope protein ectodomain sE (soluble E) performed under varied pH conditions designed to test the histidine switch hypothesis of Kampmann et al. (2006).

Modeling the mechanisms of protofibril twisting: Molecular-dynamics simulations of simian viral peptide aggregates show that β sheets of 10 to 30 chains form left-handed helical ribbons with saddlelike curvature (see picture). These structures are highly dynamic, with oscillations around an average twist angle of 9-10°, and a pitch of 15-20 nm, depending on β-sheet length. The peptides studied are key to viral entry into host cells.

To test and validate the Automated force field Topology Builder and Repository (ATB; http://compbio.biosci.uq.edu.au/atb/ ) the hydration free enthalpies for a set of 214 drug-like molecules, including 47 molecules that form part of the SAMPL4 challenge have been estimated using thermodynamic integration and compared to experiment. The calculations were performed using a fully automated protocol that incorporated a dynamic analysis of the convergence and integration error in the selection of intermediate points. The system has been designed and implemented such that hydration free enthalpies can be obtained without manual intervention following the submission of a molecule to the ATB. The overall average unsigned error (AUE) using ATB 2.0 topologies for the complete set of 214 molecules was 6.7 kJ/mol and for molecules within the SAMPL4 7.5 kJ/mol. The root mean square error (RMSE) was 9.5 and 10.0 kJ/mol respectively. However, for molecules containing functional groups that form part of the main GROMOS force field the AUE was 3.4 kJ/mol and the RMSE was 4.0 kJ/mol. This suggests it will be possible to further refine the parameters provided by the ATB based on hydration free enthalpies.

A series of three molecular dynamics simulations at 300 K in explicit solvent environments of chloroform, methanol and water has been performed on the pulmonary surfactant lipoprotein, SP-C, comprising several consecutive valine residues in order to investigate the stability of the alpha-helical conformation. Two additional simulations were performed on truncated SP-C with a five-residue N-terminal deletion at 300 K and 500 K in water, the high temperature run in order to increase the rate of peptide denaturation. Indications of destabilization appear in chloroform during 1 ns while the SP-C alpha-helix is remarkably stable during 1 ns in methanol and water. In particular the polyvalyl part comprising residues Val15 to Val21 remains intact even at elevated temperature, and the valines do not disrupt the alpha-helical conformation. The valyl-rotamer sampling is partly restricted. Unfolding takes place successively along the primary sequence starting from the C-terminal end. Factors affecting polypeptide stability in molecular dynamics simulations are addressed. The intrinsic helix-forming tendency of valine residues and its dependence on the sequence context, and the role of the solvent environment in stabilizing or destabilizing an alpha-helical fold, are discussed.

Four methods are compared to drive the unfolding of a protein: (1) high tem- perature (T-run), (2) high pressure (P-run), (3) by imposing a gradual increase in the mean ra- dius of the protein using a penalty function added to the physical interaction function (F- run, radial force driven unfolding), and (4) by weak coupling of the difference between the temperature of the radially outward moving at- oms and the radially inward moving atoms to an external temperature bath (K-run, kinetic energy driven unfolding). The characteristic features of the four unfolding pathways are an- alyzed in order to detect distortions due to the size or the type of the applied perturbation, as well as the features that are common to all of them. Hen egg white lysozyme is used as a test system. The simulations are analyzed and com- pared to experimental data like 'H-NMR amide proton exchange-folding competition, heat ca- pacity, and compressibility measurements. 0 1995 Wiley-Liss, Inc.

European Biophysics Journal, 40 (7), ISSN:0175-7571, ISSN:1432-1017

The partition coefficient of cyclosporin A (CsA) was measured in octanol/water and heptane/water by centrifugal partition chromatography. By comparison with results from model compounds, it was deduced that the hydrogen-bonding capacity of CsA changed dramatically from an apolar solvent (where it is internally H-bonded) to polar solvents (where it exposes its H-bonding groups to the solvent). Molecular dynamics simulations in water and CCl4 support the suggestion that CsA undergoes a solvent-dependent conformational changes and that the interconversion process is slow on the molecular dynamics time scale.

The continuous advance of experimental techniques is steadily increasing our knowledge of biochemical systems and processes. A detailed picture of many biomolecular processes is emerging due to the possibility of measuring atomic properties of biological macromolecules, such as proteins. X-ray diffraction has provided a three-dimensional picture of biomolecular assemblies of ever increasing size, ranging from small proteins to a complete macromolecular reaction centre. Other experimental techniques, like nuclear magnetic resonance (NMR) and other spectroscopic methods, yield less complete information at the atomic level. However, the development of multi-dimensional NMR has made it possible to determine the spatial structure of small proteins, albeit at low resolution. The advantage of spectroscopic measuring techniques over X-ray diffraction methods is that the former can be used to obtain information on the dynamics of specific atoms or groups of atoms in the biomolecule, while the latter only yield an indication of the mobility of atoms, not of the time scale of their motion. The other major step forward is the advent of the possibility to change the amino acid composition of a protein at will. Instead of studying the proteins with which nature has provided us, we may now make specific mutations and study their effect on protein properties.

Resorcinolic lipids, or resorcinols, are commonly found in plant membranes. They consist of a substituted benzene ring forming the hydrophilic lipid head, attached to an alkyl chain forming the hydrophobic tail. Experimental results show alternative effects of resorcinols on lipid membranes. Depending on whether they are added to lipid solutions before or after the formation of the liposomes, they either stabilize or destabilize these liposomes. Here we use atomistic molecular dynamics simulations to elucidate the molecular nature of this dual effect. Systems composed of either one of three resorcinol homologs, differing in the alkyl tail length, interacting with dimyristoylphosphatidylcholine lipid bilayers were studied. It is shown that resorcinols preincorporated into bilayers induce order within the lipid acyl chains, decrease the hydration of the lipid headgroups, and make the bilayers less permeable to water. In contrast, simulations in which the resorcinols are incorporated from the aqueous solution into a preformed phospholipid bilayer induce local disruption, leading to either transient pore formation or even complete rupture of the membrane. In line with the experimental data, our simulations thus demonstrate that resorcinols can either disturb or stabilize the membrane structure, and offer a detailed view of the underlying molecular mechanism.

A combination of MD simulations and free energy calculations have been used to propose a new model for the binding of amino acids to threonyl-tRNA-synthetase which not only yields a stable binding mode for l-Ser but also can explain the mechanism by which the editing domains of aminoacyl-tRNA-synthetases are enantiomeric selective preferentially binding d-amino acids.

The purple membrane (PM) is a specialized membrane patch found in halophilic archaea, containing the photoreceptor bacteriorhodopsin (bR). It is long known that calcium ions bind to the PM, but their position and role remain elusive to date. Molecular dynamics simulations in conjunction with a highly detailed model of the PM have been used to investigate the stability of calcium ions placed at three proposed cation binding sites within bR, one near the Schiff base, one in the region of the proton release group, and one near Glu9. The simulations suggest that, of the sites investigated, the binding of calcium ions was most likely at the proton release group. Binding in the region of the Schiff base, while possible, was associated with significant changes in local geometry. Calcium ions placed near Glu9 in the interior of bR (simultaneously to a Ca(2+) near the Schiff base and another one near the Glu194-Glu204 site) were not stable. The results obtained are discussed in relation to recent experimental observations and theoretical considerations.

When estimating binding affinities of a ligand, which can exists in multiple forms, for a target molecule, one must consider all possible competing equilibria. Here, a method is presented that estimates the contribution of the protonation equilibria of a ligand in solution to the measured or calculated binding affinity. The method yields a correction to binding constants that are based on the total concentration of inhibitor (the sum of all ionized forms of the inhibitor in solution) to account for the complexed form of the inhibitor only. The method is applied to the calculation of the difference in binding affinity of two inhibitors, 2-phosphoglycolate (PGA) and its phoshonate analog 3-phosphonopropionate (3PP), for the glycolytic enzyme triosephosphate isomerase. Both inhibitors have three titrating sites and exist in solution as a mixture of different forms. In this case the form that actually binds to the enzyme is present at relative low concentrations. The contributions of the alternative forms to the difference in binding energies is estimated by means of molecular dynamics simulations and corrections. The inhibitors undergo a pK(a) shift upon binding that is estimated by ab initio calculations. An interesting finding is that the affinity difference of the two inhibitors is not due to different interactions in the active site of the enzyme, but rather due to the difference in the solvation properties of the inhibitors.

We have determined the three-dimensional structure of a two-disulfide intermediate (Cys(8)-Cys(20), Cys(14)-Cys(26)) on the oxidative folding pathway of the cyclotide MCoTI-II. Cyclotides have a range of bioactivities and, because of their exceptional stability, have been proposed as potential molecular scaffolds for drug design applications. The three-dimensional structure of the stable two-disulfide intermediate shows for the most part identical secondary and tertiary structure to the native state. The only exception is a flexible loop, which is collapsed onto the protein core in the native state, whereas in the intermediate it is more loosely associated with the remainder of the protein. The results suggest that the native fold of the peptide does not represent the free energy minimum in the absence of the Cys(1)-Cys(18) disulfide bridge and that although there is not a large energy barrier, the peptide must transiently adopt an energetically unfavorable state before the final disulfide can form.

Computer-based molecular simulation techniques are increasingly used to interpret experimental data on biomolecular systems at an atomic level. Direct comparison between experiment and simulation is, however, seldom straightforward. The available experimental data are limited in scope and generally correspond to averages over both time and space. A critical analysis of the various factors that may influence the apparent degree of agreement between the results of simulations and experimentally measured quantities is presented and illustrated using examples from recent literature.

The sensitivity of molecular dynamics simulations to variations in the force field has been examined in relation to a set of 36 structures corresponding to 31 proteins simulated by using different versions of the GROMOS force field. The three parameter sets used (43a1, 53a5, and 53a6) differ significantly in regard to the nonbonded parameters for polar functional groups and their ability to reproduce the correct solvation and partitioning behavior of small molecular analogues of the amino acid side chains. Despite the differences in the force field parameters no major differences could be detected in a wide range of structural properties such as the root-mean-square deviation from the experimental structure, radii of gyration, solvent accessible surface, secondary structure, or hydrogen bond propensities on a 5 to 10 ns time scale. The small differences that were observed correlated primarily with the presence of charged residues as opposed to residues that differed most between the parameter sets. The work highlights the variation that can be observed in nanosecond simulations of protein systems and implications of this for force field validation, as well as for the analysis of protein simulations in general.

Three host-guest systems have been characterized using surface tension (sigma), calorimetry, and molecular dynamics simulations (MD). The hosts were three native cyclodextrins (CD) and the guest the non-ionic carbohydrate surfactant octyl-beta-d-glucopyranoside. It is shown that, for any host-guest system, a rough screening of the most probable complex stoichiometries can be obtained in a model free form, using only calorimetric data. The sigma data were analyzed using a model that includes a newly proposed adsorption isotherm. The equilibrium constants for several stoichiometries were simultaneously obtained through fitting the sigma data. For alpha- and beta-CD, the predominant species is 1:1 and to a lesser extent 2:1, disregarding the existence of the 1:2. For gamma-CD, the 1:2 species dominates, the other two being also present. In an attempt to confirm these results, 10 ns MD simulations for each CD were performed using seven different starting conformations. The MD stable conformations agree with the results found from the experimental data. In one case, the spontaneous dissociation-formation of a complex was observed. Analysis of the trajectories indicates that hydrophobic interactions are primarily responsible for the formation and stability of the inclusion complexes. For the 2:1 species, intermolecular H-bonds between CD molecules result in a tight packed structure where their original truncated cone shape is lost in favor of a cylindrical geometry. Together, the results clearly demonstrate that the often used assumption of considering only a 1:1 species is inappropriate.

Replica exchange methods REMs are increasingly used to improve sampling in molecular dynamics MD simulations of biomolecular systems. However, despite having been shown to be very effective on model systems, the application of REM in complex systems such as for the simulation of protein and peptide folding in explicit solvent has not been objectively tested in detail. Here we present a comparison of conventional MD and temperature replica exchange MD T-REMD simulations of a -heptapeptide in explicit solvent. This system has previously been shown to undergo reversible folding on the time scales accessible to MD simulation and thus allows a direct one-to-one comparison of efficiency. The primary properties compared are the free energy of folding and the relative populations of different conformers as a function of temperature. It is found that to achieve a similar degree of precision T-REMD simulations starting from a random set of initial configurations were approximately an order of magnitude more computationally efficient than a single 800 ns conventional MD simulation for this system at the lowest temperature investigated 275 K. However, whereas it was found that T-REMD simulations are more than four times more efficient than multiple independent MD simulations at one temperature 300 K the actual increase in conformation sampling was only twofold. The overall gain in efficiency using REMD resulted primarily from the ordering of different conformational states over temperature, as opposed to a large increase of conformational sampling. It is also shown that in this system exchanges are accepted primarily based on random fluctuations within the solvent and are not strongly correlated with the instantaneous peptide conformation raising questions in regard to the efficiency of T-REMD in larger systems. © 2007 American Institute of Physics. DOI: 10.1063/1.2404954

We report here an advanced approach for the characterization of the folding pattern of a de novo designed antiparallel coiled coil peptide by high-resolution methods. Incorporation of two fluorescence labels at the C- and N-terminus of the peptide chain as well as modi. cation of two hydrophobic core positions by Phe/[N-15,C-13]Leu enable the study of the folding characteristics and of distinct amino acid side chain interactions by fluorescence resonance energy transfer (FRET) and NMR spectroscopy. Results of both experiments reveal the antiparallel alignment of the helices and thus prove the design concept. This finding is also supported by molecular dynamics simulations. Electrospray ionization Fourier transform ion cyclotron resonance mass spectrometry (ESI-FTICR-MS) in combination with NMR experiments was used for veri. cation of the oligomerization equilibria of the coiled coil peptide.

Molecular dynamics simulations have been used to study the spontaneous aggregation of a concentrated solution of dipalmitoylphosphatidylcholine (DPPC) molecules in water into a small vesicle. The molecules were represented in atomistic detail. Starting from a DPPC solution in water, an oblong vesicle with a long axis of 15 nm and short axes of 10 nm was formed spontaneously. After 90 ns of simulation, the vesicle contained a number of water pores. Water pores were shown to facilitate exchange of lipids between inner and outer leaflets. Lipid tails were shown to be less ordered in the inner leaflet of the vesicle, as compared to those in the outer leaflet of the vesicle and an equilibrated lamellar bilayer.