The structure of water molecules near a hydrophobic solute remains elusive despite a long history of scrutiny. Here, we re-examine the subtle issue by a combination of thermodynamic analysis for Henry's constants of several nonpolar gases over a broad range of temperatures and molecular dynamic simulations for the water structure in the hydration shell using several popular semiempirical models of liquid water. Both the structural and thermodynamic data indicate that hydrophobic hydration reduces the degree of the hydrogen bonding among water molecules, and the effect becomes more prominent at high temperatures. Hydrogen-bond formation is slightly hindered near a hydrophobic solute due to the restriction of the degree of freedom for water molecules in the solvation shell, and the confinement effect becomes more significant as temperature increases. Reduction in the extent of hydrogen bonding is fully consistent with a positive contribution of a small hydrophobic solute to the solution heat capacity. As predicted by the scaled-particle theory, both Henry's constants and simulation results suggest that the hydration entropy is determined primarily by cavity formation in liquid water, with its magnitude rising with the solute size but declining with temperature.

Conventional wisdom for controlling the nanoparticle size and shape during synthesis is that particle growth favors the direction of a facet with the highest surface energy. However, the particle solvation free energy, which dictates the particle stability and growth, depends not only on the surface area and surface free energy but also on other geometric measures such as the solvent excluded volume and the surface curvature and their affiliated thermodynamic properties. In this work, we study the geometrical effects on the solvation free energies of nonspherical nanoparticles using morphometric thermodynamics and density functional theories. For idealized systems that account for only molecular excluded-volume interactions, morphometric thermodynamics yields a reliable solvation free energy when the particle size is significantly larger than the solvent correlation length. However, noticeable deviations can be identified in comparison to the microscopic theories for predicting the solvation free energies of small nanoparticles. This conclusion also holds for predicting the potential of mean force underlying the colloidal "key-and-lock" interactions. Complementary to the microscopic theories, morphometric thermodynamics requires negligible computational cost, therefore making it very appealing for a broad range of practical applications.

Hepatitis B virus (HBV) controls genome encapsidation and reverse transcription from a single-stranded RNA to a double-stranded DNA through the flexible C-terminal domain (CTD) of the capsid proteins. Although the microscopic structure of the nucleocapsid plays a critical role in the life cycle of HBV, the location of CTD residues at different stages of viral replication remains poorly understood. In this work, we report the radial distributions of individual amino-acid residues of the CTD tails for both empty and RNA-containing HBV capsids by using a coarse-grained model for the key biological components and the classical density functional theory. The density functional theory calculations reveal substantial exposure of the CTD residues outside the capsid, in particular when it is devoid of any nucleic materials. The outermost layer of the capsid surface mainly consists of residues from (170)Ser-(175)Arg of the CTD tails, i.e., the serine-arginine protein kinase binding motif. The theoretical description corroborates recent in vitro studies that show a transient CTD distribution captured by serine-arginine protein kinase binding. We have also investigated the nucleocapsid structural changes due to phosphorylation of serine residues and shown a correlation between the CTD location and the internal distribution of RNA segments.

Understanding the fundamentals of genome packaging in viral capsids is important for finding effective antiviral strategies and for utilizing benign viral particles for gene therapy. While the structure of encapsidated genomic materials has been routinely characterized with experimental techniques such as cryo-electron microscopy and x-ray diffraction, much less is known about the molecular driving forces underlying genome assembly in an intracellular environment and its in vivo interactions with the capsid proteins. Here we study the thermodynamic basis of the pregenomic RNA encapsidation in human Hepatitis B virus in vivo using a coarse-grained molecular model that captures the essential components of nonspecific intermolecular interactions. The thermodynamic model is used to examine how the electrostatic interaction between the packaged RNA and the highly charged C-terminal domains (CTD) of capsid proteins regulate the nucleocapsid formation. The theoretical model predicts optimal RNA content in Hepatitis B virus nucleocapsids with different CTD lengths in good agreement with mutagenesis measurements, confirming the predominant role of electrostatic interactions and molecular excluded-volume effects in genome packaging. We find that the amount of encapsidated RNA is not linearly correlated with the net charge of CTD tails as suggested by earlier theoretical studies. Our thermodynamic analysis of the nucleocapsid structure and stability indicates that ∼10% of the CTD residues are free from complexation with RNA, resulting in partially exposed CTD tails. The thermodynamic model also predicts the free energy of complex formation between macromolecules, which corroborates experimental results for the impact of CTD truncation on the nucleocapsid stability.

No abstract