The hydroxyl orientations in 31 dioctahedral and trioctahedral 2:1 phyllosilicate structures have been determined by electrostatic energy calculations. These structures included micas, brittle micas, and other related minerals exhibiting ordered as well as disordered cation distributions. The dioctahedral micas and brittle micas were examined with and without the interlayer cation. A range of orientations from 1.3° to 183.3° (the angle ρ between the O-H and (001) measured with respect to the M1 site) were found. The orientations for the dioctahedral structures represent a continuum of values whereas the trioctahedral species exhibit two possible orientations separated by an energy barrier. One orientation is near 90° the other is near 180°. The latter orientation results from a concentration of charge on the interlayer (IC) and tetrahedral (T) sites at the expense of the octahedral (M) sites. A multiple regression analysis of all 31 structures, using as predictors the a and b cell parameters, d001, and the charges for T, IC, M1, and M2 sites, was performed. This analysis indicated that the important factors are the charges for IC, T, and M2 sites. When treated as a separate group, one finds the same factors for the dioctahedral structures. The trioctahedral orientations are determined by the charge on the M2 site and the amount of tetrahedral rotation. Using these two predictor equations, the value of ρ can be estimated with a standard deviation of 4.7° and 2.9° for the dioctahedral and trioctahedral cases, respectively.

Many molecular events are associated with small or large conformational changes occurring in the corresponding proteins. Modeling such changes is a challenge and requires significant amount of computing time. From point of view of electrostatics, these changes can be viewed as a reorganization of local charges and dipoles in response to the changes of the electrostatic field, if the cause is insertion or deletion of a charged amino acid. Here we report a large scale investigation of modeling the changes of the folding energy due to single mutations involving charged group. This allows the changes of the folding energy to be considered mostly electrostatics in origin and to be calculated with DelPhi assigning residue-specific value of the internal dielectric constant of protein. The predicted energy changes are benchmarked against experimentally measured changes of the folding energy on a set of 257 single mutations. The best fit between experimental values and predicted changes is used to find out the effective value of the internal dielectric constant for each type of amino acid. The predicted folding free energy changes with the optimal, amino acid specific, dielectric constants are within RMSD=0.86 kcal/mol from experimentally measured changes.