A comparison of optimal and suboptimal RNA secondary structures predicted by free energy minimization with structures determined by phylogenetic comparison.

This article describes the latest version of an RNA folding algorithm that predicts both optimal and suboptimal solutions based on free energy minimization. A number of RNA's with known structures deduced from comparative sequence analysis are folded to test program performance. The group of solutions obtained for each molecule is analysed to determine how many of the known helixes occur in the optimal solution and in the best suboptimal solution. In most cases, a structure about 80% correct is found with a free energy within 2% of the predicted lowest free energy structure.     

Dynamics of phosphatidylcholine-cholesterol mixed model membranes in the liquid crystalline state.

The effects of cholesterol on the dynamics of cholestane spin probe (CSL) in various phosphatidylcholine-cholesterol mixed model membranes are examined. The lateral diffusion, D of CSL in DMPC/POPC/cholesterol ternary mixtures, is measured utilizing an improved dynamic imaging electron spin resonance method. It shows a factor of two decrease at 10 mol % and 25 degrees C, whereas it shows only a 40% decrease at 50 mol % and 50 degrees C. A comparison with results in POPC/cholesterol mixtures, which show a stronger effect of cholesterol on D, indicates that acyl chain unsaturation leads to stronger self association of cholesterol in PC model membranes. An S2CSL dependence of the activation energy for D, has been confirmed for the DMPC/POPC/cholesterol mixtures. Here SCSL is the order parameter for CSL. A similar correlation of R perpendicular, the perpendicular component of the rotational diffusion coefficient, with SCSL, which is true for all three mixtures (DMPC/cholesterol, POPC/cholesterol, and DMPC/POPC/cholesterol) we have studied, is also found. These are associated with the effects of enhanced local ordering on the free volume needed for translation and reorientation. Such correlations of dynamic properties D and R perpendicular with the thermodynamic quantity S, as well as the consistent interpretations of the effect of acyl chain unsaturation on the dynamics in terms of the activity coefficients, strongly emphasize the interrelation between the dynamic structure and the thermodynamics of the PC/cholesterol mixtures.     

Accurate NMR structures through minimization of an extended hybrid energy

The use of generous distance bounds has been the hallmark of NMR structure determination. However, bounds necessitate the estimation of data quality before the calculation, reduce the information content, introduce human bias, and allow for major errors in the structures. Here, we propose a new rapid structure calculation scheme based on Bayesian analysis. The minimization of an extended energy function, including a new type of distance restraint and a term depending on the data quality, results in an estimation of the data quality in addition to coordinates. This allows for the determination of the optimal weight on the experimental information. The resulting structures are of better quality and closer to the X-ray crystal structure of the same molecule. With the new calculation approach, the analysis of discrepancies from the target distances becomes meaningful. The strategy may be useful in other applications-for example, in homology modeling.     

Configurations of fatty acyl chains in egg phosphatidylcholine-cholesterol mixed bilayers

Based on the structural properties of cholesterol and egg phosphatidylcholine, and on the assumption that the van der Waals' type attactive interaction between the steroid nucleus and the fatty acyl chains provides a stabilizing force for the cholesterol-egg phosphatidylcholine complex, some specific orientation and configurations of the fatty acyl chains around the steroid nucleus in the interacting system are proposed in terms of an optimal packing. The proposed model suggests that the saturated chains are largely facing the flattened (α) surface of the steroid nucleus of cholesterol, while the unsaturated chains can interact with both the α and β surfaces of the steroid nucleus. It is also suggested that the angular methyl groups on the β surface of the steroid nucleus lock the unsaturated fatty acyl chain in a relatively immobile configuration. Experimental evidence which provides support for the proposed stereochemical model is presented.     

Thermodynamics of phosphatidylcholine-cholesterol mixed model membranes in the liquid crystalline state studied by the orientational order parameter.

It is shown that good estimates of the activity of cholesterol in phosphatidylcholine-cholesterol mixed model membranes are obtained by examining the orientational order parameter S of cholestane spin probe (CSL) that is obtained from electron spin resonance by spectral simulation. By introducing thermodynamic stability conditions of liquid mixtures, the variation of activity (or S) as a function of cholesterol mole fraction is utilized to predict the concentration at which the phase separation occurs. These results for DMPC and cholesterol binary mixtures agree very well with those of Tempo-partitioning experiments. The comparison of activity coefficients and the phase boundary in DMPC/cholesterol mixtures with those of POPC/cholesterol mixtures suggests that acyl chain unsaturation leads to poorer mixing of cholesterol in phosphatidylcholine model membranes at higher temperatures (i.e., greater than 35 degrees C). In ternary solutions of DMPC, POPC, and cholesterol, it is found that cholesterol shows less deviation from ideality than in either of the two binary mixtures, and this implies that the phase separation occurs at higher cholesterol concentration than in either of the two binary mixtures. The present analysis suggests that there may not be a critical point in DMPC/cholesterol mixtures, even though phase separation does occur.     

Conformational characteristics of mixed sugar puckered deoxytetranucleoside triphosphate d-GpCpGpC from energy minimization studies

As a continuation of our theoretical studies on nucleic acid subunit systems, in this article we consider the case of the tetranucleoside d-GpCpGpC, the minimally ideal representative unit for analyzing the relative stabilities of different forms of homo- and mixed helical conformation of polynucleotides. The four sugar rings are kept so as to generate B-genus, B+A genus and Z-genus conformations. Twenty five helical conformational states which resulted from judicious mixing of A-, B-, C-, W-, and Z-, states locally are subjected to energy minimization permitting the 19 dihedral angles to vary simultaneously. Conformational states corresponding to regular helical forms and mixed helical forms, when analyzed provide valuable information as to the local conformational flexibility and transitions available to polynucleotides.     

Prediction of RNA secondary structure by free energy minimization

RNA secondary structure is often predicted from sequence by free energy minimization. Over the past two years, advances have been made in the estimation of folding free energy change, the mapping of secondary structure and the implementation of computer programs for structure prediction. The trends in computer program development are: efficient use of experimental mapping of structures to constrain structure prediction; use of statistical mechanics to improve the fidelity of structure prediction; inclusion of pseudoknots in secondary structure prediction; and use of two or more homologous sequences to find a common structure.     

The solubilisation of some steroids by phosphatidylcholine and phosphatidylcholine-cholesterol vesicles.

The solubility of the three steroid hormones, progesterone, testosterone, and estradiol-17 beta in water and phosphatidylcholine vesicles was measured after shaking and ultrasonication. All three steroids have low water solubility, which increases considerably at sonication for testosterone and estradiol-17 beta. The phosphatidylcholine vesicles have a very small solubilising capacity for the steroids; about 20 mumol/mol. This increases at sonication for estradiol-17 beta and decreases for testosterone. The capacity for progesterone is almost unaltered. The incorporation of cholesterol in the vesicles decreased the solubilisation capacity for testosterone and estradiol-17 beta but increased that for progesterone of shaked preparations. For the sonicated systems the cholesterol decreased the solubilising capacity for estradiol-17 beta but increased that for testosterone. The solubilisation experiments indicate that the steroid hormones are solubilised in the hydrocarbon part of the phosphatidylcholine bilayer and also 13CNMR results support this conclusion.     

Prediction of the structure of proteins using related structures,energy minimization and computer graphics

Insight into the functions and interactions of proteins may be gained by correlating a variety of types of experimental data (including kinetics, spectroscopy, biophysical measurements, among others) with three-dimensional structural models displayed and manipulated using interactive computer graphics. Although tertiary structures have been determined for a large number of proteins, one limiting factor in structure-function studies is the lack of availability of the structural coordinates of specific proteins for which other types of detailed experimental data are known. However, as the data base of known structures grows, it becomes more and more likely that the structure of a closely related protein will be available. Here we present a method for predicting structures by ( 1 ) careful alteration of a known structure of a homologous, functionally analogous protein followed by (2) energy minimization to optimize the predicted structure. This method provides a rapid and effective solution to the initial problem of obtaining a working structure for modeling studies.     

Conformational characteristics of mixed sugar puckered deoxydinucleoside triphosphate units d-pCpGp and d-pGpCp from energy minimization studies

The deoxydinucleoside triphosphate units d-pCpGp and d-pGpCp were subjected to a rigorous theoretical investigation with a view to describing their distinctive conformational characteristics. For each unit 216 probable three-dimensional forms defined by the backbone-base dihedral angles and sugar pucker modes were considered for conformational energy minimization process and scrutinized with reference to properties, such as base-stacking, hydrogen-bonding, internal flexibility and base sequence-phosphate influence. The P-O bond torsions and the phosphate groups were treated with special attention. The results reveal a number of preferred conformational states other than the known helical forms, such as, A-, B-, C-, Z-, and Watson-Crick conformation. Many interesting one-step (change in only one of the dihedral angles or sugar puckers) conformational transitions which involve just about a kcal/mol of energy came to light. The two base sequences CG and GC were noted to differ strikingly in many of their conformational characteristics.     

Thermodynamics and dynamics of phosphatidylcholine-cholesterol mixed model membranes in the liquid crystalline state: effects of water.

A method for obtaining the thermodynamic activity of each membrane component in phosphatidylcholine (PC)/cholesterol mixtures, that is based upon ESR spin labeling is examined. The thermodynamic activity coefficients, gamma PC and gamma chol, for the PC and cholesterol, respectively, are obtained from the measured orientational order parameters, SPC and S(chol), as a function of cholesterol content for a spin-labeled PC and the sterol-type cholestane spin probe (CSL), respectively, and the effects of water concentration are also considered. At water content of 24 weight%, the thermodynamics of DMPC/cholesterol/water mixtures in the liquid-crystalline state may be treated as a two-component solution ignoring the water, but at lower water content the role of water is important, especially at lower cholesterol concentrations. At lower water content (17 wt%), gamma chol decreases with increasing cholesterol content which implies aggregation. However, at higher water content (24 wt%), gamma chol is found initially to increase as a function of cholesterol content before decreasing at higher cholesterol content. This implies a favorable accommodation for the cholesterol in the membrane at high water and low cholesterol content. Good thermodynamic consistency according to the Gibbs-Duhem equation was obtained for gamma PC and gamma chol at 24 wt% water. The availability of gamma chol (and gamma PC) as a function of cholesterol concentration permits the estimate of the boundary for phase separation. The rotational diffusion coefficients of the labeled PC and of CSL were also obtained from the ESR spectra. A previously proposed universal relation for the perpendicular component of the rotational diffusion tensor, R perpendicular, for CSL in PC/cholesterol mixtures (i.e., R perpendicular = R0 perpendicular exp(-AS2chol/RT)) is confirmed. A change in composition of cholesterol or of water for DMPC/cholesterol/water mixtures affects R perpendicular only through the dependence of S(chol) on the composition. In particular, the amount of water affects the membrane fluidity, monitored by R perpendicular for CSL, solely by the structural changes it induces in the membrane for the compositions studied. Rotational diffusion for the labeled PC is found to be more complex, most likely due to the combined action of the internal modes of motion of the flexible chain and of the overall molecular reorientation.     

The T<-->R structural transition of insulin; pathways suggested by targeted energy minimization.

The transition of insulin between its crystallographically defined states T and R is connected with considerable change even of backbone structure: the N-terminal B chain (residues B1-B8) refolds from extended conformation in T into helical in R, and vice versa. Although hitherto observed only in hexamers the transition of the monomer was adequate for developing and testing the method of 'targeted energy minimization' (TEM), capable of coping with conformational changes of such extent at moderate computational expenditure. The simulation is performed in a predetermined number of steps consisting of two atomic displacements each, one by force in the direction of the target structure, the second by energy minimization releasing the constraint caused in the first. The transition pathway is represented by the string of energy minimized transient structures. Due to the directedness of the algorithm the simulated pathway for R-->T is not the reversal of that for T-->R. It is, therefore, not pretended that the minimum energy pathway was identified. In the T-->R direction the N-terminal B chain first swivels while remaining largely stretched and then winds up extending the pre-existing helix B9-B19. The A chain advances into the space abandoned and withdraws from it in the R-->T simulation. In the latter the extended helix first kinks at B8/B9, and then the B1-B8 segment is unwound and stretched. The helical H-bonds of that segment are formed late in T-->R and are maintained during almost half of R-->T. The AN helix is less stable and more involved in the transitions than helix AC.(ABSTRACT TRUNCATED AT 250 WORDS)     

Multibilayer structure of dimyristoylphosphatidylcholine dihydrate as determined by energy minimization

Complete energy minimization was carried out on the multibilayer crystal structure of 1,2-dimyristoyl-sn-glycero-3-phosphocholine dihydrate (DMPC.2H2O), starting from the X-ray structure determination reported by Pearson and Pascher (1979) Nature 281, 499-501. The asymmetric unit contains two nonidentical DMPC molecules and four water molecules. Minimization removed the acyl chain disorder present in the X-ray structure and caused the carbon planes of the acyl chains to become mutually parallel. Two energy-minimized structures (structures I and II) were found which mainly differed in the hydrogen-bonding arrangement of the waters of hydration. In structure I as in the X-ray structure, one of the water molecules forms a hydrogen-bonded bridge between successive bilayers; but in structure II, all hydrogen bonds are satisfied on the same bilayer. Structure II corresponds to the global energy minimum and is also a suitable structure for single bilayers. The lattice constants and cell volume of the minimized structures are close to the experimental values. The electrostatic force between DMPC bilayers is attractive. The mean hydration energy of the water is -14.2 kcal/mol, which is 2.5 kcal/mol lower than the binding energy of ice.     

PMR studies of phosphatidylcholine-cholesterol vesicles interacting with lucensomycin.

Spin-lattice relaxation times (T1) were measured above the phase transition temperature on sonicated vesicles of egg- or dipalmitoyl-phosphatidylcholine containing cholesterol and/or the polyenic antibiotic, lucensomycin. T1 values of only the terminal methyl groups of the fatty acyl chains were significantly reduced by cholesterol. Lucensomycin caused, more markedly in cholesterol-containing vesicles, a selective reduction of the T1 values of the N-methyl groups. An even more conspicuous decrease, occurring only in cholesterol-containing vesicles, was observed for the transverse relaxation times (T2) of the N-methyl signals upon addition of lucensomycin. The polyene failed to remove the well-known broadening effect of cholesterol on phosphatidylcholine methylene signals. These results indicate that as lucensomycin binds to cholesterol-containing membranes, there is a detectable perturbation of the dynamic structure of the N-methyl groups with an increase in the degree of motional anisotropy. But the non-polar region of the bilayer seems not significantly perturbed by the polyene.     

Conformational studies of cyclic peptide structures in solution from 1H-Nmr data by distance geometry calculation and restrained energy minimization

The three-dimensional structure of a cyclic bouvardin analogue, cyclo (-Pro-MeTyr-Ala-MeTyr-MeTyr-D-Ala-) has been determined by distance geometry calculation and restrained energy minimization from nmr data. The preparation of the input for the distance geometry calculations, the modification of the amino acid library, and the analysis of the structures were done with the aid of a recently developed software package, GEOM. A great variety of different initial structures were explored to check the uniqueness of the determined solution structure. Calculations with 500 different initial structures and two different strategies led to a uniquely determined backbone conformation with a root mean square deviations value of 0.4 A. The backbone structure consists of two beta-turns, a beta-II turn at Pro1-MeTyr2, and a beta-VI turn at MeTyr4-MeTyr5. The efficiency of the two calculation strategies were compared in order to propose an optimal means for performing distance geometry calculations with cyclic structures.     

Inverted micellar structures in bilayer membranes. Formation rates and half-lives.

Two sorts of inverted micellar structures have previously been proposed to explain morphological and 31P-NMR observations of bilayer systems. These structures only form in systems with components that can adopt the inverse hexagonal (HII) phase. LIP (lipidic particles) are intrabilayer structures, whereas IMI (inverted micellar intermediates) are structures that form between apposed bilayers. Here, we calculate the formation rates and half-lives of these structures to determine which (or if either) of these proposed structures is a likely explanation of the data. Calculations for the egg phosphatidylethanolamine and the Ca+-cardiolipin systems show that IMI form orders of magnitude faster than LIP, which should form slowly, if at all. This result is probably true in general, and indicates that "lipidic particle" electron micrograph images probably represent interbilayer structures, as some have previously proposed. It is shown here that IMI are likely intermediates in the lamellar----HII phase transitions and in the process of membrane fusion in some systems. The calculated formation rates, half-lives, and vesicle-vesicle fusion rates are in agreement with this observation.     

Pseudocontact shifts as constraints for energy minimization and molecular dynamics calculations on solution structures of paramagnetic metalloproteins

The pseudocontact shifts of NMR signals, which arise from the magnetic susceptibility anisotropy of paramagnetic molecules, have been used as structural constraints under the form of a pseudopotential in the SANDER module of the AMBER 4.1 molecular dynamics software package. With this procedure, restrained energy minimization (REM) and restrained molecular dynamics (RMD) calculations can be performed on structural models by using pseudocontact shifts. The structure of the cyanide adduct of the Met80Ala mutant of the yeast iso-1-cytochrome c has been used for successfully testing the calculations. For this protein, a family of structures is available, which was obtained by using NOE and pseudocontact shifts as constraints in a distance geometry program. The structures obtained by REM and RMD calculations with the inclusion of pseudocontact shifts are analyzed. Proteins 29:68–76, 1997. © 1997 Wiley-Liss, Inc.