Understanding the molecular mechanism of enzyme dynamics of ribonuclease A through protonation/deprotonation of HIS48

Insight into binding mechanisms of inhibitors MKP56, MKP73, MKP86, and MKP97 to HIV-1 protease by using molecular dynamics simulation

HIV-1 protease (PR) has been a significant target for design of potent inhibitors curing acquired immunodeficiency syndrome. Molecular dynamics simulations coupled with molecular mechanics Poisson-Boltzmann surface area method were performed to study interaction modes of four inhibitors MKP56, MKP73, MKP86, and MKP97 with PR. The results suggest that the main force controlling interactions of inhibitors with PR should be contributed by van der Waals interactions between inhibitors and PR. The cross-correlation analyses based on MD trajectories show that inhibitor binding produces significant effect on the flap dynamics of PR. Hydrogen bond analyses indicate that inhibitors can form stable hydrogen bonding interactions with the residues from the catalytic strands of PR. The contributions of separate residues to inhibitor bindings are evaluated by using residue-based free energy decomposition method and the results demonstrate that the CH-π and CH-CH interactions between the hydrophobic groups of inhibitors with residues drive the associations of inhibitors with PR. We expect that this study can provide a significant theoretical aid for design of potent inhibitors targeting PR.

Large-scale molecular dynamics simulation: Effect of polarization on thrombin-ligand binding energy

Molecular dynamics (MD) simulations lasting 500 ns were performed in explicit water to investigate the effect of polarization on the binding of ligands to human α-thrombin based on the standard nonpolarizable AMBER force field and the quantum-derived polarized protein-specific charge (PPC). The PPC includes the electronic polarization effect of the thrombin-ligand complex, which is absent in the standard force field. A detailed analysis and comparison of the results of the MD simulation with experimental data provided strong evidence that intra-protein, protein-ligand hydrogen bonds and the root-mean-square deviation of backbone atoms were significantly stabilized through electronic polarization. Specifically, two critical hydrogen bonds between thrombin and the ligand were broken at approximately 190 ns when AMBER force field was used and the number of intra-protein backbone hydrogen bonds was higher under PPC than under AMBER. The thrombin-ligand binding energy was computed using the molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) method, and the results were consistent with the experimental value obtained using PPC. Because hydrogen bonds were unstable, it was failed to predict the binding affinity under the AMBER force field. Furthermore, the results of the present study revealed that differences in the binding free energy between AMBER and PPC almost comes from the electrostatic interaction. Thus, this study provides evidence that protein polarization is critical to accurately describe protein-ligand binding.

Drug resistance mechanisms of three mutations V32I, I47V and V82I in HIV-1 protease toward inhibitors probed by molecular dynamics simulations and binding free energy predictions

Molecular dynamics simulation and binding free energy calculations were used to probe drug resistance of HIV-1 protease mutations toward inhibitors.

The protonation site of para-dimethylaminobenzoic acid using atmospheric pressure ionization methods

The protonation site of para-dimethylaminobenzoic acid (p-DMABA) was investigated using atmospheric pressure ionization methods (ESI and APCI) coupled with collision-induced dissociation (CID), nuclear magnetic resonance (NMR), and computational chemistry. Theoretical calculations and NMR experiments indicate that the dimethyl amino group is the preferred site of protonation both in the gas phase and aqueous solution. Protonation of p-DMABA occurs at the nitrogen atom by ESI independent of the solvents and other operation conditions under typical thermodynamic control. However, APCI produces a mixture of the nitrogen- and carbonyl oxygen-protonated p-DMABA when aprotic organic solvents (acetonitrile, acetone, and tetrahydrofuran) are used, exhibiting evident kinetic characteristics of protonation. But using protic organic solvents (methanol, ethanol, and isopropanol) in APCI still leads to the formation of thermodynamically stable N-protonated p-DMABA. These structural assignments were based on the different CID behavior of the N- and O-protonated p-DMABA. The losses of methyl radical and water are the diagnostic fragmentations of the N- and O-protonated p-DMABA, respectively. In addition, the N-protonated p-DMABA is more stable than the O-protonated p-DMABA in CID revealed by energy resolved experiments and theoretical calculations.

Multiscale methods for computational RNA enzymology

RNA catalysis is of fundamental importance to biology and yet remains ill-understood due to its complex nature. The multidimensional "problem space" of RNA catalysis includes both local and global conformational rearrangements, changes in the ion atmosphere around nucleic acids and metal ion binding, dependence on potentially correlated protonation states of key residues, and bond breaking/forming in the chemical steps of the reaction. The goal of this chapter is to summarize and apply multiscale modeling methods in an effort to target the different parts of the RNA catalysis problem space while also addressing the limitations and pitfalls of these methods. Classical molecular dynamics simulations, reference interaction site model calculations, constant pH molecular dynamics (CpHMD) simulations, Hamiltonian replica exchange molecular dynamics, and quantum mechanical/molecular mechanical simulations will be discussed in the context of the study of RNA backbone cleavage transesterification. This reaction is catalyzed by both RNA and protein enzymes, and here we examine the different mechanistic strategies taken by the hepatitis delta virus ribozyme and RNase A.

Binding modes of three inhibitors 8CA, F8A and I4A to A-FABP studied based on molecular dynamics simulation

Adipocyte fatty-acid binding protein (A-FABP) is an important target of drug designs treating some diseases related to lipid-mediated biology. Molecular dynamics (MD) simulations coupled with solvated interaction energy method (SIE) were carried out to study the binding modes of three inhibitors 8CA, F8A and I4A to A-FABP. The rank of our predicted binding affinities is in accordance with experimental data. The results show that the substitution in the position 5 of N-benzyl and the seven-membered ring of N-benzyl-indole carboxylic acids strengthen the I4A binding, while the substitution in the position 2 of N-benzyl weakens the F8A binding. Computational alanine scanning and dynamics analyses were performed and the results suggest that the polar interactions of the positively charged residue R126 with the three inhibitors provide a significant contribution to inhibitor bindings. This polar interaction induces the disappearance of the correlated motion of the C terminus of A-FABP relative to the N terminus and favors the stability of the binding complex. This study is helpful for the rational design of potent inhibitors within the fields of metabolic disease, inflammation and atherosclerosis.

Protein's native structure is dynamically stabilized by electronic polarization

In this paper, molecular dynamics (MD) simulations were performed for a number of benchmark proteins using both the standard assisted model building with energy refinement (AMBER) charge and the dynamically adjusted polarized protein-specific charge (DPPC) from quantum fragment calculations to provide accurate electrostatic interactions. Our result shows that proteins' dynamic structures drifted away from the native structures in simulations under standard (nonpolarizable) AMBER force field. For comparison, proteins' native structures were dynamically stable after a long time simulation under DPPC. The free energy landscape reveals that the native structure is the lowest energy conformation under DPPC, while it is not under standard AMBER charge. To further investigate the polarization effect on the stability of native structures of proteins, we restarted from some decoy structures generated from simulations using standard AMBER charges and then carried out further MD simulation using DPPC to refine those structures. Our study shows that the native structures from these decoy structures can be mostly recovered using DPPC and that the dynamic structures with the highest population in cluster analysis are in close agreement with the corresponding native structures. The current study demonstrates the importance of electronic polarization of protein in stabilizing the native structure.

Insight into conformational change for 14-3-3σ protein by molecular dynamics simulation

14-3-3σ is a member of a highly conserved family of 14-3-3 proteins that has a double-edged sword role in human cancers. Former reports have indicated that the 14-3-3 protein may be in an open or closed state. In this work, we found that the apo-14-3-3σ is in an open state compared with the phosphopeptide bound 14-3-3σ complex which is in a more closed state based on our 80 ns molecular dynamics (MD) simulations. The interaction between the two monomers of 14-3-3σ in the open state is the same as that in the closed state. In both open and closed states, helices A to D, which are involved in dimerization, are stable. However, large differences are found in helices E and F. The hydrophobic contacts and hydrogen bonds between helices E and G in apo-14-3-3σ are different from those in the bound 14-3-3σ complex. The restrained and the mutated (Arg56 or Arg129 to alanine) MD simulations indicate that the conformation of four residues (Lys49, Arg56, Arg129 and Tyr130) may play an important role to keep the 14-3-3σ protein in an open or closed state. These results would be useful to evaluate the 14-3-3σ protein structure-function relationship.

Considering protonation as a posttranslational modification regulating protein structure and function

Posttranslational modification is an evolutionarily conserved mechanism for regulating protein activity, binding affinity, and stability. Compared with established posttranslational modifications such as phosphorylation or ubiquitination, posttranslational modification by protons within physiological pH ranges is a less recognized mechanism for regulating protein function. By changing the charge of amino acid side chains, posttranslational modification by protons can drive dynamic changes in protein conformation and function. Addition and removal of a proton is rapid and reversible and, in contrast to most other posttranslational modifications, does not require an enzyme. Signaling specificity is achieved by only a minority of sites in proteins titrating within the physiological pH range. Here, we examine the structural mechanisms and functional consequences of proton posttranslational modification of pH-sensing proteins regulating different cellular processes.