Diffusional interaction behavior of NSAIDs in lipid bilayer membrane using molecular dynamics (MD) simulation: Aspirin and Ibuprofen

In this research, for the first time, molecular dynamics (MD) method was used to simulate aspirin and ibuprofen at various concentrations and in neutral and charged states. Effects of the concentration (dosage), charge state, and existence of an integral protein in the membrane on the diffusion rate of drug molecules into lipid bilayer membrane were investigated on 11 systems, for which the parameters indicating diffusion rate and those affecting the rate were evaluated. Considering the diffusion rate, a suitable score was assigned to each system, based on which, analysis of variance (ANOVA) was performed. By calculating the effect size of the indicative parameters and total scores, an optimum system with the highest diffusion rate was determined. Consequently, diffusion rate controlling parameters were obtained: the drug–water hydrogen bond in protein-free systems and protein–drug hydrogen bond in the systems containing protein.     

Deprotonation states of the two active site water molecules regulate the binding of protein phosphatase 5 with its substrate: A molecular dynamics study

Protein phosphatase 5 (PP5), mainly localized in human brain, can dephosphorylate tau protein whose high level of phosphorylation is related to Alzheimer's disease. Similar to other protein phosphatases, PP5 has a conserved motif in the catalytic domain that contains two binding sites for manganese (Mn²⁺) ions. Structural data indicate that two active site water molecules, one bridging the two Mn²⁺ ions and the other terminally coordinated with one of the Mn²⁺ ions (Mn1), are involved in catalysis. Recently, a density functional theory study revealed that the two water molecules can be both deprotonated to keep a neutral active site for catalysis. The theoretical study gives us an insight into the catalytic mechanism of PP5, but the knowledge of how the deprotonation states of the two water molecules affect the binding of PP5 with its substrate is still lacking. To approach this problem, molecular dynamics simulations were performed to model the four possible deprotonation states. Through structural, dynamical and energetic analyses, the results demonstrate that the deprotonation states of the two water molecules affect the structure of the active site including the distance between the two Mn²⁺ ions and their coordination, impact the interaction energy of residues R275, R400 and H304 which directly interact with the substrate phosphoserine, and mediate the dynamics of helix αJ which is involved in regulation of the enzyme's activity. Furthermore, the deprotonation state that is preferable for PP5 binding of its substrate has been identified. These findings could provide new design strategy for PP5 inhibitor.     

Investigation on the binding mechanism of loratinib with the c-ros oncogene 1 (ROS1) receptor tyrosine kinase via molecular dynamics simulation and binding free energy calculations

The c-ros oncogene 1 (ROS1) has proven to be an important cancer target for the treatment of various human cancers. The anaplastic lymphoma kinase inhibitor crizotinib has been granted approval for the treatment of patients with ROS1 positive metastatic non-small-cell lung cancer by the Food and Drug Administration on 2016. However, serious resistance due to the secondary mutation of glycine 2032 to arginine (G2032R) was developed in clinical studies. Loratinib (PF-06463922), a macrocyclic analog of crizotinib, showed significantly improved inhibitory activity against wild–type (WT) ROS1 and ROS1^G2032R mutant. To provide insights into the inhibition mechanism, molecular dynamics simulations and free energy calculations were carried out for the complexes of loratinib with WT and G2032R mutated ROS1. The apo-ROS1^WT and apo-ROS1^G2032R systems showed similar RMSF distributions, while ROS1^G2032R-loratinib showed significantly higher than that of WT ROS1-loratinib, which revealed that the binding of loratinib to ROS1^G2032R significantly interfered the ?uctuation of protein. Calculations of binding free energies indicate that G2032R mutation significantly reduces the binding affinity of loratinib for ROS1, which arose mostly from the increase of conformation entropy and the decrease of solvation energy. Furthermore, detailed per-residue binding free energies highlighted the increased and decreased contributions of some residues in the G2032R mutated systems. The present study revealed the detailed inhibitory mechanism of loratinib as potent WT and G2032R mutated ROS1 inhibitor, which was expected to provide a basis for rational drug design.     

Immunophysical exploration of C3d-CR2(CCP1-2) interaction using molecular dynamics and electrostatics

The formation of the complex between the d-fragment of the complement component C3 (C3d) and the modular complement receptor-2 (CR2) is important for cross-linking foreign antigens with surface-bound antibodies and C3d on the surface of B cells. The first two modules of CR2, complement control protein modules (CCPs), participate in non-bonded interactions with C3d. We have used computational methods to analyze the dynamic and electrostatic properties of the C3d-CR2(CCP1-2) complex. The interaction between C3d and CR2 is known to depend on pH and ionic strength. Also, the intermodular mobility of the CR2 modules has been questioned before. We performed a 10 ns molecular dynamics simulation to generate a relaxed structure from crystal packing effects for the C3d-CR2(CCP1-2) complex and to study the energetics of the C3d-CR2(CCP1-2) association. The MD simulation suggests a tendency for intermodular twisting in CR2(CCP1-2). We propose a two-step model for recognition and binding of C3d with CR2(CCP1-2), driven by long and short/medium-range electrostatic interactions. We have calculated the matrix of specific short/medium-range pairwise electrostatic free energies of interaction involved in binding and in intermodular communications. Electrostatic interactions may mediate allosteric effects important for C3d-CR2(CCP1-2) association. We present calculations for the pH and ionic strength-dependence of C3d-CR2(CCP1-2) ionization free energies, which are in overall agreement with experimental binding data. We show how comparison of the calculated and experimental data allows for the decomposition of the contributions of electrostatic from other effects in association. We critically compare predicted stabilities for several mutants of the C3d-CR2(CCP1-2) complex with the available experimental data for binding ability. Finally, we propose that CR2(CCP1-2) is capable of assuming a large array of intermodular topologies, ranging from closed V-shaped to open linear states, with similar recognition properties for C3d, but we cannot exclude an additional contact site with C3d.