Simple estimation of absolute free energies for biomolecules

One reason that free energy difference calculations are notoriously difficult in molecular systems is due to insufficient conformational overlap, or similarity, between the two states or systems of interest. The degree of overlap is irrelevant, however, if the absolute free energy of each state can be computed. We present a method for calculating the absolute free energy that employs a simple construction of an exactly computable reference system which possesses high overlap with the state of interest. The approach requires only a physical ensemble of conformations generated via simulation, and an auxiliary calculation of approximately equal central-processing-unit (CPU) cost. Moreover, the calculations can converge to the correct free energy value even when the physical ensemble is incomplete or improperly distributed. As a "proof of principle," we use the approach to correctly predict free energies for test systems where the absolute values can be calculated exactly, and also to predict the conformational equilibrium for leucine dipeptide in implicit solvent.Comment: To appear in J. Chem. Phys., 10 pages, 6 figure

A lattice model with a theta term in three dimensions

We study a three-dimensional abelian lattice model in which the analogue of a theta term can be defined. This term is defined by introducing a neutral scalar field and its effect is to couple magnetic monopoles to the scalar field and vortices to the gauge field. An interesting feature of this model is the presence of an exact duality symmetry that acts on a three parameter space. It is shown that this model has an interesting phase structure for non-zero values of theta. In addition to the usual confinement and vortex phases there are phases in which loops with composite charges condense. The presence of novel point like excitations also alters the physical properties of the system.Comment: 32 pages in latex and three figure

Hydrogen bond networks determine emergent mechanical and thermodynamic properties across a protein family

<p>Abstract</p> <p>Background</p> <p>Gram-negative bacteria use periplasmic-binding proteins (bPBP) to transport nutrients through the periplasm. Despite immense diversity within the recognized substrates, all members of the family share a common fold that includes two domains that are separated by a conserved hinge. The hinge allows the protein to cycle between open (apo) and closed (ligated) conformations. Conformational changes within the proteins depend on a complex interplay of mechanical and thermodynamic response, which is manifested as an increase in thermal stability and decrease of flexibility upon ligand binding.</p> <p>Results</p> <p>We use a distance constraint model (DCM) to quantify the give and take between thermodynamic stability and mechanical flexibility across the bPBP family. Quantitative stability/flexibility relationships (QSFR) are readily evaluated because the DCM links mechanical and thermodynamic properties. We have previously demonstrated that QSFR is moderately conserved across a mesophilic/thermophilic RNase H pair, whereas the observed variance indicated that different enthalpy-entropy mechanisms allow similar mechanical response at their respective melting temperatures. Our predictions of heat capacity and free energy show marked diversity across the bPBP family. While backbone flexibility metrics are mostly conserved, cooperativity correlation (long-range couplings) also demonstrate considerable amount of variation. Upon ligand removal, heat capacity, melting point, and mechanical rigidity are, as expected, lowered. Nevertheless, significant differences are found in molecular cooperativity correlations that can be explained by the detailed nature of the hydrogen bond network.</p> <p>Conclusion</p> <p>Non-trivial mechanical and thermodynamic variation across the family is explained by differences within the underlying H-bond networks. The mechanism is simple; variation within the H-bond networks result in altered mechanical linkage properties that directly affect intrinsic flexibility. Moreover, varying numbers of H-bonds and their strengths control the likelihood for energetic fluctuations as H-bonds break and reform, thus directly affecting thermodynamic properties. Consequently, these results demonstrate how unexpected large differences, especially within cooperativity correlation, emerge from subtle differences within the underlying H-bond network. This inference is consistent with well-known results that show allosteric response within a family generally varies significantly. Identifying the hydrogen bond network as a critical determining factor for these large variances may lead to new methods that can predict such effects.</p