Entropy calculations on the molten globule state of a protein: side-chain entropies of alpha-lactalbumin.

We present entropy estimates based on molecular dynamics simulations of models of the molten globule state of the protein alpha-lactalbumin at low pH. The entropy calculations use the covariance matrix of atom-positional fluctuations and yield the complete configurational entropy. The configurational entropy of the entire protein and of each of its side chains is calculated.

Fibrillation power, an alternative method of ECG spectral analysis for prediction of countershock success in a porcine model of ventricular fibrillation.

Background: Noninvasive prediction of defibrillation success after cardiac arrest and cardiopulmonary resuscitation (CPR) may help in determining the optimal time for a countershock, and thus increase the chance for survival.

Methods: In a porcine model (n=25) of prolonged cardiac arrest, advanced cardiac life support was provided by administration of two or three doses of either vasopressin or epinephrine after 3 or 8 min of basic life support. After 4 min of ventricular fibrillation and 18 min of life support, defibrillation was attempted.

Molecular dynamics simulation of n-dodecyl phosphate aggregate structures.

Aggregates of n-dodecyl phosphate present an attractive model system of simple phospholipid amphiphile supramolecular structures for study by molecular dynamics simulation, since these systems have previously been studied experimentally under various conditions. A detailed molecular dynamics description of the properties of planar bilayer membranes (as a model for unilamellar vesicular membranes) and spherical micelles under various simulated conditions is presented. It is shown that the united-atom model of GROMOS96 applying the force-field parameter set 43A2 for biomolecular systems yields properties in agreement with experimental ones in most cases.

Calculation of NMR-relaxation parameters for flexible molecules from molecular dynamics simulations.

Comparatively small molecules such as peptides can show a high internal mobility with transitions between several conformational minima and sometimes coupling between rotational and internal degrees of freedom. In those cases the interpretation of NMR relaxation data is difficult and the use of standard methods for structure determination is questionable. On the other hand, in the case of those system sizes, the timescale of both rotational and internal motions is accessible by molecular dynamics (MD) simulations using explicit solvent.

The beta-peptide hairpin in solution: conformational study of a beta-hexapeptide in methanol by NMR spectroscopy and MD simulation.

The structural and thermodynamic properties of a 6-residue beta-peptide that was designed to form a hairpin conformation have been studied by NMR spectroscopy and MD simulation in methanol solution. The predicted hairpin would be characterized by a 10-membered hydrogen-bonded turn involving residues 3 and 4, and two extended antiparallel strands. The interproton distances and backbone torsional dihedral angles derived from the NMR experiments at room temperature are in general terms compatible with the hairpin conformation.

Dielectric properties of proteins from simulation: the effects of solvent, ligands, pH, and temperature.

We have used a standard Fröhlich-Kirkwood dipole moment fluctuation model to calculate the static dielectric permittivity, epsilon(0), for four different proteins, each of which was simulated under at least two different conditions of pH, temperature, solvation, or ligand binding. For the range of proteins and conditions studied, we calculate values for epsilon(0) between 15 and 40. Our results show, in agreement with prior work, that the behavior of charged residues is the primary determinant of the effective permittivity.

Assessing the effect of conformational averaging on the measured values of observables.

Experiment and computer simulation are two complementary tools to understand the dynamics and behavior of biopolymers in solution. One particular area of interest is the ensemble of conformations populated by a particular molecule in solution. For example, what fraction of a protein sample exists in its folded conformation? How often does a particular peptide form an alpha helix versus a beta hairpin? To address these questions, it is important to determine the sensitivity of a particular experiment to changes in the distribution of molecular conformations.

Comparison of different schemes to treat long-range electrostatic interactions in molecular dynamics simulations of a protein crystal.

Eight molecular dynamics simulations of a ubiquitin crystal unit cell were performed to investigate the effect of different schemes to treat the long-range electrostatic interactions as well as the need to include counter ions. A crystal system was chosen as the test system, because the higher charge density compared with a protein in solution makes it more sensitive to the way of treating the electrostatic interactions. Three different schemes of treating the long-range interactions were compared: straight cutoff, reaction-field approximation, and a lattice-sum method (P3M).

Free energy barrier estimation of unfolding the alpha-helical surfactant-associated polypeptide C.

Molecular dynamics simulations were conducted to estimate the free energy barrier of unfolding surfactant-associated polypeptide C (SP-C) from an alpha-helical conformation. Experimental studies indicate that while the helical fold of SP-C is thermodynamically stable in phospholipid micelles, it is metastable in a mixed organic solvent of CHCl3/CH3OH/0.1 M HCl at 32:64:5 (v/v/v), in which it undergoes an irreversible transformation to an insoluble aggregate that contains beta-sheet.

Entropy calculations on a reversibly folding peptide: changes in solute free energy cannot explain folding behavior.

The configurational entropy of a beta-heptapeptide in solution at four different temperatures is calculated. The contributions of the backbone and of the side-chain atoms to the total peptide entropy are analyzed separately and the effective contribution to the entropy arising from correlations between these terms determined. The correlation between the backbone and side-chain atoms amounts to about 17% and is rather insensitive to the temperature.

Folding study of an Aib-rich peptide in DMSO by molecular dynamics simulations.

To evaluate the ability of molecular dynamics (MD) simulations using atomic force-fields to correctly predict stable folded conformations of a peptide in solution, we show results from MD simulations of the reversible folding of an octapeptide rich in alpha-aminoisobutyric acid (2-amino-2-methyl-propanoic acid, Aib) solvated in di-methyl-sulfoxide (DMSO). This solvent generally prevents the formation of secondary structure, whereas Aib-rich peptides show a high propensity to form secondary structural elements, in particular 3(10)- and alpha-helical structures. Aib is, moreover, achiral, so that Aib-rich peptides can form left- or right-handed helices depending on the overall composition of the peptide, the temperature, and the solvation conditions.

Conformational study of an Aib-rich peptide in DMSO by NMR.

The strong propensity of 2-amino-2-methyl propanoic acid (Aib)-rich peptides to form stable helical structures is well documented. NMR analysis of the short peptide Z-(Aib)5-L-Leu-(Aib)2-OMe indicates the presence of a well-characterized 3(10)-helix even in dimethylsulfoxide (DMSO), a solvent known to disrupt helical structures. The structure remains stable at least up to 348 K.

The Key to Solving the Protein-Folding Problem Lies in an Accurate Description of the Denatured State.

Accurate simulation at the atomic level of the folding process of a variety of peptides into different native folds can be achieved with a general purpose force field and Newton's equations of motion. The key to understanding this peptide folding lies in the unexpectedly small size of the denatured state and an accurate description thereof.

Viscosity dependence of protein dynamics.

The influence of solvent viscosity on protein dynamics was investigated with molecular dynamics simulations of factor Xa in two solvents differing only in viscosity, by a factor of 10. We obtained this viscosity change by changing the masses of the solvent atoms by a factor of 100. Equilibrium properties of the protein, that is, the average structure, its fluctuations, and the secondary structure, show no significant dependence on the solvent viscosity.

On the similarity of properties in solution or in the crystalline state: a molecular dynamics study of hen lysozyme.

As protein crystals generally possess a high water content, it is assumed that the behaviour of a protein in solution and in crystal environment is very similar. This assumption can be investigated by molecular dynamics (MD) simulation of proteins in the different environments. Two 2ns simulations of hen egg white lysozyme (HEWL) in crystal and solution environment are compared to one another and to experimental data derived from both X-ray and NMR experiments, such as crystallographic B-factors, NOE atom-atom distance bounds, 3J(H N alpha)-coupling constants, and 1H-15N bond vector order parameters.

Molecular-dynamics simulation of the beta domain of metallothionein with a semi-empirical treatment of the metal core.

The three-metal-containing beta domain of rat liver metallothionein-2 in aqueous solution was simulated with different metal contents. The Cd(3), the CdZn(2), and the Zn(3) variant were investigated using a conventional molecular dynamics simulation, as well as a simulation with a semi-empirical quantum-chemical description (MNDO and MNDO/d) of the metal core embedded in a classical environment. For the purely classical simulations, the standard GROMOS96 force-field parameters were used, and parameters were estimated for cadmium.

Factor Xa: simulation studies with an eye to inhibitor design.

Factor Xa is a serine protease which activates thrombin and plays a key regulatory role in the blood-coagulation cascade. Factor Xa is at the crossroads of the extrinsic and intrinsic pathways of coagulation and, hence, has become an important target for the design of anti-thrombotics (inhibitors). It is not known to be involved in other processes than hemostasis and its binding site is different to that of other serine proteases, thus facilitating selective inhibition.

Molecular dynamics simulations highlight mobile regions in proteins: A novel suggestion for converting a murine V(H) domain into a more tractable species.

The V(H) region of the murine antibody 1F7 has been identified as a single-domain chorismate mutase, but a tendency to denature and aggregate has hampered its biochemical characterization. Standard mutagenesis approaches targeting antibody chain dimerization areas have been exhausted. We describe a new approach to the problem, where we use molecular dynamics (MD) simulations to find the differences between the untractable protein and the known soluble V(H) domain from a llama antibody.

On the temperature and pressure dependence of a range of properties of a type of water model commonly used in high-temperature protein unfolding simulations.

Molecular dynamics simulations of protein folding and unfolding are often carried out at temperatures (400-600 K) that are much higher than physiological or room temperature to speed up the (un)folding process. Use of such high temperatures changes both the protein and solvent properties considerably, compared to physiological or room temperature. Water models designed for use in conjunction with biomolecules, such as the simple point charge (SPC) model, have generally been calibrated at room temperature and pressure.

Molecular dynamics simulation of hen egg white lysozyme: a test of the GROMOS96 force field against nuclear magnetic resonance data.

Biomolecular force fields for use in molecular dynamics (MD) simulations of proteins, DNA, or membranes are generally parametrized against ab initio quantum-chemical and experimental data for small molecules. The application of a force field in a simulation of a biomolecular system, such as a protein in solution, may then serve as a test of the quality and transferability of the force field. Here, we compare various properties obtained from two MD simulations of the protein hen egg white lysozyme (HEWL) in aqueous solution using the latest version, GROMOS96, of the GROMOS force field and an earlier version, GROMOS87+, with data derived from nuclear magnetic resonance (NMR) experiments: NOE atom-atom distance bounds, (3)J(HNalpha)-coupling constants, and backbone and side-chain order parameters.

beta-hairpin stability and folding: molecular dynamics studies of the first beta-hairpin of tendamistat.

The stability and (un)folding of the 19-residue peptide, SCVTLYQSWRYSQADNGCA, corresponding to the first beta-hairpin (residues 10 to 28) of the alpha-amylase inhibitor tendamistat (PDB entry 3AIT) has been studied by molecular dynamics simulations in explicit water under periodic boundary conditions at several temperatures (300 K, 360 K and 400 K), starting from various conformations for simulation lengths, ranging from 10 to 30 ns. Comparison of trajectories of the reduced and oxidized native peptides reveals the importance of the disulphide bridge closing the beta-hairpin in maintaining a proper turn conformation, thereby insuring a proper side-chain arrangement of the conserved turn residues. This allows rationalization of the conservation of those cysteine residues among the family of alpha-amylase inhibitors.

The effect of motional averaging on the calculation of NMR-derived structural properties.

The effect of motional averaging when relating structural properties inferred from nuclear magnetic resonance (NMR) experiments to molecular dynamics simulations of peptides is considered. In particular, the effect of changing populations of conformations, the extent of sampling, and the sampling frequency on the estimation of nuclear Overhauser effect (NOE) inter-proton distances, vicinal (3)J-coupling constants, and chemical shifts are investigated. The analysis is based on 50-ns simulations of a beta-heptapeptide in methanol at 298 K, 340 K, 350 K, and 360 K.

Accessibility and order of water sites in and around proteins: A crystallographic time-averaging study.

Water plays an essential role in most biological processes. Water molecules solvating biomolecules are generally in fast exchange with the environment. Nevertheless, well-defined electron density is seen for water associated with proteins whose crystal structure is determined to high resolution.