It is a central issue to elucidate the new type of molecular recognition accompanied by a global structural change of a molecule upon binding to its targets. Here we investigate the driving force for the binding of R12 (a ribonucleic acid aptamer) and P16 (a partial peptide of a prion protein) during which P16 exhibits the global structural change. We calculate changes in thermodynamic quantities upon the R12–P16 binding using a statistical-mechanical approach combined with molecular models for water which is currently best suited to studies on hydration of biomolecules. The binding is driven by a water-entropy gain originating primarily from an increase in the total volume available to the translational displacement of water molecules in the system. The energy decrease due to the gain of R12–P16 attractive (van der Waals and electrostatic) interactions is almost canceled out by the energy increase related to the loss of R12–water and P16–water attractive interactions. We can explain the general experimental result that stacking of flat moieties, hydrogen bonding and molecular-shape and electrostatic complementarities are frequently observed in the complexes. It is argued that the water-entropy gain is largely influenced by the geometric characteristics (overall shapes, sizes and detailed polyatomic structures) of the biomolecules.

Although the two membrane proteins, thermophilic rhodopsin (TR) and xanthorhodopsin (XR), share a high similarity in amino-acid sequence and an almost indistinguishable three-dimensional structure, TR is much more thermostable than XR. This is counterintuitive also because TR possesses only a smaller number of intramolecular hydrogen bonds (HBs) than XR. Here we investigate physical origins of the remarkable difference between XR and TR in the stability. Our free-energy function (FEF) is improved so that not only the portion within the transmembrane (TM) region but also the extracellular and intracellular portions within the water-immersed (WI) regions can be considered in assessing the stability. The assessment is performed on the basis of the FEF change upon protein folding, which is calculated for the crystal structure of XR or TR. Since the energetics within the TM region is substantially different from that within the WI regions, we determine the TM and WI portions of XR or TR by analyzing the distribution of water molecules using all-atom molecular dynamics simulations. The energetic component of the FEF change consists of a decrease in energy arising from the formation of intramolecular HBs and an increase in energy caused by the break of protein-water HBs referred to as “energetic dehydration penalty.” The entropic component is a gain of the translational, configurational entropies of hydrocarbon groups within the lipid bilayer and of water molecules. The entropic component is calculated using the integral equation theory combined with our morphometric approach. The energetic one is estimated by a simple but physically reasonable method. We show that TR is much more stable than XR for the following reasons: The decrease in energy within the TM region is larger, and the energetic dehydration penalty within the WI regions is smaller, leading to higher energetic stabilization, and tighter packing of side chains accompanying the association of seven helices confers higher entropic stabilization on TR.

Taking protein G with 56 residues for a case study, we investigate the mechanism of protein folding. In addition to its native structure possessing α-helix and β-sheet contents of 27% and 39%, respectively, we construct a number of misfolded decoys with a wide variety of α-helix and β-sheet contents. We then consider a hierarchy of 8 different models with increasing level of detail in terms of the number of entropic and energetic physical factors incorporated. The polyatomic structure is always taken into account, but the side chains are removed in half of the models. The solvent is formed by either neutral hard spheres or water molecules. Protein intramolecular hydrogen bonds (H-bonds) and protein-solvent H-bonds (the latter is present only in water) are accounted for or not, depending on the model considered. We then apply a physics-based free-energy function (FEF) corresponding to each model and investigate which structures are most stabilized. This special approach taken on a step-by-step basis enables us to clarify the role of each physical factor in contributing to the structural stability and separately elucidate its effect. Depending on the model employed, significantly different structures such as very compact configurations with no secondary structures and configurations of associated α-helices are optimally stabilized. The native structure can be identified as that with lowest FEF only when the most detailed model is employed. This result is significant for at least the two reasons: The most detailed model considered here is able to capture the fundamental aspects of protein folding notwithstanding its simplicity; and it is shown that the native structure is stabilized by a complex interplay of minimal multiple factors that must be all included in the description. In the absence of even a single of these factors, the protein is likely to be driven towards a different, more stable state.

Molecularly imprinted polymers (MIPs) are functional materials capable of molecular recognition and other biorelevant functions and can be prepared with ease in a tailor-made fashion by copolymerization in the presence of a crosslinker and a template molecule (the target molecule or its analog) conjugated with a functional monomer(s). These robust and affordable synthetic polymer receptors are highly accepted as favorable alternatives to biomacromolecules such as antibodies and enzymes. This article provides a critical review of the present state of MIPs to establish perspectives on this technique via a survey of early to contemporary work, mainly conducted by the authors, covering new principles and methodology to generate MIPs. These include the design and synthesis of new functional monomers and crosslinkers to develop and/or introduce new functionalities, new polymerization methods to improve the imprint effect, and highly sensitive MIP-based assays and sensors.

We develop a new methodology best suited to the identification of thermostabilizing mutations for an intrinsically stable membrane protein. The recently discovered thermophilic rhodopsin, whose apparent midpoint temperature of thermal denaturation T m is measured to be ∼91.8 °C, is chosen as a paradigmatic target. In the methodology, we first regard the residues whose side chains are missing in the crystal structure of the wild type (WT) as the “residues with disordered side chains,” which make no significant contributions to the stability, unlike the other essential residues. We then undertake mutating each of the residues with disordered side chains to another residue except Ala and Pro, and the resultant mutant structure is constructed by modifying only the local structure around the mutated residue. This construction is based on the postulation that the structure formed by the other essential residues, which is nearly optimized in such a highly stable protein, should not be modified. The stability changes arising from the mutations are then evaluated using our physics-based free-energy function (FEF). We choose the mutations for which the FEF is much lower than for the WT and test them by experiments. We successfully find three mutants that are significantly more stable than the WT. A double mutant whose T m reaches ∼100 °C is also discovered.