Molecular dynamic simulations and structure-based pharmacophore development for farnesyltransferase inhibitors discovery

Farnesyltransferase is one of the enzyme targets for the development of drugs for diseases, including cancer, malaria, progeria, etc. In the present study, the structure-based pharmacophore models have been developed from five complex structures (1LD7, 1NI1, 2IEJ, 2ZIR and 2ZIS) obtained from the protein data bank. Initially, molecular dynamic (MD) simulations were performed for the complexes for 10 ns using AMBER 12 software. The conformers of the complexes (75) generated from the equilibrated protein were undergone protein-ligand interaction fingerprint (PLIF) analysis. The results showed that some important residues, such as LeuB96, TrpB102, TrpB106, ArgB202, TyrB300, AspB359 and TyrB361, are predominantly present in most of the complexes for interactions. These residues form side chain acceptor and surface (hydrophobic or π-π) kind of interactions with the ligands present in the complexes. The structure-based pharmacophore models were generated from the fingerprint bits obtained from PLIF analysis. The pharmacophore models have 3-4 pharmacophore contours consist of acceptor and metal ligation (Acc & ML), hydrophobic (HydA) and extended acceptor (Acc2) features with the radius ranging between 1-3 Å for Acc & ML and 1-2 Å for HydA. The excluded volumes of the pharmacophore contours radius are between 1-2 Å. Further, the distance between the interacting groups, root mean square deviation (RMSD), root mean square fluctuation (RMSF) and radial distribution function (RDF) analysis were performed for the MD-simulated proteins using PTRAJ module. The generated pharmacophore models were used to screen a set of natural compounds and database compounds to select significant HITs. We conclude that the developed pharmacophore model can be a significant model for the identification of HITs as FTase inhibitors.

Melatonergic system-based two-gene index is prognostic in human gliomas

Gliomas, the most common primary brain tumors in adults, are classified into four malignancy grades according to morphological features. Recent studies have shown that melatonin treatment induces cytotoxicity in glioma-initiating cells and reduces the invasion and migration of glioma cell lines, inhibiting the nuclear factor κB (NFκB) oncopathway. Given that C6 rat glioma cells produce melatonin, we investigated the correlation between the capacity of gliomas to synthesize/metabolize melatonin and their overall malignancy. We first characterized the melatonergic system of human gliomas cell lines with different grades of aggressiveness (HOG, T98G, and U87MG) and demonstrated that glioma-synthesized melatonin exerts an autocrine antiproliferative effect. Accordingly, the sensitivity to exogenous melatonin was higher for the most aggressive cell line, U87MG, which synthesized/accumulated less melatonin. Using The Cancer Genome Atlas RNAseq data of 351 glioma patients, we designed a predictive model of the content of melatonin in the tumor microenvironment, the ASMT:CYP1B1 index, combining the gene expression levels of melatonin synthesis and metabolism enzymes. The ASMT:CYP1B1 index negatively correlated with tumor grade, as well as with the expression of pro-proliferation and anti-apoptotic NFκB target genes. More importantly, the index was a grade- and histological type-independent prognostic factor. Even when considering only high-grade glioma patients, a low ASMT:CYP1B1 value, which suggests decreased melatonin and enhanced aggressiveness, was strongly associated with poor survival. Overall, our data reveal the prognostic value of the melatonergic system of gliomas and provide insights into the therapeutic role of melatonin.

QM/MM study of the mechanism of reduction of 3-hydroxy-3-methylglutaryl coenzyme A catalyzed by human HMG-CoA reductase

Detailing with atomistic resolution the reaction mechanism of human HMG-CoA reductase (HMG-CoA-R) might provide valuable insights for the development of new cholesterol-lowering drugs.

Prediction of Solvation Free Energies with Thermodynamic Integration Using the General Amber Force Field

Computer-aided drug design (CADD) techniques can be very effective in reducing costs and speeding up drug discovery. The determination of binding and solvation free energies is pivotal for this process and is, therefore, the subject of many studies. In this work, the solvation free energy change (ΔΔGsolv) for a total of 92 transformations in small molecules was predicted using Thermodynamic Integration (TI). It was our aim to compare experimental and calculated solvation free energies for typical and prime additions considered in drug optimizations, analyzing trends, and optimizing a TI protocol. The results showed a good agreement between experimental and predicted values, with an overestimation of the predicted values for CH3, halogens, and NH2, as well as an underestimation for CONH2, but all fall within ±1 kcal/mol. NO2 addition showed a larger and systematic underestimation of the predicted ΔΔGsolv, indicating the need for special attention in these cases. For small molecules, if no experimental data is available, using TI as a theoretical strategy thus appears to be a suitable choice in CADD. It provides a good compromise between time and accuracy.

Benchmarking of Density Functionals for the Accurate Description of Thiol-Disulfide Exchange

A set of 92 density functionals was employed to accurately characterize thiol-disulfide exchange. The properties we have benchmarked throughout the study include the geometry of a 15 atoms model system, the potential energy surface, the activation barrier, and the energy of reaction for thiol-disulfide exchange. Reference energies were determined at the CCSD(T)/CBS//MP2/aug-cc-pVDZ level of theory, and reference geometries were calculated at the MP2/aug-cc-pVTZ level. M11-L, M06-2X, M06-HF, N12-SX, PBE1PBE, PBEh1PBE, and OHSE2PBE described better the geometry of the model system, with average deviations of 0.06 Å in bond lengths (0.06 Å in bond-breaking lengths) and 1.9° in bond angles. On the other hand, the potential energy surface and its gradient were more accurately described by the hybrid density functional BHandH, closely followed by mPW1N, mPW1K, and mPWB1K. The barrier height and energy of reaction were better reproduced by the BMK and M06-2X functionals (deviations of 0.17 and 0.07 kcal·mol(-1), respectively) for a set of 10 Pople's basis sets. MN12-SX and M11-L showed very good results for the widely used 6-311++G(2d,2p) basis set, with deviations of 0.02 and 0.05 kcal·mol(-1), respectively. We studied the effect of the split-valence, diffuse, and polarized functions in the activation barrier of thiol-disulfide exchange, for a set of 10 Pople's basis sets. While increasing the splitting and polarization may increase the activation barrier in approximately 1 kcal·mol(-1), diffuse functions generally contribute to decreasing it no more than 0.10 kcal·mol(-1). In general, 13 functionals provided energies within 1 kcal·mol(-1) of the reference value. The BB1K density functional is one of the best density functionals to characterize thiol-disulfide exchange reactions; however, several density functionals with modified Perdew-Wang exchange and about 40% Hartree-Fock exchange, such as mPW1K, mPW1N, and mPWB1K, show a good performance, too.

Periplasmic nitrate reductase and formate dehydrogenase: similar molecular architectures with very different enzymatic activities

It is remarkable how nature has been able to construct enzymes that, despite sharing many similarities, have simple but key differences that tune them for completely different functions in living cells. Periplasmic nitrate reductase (Nap) and formate dehydrogenase (Fdh) from the DMSOr family are representative examples of this. Both enzymes share almost identical three-dimensional protein foldings and active sites, in terms of coordination number, geometry and nature of the ligands. The substrates of both enzymes (nitrate and formate) are polyatomic anions that also share similar charge and stereochemistry. In terms of the catalytic mechanism, both enzymes have a common activation mechanism (the sulfur-shift mechanism) that ensures a constant coordination number around the metal ion during the catalytic cycle. In spite of these similarities, they catalyze very different reactions: Nap abstracts an oxygen atom from nitrate releasing nitrite, whereas FdH catalyzes a hydrogen atom transfer from formate and releases carbon dioxide. In this Account, a critical analysis of structure, function, and catalytic mechanism of the molybdenum enzymes periplasmic nitrate reductase (Nap) and formate dehydrogenase (Fdh) is presented. We conclude that the main structural driving force that dictates the type of reaction, catalyzed by each enzyme, is a key difference on one active site residue that is located in the top region of the active sites of both enzymes. In both enzymes, the active site is centered on the metal ion of the cofactor (Mo in Nap and Mo or W in Fdh) that is coordinated by four sulfur atoms from two pyranopterin guanosine dinucleotide (PGD) molecules and by a sulfido. However, while in Nap there is a Cys directly coordinated to the Mo ion, in FdH there is a SeCys instead. In Fdh there is also an important His that interacts very closely with the SeCys, whereas in Nap the same position is occupied by a Met. The role of Cys in Nap and SeCys in FdH is similar in both enzymes; however, Met and His have different roles. His participates directly on catalysis, and it is therefore detrimental for the catalytic cycle of FdH. Met only participates in substrate binding. We concluded that this small but key difference dictates the type of reaction that is catalyzed by each enzyme. In addition, it allows explaining why formate can bind in the Nap active site in the same way as the natural substrate (nitrate), but the reaction becomes stalled afterward.

Receptor-based virtual screening protocol for drug discovery

Computational aided drug design (CADD) is presently a key component in the process of drug discovery and development as it offers great promise to drastically reduce cost and time requirements. In the pharmaceutical arena, virtual screening is normally regarded as the top CADD tool to screen large libraries of chemical structures and reduce them to a key set of likely drug candidates regarding a specific protein target. This chapter provides a comprehensive overview of the receptor-based virtual screening process and of its importance in the present drug discovery and development paradigm. Following a focused contextualization on the subject, the main stages of a virtual screening campaign, including its strengths and limitations, are the subject of particular attention in this review. In all of these stages special consideration will be given to practical issues that are normally the Achilles heel of the virtual screening process.

Re(I) and Tc(I) complexes for targeting nitric oxide synthase: influence of the chelator in the affinity for the enzyme

Aiming to design (99m) Tc complexes for probing nitric oxide synthase (NOS) by SPECT, we synthesized conjugates (L4-L6) comprising a NOS-recognizing moiety connected to a diamino-propionic acid (dap) chelating unit. The conjugates led to complexes of the type fac-[M(CO)3 (ĸ(3) -L)] (M = Re/(99m) Tc; Re4/Tc4: L = L4; Re5/Tc5: L = L5; Re6/Tc6: L = L6). Enzymatic studies showed that L4 and L5, but not L6, gave complexes (Re4 and Re5) that are less potent than the conjugates. To rationalize these results, we performed docking and molecular dynamics simulations. The high affinity of L4 and L5 is due to the strong interactions between the dap chelator and polar residues of the binding cavity. These interactions are hampered by metallation resulting in complexes with lower affinity. The higher potency of Re5 compared to Re4 was assigned to the increased bulkiness of Re5 and the presence of additional anchoring groups that better fit the active site and provide more extensive contacts. In turn, Re6 is too bulky and its organometallic tail is oriented toward the peripheral pocket of iNOS, leading to loss of contacts and a lower affinity. These results were compared with our previous results obtained with analogue complexes stabilized by a pyrazolyl-diamine chelating unit.

Relevant interactions of antimicrobial iron chelators and membrane models revealed by nuclear magnetic resonance and molecular dynamics simulations

The dynamics and interaction of 3-hydroxy-4-pyridinone fluorescent iron chelators, exhibiting antimicrobial properties, with biological membranes were evaluated through NMR and molecular dynamics simulations. Both NMR and MD simulation results support a strong interaction of the chelators with the lipid bilayers that seems to be strengthened for the rhodamine containing compounds, in particular for compounds that include ethyl groups and a thiourea link. For the latter type of compounds the interaction reaches the hydrophobic core of the lipid bilayer. The molecular docking and MD simulations performed for the potential interaction of the chelators with DC-SIGN receptors provide valuable information regarding the cellular uptake of these compounds since the results show that the fluorophore fragment of the molecular framework is essential for an efficient binding. Putting together our previous and present results, we put forward the hypothesis that all the studied fluorescent chelators have access to the cell, their uptake occurs through different pathways and their permeation properties correlate with a better access to the cell and its compartments and, consequently, with the chelators antimicrobial properties.

A new scoring function for protein-protein docking that identifies native structures with unprecedented accuracy

Protein-protein (P-P) 3D structures are fundamental to structural biology and drug discovery. However, most of them have never been determined. Many docking algorithms were developed for that purpose, but they have a very limited accuracy in generating native-like structures and identifying the most correct one, in particular when a single answer is asked for. With such a low success rate it is difficult to point out one docked structure as being native-like. Here we present a new, high accuracy, scoring method to identify the 3D structure of P-P complexes among a set of trial poses. It incorporates alanine scanning mutagenesis experimental data that need to be obtained a priori. The scoring scheme works by matching the computational and the experimental alanine scanning mutagenesis results. The size of the trial P-P interface area is also taken into account. We show that the method ranks the trial structures and identifies the native-like structures with unprecedented accuracy (∼94%), providing the correct P-P 3D structures that biochemists and molecular biologists need to pursue their studies. With such a success rate, the bottleneck of protein-protein docking moves from the scoring to searching algorithms.

Insights into the structural determinants of substrate specificity and activity in mouse aldehyde oxidases

In this work, a combination of homology modeling and molecular dynamics (MD) simulations was used to investigate the factors that modulate substrate specificity and activity of the mouse AOX isoforms: mAOX1, mAOX2 (previously mAOX3l1), mAOX3 and mAOX4. The results indicate that the AOX isoform structures are highly preserved and even more conserved than the corresponding amino acid sequences. The only differences are at the protein surface and substrate-binding site region. The substrate-binding site of all isoforms consists of two regions: the active site, which is highly conserved among all isoforms, and a isoform-specific region located above. We predict that mAOX1 accepts a broader range of substrates of different shape, size and nature relative to the other isoforms. In contrast, mAOX4 appears to accept a more restricted range of substrates. Its narrow and hydrophobic binding site indicates that it only accepts small hydrophobic substrates. Although mAOX2 and mAOX3 are very similar to each other, we propose the following pairs of overlapping substrate specificities: mAOX2/mAOX4 and mAOX3/mAXO1. Based on these considerations, we propose that the catalytic activity between all isoforms should be similar but the differences observed in the binding site might influence the substrate specificity of each enzyme. These results also suggest that the presence of several AOX isoforms in mouse allows them to oxidize more efficiently a wider range of substrates. This contrasts with the same or other organisms that only express one isoform and are less efficient or incapable of oxidizing the same type of substrates.

The catalytic mechanism of carboxylesterases: a computational study

The catalytic mechanism of carboxylesterases (CEs, EC 3.1.1.1) is explored by computational means. CEs hydrolyze ester, amide, and carbamate bonds found in xenobiotics and endobiotics. They can also perform transesterification, a reaction important, for instance, in cholesterol homeostasis. The catalytic mechanisms with three different substrates (ester, thioester, and amide) have been established at the M06-2X/6-311++G**//B3LYP/6-31G* level of theory. It was found that the reactions proceed through a mechanism involving four steps instead of two as is generally proposed: (i) nucleophilic attack of serine to the substrate, forming the first tetrahedral intermediate, (ii) formation of the acyl-enzyme complex concomitant with the release of the alcohol product, (iii) nucleophilic attack of a water or alcohol molecule forming the second tetrahedral intermediate, and (iv) the release of the second product of the reaction. The results agree very well with the available experimental data and show that the hydrolytic and the transesterification reactions are competitive processes when the substrate is an ester. In all the other studied substrates (thioester or amide), the hydrolytic and transesterification process are less favorable and some of them might not even take place under in vivo conditions.

Farnesyltransferase inhibitors: a comprehensive review based on quantitative structural analysis

Farnesyltransferase inhibitors (FTIs) have mainly been used in cancer therapy. However, more recently, investigations on these inhibitors revealed that FTIs can be used for the treatment of other diseases such as Progeria, P. falciparum resistant malaria, Trypnosomatid, etc. Hence the development of novel FTIs is an important task for the drug discovery program. Initially, numerous peptidomimetic FTIs were developed from the template of CAAX (CVIM was the first pharmacophore model used as a peptidomimetic). Later, many non-peptidomimetic FTIs have been discovered with the structural modification of the peptidomimetics. The structural analysis of those developed FTIs by various researchers suggested that the presence of a heterocycle or a polar group in place of the thiol group is required for interaction with the Zn(2+) ion. The bulky naphthyl, quinolinyl, phenyl, phenothazine, etc in this position provide better hydrophobicity to the molecules which interact with the aromatic amino acid moieties in the hydrophobic pocket. A hydrophilic region with polar groups is necessary for the polar or hydrogen bonding interactions with the amino acids or water molecules in the active site. Many FTIs have been isolated from natural products, which possessed inhibitory activity against farnesyltransferase (FTase). Among them, pepticinnamin E (9R), fusidienol (9T), gliotoxin (9V), cylindrol A (9X), etc possessed potential FTase inhibitory activities and their structural features are comparable to those of the synthetic molecules. The clinical studies progressing on FTIs showed that tipifarnib in combination with bortezomib is used for the treatment of patients with advanced acute leukemias. Successful phase I and II studies are undergoing for tipifarnib alone or in combination with other drugs/radiation for the treatment of multiple myeloma, AML, breast cancer, mantle cell lymphoma, solid tumors, non-small cell lung cancer (NSCLC), pancreatic cancer, glioblastoma, etc. Phase I pharmacokinetic (maximum tolerated dose, toxicity) and pharmacodynamic studies of AZD3409 (an orally active double prodrug) is progressing on patients with solid malignancies taking 500 mg once a day. A phase II study is undergoing on lonafarnib alone and in combination with zoledronic acid and pravastatin for the treatment of Hutchinson-Gilford Progeria syndrome (HGPS) and progeroid laminopathies. Lonafarnib therapy improved cardiovascular status of children with HGPS, by improved peripheral arterial stiffness, bone structure and audiological status in the patients. Other important FTIs such as BMS-214662, LB42908, LB42708, etc are under clinical studies for the treatment of various cancers. This review concluded that the quantitative structural analysis report with an elaborative study on the natural product compounds provides ideas for development of novel molecules for the FTase inhibitory activity. The fragment based analysis is also needed to select the substituents, which provides significant inhibitory activities and can also have good pharmacokinetic properties in the clinical studies.

Development of ribonucleotide reductase inhibitors: a review on structure activity relationships

Ribonucleotide reductase (RNR, E.C. 1.17.4.1), which is composed of two dissimilar proteins (subunits), often referred as R1 (containing polythiols) and R2 (containing non-heme iron and a free tyrosyl radical), which contribute to the role played by the enzyme. RNRs are one of the important targets in anticancer and antiviral drug development and many RNR inhibitors have been discovered at the end of the 20(th) century; many of them are already in clinical use. Triapine (3-AP) is one of the important RNR inhibitors belonging to the class of thiosemicarbazone derivatives, used in the treatment of various cancers. The structure activity relationship (SAR) studies on the investigated RNR inhibitors showed that the nitrogen atom in the pyridine (or other heterocycles) forms coordination complexes with the metal ions along with the imine, oxo and thio atoms of the thiosemicarbazone or semicarbazone pharmacophores. The computational analyses results in the adenine and purine derivatives suggest that the nitrogen atoms in the adenine rings make several hydrogen bonds with the water molecules present in the active site, as well as Gly249 and Glu288 residues. The OH group in third position of the sugar moiety interacts with the Ser217 (C=O) and the water molecules through hydrogen bonds. The aromatic rings in the molecules interact with the tyrosine residues. The thiosemicarbazone or semicarbazone derivatives explain that the flexibility and polar properties in the thiosemicarbazone or semicarbazone pharmacophoric regions allow the molecules to coordinate with the metal ion (especially iron) present in the RNR enzymes. This review concluded that RNR inhibitors composed of different fragments such as aryl, heteroaryl, sugar moiety, polar groups, flexible bonds, etc which are required for the binding of the molecules to the RNR enzymes. Further, the fragmental analysis of the RNR inhibitors on different toxicological and metabolic targets can provide significant novel molecules with acceptable pharmacokinetic properties.

Discovery of new sites for drug binding to the hypertension-related renin-angiotensinogen complex

Renin (REN) is a key drug target to stop the hypertension cascade, but thus far only one direct inhibitor has been made commercially available. In this study, we assess an innovative REN inhibition strategy, by targeting the interface of the renin:angiotensinogen (REN:ANG) complex. We characterized the energetic role of interfacial residues of REN:ANG and identified the ones responsible for protein:protein binding, which can serve as drug targets for disruption of the REN:ANG association. For this purpose, we applied a computational alanine scanning mutagenesis protocol, which measures the contribution of each side chain for the protein:protein binding free energy with an accuracy of ≈ 1 kcal/mol. As a result, in REN and ANG, six and eight residues were found to be critical for binding, respectively. The leading force behind REN:ANG complexation was found to be the hydrophobic effect. The binding free energy per residue was found to be proportional to the buried area. Residues responsible for binding were occluded from water at the complex, which promotes an efficient pairing between the two proteins. Two druggable pockets involving critical residues for binding were found on the surface of REN, where small druglike molecules can bind and disrupt the ANG:REN association that may provide an efficient way to achieve REN inhibition and control hypertension.

Biomolecular structure manipulation using tailored electromagnetic radiation: a proof of concept on a simplified model of the active site of bacterial DNA topoisomerase

We report on the viability of breaking selected bonds in biological systems using tailored electromagnetic radiation.

Classification study of solvation free energies of organic molecules using machine learning techniques

Classification models to predict the solvation free energies of organic molecules were developed using decision tree, random forest and support vector machine approaches and with MACCS fingerprints, MOE and PaDEL descriptors.

Binding mode prediction and identification of new lead compounds from natural products as renin and angiotensin converting enzyme inhibitors

In this study a novel renin and ACE inhibitor was developed from natural products using computational techniques. Molecular dynamic simulations showed that the new lead compound has significant binding to the targets.

The sulfur shift: an activation mechanism for periplasmic nitrate reductase and formate dehydrogenase

A structural rearrangement known as sulfur shift occurs in some Mo-containing enzymes of the DMSO reductase family. This mechanism is characterized by the displacement of a coordinating cysteine thiol (or SeCys in Fdh) from the first to the second shell of the Mo-coordination sphere metal. The hexa-coordinated Mo ion found in the as-isolated state cannot bind directly any exogenous ligand (substrate or inhibitors), while the penta-coordinated ion, attained upon sulfur shift, has a free binding site for direct coordination of the substrate. This rearrangement provides an efficient mechanism to keep a constant coordination number throughout an entire catalytic pathway. This mechanism is very similar to the carboxylate shift observed in Zn-dependent enzymes, and it has been recently detected by experimental means. In the present paper, we calculated the geometries and energies involved in the sulfur-shift mechanism using QM-methods (M06/(6-311++G(3df,2pd),SDD)//B3LYP/(6-31G(d),SDD)). The results indicated that the sulfur-shift mechanism provides an efficient way to enable the metal ion for substrate coordination.

Comparative analysis of the performance of commonly available density functionals in the determination of geometrical parameters for copper complexes

In this study, a set of 50 transition-metal complexes of Cu(I) and Cu(II), were used in the evaluation of 18 density functionals in geometry determination. In addition, 14 different basis sets were considered, including four commonly used Pople's all-electron basis sets; four basis sets including popular types of effective-core potentials: Los Alamos, Steven-Basch-Krauss, and Stuttgart-Dresden; and six triple-ζ basis sets. The results illustrate the performance of different methodological alternatives for the treatment of geometrical properties in relevant copper complexes, pointing out Double-Hybrid (DH) and Long-range Correction (LC) Generalized Gradient Approximation (GGA) methods as better descriptors of the geometry of the evaluated systems. These however, are associated with a computational cost several times higher than some of the other methods employed, such as the M06 functional, which has also demonstrated a comparable performance. Regarding the basis sets, 6-31+G(d) and 6-31+G(d,p) were the best performing approaches. In addition, the results show that the use of effective-core potentials has a limited impact, in terms of the accuracy in the determination of metal-ligand bond-lengths and angles in our dataset of copper complexes. Hence, these could become a good alternative for the geometrical description of these systems, particularly CEP-121G and SDD basis sets, if one is considering larger copper complexes where the computational cost could be an issue.

Molecular dynamics studies on both bound and unbound renin protease

The aspartic protease renin (REN) catalyses the rate-limiting step in the Renin-Angiotensin-Aldosterone System (RAAS), which regulates cardiovascular and renal homoeostasis in living organisms. Renin blockage is therefore an attractive therapeutic strategy for the treatment of hypertension. Herein, computational approaches were used to provide a structural characterization of the binding site, flap opening and dynamic rearrangements of REN in the key conserved residues and water molecules, with the binding of a dodecapeptide substrate or different inhibitors. All these structural insights during catalysis may assist future studies in developing novel strategies for REN inactivation. Our molecular dynamics simulations of several unbound-REN and bound-REN systems indicate similar flexible-segments plasticity with larger fluctuations in those belonging to the C-domain (exposed to the solvent). These segments are thought to assist the flap opening and closure to allow the binding of the substrate and catalytic water molecules. The unbound-REN simulation suggests that the flap can acquire three different conformations: closed, semi-open and open. Our results indicate that the semi-open conformation is already sufficient and appropriate for the binding of the angiotensinogen (Ang) tail, thus contributing to the high specificity of REN, and that both semi-open and open flap conformations are present in free and complexed enzymes. We additionally observed that the Tyr75-Trp39 H-bond has an important role in assisting flap movement, and we highlight several conserved water molecules and amino acids that are essential for the proper catalytic activity of REN.

Computational Alanine Scanning Mutagenesis: MM-PBSA vs TI

Understanding protein-protein association and being able to determine the crucial residues responsible for their association (hot-spots) is a key issue with huge practical applications such as rational drug design and protein engineering. A variety of computational methods exist to detect hot-spots residues, but the development of a fast and accurate quantitative alanine scanning mutagenesis (ASM) continues to be crucial. Using four protein-protein complexes, we have compared a variation of the standard computational ASM protocol developed at our group, based on the Molecular Mechanics/Poisson-Boltzmann Surface Area (MM-PBSA) approach, against Thermodynamic Integration (TI), a well-known and accurate but computationally expensive method. To compare the efficiency and the accuracy of the two methods, we have calculated the protein-protein binding free energy differences upon alanine mutation of interfacial residues (ΔΔGbind). In relation to the experimental ΔΔGbind values, the average error obtained with TI was 1.53 kcal/mol, while the ASM protocol resulted in an average error of 1.18 kcal/mol. The results demonstrate that the much faster ASM protocol gives results at the same level of accuracy as the TI method but at a fraction of the computational time required to run TI. This ASM protocol is therefore a strong and efficient alternative to the systematic evaluation of protein-protein interfaces, involving hundreds of amino acid residues in search of hot-spots.