Free energy and kinetic rate calculation via non-equilibrium molecular simulation: application to biomolecules

Calculating Binding Free Energies in Model Host-Guest Systems with Unrestrained Advanced Sampling

Host-guest interactions are important to the design of pharmaceuticals and, more broadly, to soft materials as they can enable targeted, strong, and specific interactions between molecules. The binding process between the host and guest may be classified as a "rare event" when viewing the system at atomic scales, such as those explored in molecular dynamics simulations. To obtain equilibrium binding conformations and dissociation constants from these simulations, it is essential to resolve these rare events. Advanced sampling methods such as the adaptive biasing force (ABF) promote the occurrence of less probable configurations in a system, therefore facilitating the sampling of essential collective variables that characterize the host-guest interactions. Here, we present the application of ABF to a rod-cavitand coarse-grained model of host-guest systems to acquire the potential of mean force. We show that the employment of ABF enables the computation of the configurational and thermodynamic properties of bound and unbound states, including the free energy landscape. Moreover, we identify important dynamic bottlenecks that limit sampling and discuss how these may be addressed in more general systems.

Treatment of flexibility of protein backbone in simulations of protein-ligand interactions using steered molecular dynamics

To ensure that an external force can break the interaction between a protein and a ligand, the steered molecular dynamics simulation requires a harmonic restrained potential applied to the protein backbone. A usual practice is that all or a certain number of protein's heavy atoms or Cα atoms are fixed, being restrained by a small force. This present study reveals that while fixing both either all heavy atoms and or all Cα atoms is not a good approach, while fixing a too small number of few atoms sometimes cannot prevent the protein from rotating under the influence of the bulk water layer, and the pulled molecule may smack into the wall of the active site. We found that restraining the Cα atoms under certain conditions is more relevant. Thus, we would propose an alternative solution in which only the Cα atoms of the protein at a distance larger than 1.2 nm from the ligand are restrained. A more flexible, but not too flexible, protein will be expected to lead to a more natural release of the ligand.

Loop Evolutionary Patterns Shape Catalytic Efficiency of TRI101/201 for Trichothecenes: Insights into Protein-Substrate Interactions

Trichothecenes are highly toxic mycotoxins produced by Fusarium fungi, while TRI101/201 family enzymes play a crucial role in detoxification through acetylation. Studies on the substrate specificity and catalytic kinetics of TRI101/201 have revealed distinct kinetic characteristics, with significant differences observed in catalytic efficiency toward deoxynivalenol, while the catalytic efficiency for T-2 toxin remains relatively consistent. In this study, we used structural bioinformatics analysis and a molecular dynamics simulation workflow to investigate the mechanism underlying the differential catalytic activity of TRI101/201. The findings revealed that the binding stability between trichothecenes and TRI101/201 hinges primarily on a hydrophobic cage structure within the binding site. An intrinsic disordered loop, termed loop cover, defined the evolutionary patterns of the TRI101/201 protein family that are categorized into four subfamilies (V1/V2/V3/M). Furthermore, the unique loop displayed different conformations among these subfamilies' structures, which served to disrupt (V1/V2/V3) or reinforce (M) the hydrophobic cages. The disrupted cages enhanced the water exposure of the hydrophilic moieties of substrates like deoxynivalenol and thereby hindered their binding to the catalytic sites of V-type enzymes. In contrast, this water exposure does not affect substrates like T-2 toxin, which have more hydrophobic substituents, resulting in a comparable catalytic efficiency of both V- and M-type enzymes. Overall, our studies provide theoretical support for understanding the catalytic mechanism of TRI101/201, which shows how an intrinsic disordered loop could impact the protein-ligand binding and suggests a direction for rational protein design in the future.

Editors' Roundup: June 2023

This commentary article represents the latest edition of the Biophysical Reviews 'Editors' Roundup' Series - a platform made available to the editorial board members of any journal with a genuine interest in promoting biophysical content. Each journal associated editor is able to submit a short description of up to five articles recently appearing in their journals with an explanation of why these articles are of interest. This edition (Vol. 15 Issue 3 June 2023) carries contributions from editorial members associated with Biophysics and Physicobiology (Biophysical Society of Japan), Biophysics (Russian Academy of Sciences), Cell Biochemistry and Biophysics (Springer), and Biophysical Reviews (IUPAB-International Union for Pure and Applied Biophysics).

Computational biophysics and structural biology of proteins-a Special Issue in honor of Prof. Haruki Nakamura's 70th birthday

Receiving his initial training jointly in theoretical and applied physics at the University of Tokyo, Professor Haruki Nakamura has had a long and eventful scientific career, along the way helping to shape the way that biophysics is carried out in Japan. Concentrating his research efforts on the simulation of protein structure and function, he has, over his career arc, acted as director of the Institute for Protein Research (Osaka, Japan), director of the Protein Data Bank of Japan (PDBj), president of the Biophysical Society of Japan (BSJ), president of the Protein Science Society of Japan (PSSJ), and group leader and professor of Bioinformatics and Computational Structural Biology at Osaka University. In 2022, Prof. Haruki Nakamura turned 70 years old, and to mark this occasion, his scientific colleagues from around the world have combined their efforts to produce this Festschrift Issue of the IUPAB Biophysical Reviews journal around the theme of the computational biophysics and structural biology of proteins.