Free energy simulations for protein ligand binding and stability

Discovery of novel macrocyclic derivatives as potent and selective cyclin-dependent kinase 2 inhibitors

Cyclin-dependent kinase 2 (CDK2) is a member of CDK family of kinases (CDKs) that regulate the cell cycle. Its inopportune or over-activation leads to uncontrolled cell cycle progression and drives numerous types of cancers, especially ovarian, uterine, gastric cancer, as well as those associated with amplified CCNE1 gene. However, developing selective lead compound as CDK2 inhibitors remains challenging owing to similarities in the ATP pockets among different CDKs. Herein, we described the optimization of compound 1, a novel macrocyclic inhibitor targeting CDK2/5/7/9, aiming to discover more selective and metabolically stable lead compound as CDK2 inhibitor. Molecular dynamic (MD) simulations were performed for compound 1 and 9 to gain insights into the improved selectivity against CDK5. Further optimization efforts led to compound 22, exhibiting excellent CDK2 inhibitory activity, good selectivity over other CDKs and potent cellular effects. Based on these characterizations, we propose that compound 22 holds great promise as a potential lead candidate for drug development.

Lysozyme and Human Serum Albumin Proteins as Potential Nitric Oxide Cardiovascular Drug Carriers: Theoretical and Experimental Investigation

Nitric oxide-containing drugs present a critical remedy for cardiovascular diseases. Nitroglycerin (NG, O-NO) and S-nitrosoglutathione (SNG, S-NO) are the most common nitric oxide drugs for cardiovascular diseases. Insights regarding the binding affinity of NO drugs with lysozyme and human serum albumin (HSA) proteins and their dissociation mechanism will provide inquisitive information regarding the potential of the proteins as drug carriers. For the first time, the binding interactions and affinities are investigated using molecular docking, conventional molecular dynamics, steered molecular dynamics, and umbrella sampling to explore the ability of both proteins to act as nitric oxide drug carriers. The molecular dynamics simulation results showed higher stability of lysozyme-drug complexes compared to HSA. For lysozyme, cardiovascular drugs were bound in the protein cavity mainly by the electrostatic and hydrogen bond interactions with residues ASP53, GLN58, ILE59, ARG62, TRP64, ASP102, and TRP109. For HSA, key binding residues were ARG410, TYR411, LYS414, ARG485, GLU450, ARG486, and SER489. The free energy profiles produced from umbrella sampling also suggest that lysozyme-drug complexes had better binding affinity than HSA-drug. Binding characteristics of nitric oxide-containing drugs NG and SNG to lysozyme and HSA proteins were studied using fluorescence and UV-vis absorption spectroscopy. The relative change in the fluorescence intensity as a function of drug concentrations was analyzed using Stern-Volmer calculations. This was also confirmed by the change in the UV-vis spectra. Fluorescence quenching results of both proteins with the drugs, based on the binding constant values, demonstrated significantly weak binding affinity to NG and strong binding affinity to SNG. Both computational and experimental studies provided important data for understanding protein-drug interactions and will aid in developing potential drug carrier systems in cardiovascular diseases.

Discovery of human coronaviruses pan-papain-like protease inhibitors using computational approaches

The papain-like protease (PLpro) is vital for the replication of coronaviruses (CoVs), as well as for escaping innate-immune responses of the host. Hence, it has emerged as an attractive antiviral drug-target. In this study, computational approaches were employed, mainly the structure-based virtual screening coupled with all-atom molecular dynamics (MD) simulations to computationally identify specific inhibitors of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) PLpro, which can be further developed as potential pan-PLpro based broad-spectrum antiviral drugs. The sequence, structure, and functional conserveness of most deadly human CoVs PLpro were explored, and it was revealed that functionally important catalytic triad residues are well conserved among SARS-CoV, SARS-CoV-2, and middle east respiratory syndrome coronavirus (MERS-CoV). The subsequent screening of a focused protease inhibitors database composed of ∼7,000 compounds resulted in the identification of three candidate compounds, ADM_13083841, LMG_15521745, and SYN_15517940. These three compounds established conserved interactions which were further explored through MD simulations, free energy calculations, and residual energy contribution estimated by MM-PB(GB)SA method. All these compounds showed stable conformation and interacted well with the active residues of SARS-CoV-2 PLpro, and showed consistent interaction profile with SARS-CoV PLpro and MERS-CoV PLpro as well. Conclusively, the reported SARS-CoV-2 PLpro specific compounds could serve as seeds for developing potent pan-PLpro based broad-spectrum antiviral drugs against deadly human coronaviruses. Moreover, the presented information related to binding site residual energy contribution could lead to further optimization of these compounds.

Water Dynamics Around Proteins: T- and R-States of Hemoglobin and Melittin

The water dynamics, as characterized by the local hydrophobicity (LH), is investigated for tetrameric hemoglobin (Hb) and dimeric melittin. For the T₀ to R₀ transition in Hb, it is found that LH provides additional molecular-level insight into the Perutz mechanism, i.e., the breaking and formation of salt bridges at the α₁/β₂ and α₂/β₁ interface is accompanied by changes in LH. For Hb in cubic water boxes with 90 and 120 Å edge length it is observed that following a decrease in LH as a consequence of reduced water density or change of water orientation at the protein/water interface the α/β interfaces are destabilized; this is a hallmark of the Perutz stereochemical model for the T to R transition in Hb. The present work thus provides a dynamical view of the classical structural model relevant to the molecular foundations of Hb function. For dimeric melittin, earlier results by Cheng and Rossky [ Nature 1998, 392, 696-699] are confirmed and interpreted on the basis of LH from simulations in which the protein structure is frozen. For the flexible melittin dimer, the changes in the local hydration can be as much as 30% greater than for the rigid dimer, reflecting the fact that protein and water dynamics are coupled.

Effect of Protein Flexibility from Coarse-Grained Elastic Network Parameterizations on the Calculation of Free Energy Profiles of Ligand Binding

The characterization of the affinity and binding mechanism of specific molecules to a protein active site is scientifically and industrially relevant for many applications. In principle, this information can be obtained using molecular dynamics (MD) simulations by calculating the free energy profile of the process. However, this is a computationally demanding calculation. Currently, coarse-grained (CG) force fields are very well implemented for MD simulations of biomolecular systems. These computationally efficient force fields are a major advantage to the study of large model systems and/or those requiring long simulation times. The Martini model is currently one of the most popular CG force fields for these systems. For the specific case of protein simulations, to correctly maintain the macromolecular three-dimensional structure, the Martini model needs to include an elastic network (EN). In this work, the effect of protein flexibility, as induced by three EN models compatible with the Martini force field, was tested on the calculation of free energy profiles for protein-ligand binding. The EN models used were ElNeDyn, GoMartini, and GEN. The binding of triolein (TOG) and triacetin (TAG) to a lipase protein (thermomyces lanuginosa lipase-TLL) was used as a case study. The results show that inclusion of greater flexibility in the CG parameterization of proteins is of high importance in the calculation of the free energy profiles of protein-ligand systems. However, care must be taken in order to avoid unjustified large protein deformations. In addition, due to molecular flexibility there may be no absolute need for the center of the ligand to reach the center of the protein-binding site. The calculation of the energy profile to a distance of about 0.5 nm from the active site center can be sufficient to differentiate the affinity of different ligands to a protein.

Calculating the absolute binding free energy of the insulin dimer in an explicit solvent

Insulin is a significant hormone in the regulation of glucose level in the blood. Its monomers bind to each other to form dimers or hexamers through a complex process. To study the binding of the insulin dimer, we first calculate its absolute binding free energy by the steered molecular dynamics method and the confinement method based on a fictitious thermodynamic cycle. After considering some special correction terms, the final calculated binding free energy at 298 K is −8.97 ± 1.41 kcal mol⁻¹, which is close to the experimental value of −7.2 ± 0.8 kcal mol⁻¹. Furthermore, we discuss the important residue–residue interactions between the insulin monomers, including hydrophobic interactions, π–π interactions and hydrogen bond interactions. The analysis reveals five key residues, Vla^B12, Tyr^B16, Phe^B24, Phe^B25, and Tyr^B26, for the dimerization of the insulin. We also perform MM-PBSA calculations for the wild-type dimer and some mutants and study the roles of the key residues by the change of the binding energy of the insulin dimer.

In this paper, we calculate the absolute binding free energy of an insulin dimer by steered MD method. The result of −8.97 kcal mol⁻¹ is close to the experimental value −7.2 kcal mol⁻¹. We also analyze the residue–residue interactions.

Exploiting Sparse Statistics for a Sequence-Based Prediction of the Effect of Mutations

Recent work showed that there is a significant difference between the statistics of amino acid triplets and quadruplets in sequences of folded proteins and randomly generated sequences. These statistics were used to assign a score to each sequence and make a prediction whether a sequence is likely to fold. The present paper extends the statistics to higher multiplets and suggests a way to handle the treatment of multiplets that were not found in the set of folded proteins. In particular, foldability predictions were done along the line of the previous work using pentuplet statistics and a way was found to combine the quadruplet and pentuplets statistics to improve the foldability predictions. A different, simpler, score was defined for hextuplets and heptuplets and were used to predict the direction of stability change of a protein upon mutation. With the best score combination the accuracy of the prediction was 73.4%.