An analysis of protein folding pathways

We have developed a model of the protein folding process based on three primary assumptions: that burying of hydrophobic area is the dominant contribution to the relative free energy of a conformation, that a record of the folding process is largely preserved in the final structure, and that the denatured state is a random coil. Detailed folding pathways are identified for 19 protein structures. The picture of the folding process that emerges from this analysis is one of nucleation by regions of 8-16 residues. Nucleation sites then lead to larger structures by two mechanisms: propagation and diffusion/collision. A Monte Carlo simulation is used to follow the folding pathway when propagation is the dominant mechanism. Because detailed pathways are derived for each protein, the models are susceptible to experimental verification.

Molecular dynamics study of the structure and dynamics of a protein molecule in a crystalline ionic environment, Streptomyces griseus protease A

A large-scale molecular dynamics simulation of the behavior of a serine protease (Streptomyces griseus protease A) in a crystalline environment has been performed. All atoms (including hydrogens) of two protein molecules and the surrounding solvent of crystallization, consisting of both water and salt ions, were explicitly represented, and a relatively long range of interactions (up to 15 A) were included. The simulation is the longest so far reported for a protein in such an environment (60 ps). The use of the full crystalline environment allows a direct comparison of the structure and dynamic properties of the protein and surrounding solvent to be made with the experimental X-ray structure. Here we report the comparison of the protein structures and analyze the energetics of the system, including interaction with the aqueous environment. Subsequent papers will deal with other aspects of the simulation. The overall root mean square differences between the time-averaged molecular dynamics structure and that from crystallography, for all well-ordered, non-hydrogen atoms, are 1.67 and 1.25 A for the two molecules taken as the asymmetric unit. An extensive analysis of the conformation of substructural elements and individual residues and their deviation from experiment has revealed a strong influence of the ionic medium on their behavior. Implications of the results for free energy calculations and for future directions are also discussed.

Probing the calcium-induced conformational transition of troponin C with site-directed mutants

The contraction of skeletal muscle is regulated by calcium binding to troponin C (TnC). TnC consists of two spatially independent domains, each of which contains two metal ion binding sites. Calcium binding to the regulatory sites of the N-terminal domain triggers muscle contraction by a series of conformational changes. Site-directed mutagenesis offers a means of elucidating the links in this signal path between TnC and actin-myosin crossbridges. Such mapping is possible if the mutants shift the equilibrium between 'on' and 'off' states of the regulatory complex while maintaining the coupling between calcium binding and tension development. Candidate amino-acid residues for yielding this information would be in positions remote from the calcium-binding sites and from the site of development of tension. Analysis of the crystal structure of TnC and of the model of the calcium-activated molecule has enabled us to identify two such residues: Glu 57 and Glu 88. In separate experiments we have replaced each of these residues by lysines. The resulting reduction in calcium affinity indicates that these residues have a long-range effect on calcium binding. This result may reflect the formation of a salt bridge between positions 57 and 88 that is not present in the native molecule. Moreover, the level of tension recovery when the mutants are incorporated into muscle suggests that the interaction between TnC and other muscle components has also been altered. Thus, these residues may participate in the contraction signal transmission.

Antibody-dependent cellular cytotoxicity of vascular endothelium: characterization and pathogenic associations in systemic sclerosis

Ten sera from 48 patients with systemic sclerosis were found to be capable of producing cytotoxicity of human umbilical venous and arterial endothelium when co-cultured with peripheral blood mononuclear cells. Fractionation of sera on Ultrogel and the preparation of monomeric IgG by ion exchange and affinity chromatography suggested that the cytotoxicity was mediated by anti-endothelial antibodies capable of pre-sensitizing target cells in a mechanism that resembled antibody-dependent cellular cytotoxicity. These anti-endothelial antibodies together with C1q-binding immune complexes and anti-cardiolipin antibodies were found in 18 of 28 patients so investigated, suggesting that multiple immunological mechanisms may be involved in the pathogenesis of the vascular lesion of systemic sclerosis.

Prostacyclin production by human umbilical vein endothelium in response to serum from patients with systemic sclerosis

Sera from 29 patients with systemic sclerosis and 30 normal controls were examined for their effect on prostacyclin release by human umbilical vein endothelial cells during periods of response of 15 min and 72 h and for their effect on endothelial growth and 3H-thymidine uptake during a 72-h culture period. In contrast to previous reports, no significant differences were detected between patient and control sera in their effect on endothelial cell prostacyclin release, growth or 3H-thymidine uptake.

Elevated von Willebrand factor antigen in systemic sclerosis: relationship to visceral disease

Plasma levels of the factor VIII complex (von Willebrand factor antigen, factor VIII coagulant and ristocetin co-factor) were measured in 28 patients with systemic sclerosis. Elevated von Willebrand factor antigen was found in 12 patients overall and in 10 of 16 patients characterized by severe extensive visceral disease, with a resulting positive correlation between the extent of visceral involvement and the plasma level of von Willebrand factor antigen (r = 0.60, p less than 0.001). Factor VIII coagulant and ristocetin co-factor levels, however, frequently failed to parallel the increases of von Willebrand factor antigen, supporting the view that these increases were due to in vivo endothelial damage. The findings suggest that vascular damage is an important aspect of the visceral lesions of systemic sclerosis.

Comparative molecular model building of two serine proteinases from cytotoxic T lymphocytes

Two genes that are expressed when precursor cytotoxic T lymphocytes are transformed to T killer cells have been cloned and sequenced. The derived amino acid sequences, coding for cytotoxic cell protease 1 (CCP1) and Hannuka factor (HF) are highly homologous to members of the serine proteinase family. Comparative molecular model building using the known three-dimensional structures and the derived amino acid sequences of the lymphocyte enzymes has provided useful structural information, especially in predicting the conformations of the substrate binding sites. In applying this modelling procedure, we used the X-ray structures of four serine proteinases to provide a structurally based sequence alignment: alpha-chymotrypsin (CHT), bovine trypsin (BT), Streptomyces griseus trypsin (SGT), and rat mast cell protease 2 (RMCP2). The root mean square differences in alpha-carbon atom positions among these four structures when compared in a pairwise fashion range from 0.79 to 0.97 A for structurally equivalent residues. The sequences of the two lymphocyte enzymes were then aligned to these proteinases using chemical criteria and the superimposed X-ray structures as guides. The alignment showed that the sequence of CCP1 was most similar to RMCP2, whereas HF has regions of homology with both RMCP2 and BT. With RMCP2 as a template for CCP1 and the two enzymes RMCP2 and BT as templates for HF, the molecular models were constructed. Intramolecular steric clashes that resulted from the replacement of amino acid side chains of the templates by the aligned residues of CCP1 and HF were relieved by adjustment of the side chain conformational angles in an interactive computer graphics device. This process was followed by energy minimization of the enzyme model to optimize the stereochemical geometry and to relieve any remaining unacceptably close nonbonded contacts. The resulting model of CCP1 has an arginine residue at position 226 in the specificity pocket, thereby predicting a substrate preference for P1 aspartate or glutamate residues. The model also predicts favorable binding for a small hydrophobic residue at the P2 position of the substrate. The primary specificity pocket of HF resembles that of BT and therefore predicts a lysine or arginine preference for the P1 residue. The arginine at position 99 in the model of HF suggests a preference for aspartate or glutamate side chains in the P2 position of the substrate. Both CCP1 and HF have a free cysteine in the segment of polypeptide 88 to 93.(ABSTRACT TRUNCATED AT 400 WORDS)

Bacterial resistance to beta-lactam antibiotics: crystal structure of beta-lactamase from Staphylococcus aureus PC1 at 2.5 A resolution

beta-lactamases are enzymes that protect bacteria from the lethal effects of beta-lactam antibiotics, and are therefore of considerable clinical importance. The crystal structure of beta-lactamase from the Gram-positive bacterium Staphylococcus aureus PC1 has been determined at 2.5 angstrom resolution. It reveals a molecule of novel topology, made up of two closely associated domains. The active site is located at the interface between the domains, with the key catalytic residue Ser70 at the amino terminus of a buried helix. Examination of the disposition of the functionally important residues within the active site depression leads to a model for the binding of a substrate and a functional analogy to the serine proteases. The unusual topology of the secondary structure units is relevant to questions concerning the evolutionary relation to the beta-lactam target enzymes of the bacterial cell wall.

An algorithm for determining the conformation of polypeptide segments in proteins by systematic search

The feasibility of determining the conformation of segments of a polypeptide chain up to six residues in length in globular proteins by means of a systematic search through the possible conformations has been investigated. Trial conformations are generated by using representative sets of phi, psi, and chi angles that have been derived from an examination of the distributions of these angles in refined protein structures. A set of filters based on simple rules that protein structures obey is used to reduce the number of conformations to a manageable total. The most important filters are the maintenance of chain integrity and the avoidance of too-short van der Waals contacts with the rest of the protein and with other portions of the segment under construction. The procedure is intended to be used with approximate models so that allowance is made throughout for errors in the rest of the structure. All possible main chains are first constructed and then all possible side-chain conformations are built onto each of these. The electrostatic energy, including a solvent screening term, and the exposed hydrophobic area are evaluated for each accepted conformation. The method has been tested on two segments of chain in the trypsin like enzyme from Streptomyces griseus. It is found that there is a wide spread of energies among the accepted conformations, and the lowest energy ones have satisfactorily small root mean square deviations from the X-ray structure.

A model for the Ca2+-induced conformational transition of troponin C. A trigger for muscle contraction

The initial contractile event in muscle is the binding of Ca2+ ions to troponin C of the troponin complex, leading to a series of conformational changes in the members of the thin and thick filaments. Knowledge of the crystal structure of turkey skeletal muscle troponin C has provided a structural basis for the modeling of the first stage of this process in atomic detail. This crystal structure probably represents the molecule in the relaxed state of muscle, with two of the maximum of 4 Ca2+ ions bound. The basis for the model presented here is that upon binding of the additional two Ca2+ ions, the regulatory domain of the molecule undergoes a conformational transition to become closely similar in structure to the domain which always binds Ca2+ or Mg2+ under physiological conditions. The root mean square discrepancy in atomic coordinates between the apo and the modeled Ca2+-bound states of the regulatory domain is 4.8 A, with some shifts as large as 10-15 A in the region near the linker between the two Ca2+ binding sites. It is demonstrated that this Ca2+-bound conformation of the regulatory domain conforms to accepted protein structure rules and that the change in conformation can be accomplished without encountering any barriers too high to be surmounted on the physiological time scale.

Electron density calculations as an extension of protein structure refinement. Streptomyces griseus protease A at 1.5 A resolution

Ab initio quantum mechanical calculations have been used to obtain details of the electron density distribution in a high-resolution refined protein structure. It is shown that with accurate atomic co-ordinates, electron density may be calculated with a quality similar to that which can be obtained directly from crystallographic studies of small organic molecules, and that this density contains information relevant to the understanding of catalysis. Atomic co-ordinates from the 1.8 A and 1.5 A resolution refinements of the crystal structure of protease A from Streptomyces griseus have been used to examine the influence of the environment on the electron density in the side-chain of the active site histidine (His57). The neighbouring aspartic acid 102 is the dominant factor in the environment, and quantum mechanical calculations have been performed on these two residues. Most interesting from the point of view of understanding the catalytic process is the effect that Asp102 has on the electron density in the region of the imidazole nitrogen (N epsilon 2) adjacent to the active site serine 195. In the positively charged imidazolium species, there is a polarization of the N epsilon 2-H bond, reducing the bonding density in a manner that may lower the height of the energy barrier for proton transfer. In the uncharged imidazole species, the proximity of Asp102 causes a movement of density from the lone pair region of the N epsilon 2 into the pi bonding region above and below the plane of the ring. Although it is shown that the primary effect of the aspartic acid is electrostatic, this movement is perpendicular to the direction of the electric field inducing it.

Computer simulation of the solvent structure around biological macromolecules

The structure and energetics of the solvent, both ordered and disordered, in a small peptide crystal and in the triclinic lysozyme crystal have been simulated using the Monte Carlo method. The results are in good agreement with the experimental data and provide insight into the role of solvent structure around biological macromolecules in solution.

Structural analysis of denaturant-protein interactions: comparison between the effects of bromoethanol and SDS on denaturation and renaturation of triclinic lysozyme

This paper summarizes our crystallographic studies of the interaction of denaturants with cross-linked triclinic lysozyme. Electron density maps of various bromoethanol-lysozyme complexes are analyzed and compared to those reported earlier for SDS-lysozyme complexes. Despite differences in the chemical nature and size of the two denaturants their mode of interaction with the protein is quite similar, suggesting the existence of a general mechanism for binding of hydrophobic-hydrophilic denaturants to proteins. Our results are consistent with the conclusion that lysozyme consists of two domains connected by a flexible segment and that this segment represents an internal degree of freedom of the protein.

Crystallographic studies of protein denaturation and renaturation. 1. Effects of denaturants on volume and X-ray pattern of cross-linked triclinic lysozyme crystals

Triclinic crystals of hen egg-white lysozyme cross-linked with glutaraldehyde have been treated with various denaturants and found to be susceptible to x-ray structure analysis even after major conformational changes in the protein. Cross-linked crystals were isomorphous with the native form, and electron density difference maps indicated the locations of intermolecular corss-links, but showed no appreciable differences in the protein conformation. Soaking of the cross-linked crystals in danaturant solutions of increasing concentrations caused corresponding increases in crystal volume and decreases in minimum observable x-ray spacings. These changes proved partly reversible on diluting the solutions, and measurements of crystal volume and minimums x-ray spacing were used to follow denaturation and renaturation as a function of concentration for several denaturants. Some of these, including bromoethanol and sodium dodecyl sulfate, had little effect on the crystals below critical concentrations at which there was a sharp volume increase and loss of x-ray pattern, which could, however, be regenerated to about 3.2-A resolution. Others, including KCNS and urea, caused more gradual changes, but with a smaller degree of recovery. It is suggested that at least two different denaturation mechanisms are involved with detergent-like reagents disrupting the hydrophobic interactions joining the two wings of the lysozyme molecule and hydrophilic denaturants interacting primarily with polar groups on the molecular surface.