A critical analysis of continuum electrostatics: the screened Coulomb potential--implicit solvent model and the study of the alanine dipeptide and discrimination of misfolded structures of proteins

An analysis of the screened Coulomb potential--implicit solvent model (SCP--ISM) is presented showing that general equations for both the electrostatic and solvation free energy can be derived in a continuum approach, using statistical averaging of the polarization field created by the solvent around the molecule. The derivation clearly shows how the concept of boundary, usually found in macroscopic approaches, is eliminated when the continuum model is obtained from a microscopic treatment using appropriate averaging techniques. The model is used to study the alanine dipeptide in aqueous solution, as well as the discrimination of native protein structures from misfolded conformations. For the alanine dipeptide the free energy surface in the phi--psi space is calculated and compared with recently reported results of a detailed molecular dynamics simulation using an explicit representation of the solvent, and with other available data. The study showed that the results obtained using the SCP--ISM are comparable to those of the explicit water calculation and compares favorably to the FDPB approach. Both transition states and energy minima show a high correlation (r > 0.98) with the results obtained in the explicit water analysis. The study of the misfolded structures of proteins comprised the analysis of three standard decoy sets, namely, the EMBL, Park and Levitt, and Baker's CASP3 sets. In all cases the SCP--ISM discriminated well the native structures of the proteins, and the best-predicted structures were always near-native (cRMSD approximately 2 A).     

Influence of local interactions on protein structure. I. Conformational energy studies of N-acetyl-N′-methylamides of pro-X and X-pro dipeptides

Conformational energy calculations using an Empirical Conformational Energy Program for Peptides (ECEPP) were carried out on the N-acetyl-N′-methylamides of Pro-X, where X = Ala, Asn, Asp, Gly, Leu, Phe, Ser, and Val, and of X-Pro, where X = Ala, Asn, Gly, and Pro. The conformational energy was minimized from starting conformations which included all combinations of low-energy single-residue minima and several standard bend structures. It was found that almost all resulting minima are combinations of low-energy single-residue minima, suggesting that intra residue interactions predominate in determining conformation. The calculations also indicate, however, that inter residue interactions can be important. In addition, librational entropy was found to influence the relative stabilities of some minima. Because of the existence of 10–100 low-energy minima for each dipeptide, the normalized statistical weight of an individual minimum rarely exceeds 0.3, suggesting that these dipeptides have considerable conformational flexibility and exist as statistical ensembles of low-energy structures. The propensity of each dipeptide to form bend conformations was calculated, and the results were compared with available experimental data. It was found that bends are favored in Pro-X dipeptides because ?_Pro is fixed by the pyrrolidine ring in a conformation which is frequently found in bends, but that bends are not favored in X-Pro dipeptides because interactions between the X residue and the pyrrolidine ring restrict the X residue to conformations which are not usually found in bends.     

Free energy distributions in proteins.

Protein molecules in solution have a broad distribution of enthalpy states. A good approximation to the distribution function for enthalpy states can be calculated, using the maximum-entropy method, from the moments of the distribution that, in turn, are obtained from the experimental temperature dependence of the heat capacity. In the present paper, we show that the enthalpy probability distribution can then be formulated in terms of a free energy function that gives the free energy of the protein corresponding to a particular value of the enthalpy. By the location of the minima in this function, the free energy distribution graphically indicates the most probable values of the enthalpy for the protein. We find that the behavior of the free energy functions for proteins falls somewhere between two different cases: a two-state like function with two minima, the relative levels of the two states changing with temperature; and, a single-minimum function where the position of the minimum shifts to higher enthalpy values as the temperature is increased. We show that the temperature dependence of the free energy function can be expressed in terms of a central free energy distribution for a given, fixed temperature (which is most conveniently chosen as the temperature of the maximum in the heat capacity). The nature of this central free energy function for a given protein thus yields all of the thermodynamic behavior of that protein over the temperature range of the denaturation process.     

Influence of local interactions on protein structure. II. Conformational energy studies of N-acetyl-N′-methylamides of Ala-X and X-Ala dipeptides

Conformational energy calculations using an empirical conformational energy program for peptides (ECEPP) were carried out on 17 N-acetyl-N′-methylamides of Ala-X and X-Ala dipeptides, Where X = Ala, Asn, Asp, Gly, Phe, Ser, Tyr, Val, and Pro. Each dipeptide was found to have many low-energly minima, some of which corresponded to bend structures. The stability of bends was found to depend on the amino acid composition and sequence, with the Ala-X dipeptide generally favoring bends more than the X-Ala dipeptide for a particular X. In bends and nonbends alike, intraresidue interactions dominate over interresidue interactions in determining conformational propeties. The calcutions were shown to be in good agreement with available experimental data.     

Inclusion of solvation free energy with molecular mechanics energy: alanyl dipeptide as a test case.

A combined force field of molecular mechanics and solvation free energy is tested by carrying out energy minimization and molecular dynamics on several conformations of the alanyl dipeptide. Our results are qualitatively consistent with previous experimental and computational studies, in that the addition of solvation energy stabilizes the C5 conformation of the alanyl dipeptide relative to the C7.     

Stability of a model beta-sheet in water.

We have used molecular dynamics simulations to determine the stability in water of a model beta-sheet formed by two alanine dipeptide molecules with two intermolecular hydrogen bonds in the closely spaced antiparallel arrangement. In this paper we describe our computations of the binding free energy of the model sheet and a portion of the free energy surface as a function of a reaction co-ordinate for sheet formation. We used the free energy surface to identify stable conformations along the reaction co-ordinate. To determine whether or not the model sheet with two hydrogen bonds is more stable than a single amide hydrogen bond in water, we compared the results of the present calculations to results from our earlier study of linear hydrogen bond formation between two formamide molecules (the formamide "dimer"). The free energy surfaces for the sheet and formamide dimer each have two minima corresponding to locally stable hydrogen-bonded and solvent-separated configurations. The binding free energies of the model sheet and the formamide dimer are -5.5 and -0.34 kcal/mol, respectively. Thus, the model sheet with two hydrogen bonds is quite stable while the simple amide hydrogen bond is only marginally stable. To understand the relative stabilities of the model sheet and formamide dimer in terms of solute-solute and solute-water interactions, we decomposed the free energy differences between hydrogen-bonded and solvent-separated conformations into energetic and entropic contributions. The changes in the peptide-peptide energy and the entropy are roughly twice as large for the sheet as they are for the formamide dimer. The magnitude of the peptide-water energy difference for the sheet is less than twice (by about 3.5 kcal/mol) that for the formamide dimer, and this accounts for the stability of the sheet. The presence of the side-chains and/or blocking groups apparently prevents the amide groups in the sheet from being solvated as favorably in the separated arrangement as in the formamide dimer, where the amide groups are completely exposed to the solvent.     

Conformations of the glycine dipeptide

Integral equation theory is applied to the determination of the intramolecular potential of mean force for the glycine dipeptide, N-acetyl glycyl-N-methylamide, in aqueous solution. The solvated free energy for the dipeptide as a function of the dihedral angles ? and ψ (Ramachandran plot) is determined and compared with the vacuum surface. Conformations forbidden in vacuum are found to be populated in aqueous solution. The results of the glycine dipeptide are compared to a parallel study on the alanine dipeptide. Solvent effects are found to be responsible for the extent of many of glycines properties related to flexibility.     

Applications of simulated annealing to peptides

We report the application of a new conformation searching algorithm called simulated annealing to the location of the global minimum energy conformation of peptides. Simulated annealing is a Metropolis Monte Carlo approach to conformation generation in which both the energy and temperature dependence of the Boltzmann distribution guides the search for the global minimum. Both uphill and downhill moves are possible, which allows the molecule to escape from local minima. Applications to the 20 natural amino acid "dipeptide models" as well as to polyalanines up to Ala80 are very successful in finding the lowest energy conformation. A history file of the simulated annealing process allows reconstruction and examination of the random walk around conformation space. A separate program, Conf-Gen, reads the history file and extracts all low-energy conformations visited during the run.     

Continuum electrostatic analysis of preferred solvation sites around proteins in solution

To understand water-protein interactions in solution, the electrostatic field is calculated by solving the Poisson-Boltzmann equation, and the free energy surface of water is mapped by translating and rotating an explicit water molecule around the protein. The calculation is applied to T4 lysozyme with data available on the conservation of solvent binding sites in 18 crystallographically independent molecules. The free energy maps around the ordered water sites provide information on the relationship between water positions in crystal structure and in solution. Results show that almost all conserved sites and the majority of nonconserved sites are within 1.3 A of local free energy minima. This finding is in sharp contrast to the behavior of randomly placed water molecules in the boundary layer, which, on the average, must travel more than 3 A to the nearest free energy minimum. Thus, the solvation sites are at least partially determined by protein-water interactions rather than by crystal packing alone. The characteristic water residence times, obtained from the free energies at the local minima, are in good agreement with nuclear magnetic resonance experiments. Only about half of the potential sites show up as ordered water in the 1.7 A resolution X-ray structure. Crystal packing interactions can stabilize weak or mobile potential sites (in fact, some ordered water positions are not close to free energy minima) or can prevent water from occupying certain sites. Apart from a few buried water molecules that are strong binders, the free energies are not very different for conserved and nonconserved sites. We show that conservation of a water site between two crystals occurs if the positions of protein atoms, primarily contributing to the free energy at the local minimum, do not substantially change from one structure to the other. This requirement can be correlated with the nature of the side chain contacting the water molecule in the site.     

Side-chain torsional potentials: effect of dipeptide, protein, and solvent environment.

Side-chain torsional potentials in the bovine pancreatic trypsin inhibitor are calculated from empirical energy functions by use of the known X-ray structure of the protein and the rigid-geometry mapping technique. The potentials are analyzed to determine the roles and relative importance of contributions from the dipeptide backbone, the protein, and the crystalline environment of solvent and other protein molecules. The structural characteristics of the side chains determine two major patterns of energy surfaces, E(X1,X2): a gamma-branched pattern and a pattern for longer, straight side chains (Arg, Lys, Glu, and Met). Most of the dipeptide potential curves and surfaces have a local minimum corresponding to the side-chain torsional angles in the X-ray structure. Addition of the protein forces sharpens and/or selects from these minima, providing very good agreement with the experimental conformation for most side chains at the surface or in the core of the protein. Inclusion of the crystalline environment produces still better results, especially for the side chains extending away from the protein. The results are discussed in terms of the details of the interactions due to the surrounding, calculated solvent-accessibility figures and the temperature factors derived from the crystallographic refinement of the pancreatic trypsin inhibitor.     

A structural and free energy analysis of Ag complexes to five small peptides

The complexes of Ag⁺ with the peptides MetGly, ProGly, GlyPro, GlyHis and GlyProAla were investigated using hybrid density functional theory at the B3LYP/DZVP level. The silver ion binding free energies at 298 K to each of these peptides was calculated to be 60.8, 52.0, 54.3, 71.2 and 63.3 kcal mol⁻¹, respectively. Structural information and relative free energies are presented for several isomers for each of the five complexes. Each of the global minima found for the five complexes is a charge-solvated ion. An important finding is that the Ag⁺-ProGly is the only complex where a salt bridge structure is energetically favored occurring at 4.0 kcal mol⁻¹ higher in free energy than the global minimum. The Ag⁺ ion in this salt bridge structure is attached to the carboxylate anion of zwitterionic ProGly in which the terminal amino nitrogen is protonated. For all the other complexes studied, the salt bridge structure occurs at much higher energies. All the dipeptide complexes with Ag⁺, but one, exhibit a di- or tri-coordinate metal where the sites of attachment are amino and carbonyl groups. However, the highest coordination numbers are not always the global minima due to steric costs. The global minimum of the Ag⁺-GlyProAla complex is the only structure found in this study where the metal is tetra-coordinated, binding to the terminal amino nitrogen and all three carbonyl oxygen atoms. Silver binding to sulphur and imidazole nitrogen atoms of MetGly and GlyHis, respectively, are present in the three most energetically favored species in each of these cases.     

Evaluation of the Adaptive Umbrella Sampling Method

Abstract

The adaptive umbrella sampling technique, introduced recently to improve the probability ratio method and found to perform more reliably than the customary harmonic umbrella sampling, is tested and compared with other free energy methods. One of the tests applies the method to a transition involving a chemical change: calculation of the hydration free energy difference between acetone and dimethylamine and the other test calculates the conformational free energy difference between the C ₇ and α_R conformations of the alanide dipeptide. The dipeptide problem is also treated by two types of thermodynamic integrations and by the perturbation method. The result for the acetone-dimethylamine problem is compared with previous calculations on the same system using the perturbation method, overlap ratio method and finite difference thermodynamic integration. Enhancements to the adaptive umbrella sampling method are also presented.     

The protein folding network

The conformation space of a 20 residue antiparallel beta-sheet peptide, sampled by molecular dynamics simulations, is mapped to a network. Snapshots saved along the trajectory are grouped according to secondary structure into nodes of the network and the transitions between them are links. The conformation space network describes the significant free energy minima and their dynamic connectivity without requiring arbitrarily chosen reaction coordinates. As previously found for the Internet and the World-Wide Web as well as for social and biological networks, the conformation space network is scale-free and contains highly connected hubs like the native state which is the most populated free energy basin. Furthermore, the native basin exhibits a hierarchical organization, which is not found for a random heteropolymer lacking a predominant free-energy minimum. The network topology is used to identify conformations in the folding transition state (TS) ensemble, and provides a basis for understanding the heterogeneity of the TS and denatured state ensemble as well as the existence of multiple pathways.     

Free energy landscape of a biomolecule in dihedral principal component space: sampling convergence and correspondence between structures and minima

Dihedral principal component analysis (dPCA) has recently been developed and shown to display complex features of the free energy landscape of a biomolecule that may be absent in the free energy landscape plotted in principal component space due to mixing of internal and overall rotational motion that can occur in principal component analysis (PCA) [Mu et al., Proteins: Struct Funct Bioinfo 2005;58:45-52]. Another difficulty in the implementation of PCA is sampling convergence, which we address here for both dPCA and PCA using a tetrapeptide as an example. We find that for both methods the sampling convergence can be reached over a similar time. Minima in the free energy landscape in the space of the two largest dihedral principal components often correspond to unique structures, though we also find some distinct minima to correspond to the same structure.     

Empirical solvation models in the context of conformational energy searches: application to bovine pancreatic trypsin inhibitor.

Continuum solvation models that estimate free energies of solvation as a function of solvent accessible surface area are computationally simple enough to be useful for predicting protein conformation. The behavior of three such solvation models has been examined by applying them to the minimization of the conformational energy of bovine pancreatic trypsin inhibitor. The models differ only with regard to how the constants of proportionality between free energy and surface area were derived. Each model was derived by fitting to experimentally measured equilibrium solution properties. For two models, the solution property was free energy of hydration. For the third, the property was NMR coupling constants. The purpose of this study is to determine the effect of applying these solvation models to the nonequilibrium conformations of a protein arising in the course of global searches for conformational energy minima. Two approaches were used: (1) local energy minimization of an ensemble of conformations similar to the equilibrium conformation and (2) global search trajectories using Monte Carlo plus minimization starting from a single conformation similar to the equilibrium conformation. For the two models derived from free energy measurements, it was found that both the global searches and local minimizations yielded conformations more similar to the X-ray crystallographic structures than did searches or local minimizations carried out in the absence of a solvation component of the conformational energy. The model derived from NMR coupling constants behaved similarly to the other models in the context of a global search trajectory. For one of the models derived from measured free energies of hydration, it was found that minimization of an ensemble of near-equilibrium conformations yielded a new ensemble in which the conformation most similar to the X-ray determined structure PTI4 had the lowest total free energy. Despite the simplicity of the continuum solvation models, the final conformation generated in the trajectories for each of the models exhibited some of the characteristics that have been reported for conformations obtained from molecular dynamics simulations in the presence of a bath of explicit water molecules. They have smaller root mean square (rms) deviations from the experimentally determined conformation, fewer incorrect hydrogen bonds, and slightly larger radii of gyration than do conformations derived from search trajectories carried out in the absence of solvent.     

Protein folding transition states probed by loop extension

We propose a new way to characterize protein folding transition states by (1) insertion of one or more residues into an unstructured protein loop, (2) measurement of the effect on protein folding kinetics and thermodynamics, and (3) analysis of the results in terms of a rate-equilibrium free energy relationship, alpha(Loop). alpha(Loop) reports on the fraction of molecules that form the perturbed loop in the transition state. Interpretation of the changes in equilibrium free energy using standard polymer theory can help detect residual structure in the unfolded state. We illustrate our approach with data for the model proteins CI2 and the alpha spectrin SH3 domain.     

Mechanical implications of the domain structure of fiber-forming collagens: comparison of the molecular and fibrillar flexibilities of the alpha1-chains found in types I-III collagen

Fibrillar collagens store, transmit and dissipate elastic energy during tensile deformation. Results of previous studies suggest that the collagen molecule is made up of alternating rigid and flexible domains, and extension of the flexible domains is associated with elastic energy storage. In this study, we model the flexibility of the alpha1-chains found in types I-III collagen molecules and microfibrils in order to understand the molecular basis of elastic energy storage in collagen fibers by analysing the areas under conformational plots for dipeptide sequences. Results of stereochemical modeling suggest that the collagen triple helix is made up of rigid and flexible domains that alternate with periods that are multiples of three amino acid residues. The relative flexibility of dipeptide sequences found in the flexible regions is about a factor of five higher than that found for the flexibility of the rigid regions, and the flexibility of types II and III collagen molecules appears to be higher than that found for the type I collagen molecule. The different collagen alpha1-chains were compared by correlating the flexibilities. The results suggest that the flexibilities of the alpha1-chains of types I and III collagen are more closely related than the flexibilities of the alpha1-chains in types I and II and II and III collagen. The flexible domains found in the alpha1-chains of types I-III collagen were found to be conserved in the microfibril and had periods of about 15 amino acid residues and multiples thereof. The flexibility profiles of types I and II collagen microfibrils were found to be more highly correlated than those for types I and III and II and III. These results suggest that the domain structure of the alpha1-chains found in types I-III collagen is an efficient means for storage of elastic energy during stretching while preserving the triple helical structure of the overall molecule. It is proposed that all collagens that form fibers are designed to act as storage elements for elastic energy. The function of fibers rich in type I collagen is to store and then transmit this energy while fibers rich in types II and III collagen may store and then reflect elastic energy for dissipation through viscous fibrillar slippage. Impaired elastic energy storage by extracellular matrices may lead to cellular damage and changes in signaling by mechanochemical transduction at the extracellular matrix-cell interface.