Biomedical subjects

W F van Gunsteren

Publications and source records attributed to W F van Gunsteren.

At least 19 recordsLinked to original sources

PubMed1996

A combined quantum/classical molecular dynamics study of the catalytic mechanism of HIV protease.

Based on available three-dimensional structures of enzyme-inhibitor complexes, the mechanism of the reaction catalysed by HIV protease is studied using molecular dynamics simulations with molecular mechanics and combined quantum-mechanics/molecular-mechanics potential energy functions. The results support the general acid/general base catalysis mechanism, with Asp25' protonated in the enzyme-substrate complex. In the enzyme-substrate complex, the lytic water molecule binds at a position different from the positions of the hydroxyl groups in various aspartic protease-inhibitor complexes. The carboxyl groups at the active site also adopt a different orientation. However, when the lytic water molecule approaches the scissile peptide, the reaction centre changes gradually to a conformation close to that derived from X-ray diffraction studies of various enzyme-inhibitor complexes. The proton transfer processes can take place only after the lytic water molecule has approached the scissile peptide bond to a certain degree. Qualitatively, the free-energy barrier associated with the nucleophilic attack step, which takes place at physiological pH, is comparable with the acid or base-catalysed reactions of model systems. The structure of the tetrahedral intermediate resulting from the nucleophilic attack step also indicates a straightforward pathway of the next reaction step, i.e. the breaking of the C-N bond.

Computer Simulation↗

PubMed1996

Structural stability of disulfide mutants of basic pancreatic trypsin inhibitor: a molecular dynamics study.

The structure and folding of basic pancreatic trypsin inhibitor (BPTI) has been studied extensively by experimental means. We report a computer simulation study of the structural stability of various disulfide mutants of BPTI, involving eight 250-psec molecular dynamics simulations of the proteins in water, with and without a phosphate counterion. The presence of the latter alters the relative stability of the single disulfide species [5-55] and [30-51]. This conclusion can explain results of mutational studies and the conservation of residues in homologues of BPTI, and suggests a possible role of ions in stabilizing one intermediate over another in unfolding or folding processes.

Aprotinin↗

PubMed1995

Exploring the role of the solvent in the denaturation of a protein: a molecular dynamics study of the DNA binding domain of the 434 repressor.

Molecular dynamics simulations of the DNA binding domain of 434 repressor are presented which aim at unraveling the role of solvent in protein denaturation. Four altered solvent models, each mimicking various possible aspects of the addition of a denaturant to the aqueous solvent, were used in the simulations to analyze their effects on the stability of the protein. The solvent was altered by selectively changing the Coulombic interaction between water and protein atoms and between different water molecules. The use of a modified solvent model has the advantage of mimicking the presence of denaturant without having denaturant molecules present in the simulation, which would require much longer simulations. In these simulations, only an increase in the solvent-protein Coulombic interaction causes initiation of protein unfolding in a manner consistent with NMR data. The altered solvent thus provides a model of a denaturing environment for studying protein unfolding.

Amino Acid Sequence↗

PubMed1995

Fluctuation and cross-correlation analysis of protein motions observed in nanosecond molecular dynamics simulations.

Nanosecond molecular dynamics simulations of bovine pancreatic trypsin inhibitor and lysozyme in water are analyzed in terms of backbone atomic positional fluctuations and dynamical cross-correlations. It is found that although the molecular systems are stable, B-factors calculated over a time period as long as 500 ps are not representative for the motions within the proteins. This is especially true for the most mobile residues. On a nanosecond time-scale, the B-factors calculated from the simulations of the proteins in solution are considerably larger than those obtained by structure refinement of the proteins in crystals, based on X-ray data. The time evolution of the atomic fluctuations shows that for large portions of the proteins under study, atomic positional fluctuations are not yet converged after a nanosecond. Cross-correlations do not converge faster than the fluctuations themselves. Most display very erratic behavior if the sampling covers less than about 200 ps. It is also shown that inclusion of mobile atoms into the procedure used to remove rigid-body motion from the simulation can lead to spurious correlations between the motions of the atoms at the surface of the protein.

Animals↗

PubMed1995

Comparison of MD simulations and NMR experiments for hen lysozyme. Analysis of local fluctuations, cooperative motions, and global changes.

Three 1000 ps molecular dynamics simulations of hen lysozyme have been compared with a range of experimental NMR parameters in order to gain insight into the dynamical properties of the protein and to assess the significance of the motional events observed in the simulations. The simulations, one in vacuum and two in water, were used to estimate interproton distances (for comparison with NOE data), 3JHN alpha and 3J alpha beta coupling constants and 1H-15N order parameters. Comparison of these values with experimental data, particularly NOEs, enabled force field-induced changes to the structure during the simulations to be recognized. It has been shown, however, that these changes can be largely eliminated by slight modifications to the force field. Using a simulation performed in water with this modified force field, it has been found that 1H-15N order parameters calculated for side chain groups in particular correlate well with experimental values and reflect the substantial dependence of these motional properties on the environment, particularly surface exposure, in which the side chain is found. In this case, the simulation then provides models for the motional processes giving rise to the observed experimental data. The results indicate that the order parameter values reflect primarily the number of torsion angles about which rotameric interchange occurs. In addition to local motions, the two different domains of lysozyme have been found to behave differently in the simulations. Possible implications of these differences for the interpretation of unfolding simulations and experimental observations of folding intermediates for lysozyme are discussed.

Animals↗

PubMed1995

The effect of environment on the stability of an integral membrane helix: molecular dynamics simulations of surfactant protein C in chloroform, methanol and water.

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.

Amino Acid Sequence↗

PubMed1995

Internal mobility of the basic pancreatic trypsin inhibitor in solution: a comparison of NMR spin relaxation measurements and molecular dynamics simulations.

Order parameters as well as longitudinal and transverse relaxation rates are calculated for the backbone 15N and 13C alpha nuclei of the basic pancreatic trypsin inhibitor (BPTI) from a 1000 ps molecular dynamics trajectory in explicit water at 277 K using the "model free" approach of Lipari and Szabo. New NMR relaxation data at 277 K are presented, and a comparison is made between NMR relaxation measurements and molecular dynamics relaxation data. It is found that the relaxation processes determining the longitudinal (T1) relaxation rates are inadequately sampled even during this length of simulation. In effect, the calculated relaxation rates are determined almost solely by the order parameters and the overall rotational correlation time of the protein, which appears to be in clear contrast to experimental relaxation rates.

Amino Acid Sequence↗

PubMed1995

Computational approaches to study protein unfolding: hen egg white lysozyme as a case study.

Four methods are compared to drive the unfolding of a protein: (1) high temperature (T-run), (2) high pressure (P-run), (3) by imposing a gradual increase in the mean radius 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 atoms 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 analyzed 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 compared to experimental data like 1H-NMR amide proton exchange-folding competition, heat capacity, and compressibility measurements.

Amides↗

PubMed1995

Structure and internal dynamics of the bovine pancreatic trypsin inhibitor in aqueous solution from long-time molecular dynamics simulations.

Structural and dynamic properties of bovine pancreatic trypsin inhibitor (BPTI) in aqueous solution are investigated using two molecular dynamics (MD) simulations: one of 1.4 ns length and one of 0.8 ns length in which atom-atom distance bounds derived from NMR spectroscopy are included in the potential energy function to make the trajectory satisfy these experimental data more closely. The simulated properties of BPTI are compared with crystal and solution structures of BPTI, and found to be in agreement with the available experimental data. The best agreement with experiment was obtained when atom-atom distance restraints were applied in a time-averaged manner in the simulation. The polypeptide segments found to be most flexible in the MD simulations coincide closely with those showing differences between the crystal and solution structures of BPTI.

Animals↗

PubMed1995

Structure refinement with molecular dynamics and a Boltzmann-weighted ensemble.

Time-averaging restraints in molecular dynamics simulations were introduced to account for the averaging implicit in spectroscopic data. Space- or molecule-averaging restraints have been used to overcome the fact that not all molecular conformations can be visited during the finite time of a simulation of a single molecule. In this work we address the issue of using the correct Boltzmann weighting for each member of an ensemble, both in time and in space. It is shown that the molecular- or space-averaging method is simple in theory, but requires a priori knowledge of the behaviour of a system. This is illustrated using a five-atom model system and the small cycle peptide analogue somatostatin. When different molecular conformers that are separated by energy barriers insurmountable on the time scale of a simulation contribute significantly to a measured NOE intensity, the use of space- or molecule-averaged distance restraints yields a more appropriate description of the measured data than conventional single-molecule refinement with or without application of time averaging.

Algorithms↗

PubMed1995

Rapid non-empirical approaches for estimating relative binding free energies.

Rapid non-empirical methods for estimating binding free energies are reviewed. A novel approach based on the application of the free energy perturbation formula to a biased ensemble is presented. Preliminary results demonstrating the applicability of this approach in protein systems are shown and the potential of this method in structure-based drug design is discussed.

Models, Theoretical↗

PubMed1994

A reassessment of the structure of chymotrypsin inhibitor 2 (CI-2) using time-averaged NMR restraints.

Chymotrypsin inhibitor 2 (CI-2) is one of the growing family of proteins for which well-defined solution and crystal structures have been published and for which small, but distinct differences between these were found. It presents an ideal case to address the question of whether a structural difference is physically real or due to the simplifying approximations with respect to averaging that are used in the conventional methods for structure refinement. NOE distance and 3J coupling constant restrained molecular dynamics simulations were performed using conventional and time-averaged restraints, both in vacuo and in aqueous solution, and the trajectories were compared with structural properties of published structures. The time-averaged restrained molecular dynamics simulations sampled more conformations at various times and visited states consistent with both previously published solution and crystal structures. It was found that the difference between these structures is due to the refinement methodology used. Application of time-averaged restraints in structure refinement yields a physically different picture of the molecular mobility.

Computer Simulation↗

PubMed1994

Simultaneous refinement of the structure of BPTI against NMR data measured in solution and X-ray diffraction data measured in single crystals.

The structure of the bovine pancreatic trypsin inhibitor (BPTI) has been determined to high resolution by both NMR spectroscopy in solution and X-ray diffraction in crystals. The root-mean-square difference calculated between the two structures for the polypeptide backbone is 0.9 A. Several amino acid side-chains, of which all but one are charged or polar, have different conformations. We find that by refining one structure simultaneously against both the NMR and crystallographic data sets, it can accommodate both. Different starting configurations were used, including the X-ray structure 5pti, an NMR conformer, and the X-ray structure in the full unit cell with extra solvent placed in the bulk solvent region. The X-ray structures quickly converged to accommodate the NMR data in addition to the crystallographic data. Starting from an NMR conformer, however, the convergence to accommodate the more abundant X-ray data in addition to the NMR data is much slower.

Amino Acid Sequence↗

PubMed1994

Decomposition of the free energy of a system in terms of specific interactions. Implications for theoretical and experimental studies.

Recently, a number of methods have been proposed that are designed to extract contributions to the change in free energy associated with a given perturbation or mutation of a protein originating from specific residue-residue or atom-atom interactions, both based on theoretical calculations and on experimental data. We caution here that detailed analysis based on these methods is unreliable. It is demonstrated, both from first principles using statistical mechanics and by way of example, that in a general case a meaningful decomposition of the free energy in terms of specific residue-residue or atom-atom interactions is not possible.

Azurin↗

PubMed1994

The high-resolution structure of the histidine-containing phosphocarrier protein HPr from Escherichia coli determined by restrained molecular dynamics from nuclear magnetic resonance nuclear Overhauser effect data.

The solution structure of the histidine-containing phosphocarrier protein HPr from Escherichia coli has been determined by NMR in combination with distance geometry and restrained molecular dynamics. The structure is based on 1520 experimental restraints identified from both three-dimensional 1H-1H-13C and 1H-1H-15N nuclear Overhauser effect multiple-quantum coherence spectroscopy and two-dimensional 1H-1H nuclear Overhauser effect spectra. Thirty-two four-dimensional coordinate frames were produced by metric matrix distance geometry, subjected to distance bounds driven dynamics, projected into three-dimensional space and again subjected to distance-bounds driven dynamics. These 32 distance geometry structures were refined further by restrained molecular dynamics (40 ps) in the GROMOS in vacuo force field. All 32 structures reached acceptable energy minima while satisfying the imposed restraints. Two of these structures were subjected to a further 200 ps of molecular dynamics simulation in water, using time-dependent distance restraining, followed by a 200 ps free simulation without any distance restraining. The resulting structure is very similar to the X-ray structure of Bacillus subtilis HPr, but differs mainly in the position of the two loops containing the active site histidine residue 15 and residues 53 to 57 relative to the rest of the structure. The unfavorable phi torsion angle that was found for residue 16 in the active center of unphosphorylated Streptococcus faecalis HPr was proposed to play a role in the activity of the protein. Although present at the early stages of the structure calculations, this torsion-angle strain disappeared in the final model obtained from molecular dynamics simulations in water using time-averaged distance restraining and upon releasing the distance restraints. This suggests that the strain may be an artifact of crystallization conditions instead of an essential element in the phosphorylation/dephosphorylation process.

Bacterial Proteins↗

PubMed1994

Investigation of shape variations in the antibody binding site by molecular dynamics computer simulation.

Molecular dynamics simulations have been used to investigate the flexibility and variations in the shape of the binding site of an antibody against human Rhinovirus serotype 2 (HRV2) and its complex with a 15 amino acid oligopeptide, the structure of which has been recently determined by X-ray crystallography. During the simulation of the unbound antibody the binding site, defined in terms of the hypervariable regions or complementarity determining regions (CDRs), shows significant fluctuations in shape. For the complex such variations in the shape of the binding site were reduced. The largest fluctuations in the unbound antibody occurred within the CDR-H3. The largest differences between the bound and unbound crystal structures are also associated with CDR-H3. The relative displacements of the loops have been analysed in terms of internal distortions, rigid body motions of the loops and changes with respect to the framework regions. The degree to which the motions of the loops are correlated and the variation in the volume of the binding pocket during the simulation have also been examined.

Amino Acid Sequence↗

PubMed1994

Translational and rotational diffusion of proteins.

We have investigated the translational and rotational diffusion of two proteins, bovine pancreatic trypsin inhibitor and hen egg white lysozyme, using molecular dynamics simulations in explicit solvent. The translational diffusion constants obtained from the simulations compare favourably with results from experiments. However, rotational diffusion constants indicate increased rotational diffusion compared to experiments. Further analysis, using simple hydration models, suggests that the observed rotational diffusion can be explained without invoking a fixed shell of hydration water associated with the protein. Within the assumptions of the hydration models, the results suggest that the strength of protein-water interactions in the current force field may be slightly underestimated.

Animals↗

PubMed1994

Local elevation: a method for improving the searching properties of molecular dynamics simulation.

The concept of memory has been introduced into a molecular dynamics algorithm. This was done so as to persuade a molecular system to visit new areas of conformational space rather than be confined to a small number of low-energy regions. The method is demonstrated on a simple model system and the 11-residue cyclic peptide cyclosporin A. For comparison, calculations were also performed using simulated temperature annealing and a potential energy annealing scheme. Although the method can only be applied to systems with a small number of degrees of freedom, it offers the chance to generate a multitude of different low-energy structures, where other methods only give a single one or few. This is clearly important in problems such as drug design, where one is interested in the conformational spread of a system.

Algorithms↗