Variational Methods for Biomolecular Modeling

Structure, function and dynamics of many biomolecular systems can be characterized by the energetic variational principle and the corresponding systems of partial differential equations (PDEs). This principle allows us to focus on the identification of essential energetic components, the optimal parametrization of energies, and the efficient computational implementation of energy variation or minimization. Given the fact that complex biomolecular systems are structurally non-uniform and their interactions occur through contact interfaces, their free energies are associated with various interfaces as well, such as solute-solvent interface, molecular binding interface, lipid domain interface, and membrane surfaces. This fact motivates the inclusion of interface geometry, particular its curvatures, to the parametrization of free energies. Applications of such interface geometry based energetic variational principles are illustrated through three concrete topics: the multiscale modeling of biomolecular electrostatics and solvation that includes the curvature energy of the molecular surface, the formation of microdomains on lipid membrane due to the geometric and molecular mechanics at the lipid interface, and the mean curvature driven protein localization on membrane surfaces. By further implicitly representing the interface using a phase field function over the entire domain, one can simulate the dynamics of the interface and the corresponding energy variation by evolving the phase field function, achieving significant reduction of the number of degrees of freedom and computational complexity. Strategies for improving the efficiency of computational implementations and for extending applications to coarse-graining or multiscale molecular simulations are outlined.Comment: 36 page

Structural and biochemical analysis of human pathogenic astrovirus serine protease at 2.0 A resolution

Astroviruses are single-stranded RNA viruses with a replication strategy based on the proteolytic processing of a polyprotein precursor and subsequent release of the viral enzymes of replication. So far, the catalytic properties of the astrovirus protease as well as its structure have remained uncharacterized. In this study, the three-dimensional crystal structure of the predicted protease of human pathogenic astrovirus has been solved to 2.0 A resolution. The protein displays the typical properties of trypsin-like enzymes but also several characteristic features: (i) a catalytic Asp-His-Ser triad in which the aspartate side chain is oriented away from the histidine, being replaced by a water molecule; (ii) a non-common conformation and composition of the S1 pocket; and (iii) the lack of the typical surface beta-ribbons together with a "featureless" shape of the substrate-binding site. Hydrolytic activity assays indicate that the S1 pocket recognises Glu and Asp side chains specifically, which, therefore, are predicted to occupy the P1 position on the substrate cleavage site. The positive electrostatic potential featured by the S1 region underlies this specificity. The comparative structural analysis highlights the peculiarity of the astrovirus protease, and differentiates it from the human and viral serine proteases