Conformational dynamics of chymotrypsin inhibitor 2 by coarse-grained simulations

Features of large hinge-bending conformational transitions. Prediction of closed structure from open state

We performed a detailed analysis of conformational transition pathways for a set of 10 proteins, which undergo large hinge-bending-type motions with 4-12 Å RMSD (root mean-square distance) between open and closed crystal structures. Anisotropic network model-Monte Carlo (ANM-MC) algorithm generates a targeted pathway between two conformations, where the collective modes from the ANM are used for deformation at each iteration and the conformational energy of the deformed structure is minimized via an MC algorithm. The target structure was approached successfully with an RMSD of 0.9-4.1 Å when a relatively low cutoff radius of 10 Å was used in ANM. Even though one predominant mode (first or second) directed the open-to-closed conformational transition, changes in the dominant mode character were observed for most cases along the transition. By imposing radius of gyration constraint during mode selection, it was possible to predict the closed structure for eight out of 10 proteins (with initial 4.1-7.1 Å and final 1.7-2.9 Å RMSD to target). Deforming along a single mode leads to most successful predictions. Based on the previously reported free energy surface of adenylate kinase, deformations along the first mode produced an energetically favorable path, which was interestingly facilitated by a change in mode shape (resembling second and third modes) at key points. Pathway intermediates are provided in our database of conformational transitions (http://safir.prc.boun.edu.tr/anmmc/method/1).

Conformational transition pathways explored by Monte Carlo simulation integrated with collective modes

Conformational transitions between open/closed or free/bound states in proteins possess functional importance. We propose a technique in which the collective modes obtained from an anisotropic network model (ANM) are used in conjunction with a Monte Carlo (MC) simulation approach, to investigate conformational transition pathways and pathway intermediates. The ANM-MC technique is applied to adenylate kinase (AK) and hemoglobin. The iterative method, in which normal modes are continuously updated during the simulation, proves successful in accomplishing the transition between open-closed conformations of AK and tense-relaxed forms of hemoglobin (C(alpha)-root mean square deviations between two end structures of 7.13 A and 3.55 A, respectively). Target conformations are reached by root mean-square deviations of 2.27 A and 1.90 A for AK and hemoglobin, respectively. The intermediate conformations overlap with crystal structures from the AK family within a 3.0-A root mean-square deviation. In the case of hemoglobin, the transition of tense-to-relaxed passes through the relaxed state. In both cases, the lowest-frequency modes are effective during transitions. The targeted Monte Carlo approach is used without the application of collective modes. Both the ANM-MC and targeted Monte Carlo techniques can explore sequences of events in transition pathways with an efficient yet realistic conformational search.

Mutations in alpha-helical solvent-exposed sites of eglin c have long-range effects: evidence from molecular dynamics simulations

Eglin c is a small protease inhibitor whose structural and thermodynamic properties have been well studied. Previous thermodynamic measurements on mutants at solvent-accessible positions in the protein's helix have shown the unexpected result that the data could be best fit by the inclusion of residue- and position-specific parameters to the model. To explore the origins of this surprising result, long molecular dynamics simulations in explicit solvent have been performed. These simulations indicate specific long-range interactions between the solvent-exposed residues in the eglin c alpha-helix and binding loop, an unexpected observation for such a small protein. The residues involved in the interaction are on opposite sides of the protein, about 25 A apart. Simulations of alanine substitutions at the solvent-exposed helix positions, arginine 22, glutamic acid 23, threonine 26, and leucine 27, show both small and large perturbations of eglin c dynamics. Two mutations exhibit large impacts on the long-range helix-loop interactions. Previous stability measurements (Yi et al., Biochemistry 2003;42:7594-7603) had indicated that an alanine substitution at position 27 was less stabilizing than at other solvent-exposed positions in the helix. The L27A mutation effects observed in these simulations suggest that the position-dependent loss of stability measured in wet bench experiments is derived from changes in dynamics that involve long-range interactions; thus, these simulations support the hypothesis that solvent-exposed positions in helices are not always equivalent.

Monte Carlo studies of folding, dynamics, and stability in alpha-helices

Folding simulations of polyalanine peptides were carried out using an off-lattice Monte Carlo simulation technique. The peptide was represented as a chain of residues, each of which contains two interaction sites: one corresponding to the C(alpha) atom and the other to the side chain. A statistical potential was used to describe the interaction between these sites. The preferred conformations of the peptide chain on the energy surface, starting from several initial conditions, were searched by perturbations on its generalized coordinates with the Metropolis criterion. We observed that, at low temperatures, the effective energy was low and the helix content high. The calculated helix propagation (s) and nucleation (sigma) parameters of the Zimm-Bragg model were in reasonable agreement with the empirical data. Exploration of the energy surface of the alanine-based peptides (AAQAA)(3) and AAAAA(AAARA)(3)A demonstrated that their behavior is similar to that of polyalanine, in regard to their effective energy, helix content, and the temperature-dependence of their helicity. In contrast, stable secondary structures were not observed for (Gly)(20) at similar temperatures, which is consistent with the nonfolder nature of this peptide. The fluctuations in the slowest dynamics mode, which describe the elastic behavior of the chain, showed that as the temperature decreases, the polyalanine peptides become stiffer and retain conformations with higher helix content. Clustering of conformations during the folding phase implied that polyalanine folds into a helix through fewer numbers of intermediate conformations as the temperature decreases.

Interactions of hydrophobic peptides with lipid bilayers: Monte Carlo simulations with M2delta

We introduce here a novel Monte Carlo simulation method for studying the interactions of hydrophobic peptides with lipid membranes. Each of the peptide's amino acids is represented as two interaction sites: one corresponding to the backbone alpha-carbon and the other to the side chain, with the membrane represented as a hydrophobic profile. Peptide conformations and locations in the membrane and changes in the membrane width are sampled using the Metropolis criterion, taking into account the underlying energetics. Using this method we investigate the interactions between the hydrophobic peptide M2delta and a model membrane. The simulations show that starting from an extended conformation in the aqueous phase, the peptide first adsorbs onto the membrane surface, while acquiring an ordered helical structure. This is followed by formation of a helical-hairpin and insertion into the membrane. The observed path is in agreement with contemporary understanding of peptide insertion into biological membranes. Two stable orientations of membrane-associated M2delta were obtained: transmembrane (TM) and surface, and the value of the water-to-membrane transfer free energy of each of them is in agreement with calculations and measurements on similar cases. M2delta is most stable in the TM orientation, where it assumes a helical conformation with a tilt of 14 degrees between the helix principal axis and the membrane normal. The peptide conformation agrees well with the experimental data; average root-mean-square deviations of 2.1 A compared to nuclear magnetic resonance structures obtained in detergent micelles and supported lipid bilayers. The average orientation of the peptide in the membrane in the most stable configurations reported here, and in particular the value of the tilt angle, are in excellent agreement with the ones calculated using the continuum-solvent model and the ones observed in the nuclear magnetic resonance studies. This suggests that the method may be used to predict the three-dimensional structure of TM peptides.

Structure-based prediction of potential binding and nonbinding peptides to HIV-1 protease

HIV-1 protease is a major drug target against AIDS as it permits viral maturation by processing the gag and pol polyproteins of the virus. The cleavage sites in these polyproteins do not have obvious sequence homology or a binding motif and the specificity of the protease is not easily determined. We used various threading approaches, together with the crystal structures of substrate complexes which served as template structures, to study the substrate specificity of HIV-1 protease with the aim of obtaining a better differentiation between binding and nonbinding sequences. The predictions from threading improved when distance-dependent interaction energy functions were used instead of contact matrices. To rank the peptides and properly account for the peptide's conformation in the total energy, the results from using short-range potentials on multiple template structures were averaged. Finally, a dynamic threading approach is introduced which is potentially useful for cases when there is only one template structure available. The conformational energy of the peptide-especially the term accounting for the side chains-was found to be important in differentiating between binding and nonbinding sequences. Hence, the substrate specificity, and thus the ability of the virus to mature, is affected by the compatibility of the substrate peptide to fit within the limited conformational space of the active site groove.

Coordination topology and stability for the native and binding conformers of chymotrypsin inhibitor 2

We demonstrate that the stabilization of the binding region is accomplished at the expense of a loss in the stability of the rest of the protein. A novel molecular mechanics (MM) approach is introduced to distinguish residue stabilities of proteins in a given conformation. As an example, the relative stabilities of folded chymotrypsin inhibitor 2 (CI2) in unbound form, and CI2 in complex with subtilisin novo is investigated. The conformation of the molecule in the two states is almost identical, with an approximately 0.6-A root-mean-square deviation (RMSD) of the Calpha atoms. On binding, the packing density changes only at the binding loop. However, residue fluctuations in the rest of the protein are greatly altered solely due to those contacts, indicating the effective propagation of perturbation and the presence of remotely controlling residues. To quantify the interplay between packing density, packing order, residue fluctuations, and residue stability, we adopt an MM approach whereby small displacements are inserted at selected residues, followed by energy minimization; the displacement of each residue in response to such perturbations are organized in a perturbation-response matrix L. We define residue stability lambda(i) = summation operator((j)L(ij))/ summation operator((j) L(ji)) as the ratio of the amount of change to which the residue is amenable, to the ability of a given residue to induce change. We then define the free energy associated with residue stability, DeltaG(lambda) = -RT ln lambda. DeltaG(lambda) intrinsically selects the residues that are in the folding core. Upon complexation, the binding loop becomes more resistant to perturbation, in contrast to the alpha-helix that favors change. Although the two forms of CI2 are structurally similar, residue fluctuations differ vastly, and the stability of many residues is altered upon binding. The decrease in entropy introduced by binding is thus compensated by these changes.

Conformational dynamics of subtilisin-chymotrypsin inhibitor 2 complex by coarse-grained simulations

An off-lattice dynamic Monte Carlo (MC) method is used to investigate the conformational dynamics of chymotrypsin inhibitor 2 (CI2) and subtilisin in both free and complex forms over two time windows, referring to short and long time scales. The conformational dynamics of backbone bonds analysed from several independent trajectories reveal that: Both the inhibitor and the enzyme are restricted in their bond rotations, excluding a few bonds, upon binding; the effect being greatest for the loop regions, and for the inhibitor. A cooperativity in the near-neighbor bond rotations are observed on both time scales, whereas the cooperative rotations of the bonds far along the sequence appear only in the long time window, and the latter time window is where most of the interactions between the inhibitor and the enzyme are observed. Upon binding, the cooperatively rotating parts of the inhibitor and the enzyme are readjusted compared to their free forms, and new correlations appear. The binding loop, although it is the closest contact region, is not the only part of the inhibitor involved in the interactions with the enzyme. Loops 3 and 8 and the helices F and G in bound enzyme and the binding loop of the inhibitor contribute at the most to the collective motions of whole structure on the slow time scale and are apparently important for enzyme-inhibitor interactions and function. The results in general provide evidence for the contribution of the loops with cooperative motions to the extensive communication network of the complex.