Umbrella Sampling-Based Method to Compute Ligand-Binding Affinity

An atomistic scale simulation study of structural properties in the silk-fibrohexamerin complex

The use of Bombyx mori silk fibroin in composite materials has been extensively explored in many studies, owing to its remarkable mechanical properties. Recently, the N-glycan-engineered P25 protein was utilized to improve the mechanical properties of silk. However, the mechanism by which N-glycan-engineered P25 protein enhances the mechanical properties of silk remains unclear. This study analyzed the interaction between the P25 protein and silkworm silk using quantum mechanics/molecular mechanics multiscale simulations and discovered stronger hydrogen bonding between the amorphous domain and the P25 protein. The results confirmed that glycoengineering of the mannose molecule in N-glycan in orders of three, five, and seven increased the hydrogen bonding of the amorphous structures. However, P25 has fewer binding interactions with the crystalline domain. Silk amino acids and mannose molecules were analyzed using QM simulations, and hydroxyl and charged amino acids in the amorphous domains were found to have relatively higher reactivity with mannose molecules in N-glycans than basic and aliphatic amino acids in the crystalline domain. This study demonstrates how the N-glycan-engineered P25 protein can improve the mechanical properties of silk fibroin and identifies a key factor for N-glycan-engineered proteins.

Simulation of the ligand-leaving process of the human heat shock protein

The human heat shock protein plays a critical role in various diseases and is an important target for pharmacological modulation. Simulation of conformational changes and free energy profiles of the human heat shock protein derived by the ligand-leaving process is a challenging issue. In this work, steered molecular dynamics simulation was adopted to simulate the ligand-leaving process. Two composite systems of heat shock protein NHSP90 and small molecules 6FJ and 6G7 are selected as research objects. The free energy during the leaving of ligand small molecules is calculated using conventional molecular dynamics simulation, steered molecular dynamics simulation (SMD), and the umbrella sampling method. We found that the a slower pulling velocity (0.001 nm ns-1) will result in 2.19 kcal mol-1, and the umbrella sampling method gives a value of 3.26 kcal mol-1 for the free energy difference for the two systems, which reasonably agrees with experimental results. A faster-pulling velocity (0.01 nm ns-1) leads to a large overestimation of free energy. At the same time, the conformational analysis indicated that the faster pulling velocity may lead to the conformational change of NHSP90, which was proved to be false by the slower pulling velocity and the umbrella sampling method.

Molecular Insights into the Binding Behavior of Imidazolium Ionic Liquids to the Receptor Binding Domain of the SARS-CoV-2 Spike Protein

The receptor-binding domain (RBD) of the SARS-CoV-2 spike protein is considered as a key target for the design and development of COVID-19 drugs and inhibitors. Due to their unique structure and properties, ionic liquids (ILs) have many special interactions with proteins, showing great potential in biomedicine. Nevertheless, few research studies have been carried out on ILs and the spike RBD protein. Here, we explore the interaction of ILs and the RBD protein through large-scale molecular dynamics simulations (4 μs in total). It was found that IL cations with long alkyl chain lengths (nchain) could spontaneously bind to the cavity region of the RBD protein. The longer the alkyl chain is, the stabler the cations bind to the protein. The binding free energy (ΔG) had the same trend, peaking at nchain = 12 with -101.19 kJ/mol. The cationic chain lengths and their fit to the pocket are decisive factors that influence the binding strength of cations and proteins. The cationic imidazole ring has a high contact frequency with phenylalanine and tryptophan, and the hydrophobic residues phenylalanine, valine, leucine, and isoleucine are the most interacting residues with side chains of cations. Meanwhile, through analysis of the interaction energy, the hydrophobic and π-π interactions are the main contributors to the high affinity between cations and the RBD protein. In addition, the long-chain ILs would also act on the protein through clustering. These studies not only provide insights into the molecular interaction between ILs and the RBD of SARS-CoV-2 but also contribute to the rational design of IL-based drugs, drug carriers, and selective inhibitors as a therapeutic for SARS-CoV-2.

Probing interactions of aminoacyl-adenylate with Mycobacterium tuberculosis methionyl-tRNA synthetase through in silico site-directed mutagenesis and free energy calculation

Methionyl-tRNA synthetase (MetRS) is an attractive molecular target for antibiotic discovery. Recently, we have developed several classes of small-molecular inhibitors of Mycobacterium tuberculosis MetRS possessing antibacterial activity. In this article, we performed in silico site-directed mutagenesis of aminoacyl-adenylate binding site of M. tuberculosis MetRS in order to identify crucial amino acid residues for substrate interaction. The umbrella sampling algorithm was used to calculate the binding free energy (ΔG) of these mutated forms with methionyl-adenylate analogue. According to the obtained results, the replacement of Glu24 and Leu293 by alanine leads to the most significant decrease in the binding free energy (ΔG) for adenylate analogue with methionyl-tRNA synthetase indicating increasing of the affinity, which in turn causes the loss of compounds inhibitory activity. Therefore, these amino acid residues can be proposed for further experimental site-directed mutagenesis to confirm binding mode of inhibitors and should be taken into account during chemical optimization to overcome resistance due to mutations.Communicated by Ramaswamy H. Sarma.