Molecular dynamics simulations and free energy calculations on the enzyme 4-hydroxyphenylpyruvate dioxygenase.

4-Hydroxyphenylpyruvate dioxygenase is a relevant target in both pharmaceutical and agricultural research. We report on molecular dynamics simulations and free energy calculations on this enzyme, in complex with 12 inhibitors for which experimental affinities were determined. We applied the thermodynamic integration approach and the more efficient one-step perturbation. Even though simulations seem well converged and both methods show excellent agreement between them, the correlation with the experimental values remains poor. We investigate the effect of slight modifications on the charge distribution of these highly conjugated systems and find that accurate models can be obtained when using improved force field parameters. This study gives insight into the applicability of free energy methods and current limitations in force field parameterization.

Membrane protein dynamics in different environments: simulation study of the outer membrane protein X in a lipid bilayer and in a micelle.

The bacterial outer membrane protein OmpX from Escherichia coli has been investigated by molecular dynamics simulations when embedded in a phospholipid bilayer and as a protein-micelle aggregate. The resulting simulation trajectories were analysed in terms of structural and dynamic properties of the membrane protein. In agreement with experimental observations, highest relative stability was found for the ß-barrel region that is embedded in the lipophilic phase, whereas an extracellular protruding ß-sheet, which is a unique structural feature of OmpX that supposedly plays an important role in cell adhesion and invasion, shows larger structure fluctuations. Additionally, we investigated water permeation into the core of the ß-barrel protein, which contains a tight salt-bridge and hydrogen-bond network, so that extensive water flux is unlikely. Differences between the bilayer and the micellar system were observed in the length of the barrel and its position inside the lipid environment, and in the protein interactions with the hydrophilic part of the lipids near the lipid/water interface. Those variations suggest that micelles and other detergent environments might not offer a wholly membrane-like milieu to promote adoption of the physiological conformational state by OmpX.

Simulation of beta-depsipeptides: the effect of missing hydrogen-bond donors on their folding equilibria.

beta-Depsipeptides are beta-peptides in which one or more peptide linkages are replaced by ester linkages, resulting in a loss of a hydrogen-bond donor (N--H) and weakening of the corresponding carbonyl hydrogen-bond acceptor moiety. The effects of three of such peptide by ester substitutions in a hepta-beta-peptide upon its (un)folding equilibrium in methanol solution are investigated using molecular dynamics simulations and compared to experimental data from NMR spectroscopy. The simulated conformational ensembles largely reproduce the experimentally measured NOE and 3J-coupling constant data for the three different hepta-beta-peptides, and confirm the relative stabilities of the 3(14)-helical conformation, which is most weakened by substitution of the 4th peptide linkage and least by substitution of the 6th peptide linkage. The simulations are complementary to the experimental data by providing detailed insight into the conformational distributions that are compatible with the experimentally measured average values of observables.

Biomolecular modeling: Goals, problems, perspectives.

Computation based on molecular models is playing an increasingly important role in biology, biological chemistry, and biophysics. Since only a very limited number of properties of biomolecular systems is actually accessible to measurement by experimental means, computer simulation can complement experiment by providing not only averages, but also distributions and time series of any definable quantity, for example, conformational distributions or interactions between parts of systems. Present day biomolecular modeling is limited in its application by four main problems: 1) the force-field problem, 2) the search (sampling) problem, 3) the ensemble (sampling) problem, and 4) the experimental problem. These four problems are discussed and illustrated by practical examples. Perspectives are also outlined for pushing forward the limitations of biomolecular modeling.

A molecular dynamics study of the bee venom melittin in aqueous solution, in methanol, and inserted in a phospholipid bilayer.

The structural properties of melittin, a small amphipathic peptide found in the bee venom, are investigated in three different environments by molecular dynamics simulation. Long simulations have been performed for monomeric melittin solvated in water, in methanol, and shorter ones for melittin inserted in a dimyristoylphosphatidylcholine bilayer. The resulting trajectories were analysed in terms of structural properties of the peptide and compared to the available NMR data. While in water and methanol solution melittin is observed to partly unfold, the peptide retains its structure when embedded in a lipid bilayer. The latter simulation shows good agreement with the experimentally derived (3)J-coupling constants. Generally, it appears that higher the stability of the helical conformation of melittin, lower is the dielectric permittivity of the environment. In addition, peptide-lipid interactions were investigated showing that the C-terminus of the peptide provides an anchor to the lipid bilayer by forming hydrogen bonds with the lipid head groups.

Helices and other secondary structures of beta- and gamma-peptides.

The principal secondary structural motifs adopted by peptides assembled from beta-amino acid units are discussed: the 14-, 12-, 10-, 12/10-, and 8-helices, as well as the hairpin turn, extended structures, stacks, and sheets. Features that promote a particular folding propensity are outlined and illustrated by structures determined in solution (NMR) and in the solid-state (x-ray). The N-C(beta)-C(alpha)-CO dihedral angles from molecular dynamics simulations, which are indicative of a particular secondary structure, are presented. A brief description of a helix and a turn of gamma-peptides is also given.

Interpreting NMR data for beta-peptides using molecular dynamics simulations.

NMR is one of the most used techniques to resolve structure of proteins and peptides in solution. However, inconsistencies may occur due to the fact that a polypeptide may adopt more than one conformation. Since the NOE distance bounds and (3)J-values used in such structure determination represent a nonlinear average over the total ensemble of conformers, imposition of NOE or (3)J-value restraints to obtain one unique conformation is not an appropriate procedure in such cases. Here, we show that unrestrained MD simulation of a solute in solution using a high-quality force field yields a conformational ensemble that is largely compatible with the experimental NMR data on the solute. Four 100 ns MD simulations of two forms of a nine-residue beta-peptide in methanol at two temperatures produced conformational ensembles that were used to interpret the NMR data on this molecule and resolve inconsistencies between the experimental NOEs. The protected and unprotected forms of the beta-peptide adopt predominantly a 12/10-helix in agreement with the qualitative interpretation of the NMR data. However, a particular NOE was not compatible with this helix indicating the presence of other conformations. The simulations showed that 3(14)()-helical structures were present in the ensemble of the unprotected form and that their presence correlates with the fulfillment of the particular NOE. Additionally, all inter-hydrogen distances were calculated to compare NOEs predicted by the simulations to the ones observed experimentally. The MD conformational ensembles allowed for a detailed and consistent interpretation of the experimental data and showed the small but specific conformational differences between the protected and unprotected forms of the peptide.

On the influence of charged side chains on the folding-unfolding equilibrium of beta-peptides: a molecular dynamics simulation study.

The influence of charged side chains on the folding-unfolding equilibrium of beta-peptides was investigated by means of molecular dynamics simulations. Four different peptides containing only negatively charged side chains, positively charged side chains, both types of charged side chains (with the ability to form stabilizing salt bridges) or no charged side chains were studied under various conditions (different simulation temperatures, starting structures and solvent environment). The NMR solution structure in methanol of one of the peptides (A) has already been published; the synthesis and NMR analysis of another peptide (B) is described here. The other peptides (C and D) studied herein have hitherto not been synthesized. All four peptides A-D are expected to adopt a left-handed 3(14)-helix in solution as well as in the simulations. The resulting ensembles of structures were analyzed in terms of conformational space sampled by the peptides, folding behavior, structural properties such as hydrogen bonding, side chain-side chain and side chain-backbone interactions and in terms of the level of agreement with the NMR data available for two of the peptides. It was found that the presence of charged side chains significantly slows down the folding process in methanol solution due to the stabilization of intermediate conformers with side chain-backbone interactions. In water, where the solvent competes with the solute-solute polar interactions, the folding process to the 3(14)-helix is faster in the simulations.

A novel approach for designing simple point charge models for liquid water with three interaction sites.

A simultaneous improvement of the diffusion and dielectric properties of the simple point charge (SPC) model for liquid water appears to be very difficult with conventional reparametrization of the commonly used Lennard-Jones and Coulomb interaction functions and without including a self-energy correction in the effective pair-potential as is done in the SPC/E model. Here, a different approach to circumvent this problem is presented. A short-range interaction term, which corrects the oxygen-oxygen energy at small distances by small amounts of energy, was introduced in the nonbonded interaction function. This additional force-field term allows to derive new parameter sets for SPC-like water models that yield better agreement with experimental data on liquid water. Based on previous investigations of the force-field parameter dependence of the water properties of SPC-like models, the necessary parameter changes to obtain a lower diffusion coefficient and a larger dielectric permittivity were specified and accordingly six new models were developed. They all represent an improvement over SPC in terms of structural and diffusional properties, four of them show better dielectric properties also. One model, SPC/S, has been characterized in more detail, and represents most properties of liquid water better than SPC while avoiding the larger discrepancies with experimental values regarding density, thermal compressibility, energy, and free energy of the SPC/E model. We conclude that the use of a simple, short-ranged additional oxygen-oxygen interaction term makes a simultaneous improvement of the diffusion coefficient and the dielectric properties of water feasible.

Can one derive the conformational preference of a beta-peptide from its CD spectrum?

CD spectroscopy is often used to elucidate the secondary structure of peptides built from non-natural amino acids such as beta-amino acids. The interpretation of such CD spectra is not always unambiguous. Here, we present a case where two beta-hexapeptides, a dimethyl-beta-hexapeptide indicated as DM-BHP (A) and its nonmethylated analogue indicated as BHP (B), exhibit similar CD spectra, whereas they are expected to differ in secondary structure. The structural properties of both peptides were studied by molecular dynamics simulation, and from the resulting trajectories, the corresponding CD spectra were calculated. Starting from a fully extended conformation, BHP is observed to form a 3(14)-helix, while DM-BHP remains unfolded. However, even though these two peptides hardly share any conformations, their calculated CD spectra are alike and show the same features as the experimentally measured ones. Our results imply that a particular CD pattern can be induced by spatially different structures, which makes it difficult to derive the conformational preference of a peptide from its CD spectrum alone. To gain more insight into the relationship between the preferred conformation of a peptide and its CD spectrum, more accurate methods to calculate the CD spectrum for a given conformation are required.