Classical harmonic model for the behavior of pure fluids at the critical point.

We present a simple model to explain the limiting behavior of many thermodynamic properties of pure classical fluids as one approaches the critical point. The model consists of four main aspects: (i) a focus on the thermodynamic properties as described by Fluctuation Solution Theory (FST); (ii) the removal of all intramolecular energy contributions from the FST expressions; (iii) particle fluctuations (and not the energy fluctuations) dominate the thermodynamics close to the critical point; and (iv) the equipartition theorem applies to the resulting low frequency collective modes of the fluid at the critical point. Using this approach, we predict and compare the thermodynamic results for 121 pure fluids as provided by accurate equations of state. The model explains why some ratios of diverging thermodynamic properties remain finite at the critical point, and reliably predicts the value of these ratios for equations of state that mimic both mean field and Ising system types, suggesting universal behavior.

A Kirkwood-Buff derived force field for alkaline earth halide salts.

The activity and function of many macromolecules in cellular environments are coupled with the binding of divalent ions such as calcium or magnesium. In principle, computer simulations can be used to understand the molecular level aspects of how many important macromolecules interact with ions. However, most of the force fields currently available often fail to accurately reproduce the properties of divalent ions in aqueous environments. Here we develop classical non-polarizable force fields for the aqueous alkaline earth metal halides (MX2), where M = Mg2+, Ca2+, Sr2+, Ba2+ and X = Cl-, Br-, I-, which can be used in bimolecular simulations and which are compatible with the Simple Point Charge/Extended (SPC/E) water model. The force field parameters are specifically developed to reproduce the experimental Kirkwood-Buff integrals for aqueous solutions and thereby the experimental activity derivatives, partial molar volumes, and excess coordination numbers. This ensures that a reasonable balance between ion-ion, ion-water, and water-water distributions is obtained. However, this requires a scaling of the cation to water oxygen interaction strength in order to accurately reproduce the integrals. The scaling factors developed for chloride salts are successfully transferable to the bromide and iodide salts. Use of these new models leads to reasonable diffusion constants and dielectric decrements. However, the performance of the models decreases with increasing salt concentration (>4m), and simulations of the pure crystals exhibited unstable behavior.

Experimental Investigation of Triplet Correlation Approximations for Fluid Water.

Triplet correlations play a central role in our understanding of fluids and their properties. Of particular interest is the relationship between the pair and triplet correlations. Here we use a combination of Fluctuation Solution Theory and experimental pair radial distribution functions to investigate the accuracy of the Kirkwood Superposition Approximation (KSA), as given by integrals over the relevant pair and triplet correlation functions, at a series of state points for pure water using only experimental quantities. The KSA performs poorly, in agreement with a variety of other studies. Several additional approximate relationships between the pair and triplet correlations in fluids are also investigated and generally provide good agreement for the fluid thermodynamics for regions of the phase diagram where the compressibility is small. A simple power law relationship between the pair and triplet fluctuations is particularly successful for state points displaying low to moderately high compressibilities.

How Osmolytes Counteract Pressure Denaturation on a Molecular Scale.

Life in the deep sea exposes enzymes to high hydrostatic pressure, which decreases their stability. For survival, deep sea organisms tend to accumulate various osmolytes, most notably trimethylamine N-oxide used by fish, to counteract pressure denaturation. However, exactly how these osmolytes work remains unclear. Here, a rigorous statistical thermodynamics approach is used to clarify the mechanism of osmoprotection. It is shown that the weak, nonspecific, and dynamic interactions of water and osmolytes with proteins can be characterized only statistically, and that the competition between protein-osmolyte and protein-water interactions is crucial in determining conformational stability. Osmoprotection is driven by a stronger exclusion of osmolytes from the denatured protein than from the native conformation, and water distribution has no significant effect on these changes at low osmolyte concentrations.

Rationally designed peptide nanosponges for cell-based cancer therapy.

A novel type of supramolecular aggregate, named a "nanosponge" was synthesized through the interaction of novel supramolecular building blocks with trigonal geometry. The cholesterol-(K/D)nDEVDGC)3-trimaleimide unit consists of a trigonal maleimide linker to which homopeptides (either K or D) of variable lengths (n=5, 10, 15, 20) and a consensus sequence for executioner caspases (DEVDGC) are added via Michael addition. Upon mixing in aqueous buffer cholesterol-(K)nDEVDGC)3-trimaleimides and a 1:1 mixture of cholesterol-(K/D)nDEVDGC)3-trimaleimides form stable nanosponges, whereas cholesterol-(D)nDEVDGC)3-trimaleimide is unable to form supramolecular aggregates with itself. The structure of the novel nanosponges was investigated through explicit solvent and then coarse-grained molecular dynamics (MD) simulations. The nanosponges are between 80 nm and several micrometers in diameters and virtually non-toxic to monocyte/macrophage-like cells.

Experimental triplet and quadruplet fluctuation densities and spatial distribution function integrals for liquid mixtures.

Kirkwood-Buff or Fluctuation Solution Theory can be used to provide experimental pair fluctuations, and/or integrals over the pair distribution functions, from experimental thermodynamic data on liquid mixtures. Here, this type of approach is used to provide triplet and quadruplet fluctuations, and the corresponding integrals over the triplet and quadruplet distribution functions, in a purely thermodynamic manner that avoids the use of structure factors. The approach is then applied to binary mixtures of water + methanol and benzene + methanol over the full composition range under ambient conditions. The observed correlations between the different species vary significantly with composition. The magnitude of the fluctuations and integrals appears to increase as the number of the most polar molecule involved in the fluctuation or integral also increases. A simple physical picture of the fluctuations is provided to help rationalize some of these variations.

Experimental triplet and quadruplet fluctuation densities and spatial distribution function integrals for pure liquids.

Fluctuation solution theory has provided an alternative view of many liquid mixture properties in terms of particle number fluctuations. The particle number fluctuations can also be related to integrals of the corresponding two body distribution functions between molecular pairs in order to provide a more physical picture of solution behavior and molecule affinities. Here, we extend this type of approach to provide expressions for higher order triplet and quadruplet fluctuations, and thereby integrals over the corresponding distribution functions, all of which can be obtained from available experimental thermodynamic data. The fluctuations and integrals are then determined using the International Association for the Properties of Water and Steam Formulation 1995 (IAPWS-95) equation of state for the liquid phase of pure water. The results indicate small, but significant, deviations from a Gaussian distribution for the molecules in this system. The pressure and temperature dependence of the fluctuations and integrals, as well as the limiting behavior as one approaches both the triple point and the critical point, are also examined.

Simulated pressure changes in LacI suggest a link between hydration and functional conformational changes.

The functions of many proteins are associated with interconversions among conformational substates. However, these substates can be difficult to measure experimentally, and determining contributions from hydration changes can be especially difficult. Here, we assessed the use of pressure perturbations to sample the substates accessible to the Escherichia coli lactose repressor protein (LacI) in various liganded forms. In the presence of DNA, the regulatory domain of LacI adopts an Open conformation that, in the absence of DNA, changes to a Closed conformation. Increasing the simulation pressure prevented the transition from an Open to a Closed conformation, in a similar manner to the binding of DNA and anti-inducer, ONPF. The results suggest the hydration of specific residues play a significant role in determining the population of different LacI substates and that simulating pressure perturbation could be useful for assessing the role of hydration changes that accompany functionally-relevant amino acid substitutions.

Health-associated changes of the fecal microbiota in dairy heifer calves during the pre-weaning period.

Introduction: Neonatal calf diarrhea is a multifactorial condition that occurs in early life when calves are particularly susceptible to enteric infection and dysbiosis of the gut microbiome. Good calf health is dependent on successful passive transfer of immunity from the dam through colostrum. There are limited studies on the developing gut microbiota from birth to weaning in calves. Methodology: Therefore, the objective of this study was to examine the effect of immune status and diarrheal incidence on the development of the fecal microbiota in Jersey (n = 22) and Holstein (n = 29) heifer calves throughout the pre-weaning period. Calves were hand-fed a colostrum volume equivalent to 8.5% of their birthweight, from either the calf's dam (n = 28) or re-heated mixed colostrum (≤2 cows, ≤1d; n = 23) within 2 h of birth. All calves were clinically assessed using a modified Wisconsin-Madison calf health scoring system and rectal temperature at day (d) 0, d7, d21, or disease manifestation (DM) and weaning (d83). Weights were recorded at d0, d21, and d83. Calf blood samples were collected at d7 for the determination of calf serum IgG (sIgG). Fecal samples were obtained at d7, d21/DM [mean d22 (SE 0.70)], and at weaning for 16S rRNA amplicon sequencing of the fecal microbiota. Data were processed in R using DADA2; taxonomy was assigned using the SILVA database and further analyzed using Phyloseq and MaAsLin 2. Results and discussion: Significant amplicon sequence variants (ASVs) and calf performance data underwent a Spearman rank-order correlation test. There was no effect (p > 0.05) of colostrum source or calf breed on serum total protein. An effect of calf breed (p < 0.05) was observed on sIgG concentrations such that Holstein calves had 6.49 (SE 2.99) mg/ml higher sIgG than Jersey calves. Colostrum source and calf breed had no effect (p > 0.05) on health status or the alpha diversity of the fecal microbiota. There was a relationship between health status and time interaction (p < 0.001), whereby alpha diversity increased with time; however, diarrheic calves had reduced microbial diversity at DM. No difference (p > 0.05) in beta diversity of the microbiota was detected at d7 or d83. At the genus level, 33 ASVs were associated (adj.p < 0.05) with health status over the pre-weaning period.

The intrinsically disordered transcriptional activation domain of CIITA is functionally tuneable by single substitutions: An exception or a new paradigm?

During protein evolution, some amino acid substitutions modulate protein function ("tuneability"). In most proteins, the tuneable range is wide and can be sampled by a set of protein variants that each contains multiple amino acid substitutions. In other proteins, the full tuneable range can be accessed by a set of variants that each contains a single substitution. Indeed, in some globular proteins, the full tuneable range can be accessed by the set of site-saturating substitutions at an individual "rheostat" position. However, in proteins with intrinsically disordered regions (IDRs), most functional studies-which would also detect tuneability-used multiple substitutions or small deletions. In disordered transcriptional activation domains (ADs), studies with multiple substitutions led to the "acidic exposure" model, which does not anticipate the existence of rheostat positions. In the few studies that did assess effects of single substitutions on AD function, results were mixed: the ADs of two full-length transcription factors did not show tuneability, whereas a fragment of a third AD was tuneable by single substitutions. In this study, we tested tuneability in the AD of full-length human class II transactivator (CIITA). Sequence analyses and experiments showed that CIITA's AD is an IDR. Functional assays of singly-substituted AD variants showed that CIITA's function was highly tuneable, with outcomes not predicted by the acidic exposure model. Four tested positions showed rheostat behavior for transcriptional activation. Thus, tuneability of different IDRs can vary widely. Future studies are needed to illuminate the biophysical features that govern whether an IDR is tuneable by single substitutions.

Breed and ruminal fraction effects on bacterial and archaeal community composition in sheep.

While the breed of cattle can impact on the composition and structure of microbial communities in the rumen, breed-specific effects on rumen microbial communities have rarely been examined in sheep. In addition, rumen microbial composition can differ between ruminal fractions, and be associated with ruminant feed efficiency and methane emissions. In this study, 16S rRNA amplicon sequencing was used to investigate the effects of breed and ruminal fraction on bacterial and archaeal communities in sheep. Solid, liquid and epithelial rumen samples were obtained from a total of 36 lambs, across 4 different sheep breeds (Cheviot (n = 10), Connemara (n = 6), Lanark (n = 10) and Perth (n = 10)), undergoing detailed measurements of feed efficiency, who were offered a nut based cereal diet ad-libitum supplemented with grass silage. Our results demonstrate that the feed conversion ratio (FCR) was lowest for the Cheviot (most efficient), and highest for the Connemara breed (least efficient). In the solid fraction, bacterial community richness was lowest in the Cheviot breed, while Sharpea azabuensis was most abundant in the Perth breed. Lanark, Cheviot and Perth breeds exhibited a significantly higher abundance of epithelial associated Succiniclasticum compared to the Connemara breed. When comparing ruminal fractions, Campylobacter, Family XIII, Mogibacterium, and Lachnospiraceae UCG-008 were most abundant in the epithelial fraction. Our findings indicate that breed can impact the abundance of specific bacterial taxa in sheep while having little effect on the overall composition of the microbial community. This finding has implications for genetic selection breeding programs aimed at improving feed conversion efficiency of sheep. Furthermore, the variations in the distribution of bacterial species identified between ruminal fractions, notably between solid and epithelial fractions, reveals a rumen fraction bias, which has implications for sheep rumen sampling techniques.

Flexible connection of the N-terminal domain in ClpB modulates substrate binding and the aggregate reactivation efficiency.

ClpB reactivates aggregated proteins in cooperation with DnaK/J. The ClpB monomer contains two nucleotide-binding domains (D1, D2), a coiled-coil domain, and an N-terminal domain attached to D1 with a 17-residue-long unstructured linker containing a Gly-Gly motif. The ClpB-mediated protein disaggregation is linked to translocation of substrates through the central channel in the hexameric ClpB, but the events preceding the translocation are poorly understood. The N-terminal domains form a ring surrounding the entrance to the channel and contribute to the aggregate binding. It was suggested that the N-terminal domain's mobility that is maintained by the unstructured linker might control the efficiency of aggregate reactivation. We produced seven variants of ClpB with modified sequence of the N-terminal linker. To increase the linker's conformational flexibility, we inserted up to four Gly next to the GG motif. To decrease the linker's flexibility, we deleted the GG motif and converted it into GP and PP. We found that none of the linker modifications inhibited the basal ClpB ATPase activity or its capability to form oligomers. However, the modified linker ClpB variants showed lower reactivation rates for aggregated glucose-6-phosphate dehydrogenase and firefly luciferase and a lower aggregate-binding efficiency than wt ClpB. We conclude that the linker does not merely connect the N-terminal domain, but it supports the chaperone activity of ClpB by contributing to the efficiency of aggregate binding and disaggregation. Moreover, our results suggest that selective pressure on the linker sequence may be crucial for maintaining the optimal efficiency of aggregate reactivation by ClpB.

Differences in the Composition of the Rumen Microbiota of Finishing Beef Cattle Divergently Ranked for Residual Methane Emissions.

With the advent of high throughput technology, it is now feasible to study the complex relationship of the rumen microbiota with methanogenesis in large populations of ruminant livestock divergently ranked for enteric emissions. Recently, the residual methane emissions (RME) concept has been identified as the optimal phenotype for assessing the methanogenic potential of ruminant livestock due to the trait's independence from animal productivity but strong correlation with daily methane emissions. However, there is currently a dearth of data available on the bacterial and archaeal microbial communities residing in the rumens of animals divergently ranked for RME. Therefore, the objective of this study was to investigate the relationship between the rumen microbiota and RME in a population of finishing beef cattle. Methane emissions were estimated from individual animals using the GreenFeed Emissions Monitoring system for 21 days over a mean feed intake measurement period of 91 days. Residual methane emissions were calculated for 282 crossbred finishing beef cattle, following which a ∼30% difference in all expressions of methane emissions was observed between high and low RME ranked animals. Rumen fluid samples were successfully obtained from 268 animals during the final week of the methane measurement period using a trans-oesophageal sampling device. Rumen microbial DNA was extracted and subjected to 16S rRNA amplicon sequencing. Animals ranked as low RME had the highest relative abundances (P < 0.05) of lactic-acid-producing bacteria (Intestinibaculum, Sharpea, and Olsenella) and Selenomonas, and the lowest (P < 0.05) proportions of Pseudobutyrivibrio, Butyrivibrio, and Mogibacterium. Within the rumen methanogen community, an increased abundance (P < 0.05) of the genus Methanosphaera and Methanobrevibacter RO clade was observed in low RME animals. The relative abundances of both Intestinibaculum and Olsenella were negatively correlated (P < 0.05) with RME and positively correlated with ruminal propionate. A similar relationship was observed for the abundance of Methanosphaera and the Methanobrevibacter RO clade. Findings from this study highlight the ruminal abundance of bacterial genera associated with the synthesis of propionate via the acrylate pathway, as well as the methanogens Methanosphaera and members of the Methanobrevibacter RO clade as potential microbial biomarkers of the methanogenic potential of beef cattle.

Enteric methane research and mitigation strategies for pastoral-based beef cattle production systems.

Ruminant livestock play a key role in global society through the conversion of lignocellulolytic plant matter into high-quality sources of protein for human consumption. However, as a consequence of the digestive physiology of ruminant species, methane (CH4), which originates as a byproduct of enteric fermentation, is accountable for 40% of global agriculture's carbon footprint and ~6% of global greenhouse gas (GHG) emissions. Therefore, meeting the increasing demand for animal protein associated with a growing global population while reducing the GHG intensity of ruminant production will be a challenge for both the livestock industry and the research community. In recent decades, numerous strategies have been identified as having the potential to reduce the methanogenic output of livestock. Dietary supplementation with antimethanogenic compounds, targeting members of the rumen methanogen community and/or suppressing the availability of methanogenesis substrates (mainly H2 and CO2), may have the potential to reduce the methanogenic output of housed livestock. However, reducing the environmental impact of pasture-based beef cattle may be a challenge, but it can be achieved by enhancing the nutritional quality of grazed forage in an effort to improve animal growth rates and ultimately reduce lifetime emissions. In addition, the genetic selection of low-CH4-emitting and/or faster-growing animals will likely benefit all beef cattle production systems by reducing the methanogenic potential of future generations of livestock. Similarly, the development of other mitigation technologies requiring minimal intervention and labor for their application, such as anti-methanogen vaccines, would likely appeal to livestock producers, with high uptake among farmers if proven effective. Therefore, the objective of this review is to give a detailed overview of the CH4 mitigation solutions, both currently available and under development, for temperate pasture-based beef cattle production systems. A description of ruminal methanogenesis and the technologies used to estimate enteric emissions at pastures are also presented.

A Kirkwood-Buff force field for the aromatic amino acids.

In a continuation of our efforts to develop a united atom non-polarizable protein force field based upon the solution theory of Kirkwood and Buff i.e., the Kirkwood-Buff Force Field (KBFF) approach, we present KBFF models for the side chains of phenylalanine, tyrosine, tryptophan, and histidine, including both tautomers of neutral histidine and doubly-protonated histidine. The force fields were specifically designed to reproduce the thermodynamic properties of mixtures over the full composition range in an attempt to provide an improved description of intermolecular interactions. The models were developed by careful parameterization of the solution phase partial charges to reproduce the experimental Kirkwood-Buff integrals for mixtures of solutes representative of the amino acid sidechains in solution. The KBFF parameters and simulated thermodynamic and structural properties are presented for the following eleven binary mixtures: benzene + methanol, benzene + toluene, toluene + methanol, toluene + phenol, toluene + p-cresol, pyrrole + methanol, indole + methanol, pyridine + methanol, pyridine + water, histidine + water, and histidine hydrochloride + water. It is argued that the present approach and models provide a reasonable description of intermolecular interactions which ensures that the required balance between solute-solute, solute-solvent, and solvent-solvent distributions is obtained.

Fluctuation theory of molecular association and conformational equilibria.

General expressions relating the effects of pressure, temperature, and composition on solute association and conformational equilibria using the fluctuation theory of solutions are provided. The expressions are exact and can be used to interpret experimental or computer simulation data for any multicomponent mixture involving molecules of any size and character at any composition. The relationships involve particle-particle, particle-energy, and energy-energy correlations within local regions in the vicinity of each species involved in the equilibrium. In particular, it is demonstrated that the results can be used to study peptide and protein association or aggregation, protein denaturation, and protein-ligand binding. Exactly how the relevant fluctuating properties may be obtained from experimental or computer simulation data are also outlined. It is shown that the enthalpy, heat capacity, and compressibility differences associated with the equilibrium process can, in principle, be obtained from a single simulation. Fluctuation based expressions for partial molar heat capacities, thermal expansions, and isothermal compressibilities are also provided.

Local fluctuations in solution mixtures.

An extension of the traditional Kirkwood-Buff (KB) theory of solutions is outlined which provides additional fluctuating quantities that can be used to characterize and probe the behavior of solution mixtures. Particle-energy and energy-energy fluctuations for local regions of any multicomponent solution are expressed in terms of experimentally obtainable quantities, thereby supplementing the usual particle-particle fluctuations provided by the established KB inversion approach. The expressions are then used to analyze experimental data for pure water over a range of temperatures and pressures, a variety of pure liquids, and three binary solution mixtures - methanol and water, benzene and methanol, and aqueous sodium chloride. In addition to providing information on local properties of solutions it is argued that the particle-energy and energy-energy fluctuations can also be used to test and refine solute and solvent force fields for use in computer simulation studies.

Kirkwood-Buff-Derived Force Field for Peptides and Proteins: Applications of KBFF20.

Here, we perform structural, thermodynamic, and kinetics tests of the Kirkwood-Buff-derived force field, KBFF20, for peptides and proteins developed in the previous article. The physical/structural tests measure the ability of KBFF20 to capture the experimental J-couplings for small peptides, to keep globular monomeric and oligomeric proteins folded, and to produce the experimentally relevant expanded conformational ensembles of intrinsically disordered proteins. The thermodynamic-based tests probe KBFF20's ability to quantify the preferential interactions of sodium chloride around native ß-lactoglobulin and urea around native lysozyme, to reproduce the melting curves for small helix- and sheet-based peptides, and to fold the small proteins Trp-cage and Villin. The kinetics-based tests quantify how well KBFF20 can match the experimental contact formation rates of small, repeat-sequence peptides of variable lengths and the rotational diffusion coefficients of globular proteins. The results suggest that KBFF20 is naturally able to reproduce properties of both folded and disordered proteins, which we attribute to the use of the Kirkwood-Buff theory as the foundation of the force field's development. However, we show that KBFF20 tends to lose some well-defined secondary structural elements and increases the percentage of coil regions, indicating that the perfect balance of all interactions remains elusive. Nevertheless, we argue that KBFF20 is an improvement over recently modified force fields that require ad hoc interventions to prevent the collapse of intrinsically disordered proteins.

Effect of divergence in residual methane emissions on feed intake and efficiency, growth and carcass performance, and indices of rumen fermentation and methane emissions in finishing beef cattle.

Residual expressions of enteric emissions favor a more equitable identification of an animal's methanogenic potential compared with traditional measures of enteric emissions. The objective of this study was to investigate the effect of divergently ranking beef cattle for residual methane emissions (RME) on animal productivity, enteric emissions, and rumen fermentation. Dry matter intake (DMI), growth, feed efficiency, carcass output, and enteric emissions (GreenFeed emissions monitoring system) were recorded on 294 crossbred beef cattle (steers = 135 and heifers = 159; mean age 441 d (SD = 49); initial body weight (BW) of 476 kg (SD = 67)) at the Irish national beef cattle performance test center. Animals were offered a total mixed ration (77% concentrate and 23% forage; 12.6 MJ ME/kg of DM and 12% CP) ad libitum with emissions estimated for 21 d over a mean feed intake measurement period of 91 d. Animals had a mean daily methane emissions (DME) of 229.18 g/d (SD = 45.96), methane yield (MY) of 22.07 g/kg of DMI (SD = 4.06), methane intensity (MI) 0.70 g/kg of carcass weight (SD = 0.15), and RME 0.00 g/d (SD = 0.34). RME was computed as the residuals from a multiple regression model regressing DME on DMI and BW (R2 = 0.45). Animals were ranked into three groups namely high RME (>0.5 SD above the mean), medium RME (±0.5 SD above/below the mean), and low RME (>0.5 SD below the mean). Low RME animals produced 17.6% and 30.4% less (P < 0.05) DME compared with medium and high RME animals, respectively. A ~30% reduction in MY and MI was detected in low versus high RME animals. Positive correlations were apparent among all methane traits with RME most highly associated with (r = 0.86) DME. MY and MI were correlated (P < 0.05) with DMI, growth, feed efficiency, and carcass output. High RME had lower (P < 0.05) ruminal propionate compared with low RME animals and increased (P < 0.05) butyrate compared with medium and low RME animals. Propionate was negatively associated (P < 0.05) with all methane traits. Greater acetate:propionate ratio was associated with higher RME (r = 0.18; P < 0.05). Under the ad libitum feeding regime deployed here, RME was the best predictor of DME and only methane trait independent of animal productivity. Ranking animals on RME presents the opportunity to exploit interanimal variation in enteric emissions as well as providing a more equitable index of the methanogenic potential of an animal on which to investigate the underlying biological regulatory mechanisms.

Kirkwood-Buff-Derived Force Field for Peptides and Proteins: Philosophy and Development of KBFF20.

A new classical nonpolarizable force field, KBFF20, for the simulation of peptides and proteins is presented. The force field relies heavily on the use of Kirkwood-Buff theory to provide a comparison of simulated and experimental Kirkwood-Buff integrals for solutes containing the functional groups common in proteins, thus ensuring intermolecular interactions that provide a good balance between the peptide-peptide, peptide-solvent, and solvent-solvent distributions observed in solution mixtures. In this way, it differs significantly from other biomolecular force fields. Further development and testing of the intermolecular potentials are presented here. Subsequently, rotational potentials for the ϕ/ψ and χ dihedral degrees of freedom are obtained by analysis of the Protein Data Bank, followed by small modifications to provide a reasonable balance between simulated and observed α and ß percentages for small peptides. This, the first of two articles, describes in detail the philosophy and development behind KBFF20.