Monte carlo simulations of protein folding. II. Application to protein A,ROP, and crambin

The hierarchy of lattice Monte Carlo models described in the accompanying paper (Kolinski, A., Skolnick, J. Monte Carlo simulations of protein folding. I. Lattice model and interaction scheme. Proteins 18:338–352, 1994) is applied to the simulation of protein folding and the prediction of 3-dimensional structure. Using sequence information alone, three proteins have been successfully folded: the B domain of staphylococcal protein A, a 120 residue, monomeric version of ROP dimer, and crambin. Starting from a random expanded conformation, the model proteins fold along relatively well-defined folding pathways. These involve a collection of early intermediates, which are followed by the final (and rate-determining) transition from compact intermediates closely resembling the molten globule state to the native-like state. The predicted structures are rather unique, with native-like packing of the side chains. The accuracy of the predicted native conformations is better than those obtained in previous folding simulations. The best (but by no means atypical) folds of protein A have a coordinate rms of 2.25 Å from the native Cα trace, and the best coordinate rms from crambin is 3.18 Å. For ROP monomer, the lowest coordinate rms from equivalent Cαs of ROP dimer is 3.65 Å. Thus, for two simple helical proteins and a small α/β protein, the ability to predict protein structure from sequence has been demonstrated. © 1994 John Wiley & Sons, Inc.     

Exploration of disorder in protein structures by X-ray restrained molecular dynamics

Conformational disorder in crystal structures of ribonuclease-A and crambin is studied by including two independent structures in least-squares optimizations against X-ray data. The optimizations are carried out by X-ray restrained molecular dynamics (simulated annealing refinement) and by conventional least-squares optimization. Starting from two identical structures, the optimizations against X-ray data lead to significant deviations between the two, with rms backbone displacements of 0.45 A for refinement of ribonuclease at 1.53 A resolution, and 0.31 A for crambin at 0.945 A. More than 15 independent X-ray restrained molecular dynamics runs have been carried out for ribonuclease, and the displacements between the resulting structures are highly reproducible for most atoms. These include residues with two or more conformations with significant dihedral angle differences and alternative hydrogen bonding, as well as groups of residues that undergo displacements that are suggestive of rigid-body librations. The crystallographic R-values obtained are approximately 13%, as compared to 15.3% for a comparable refinement with a single structure. Least-squares optimization without an intervening restrained molecular dynamics stage is sufficient to reproduce most of the observed displacements. Similar results are obtained for crambin, where the higher resolution of the X-ray data allows for refinement of unconstrained individual anisotropic temperature factors. These are shown to be correlated with the displacements in the two-structure refinements.     

“Ensemble” iterative relaxation matrix approach: A new NMR refinement protocol applied to the solution structure of crambin

The structure in solution of crambin, a small protein of 46 residues, has been determined from 2D NMR data using an iterative relaxation matrix approach (IRMA) together with distance geometry, distance bound driven dynamics, molecular dynamics, and energy minimization. A new protocol based on an “ensemble” approach is proposed and compared to the more standard initial rate analysis approach and a “single structure” relaxation matrix approach. The effects of fast local motions are included and R-factor calculations are performed on NOE build-ups to describe the quality of agreement between theory and experiment. A new method for stereospecific assignment of prochiral groups, based on a comparison of theoretical and experimental NOE intensities, has been applied. The solution structure of crambin could be determined with a precision (rmsd from the average structure) of 0.7 Å on backbone atoms and 1.1 Å on all heavy atoms and is largely similar to the crystal structure with a small difference observed in the position of the side chain of Tyr-29 which is determined in solution by both J-coupling and NOE data. Regions of higher structural variability (suggesting higher mobility) are found hi the solution structure, in particular for the loop between the two helices (Gly-20 to Pro-22). © 1993 Wiley-Liss, Inc.     

Influence of rotational energy barriers to the conformational search of protein loops in molecular dynamics and ranking the conformations.

An adjustable-barrier dihedral angle potential was added as an extension to a novel, previously presented soft-core potential to study its contribution to the efficacy of the search of the conformational space in molecular dynamics. As opposed to the conventional soft-core potential functions, the leading principle in the design of the new soft-core potential, as well as of its extension, the soft-core and adjustable-barrier dihedral angle (SCADA) potential (referred as the SCADA potential), was to maintain the main equilibrium properties of the original force field. This qualifies the methods for a variety of a priori modeling problems without need for additional restraints typically required with the conventional soft-core potentials. In the present study, the different potential energy functions are applied to the problem of predicting loop conformations in proteins. Comparison of the performance of the soft-core and SCADA potential showed that the main hurdles for the efficient sampling of the conformational space of (loops in) proteins are related to the high-energy barriers caused by the Lennard-Jones and Coulombic energy terms, and not to the rotational barriers, although the conformational search can be further enhanced by lowering the rotational barriers of the dihedral angles. Finally, different evaluation methods were studied and a few promising criteria found to distinguish the near-native loop conformations from the wrong ones.     

Variations on a theme by Debye and Waller: From simple crystals to proteins

Debye and Waller showed how to adjust scattering intensities in diffraction experiments for harmonic motions of atoms about an average, static reference configuration. However, many motions, particularly in biological molecules as compared to simple crystals, are far from harmonic. We show how, using a variety of simple anharmonic, multiconformational models, it is possible to construct a variety of Generalized Debye-Waller Factors, and understand their meaning. A central result for these cases is that, in principle, intensity factors cannot be obtained from true total mean square displacements of the atoms. We make the distinction between the intensity factors for unimodal quasiharmonic motions and those for the anharmonic, multimodal (valley hopping) motions. Only the former affect the conventional B factors. Proteins 29:153–160, 1997. Published 1997 Wiley-Liss, Inc.

1 This article is a US Government work and, as such, is in the public domain in the United States of America.

     

Monte Carlo simulations of a protein molecule with and without hydration energy calculated by the hydration-shell model

Monte Carlo simulations of a small protein, carmbin, were carried out with and without hydration energy. The methodology presented here is characterized, as compared with the other similar simulations of proteins in solution, by two points: (1) protein conformations are treated in fixed geometry so that dihedral angles are independent variables rather than cartesian coordinates of atoms; and (2) instead of treating water molecules explicitly in the calculation, hydration energy is incorporated in the conformational energy function in the form of g _i A _i, whereA _i is the accessible surface area of an atomic groupi in a given conformation, andg _i is the free energy of hydration per unit surface area of the atomic group (i.e., hydration-shell model). Reality of this model was tested by carrying out Monte Carlo simulations for the two kinds of starting conformations, native and unfolded ones, and in the two kinds of systems,in vacuo and solution. In the simulations starting from the native conformation, the differences between the mean propertiesin vacuo and solution simulations are not very large, but their fluctuations around the mean conformation during the simulation are relatively smaller in solution thanin vacuo. On the other hand, in the simulations starting from the unfolded conformation, the molecule fluctuates much more largely in solution thanin vacuo, and the effects of taking into account the hydration energy are pronounced very much. The results suggest that the method presented in this paper is useful for the simulations of proteins in solution.