Search

Staphylococcal nuclease forms three different partially-folded intermediates at low pH in the presence of low to moderate concentration of anions, differing in the amount of secondary structure, globularity, stability, and compactness. Although these intermediates are monomeric at low protein concentration (< or =0.25 mg/mL), increasing concentrations of protein result in the formation of dimers and soluble oligomers, ultimately leading to larger insoluble aggregates. Unexpectedly, increasing protein concentration not only led to association, but also to increased structure of the intermediates. The secondary structure, stability, and globularity of the two less-ordered partially-folded intermediates (A1 and A2) were substantially increased upon association, suggesting that aggregation induces structure. An excellent correlation was found between degree of association and amount of structure measured by different techniques, including circular dichroism, fluorescence, Fourier transform infrared spectroscopy (FTIR), and small-angle X-ray scattering. The associated states were also substantially more stable toward urea denaturation than the monomeric forms. A mechanism is proposed, in which the observed association of monomeric intermediates involves intermolecular interactions which correspond to those found intramolecularly in normal folding to the native state.

Structural studies of integral membrane proteins typically rely upon detergent micelles as faithful mimics of the native lipid bilayer. Therefore, membrane protein structure determination would be greatly facilitated by biophysical techniques that are capable of evaluating and assessing the fold and oligomeric state of these proteins solubilized in detergent micelles. In this study, an approach to the characterization of detergent-solubilized integral membrane proteins is presented. Eight Thermotoga maritima membrane proteins were screened for solubility in 11 detergents, and the resulting soluble protein–detergent complexes were characterized with small angle X-ray scattering (SAXS), nuclear magnetic resonance (NMR) spectroscopy, circular dichroism (CD) spectroscopy, and chemical cross-linking to evaluate the homogeneity, oligomeric state, radius of gyration, and overall fold. A new application of SAXS is presented, which does not require density matching, and NMR methods, typically used to evaluate soluble proteins, are successfully applied to detergent-solubilized membrane proteins. Although detergents with longer alkyl chains solubilized the most proteins, further characterization indicates that some of these protein–detergent complexes are not well suited for NMR structure determination due to conformational exchange and protein oligomerization. These results emphasize the need to screen several different detergents and to characterize the protein–detergent complex in order to pursue structural studies. Finally, the physical characterization of the protein–detergent complexes indicates optimal solution conditions for further structural studies for three of the eight overexpressed membrane proteins.

We present an algorithm for simulating the long time scale dynamics of proteins and other macromolecules. Our method applies the concept of multiple time step integration to the diffusive Langevin equation, in which short time scale dynamics are replaced by friction and noise. The macromolecular force field is represented at atomic resolution. Slow motions are modeled by constrained Langevin dynamics with very large time steps, while faster degrees of freedom are kept in local thermal equilibrium. In the limit of a sufficiently large molecule, our algorithm is shown to reduce the CPU time required by two orders of magnitude. We test the algorithm on two systems, alanine dipeptide and bovine pancreatic trypsin inhibitor (BPTI), and find that it accurately calculates a variety of equilibrium and dynamical properties. In the case of BPTI, the CPU time required is reduced by nearly a factor of 60 compared to a conventional, unconstrained Langevin simulation using the same force field. Proteins 30:215–227, 1998. © 1998 Wiley-Liss, Inc.

We present a novel algorithm of constrained, overdamped dynamics to study the long-time properties of peptides, proteins, and related molecules. The constraints are applied to an all-atom model of the molecule by projecting out all components of the nonbonding interactions which tend to alter fixed bond lengths and angles. Because the overdamped dynamical equations are first order in time, the constraints are satisfied by inversion of a banded matrix at each timestep, which is computationally efficient. Thermal effects are included through a Langevin noise term in the equation of motion. Because high-frequency components of the motion have been eliminated, the timestep of the algorithm is determined by the nonbonding forces, which are two to three orders of magnitude weaker than the bonding forces. Using polyalanine as a test example, we demonstrate that trajectories simulating a microsecond of motion can be run about 103 times faster than an equivalent molecular dynamics simulation. © 1994 by John Wiley & Sons, Inc.

A Monte Carlo algorithm that searches for the optimal docking configuration of hen egg white lysozyme to an antibody is de- veloped. Both the lysozyme and the antibody are kept rigid. Unlike the work of other au- thors, our algorithm does not attempt to explic- itly maximize surface contact, but minimizes the energy computed using coarse-grained pair potentials. The final refinement of our best so- lutions using all-atom OPLS potentials (Jor- gensen and Tirado-Rives') consistently yields the native conformation as the preferred solu- tion for three different antibodies. We find that the use of an exponential distance-dependent dielectric function is an improvement over the more commonly used linear form. 0 1993 Wiley-Liss, Inc.

The current work describes a simplified representation of protein structure with uses in the simulation of protein folding. The model assumes that a protein can be represented by a freely rotating rigid chain with a single atom approximating the effect of each side chain. Potentials describing the attraction or repulsion between different types of amino acids are determined directly from the distribution of amino acids in the database of known protein structures. The optimization technique of simulated annealing has been used to dynamically sample the conformations available to this simple model, allowing the protein to evolve from an extended, random coil into a compact globular structure. Many characteristics expected of true proteins, such as the sequence-dependent formation of secondary structure, the partitioning of hydrophobic residues, and specific disulfide pairing, are reproduced by the simulation, suggesting the model may accurately simulate the folding process.