Conformational equilibria of valine studied by dynamics simulation.

The conformational probability distribution of a valine residue in the valine dipeptide and of the valine side chain in an alpha-helix, as well as the change in helix stability for replacing alanine with valine, has been calculated by molecular dynamics simulations of explicitly hydrated systems: dipeptide, tetrapeptide and 10-, 14- and 18-residue oligoalanine helices. All computed free-energy differences are means from at least eight separate slow-growth simulations, four in each direction and are reported with their root-mean-square deviations. Different values for the change in free energy of folding (delta delta G degrees) have been calculated with the use of forcefields having an all-atom and a central-atom representation of methyl groups, etc. The value obtained with the all-atom forcefield agrees well with new experimental values (3 kJ/mol = 0.7 kcal/mol). Furthermore, the most stable valine side-chain rotamer in the helix is different for these two representations. The most stable rotamer for the all atom conformation is the same one that predominates for valines in alpha-helices in proteins of known conformation. The lower conformational freedom of the valine side chain in the helix contributes 1 kJ/mol to the difference in stability computed with the all-atom potential; unfavorable interactions of the side chain with helix, even in the most stable conformation, further increase delta delta G degrees.     

Proline in alpha-helix: stability and conformation studied by dynamics simulation

Free-energy simulations have been used to estimate the change in the conformational stability of short polyalanine alpha-helices when one of the alanines is replaced by a proline residue. For substituting proline in the middle of the helix the change in free energy of folding (delta delta G degrees) was calculated as 14 kJ/mol (3.4 kcal/mol), in excellent agreement with the one available experimental value. The helix containing proline was found to be strongly kinked; the free energy for reducing the angle of the kink from 40 degrees to 15 degrees was calculated, and found to be small. A tendency to alternate hydrogen bonding schemes was observed in the proline-containing helix. These observations for the oligopeptide agree well with the observation of a range of kink angles (18-35 degrees) and variety of hydrogen bonding schemes, in the rare instances where proline occurs in helices in globular proteins. For substituting proline at the N-terminus of the helix the change in free energy of folding (delta delta G degrees) was calculated as -4 kJ/mol in the first helical position (N1) and +6 kJ/mol in the second helical position (N2). The observed frequent occurrence of proline in position N1 in alpha-helices in proteins therefore has its origin in stability differences of secondary structure. The conclusion reached here that proline may be a better helix former in position N1 than (even) alanine, and thus be a helix initiator may be testable experimentally by measurements of fraction helical conformation of individual residues in oligopeptides of appropriate sequence. The relevance of these results in regards to the frequent occurrence of proline-containing helices in certain membrane proteins is discussed.     

Destabilizing effect of proline substitutions in two helical regions of T4 lysozyme: leucine 66 to proline and leucine 91 to proline.

A class of temperature-sensitive (ts) mutants of T4 lysozyme with reduced activity at 30 degrees C and no activity at 43 degrees C has been selected. These mutants, designated "tight" ts mutants, differ from most other T4 lysozyme mutants that are active at 43 degrees C, but only manifest their ts lesion by a reduced halo size around phage plaques after exposure of the growth plates to chloroform vapors. For example, in the series of T4 lysozyme mutants at position 157, the original randomly selected mutant, T1571, is the least stable of the series, yet, apart from the halo assay and subsequent in vitro protein stability measurements, this mutant is indistinguishable from wild type (WT) even at 43 degrees C. Two mutants were identified: L91P and L66P. Both insert proline residues into alpha-helical regions of the WT protein structure. The stabilities (delta delta G) as determined by urea denaturation are 8.2 kcal/mol for L91P and 7.1 kcal/mol for L66P. CD spectra indicate that no major conformational changes have occurred in the mutant structures. The structures of the mutants were modeled with a 40-ps molecular dynamics simulation using explicit solvent. For L91P, the reduction of stability appears to be due to an unsatisfied hydrogen bond in the alpha-helix and to a new buried cavity. For L66P, the reduction of stability appears to be due to a disruption of the interdomain alpha-helix, at least two unsatisfied hydrogen bonds, and a newly formed solvent-filled pocket that protrudes into the hydrophobic core, possibly reducing the stabilizing contribution of a partially buried intrachain salt bridge.     

Conformational stability and mechanism of folding of ribonuclease T1

Urea and thermal unfolding curves for ribonuclease T1 (RNase T1) were determined by measuring several different physical properties. In all cases, steep, single-step unfolding curves were observed. When these results were analyzed by assuming a two-state folding mechanism, the plots of fraction unfolded protein versus denaturant were coincident. The dependence of the free energy of unfolding, delta G (in kcal/mol), on urea concentration is given by delta G = 5.6 - 1.21 (urea). The parameters characterizing the thermodynamics of unfolding are: midpoint of the thermal unfolding curve, Tm = 48.1 degrees C, enthalpy change at Tm, delta Hm = 97 kcal/mol, and heat capacity change, delta Cp = 1650 cal/mol deg. A single kinetic phase was observed for both the folding and unfolding of RNase T1 in the transition and post-transition regions. However, two slow kinetic phases were observed during folding in the pre-transition region. These two slow phases account for about 90% of the observed amplitude, indicating that a faster kinetic phase is also present. The slow phases probably result from cis-trans isomerization at the 2 proline residues that have a cis configuration in folded RNase T1. These results suggest that RNase T1 folds by a highly cooperative mechanism with no structural intermediates once the proline residues have assumed their correct isomeric configuration. At 25 degrees C, the folded conformation is more stable than the unfolded conformations by 5.6 kcal/mol at pH 7 and by 8.9 kcal/mol at pH 5, which is the pH of maximum stability. At pH 7, the thermodynamic data indicate that the maximum conformational stability of 8.3 kcal/mol will occur at -6 degrees C.     

Enhanced thermostability of methyl parathion hydrolase from Ochrobactrum sp. M231 by rational engineering of a glycine to proline mutation

Protein thermostability can be increased by some glycine to proline mutations in a target protein. However, not all glycine to proline mutations can improve protein thermostability, and this method is suitable only at carefully selected mutation sites that can accommodate structural stabilization. In this study, homology modeling and molecular dynamics simulations were used to select appropriate glycine to proline mutations to improve protein thermostability, and the effect of the selected mutations was proved by the experiments. The structure of methyl parathion hydrolase (MPH) from Ochrobactrum sp. M231 (Ochr-MPH) was constructed by homology modeling, and molecular dynamics simulations were performed on the modeled structure. A profile of the root mean square fluctuations of Ochr-MPH was calculated at the nanosecond timescale, and an eight-amino acid loop region (residues 186-193) was identified as having high conformational fluctuation. The two glycines nearest to this region were selected as mutation targets that might affect protein flexibility in the vicinity. The structures and conformational fluctuations of two single mutants (G194P and G198P) and one double mutant (G194P/G198P) were modeled and analyzed using molecular dynamics simulations. The results predicted that the mutant G194P had the decreased conformational fluctuation in the loop region and might increase the thermostability of Ochr-MPH. The thermostability and kinetic behavior of the wild-type and three mutant enzymes were measured. The results were consistent with the computational predictions, and the mutant G194P was found to have higher thermostability than the wild-type enzyme.     

Structural stability of beta-lactoglobulin in the presence of kosmotropic salts. A kinetic and thermodynamic study

The thiol group of beta-lactoglobulin reacted very sluggishly with dithio-bis-nitro-benzoic acid as compared to that of glutathione at pH 6.85. The pKapp value of the thiol group of the protein was 9.35. In the presence of 3 M urea, the thiol group reacted completely with dithio-bis-nitrobenzoic acid at pH 6.85. Heating (from 50 degrees to 80 degrees) increased the exposure of the thiol by dissociating the dimer unit. From the pseudo-first order rate constants of heat-exposure of thiol, thermodynamic activation parameters, delta G++, delta H++, and delta S++, for the heat-dissociation of beta-lactoglobulin dimer were estimated to be 23,290 cal/mol, 31,160 cal/mol, and 22.9 e.u. (at 70 degrees), respectively. Addition of kosmotropic salts, chloride, tartrate, sulfate, phosphate, and citrate (0.2 M) decreased the heat-induced exposure of the thiol group (at 70 degrees), probably by decreasing the dissociation of the dimer at pH 6.85. The relative change in free energy of activation for the dissociation of the dimer, delta(delta G++dimer), in the presence of the salts was positive, suggesting that these additives increase the stability of the dimer against heat. These salts also increased the conformational stability of beta-lactoglobulin as revealed by an increase in -delta(delta G0conf) values in their presence. Both delta(delta G++dimer) and -delta(delta G0conf) values followed the order, chloride less than tartrate less than sulfate less than phosphate less than citrate. These salts seem to manifest their structure-stabilizing effect by increasing both inter- and intramolecular hydrophobic interactions via changes in structure of water.     

One gene, two diseases and three conformations: molecular dynamics simulations of mutants of human prion protein at room temperature and elevated temperatures

Fatal familial insomnia (FFI) and Creutzfeldt-Jakob disease (CJD) are associated to the same mutation at codon 178 but differentiate into clinicopathologically distinct diseases determined by this mutation and a naturally occurring methionine-valine polymorphism at codon 129 of the prion protein gene. It has been suggested that the clinical and pathological difference between FFI and CJD is caused by different conformations of the prion protein. Using molecular dynamics (MD), we investigated the effect of the mutation at codon 178 and the polymorphism at codon 129 on prion protein dynamics and conformation at normal and elevated temperatures. Four model structures were examined with a focus on their dynamics and conformational changes. The results showed differences in stability and dynamics between polymorphic variants. Methionine variants demonstrated a higher stability than valine variants. Elongation of existing beta-sheets and formation of new beta-sheets was found to occur more readily in valine polymorphic variants. We also discovered the inhibitory effect of proline residue on existing beta-sheet elongation.     

Thermodynamics of ribonuclease T1 denaturation.

Differential scanning calorimetry has been used to investigate the thermodynamics of denaturation of ribonuclease T1 as a function of pH over the pH range 2-10, and as a function of NaCl and MgCl2 concentration. At pH 7 in 30 mM PIPES buffer, the thermodynamic parameters are as follows: melting temperature, T1/2 = 48.9 +/- 0.1 degrees C; enthalpy change, delta H = 95.5 +/- 0.9 kcal mol-1; heat capacity change, delta Cp = 1.59 kcal mol-1 K-1; free energy change at 25 degrees C, delta G degrees (25 degrees C) = 5.6 kcal mol-1. Both T1/2 = 56.5 degrees C and delta H = 106.1 kcal mol-1 are maximal near pH 5. The conformational stability of ribonuclease T1 is increased by 3.0 kcal/mol in the presence of 0.6 M NaCl or 0.3 M MgCl2. This stabilization results mainly from the preferential binding of cations to the folded conformation of the protein. The estimates of the conformational stability of ribonuclease T1 from differential scanning calorimetry are shown to be in remarkably good agreement with estimates derived from an analysis of urea denaturation curves.     

Molecular dynamics simulations of a bacterial autotransporter: NalP from Neisseria meningitidis

NalP is an autotransporter secretory protein found in the outer membrane of Neisseria meningitidis. The crystal structure of the NalP translocator domain revealed a transmembrane beta-barrel containing a central alpha-helix. The role of this alpha-helix, and of the conformational dynamics of the beta-barrel pore have been studied via atomistic molecular dynamics simulations. Three simulations, each of 10 ns duration, of NalP embedded within a solvated DMPC bilayer were performed. The helix was removed from the barrel interior in one simulation. The conformational stability of the protein is similar to that of other outer membrane proteins, e.g., OmpA, in comparable simulations. The transmembrane beta-barrel is stable even in the absence of the alpha-helix. Removal of the helix results in an influx of water into the pore region, suggesting the helix acts as a 'plug'. Water molecules entering the resultant pore form hydrogen bonds with the barrel lining that compensate for the loss of helix-barrel hydrogen bonds. The dimensions of the pore fluctuate over the course of the simulation revealing it to be flexible, but only wide enough to allow transport of the passenger domain in an unfolded or extended conformation. The simulations help us to understand the role of the central helix in plugging the pore and in maintaining the width of the barrel, and show that the NalP monomer is sufficient for the transport of the passenger domain in an unfolded or extended conformation.     

Cis proline mutants of ribonuclease A. I. Thermal stability.

A chemically synthesized gene for ribonuclease A has been expressed in Escherichia coli using a T7 expression system (Studier, F.W., Rosenberg, A.H., Dunn, J.J., & Dubendorff, J.W., 1990, Methods Enzymol. 185, 60-89). The expressed protein, which contains an additional N-terminal methionine residue, has physical and catalytic properties close to those of bovine ribonuclease A. The expressed protein accumulates in inclusion bodies and has scrambled disulfide bonds; the native disulfide bonds are regenerated during purification. Site-directed mutations have been made at each of the two cis proline residues, 93 and 114, and a double mutant has been made. In contrast to results reported for replacement of trans proline residues, replacement of either cis proline is strongly destabilizing. Thermal unfolding experiments on four single mutants give delta Tm approximately equal to 10 degrees C and delta delta G0 (apparent) = 2-3 kcal/mol. The reason is that either the substituted amino acid goes in cis, and cis<==>trans isomerization after unfolding pulls the unfolding equilibrium toward the unfolded state, or else there is a conformational change, which by itself is destabilizing relative to the wild-type conformation, that allows the substituted amino acid to form a trans peptide bond.     

Mutational analysis of invariant valine B12 in insulin: implications for receptor binding

An invariant residue, valine B12, is part of the insulin B-chain central alpha-helix (B9-B19), and its aliphatic side chain lies at the surface of the hydrophobic core of the insulin monomer in close contact with the neighboring aromatic side chains of phenylalanines (B24 and B25) and tyrosines (B26 and B16). This surface contributes to the dimerization of insulin, maintains the active conformation of the insulin monomer, and has been suspected to be directly involved in receptor recognition. To investigate in detail the role of the B12 residue in insulin-receptor interactions, we have synthesized nine analogues bearing natural or unnatural amino acid replacements for valine B12 by chemical synthesis of modified insulin B-chains and the subsequent combination of each synthetic B-chain with natural insulin A-chain. The receptor binding potencies of the synthetic B12 analogues relative to porcine insulin were determined by use of isolated canine hepatocytes, and the following results were obtained: isoleucine, 13%; allo-isoleucine, 77%; tert-leucine, 107%; cyclopropylglycine, 43%; threonine, 5.4%; D-valine, 3.4%; alpha-amino-n-butyric acid, 14%; alanine, 1.0%; and glycine, 0.32%. Selected analogues were also analyzed by far-UV circular dichroic spectroscopy and by absorption spectroscopy of their complexes with Co(2+). Our results indicate that beta-branched aliphatic amino acids are generally tolerated at the B12 position with specific steric preferences and that the receptor binding potencies of these analogues correlate with their abilities to form dimers. Furthermore, the structure-activity relationships of valine B12 are quite similar to those of valine A3, suggesting that valine residues at both A3 and B12 contribute to the insulin-receptor interactions in a similar manner.     

Conformational behavior of temporin A and temporin L in aqueous solution: a computational/experimental study

Molecular dynamics (MD) simulations and circular dichroism (CD) experiments were carried out on aqueous temporin A and L, two short peptides belonging to an interesting class of natural substances known to be active mainly against Gram-positive/negative bacteria and fungi. Experimental results indicate the higher propensity of temporin L, with respect to temporin A, in forming alpha-helical structures. These results were revisited by long-timescale MD simulations, in which their alpha-helical propensity was investigated in the absence of trifluoroethanol. Results clearly show the higher stability of alpha-helix conformations in temporin L; moreover, an interestingly strong mechanical analogy emerges since both temporins show the same residue interval (from 7 to 10) as the most energetically accessible for alpha-helix formation. Such studies provide some intriguing structural and mechanical evidence that may help in better understanding and rationalizing the conformational behaviour of temporins in water solution and, ultimately, the inner principles of their microbial targets selectivity and mechanism of action at the level of cell membranes.     

Molecular dynamics simulations of the whey protein beta-lactoglobulin.

Molecular dynamics simulations have been used to model the motions and conformational behavior of the whey protein bovine beta-lactoglobulin. Simulations were performed for the protein by itself and complexed to a single retinol ligand located in a putative interior binding pocket. In the absence of the retinol ligand, the backbone loops around the opening of this interior pocket shifted inward to partially close off this cavity, similar to the shifts observed in previously reported molecular dynamics simulations of the uncomplexed form of the homologous retinol binding protein. The protein complexed with retinol does not exhibit the same conformational shifts. Conformational changes of this type could serve as a recognition signal allowing in vivo discrimination between the free and retinol complexed forms of the beta-lactoglobulin molecule. The unusual bending of the single alpha-helix observed in the simulations of retinol binding protein were not observed in the present calculations.     

Tolerance of T4 lysozyme to proline substitutions within the long interdomain alpha-helix illustrates the adaptability of proteins to potentially destabilizing lesions.

To investigate the ability of a protein to accommodate potentially destabilizing amino acid substitutions, and also to investigate the steric requirements for catalysis, proline was substituted at different sites within the long alpha-helix that connects the amino-terminal and carboxyl-terminal domains of T4 lysozyme. Of the four substitutions attempted, three yielded folded, functional proteins. The catalytic activities of these three mutant proteins (Q69P, D72P, and A74P) were 60-90% that of wild-type. Their melting temperatures were 7-12 degrees C less than that of wild-type at pH 6.5. Mutant D72P formed crystals isomorphous with wild-type allowing the structure to be determined at high resolution. In the crystal structure of wild-type lysozyme the interdomain alpha-helix has an overall bend angle of 8.5 degrees. In the mutant structure the introduction of the proline causes this bend angle to increase to 14 degrees and also causes a corresponding rotation of 5.5 degrees of carboxyl-terminal domain relative to the amino-terminal one. Except for the immediate location of the proline substitution there is very little change in the geometry of the interdomain alpha-helix. The results support the view that protein structures are adaptable and can compensate for potentially destabilizing amino acid substitutions. The results also suggest that the precise shape of the active site cleft of T4 lysozyme is not critical for catalysis.     

13C multiplet nuclear magnetic resonance relaxation-derived ring puckering and backbone dynamics in proline-containing glycine-based peptides.

13CH2-multiplet nuclear magnetic resonance relaxation studies on proline (P)-containing glycine (G)-based peptides, GP, PG, GPG, PGG, and GPGG, provided numerous dipolar auto- and cross-correlation times for various motional model analyses of backbone and proline-ring bond rotations. Molecular dynamics simulations and bond rotation energy profiles were calculated to assess which motions could contribute most to observed relaxation phenomena. Results indicate that proline restricts backbone psi 1, psi 2, and phi 2 motions by 50% relative to those found for a polyglycine control peptide. psi 1 rotations are more restricted in the trans-proline isomer state than in the cis form. A two-state jump model best approximates proline ring puckering which in water could occur either by the C gamma endo-exo or by the C2 interconversion mechanism. The temperature dependence (5 degrees to 75 degrees C) of C beta, and C gamma, and C delta angular changes is rather flat, suggesting a near zero enthalpic contribution to the ring puckering process. In lower dielectric solvents, dimethylsulfoxide and methanol, which may mimic the hydrophobic environment within a protein, the endo-exo mechanism is preferred.     

Florida citrus molasses as a fermentation substrate: I. Free amino acids in molasses from early and midseason varieties of citrus fruits during 2 years of production

Samples of midseason Florida citrus molasses from 10 different manufacturing plants were studied during 2 years of production for suitability as fermentation substrates. They were analyzed for degrees Brix, pH, percentage of protein, and free amino acid content. The following ranges of values were obtained: degrees Brix, 65.4 to 73.0; pH, 4.2 to 6.4; and protein, 3.3 to 5.1%. Although the general characteristics of all samples were quite similar, the total free amino acid content and the number of amino acids differed from year to year and from plant to plant. The free amino acids consistently present were asparagine, aspartic acid, alanine, γ-amino-butyric, glycine, glutamic acid, proline, serine, and threonine. Occasionally found were arginine, histidine, isoleucine, leucine, lysine, phenylalanine, tyrosine, and valine. An extremely high level of proline was detected in all samples.     

Molecular dynamics simulations of the unfolding of beta(2)-microglobulin and its variants

In this study, we examined the unfolding processes of native beta(2)-microglobulin and two related variants, one with an N-terminal hexapeptide deletion DeltaN6 and another with Lys57-Asp58 cleavage, by high-temperature molecular dynamics simulations. Three simulation models were used, molecular dynamics (MD) simulations with explicit water solvation, MD simulations with the CHARMM EEF1 force field and Langevin dynamics with the CHARMM EEF1 force field. Our simulations reproduce many of the experimentally observed structural changes. The most striking agreement is in the beta-strands to alpha-helix transition. In our simulations, strands beta(3), beta(4) and beta(5) consistently change to alpha-helix, whereas beta(8) changes to an alpha-helix only briefly. Through comparisons of the conformational behavior of the native, the DeltaN6 and the Lys-cut beta(2)-m, using the three simulation methods, we identified the consensus conformational changes that differentiate between the native beta(2)-m and its two variants. We found that the main effect of the removal of the N-terminal hexapeptide is to increase the separation between strands beta(2) and beta(6) and to facilitate the beta to alpha transition. On the other hand, the lysine cleavage only increases the flexibility of strand beta(5) and does not affect the interactions between strands beta(2) and beta(6). These conformational changes may relate to polymerization tendencies of these variants.     

The presence of high concentrations of free D-amino acids in human saliva

Free neutral D-amino acids have previously been detected in human plasma, usually accounting for less than 2% of the total free amino acid concentration (D-amino acid ratio) [Nagata, Y., Masui, R., Akino, T., 1992a. The presence of free D-serine, D-alanine and D-proline in human plasma. Experientia 48, 986-988. Nagata, Y., Yamamoto, K., Shimojo, T., 1992b. Determination of D- and L-amino acids in mouse kidney by high-performance liquid chromatography. Journal of Chromatography 575, 147-152. Nagata, Y., Yamamoto, K., Shimojo, T., Konno, R., Yasumura, Y., Akino, T., 1992c. The presence of free D-alanine, D-proline and D-serine in mice. Biochimca et Biiophysica Acta 1115, 208-211]. In the present study to search for the source of free D-amino acids, D- and L-enantiomers of the major non-essential amino acids, i.e., the free form of serine, alanine, proline, aspartate and glutamate were analyzed by HPLC in human saliva, submandibular glands and oral epithelial cells. The D-enantiomer ratios to total of free alanine or proline were 35% and 20%, respectively, in saliva. The ratios of the other D-amino acids were less than 11%. The effect of ingested food and oral bacteria on the saliva amino acid levels was suggested to be insignificant. D-Alanine and d-aspartate were also detected in the submandibular gland in ratios up to 5%, and D-alanine and d-proline were found in oral epithelial cells in ratios of 18% and 5%, respectively. The submandibular gland and oral epithelial cells are suggested to be possible sources of the saliva D-alanine and D-aspartate.     

A statistical approach to the interpretation of molecular dynamics simulations of calmodulin equilibrium dynamics

A sample of 35 independent molecular dynamics (MD) simulations of calmodulin (CaM) equilibrium dynamics was prepared from different but equally plausible initial conditions (20 simulations of the wild-type protein and 15 simulations of the D129N mutant). CaM's radius of gyration and backbone mean-square fluctuations were analyzed for the effect of the D129N mutation, and simulations were compared with experiments. Statistical tests were employed for quantitative comparisons at the desired error level. The computational model predicted statistically significant compaction of CaM relative to the crystal structure, consistent with the results of small-angle X-ray scattering (SAXS) experiments. This effect was not observed in several previously reported studies of (Ca2+)(4)-CaM, which relied on a single MD run. In contrast to radius of gyration, backbone mean-square fluctuations showed a distinctly non-normal and positively skewed distribution for nearly all residues. Furthermore, the D129N mutation affected the backbone dynamics in a complex manner and reduced the mobility of Glu123, Met124, Ile125, Arg126, and Glu127 located in the adjacent alpha-helix G. The implications of these observations for the comparisons of MD simulations with experiments are discussed. The proposed approach may be useful in studies of protein equilibrium dynamics where MD simulations fall short of properly sampling the conformational space, and when the comparison with experiments is affected by the reproducibility of the computational model.     

Effects of base substitutions in an RNA hairpin from molecular dynamics and free energy simulations

Contributions of individual interactions in the GGCGCAAGCC hairpin containing a GCAA tetraloop were studied by computer simulations using base substitutions. The G in the first tetraloop position was replaced by inosine (I) or adenosine (A), and the G in the C-G basepair closing the tetraloop was replaced by I. These substitutions eliminate particular hydrogen bonds proposed in the nuclear magnetic resonance model of the GCAA tetraloop. Molecular dynamics simulations of the GCAA tetraloop in aqueous solvent displayed a well-defined hydrogen pattern between the first and last loop nucleotides (G and A) stabilized by a bridging water molecule. Substitution of G-->I in the basepair closing the tetraloop did not significantly influence the loop structure and dynamics. The ICAA loop maintained the overall structure, but displayed variation in the hydrogen-bond network within the tetraloop itself. Molecular dynamics simulations of the ACAA loop led to conformational heterogeneity of the resulting structures. Changes of hairpin formation free energy associated with substitutions of individual bases were calculated by the free energy perturbation method. The calculated decrease of the hairpin stability upon G-->I substitution in the C-G basepair closing the tetraloop was in good agreement with experimental thermodynamic data. Our theoretical estimates for G-->I and G-->A mutations located in the tetraloop suggest larger loop destabilization than corresponding experimental results. The extent of conformational sampling of the structures resulting from base substitutions and its impact on the calculated free energy was discussed.