Molecular Mechanics Calculations of Proteins. Comparison of Different Energy Minimization Strategies

Abstract

A general strategy for performing energy minimization of proteins using the SYBYL molecular modelling program has been developed. The influence of several variables including energy minimization procedure, solvation, dielectric function and dielectric constant have been investigated in order to develop a general method, which is capable of producing high quality protein structures. Avian pancreatic polypeptide (APP) and bovine pancreatic phospholipase A2 (BP PLA2) were selected for the calculations, because high quality X-ray structures exist and because all classes of secondary structure are represented in the structures. The energy minimized structures were evaluated relative to the corresponding X-ray structures. The overall similarity was checked by calculating RMS distances for all atom positions. Backbone conformation was checked by Ramachandran plots and secondary structure elements evaluated by the length on hydrogen bonds. The dimensions of active site in BP PLA2 is very dependent on electrostatic interactions, due to the presence of the positively charged calcium ion. Thus, the distances between calcium and the calcium-coordinating groups were used as a quality index for this protein. Energy minimized structures of the trimeric PLA2 from Indian cobra (N.n.n. PLA2) were used for assessing the impact of protein-protein interactions. Based on the above mentioned criteria, it could be concluded that using the following conditions: Dielectric constant ? = 4 or 20; a distance dependent dielectric function and stepwise energy minimization, it is possible to reproduce X-ray structures very accurately without including explicit solvent molecules.     

Conformation of azidothymidine: an anti-AIDS drug

The nucleoside antibiotic, 3′-azido-3′-deoxythymidine, or simply, azidothymidine has shown great promise in inhibiting the human immuno deficiency virus and in reducing mortality among AIDS patients. Conformational properties of azidothymidine have been investigated by quantum-mechanical PCILO method and compared with those of the parent nucleoside, thymidine. The results indicate great similarity between them and this similarity is remarkably striking in the situations that prevail in aqueous solution. This result has important biological significance in explaining the drug action of azidothymidine.     

Conformational analysis of cytokinins and analogs

The conformational energy surfaces of 12 active cytokinins and analogs are studied with the aid of PCILO quantum mechanical calculations. The resulting conformational energy maps indicate that cytokinin activity is associated with the ability of the above molecules to attain a specific conformation, presumably related to their conformation at the active site of cytokinin receptor(s). The calculations locate the conformational energy minima and describe the flexibility of the studied molecules in terms of conformational barriers and transition paths. An approximate relation is found between cytokinin activity and the values of energy barriers to transitions between certain local minima. According to this relation, active compounds should have rotational barriers within 4–12 kcal/mol, besides the known hitherto constitutional requirements for high physiological activity.     

Quantum chemical studies of methylbenz[a]anthracenes: Metabolism and correlations with carcinogenicity

In the context of the bay region, K-region and radical cation hypotheses for polycyclic aromatic carcinogens, molecular properties were calculated for fourteen methyl derivatives of benz[a]anthracene (BA) related to (1) intrinsic substrate reactivities towards activating and detoxifying metabolism and (2)-the stabilities of the putative carbocation ultimate carcinogens. Allvalence electron methods were used, avoiding the inherent difficulties found in the π-electron methods. The calculated substrate reactivities were found to predict major metabolites sucessfully, supporting the validity of their use in attempted correlations with observed carcinogenic potencies. Positive correlations were found between observed carcinogenic potencies and (1) the reactivities of the parent polycyclic aromatic hydrocarbons (PAH) towards the initial distal bay region epoxidation and (2) the stabilities of the diol epoxide carbocations. This latter correlation holds when both the methyl derivatives of BA and previously studied unsubstituted PAH are considered together, indicating its potential usefulness as a screening parameter for carcinogenic activity.