Effect of allosteric molecules on structure and drug affinity of HIV-1 protease by molecular dynamics simulations

Recent experiments show that small molecules can bind onto the allosteric sites of HIV-1 protease (PR), which provides a starting point for developing allosteric inhibitors. However, the knowledge of the effect of such binding on the structural dynamics and binding free energy of the active site inhibitor and PR is still lacking. Here, we report 200 ns long molecular dynamics simulation results to gain insight into the influences of two allosteric molecules (1H-indole-6-carboxylic acid, 1F1 and 2-methylcyclohexano, 4D9). The simulations demonstrate that both allosteric molecules change the PR conformation and stabilize the structures of PR and the inhibitor; the residues of the flaps are sensitive to the allosteric molecules and the flexibility of the residues is pronouncedly suppressed; the additions of the small molecules to the allosteric sites strengthen the binding affinities of 3TL-PR by about 12–15 kal/mol in the binding free energy, which mainly arises from electrostatic term. Interestingly, it is found that the action mechanisms of 1F1 and 4D9 are different, the former behaviors like a doorman that keeps the inhibitor from escape and makes the flaps (door) partially open; the latter is like a wedge that expands the allosteric space and meanwhile closes the flaps. Our data provide a theoretical support for designing the allosteric inhibitor.     

Identification of phenoxyacetamide derivatives as novel DOT1L inhibitors via docking screening and molecular dynamics simulation

Dot1-like protein (DOT1L) is a histone methyltransferase that has become a novel and promising target for acute leukemias bearing mixed lineage leukemia (MLL) gene rearrangements. In this study, a hierarchical docking-based virtual screening combined with molecular dynamic (MD) simulation was performed to identify DOT1L inhibitors with novel scaffolds. Consequently, 8 top-ranked hits were eventually identified and were further subjected to MD simulation. It was indicated that all hits could reach equilibrium with DOT1L in the MD simulation and further binding free energy calculations suggested that phenoxyacetamide-derived hits such as L01, L03, L04 and L05 exhibited remarkably higher binding affinity compared to other hits. Among them, L03 showed both the lowest glide score (−12.281) and the most favorable binding free energy (−303.9 +/− 16.5 kJ/mol), thereby making it a promising lead for further optimization.     

Nanosecond molecular dynamics simulations of Cdc25B and its complex with a 1,4-naphthoquinone inhibitor: implications for rational inhibitor design

Cdc25 phosphatases have been considered as attractive drug targets for anticancer therapies due to the correlation of their overexpression with a wide variety of cancers. To gain insight into designing new potent inhibitors, we investigate the dynamic properties of Cdc25B and its complex with a 1,4-naphtoquinone inhibitor NSC 95397 by means of molecular dynamics simulations in aqueous solution. It is shown from the calculated dynamic properties that the malleability of the residues 530–532 residing at the start of C-terminal region around the active site should be responsible for the catalytic action of Cdc25B. However, binding of the inhibitor in the active site leads to a substantial decrease in the motional amplitude of the flexible residues, due to the hydrophobic interactions with the side chain of Met531. The simulation results also indicate that at least four hydrogen bonds are involved in the enzyme-inhibitor complex. Among them, the hydrogen bond between the side chain carboxylate group of Glu478 and one of the hydroxyl groups of the inhibitor is found to be the most significant binding force stabilizing the inhibitor in the active site. This result supports the previous experimental implication that the possession of a single hydroxyl group is sufficient for the inhibitory activity of 1,4-naphthoquinone inhibitors.     

Binding mechanism of CDK5 with roscovitine derivatives based on molecular dynamics simulations and MM/PBSA methods

Roscovitine derivatives are potent inhibitors of cyclin-dependent kinase 5 (CDK5), but they exhibit different activities, which has not been understood clearly up to now. On the other hand, the task of drug design is difficult because of the fuzzy binding mechanism. In this context, the methods of molecular docking, molecular dynamics (MD) simulation, and binding free energy analysis are applied to investigate and reveal the detailed binding mechanism of four roscovitine derivatives with CDK5. The electrostatic and van der Waals interactions of the four inhibitors with CDK5 are analyzed and discussed. The calculated binding free energies in terms of MM-PBSA method are consistent with experimental ranking of inhibitor effectiveness for the four inhibitors. The hydrogen bonds of the inhibitors with Cys83 and Lys33 can stabilize the inhibitors in binding sites. The van der Waals interactions, especially the pivotal contacts with Ile10 and Leu133 have larger contributions to the binding free energy and play critical roles in distinguishing the variant bioactivity of four inhibitors. In terms of binding mechanism of the four inhibitors with CDK5 and energy contribution of fragments of each inhibitor, two new CDK5 inhibitors are designed and have stronger inhibitory potency.     

Identification of novel allosteric modulator binding sites in NMDA receptors: A molecular modeling study

The dysfunction of N-methyl-d-Aspartate receptors (NMDARs), a subtype of glutamate receptors, is correlated with schizophrenia, stroke, and many other neuropathological disorders. However, not all NMDAR subtypes equally contribute towards these disorders. Since NMDARs composed of different GluN2 subunits (GluN2A-D) confer varied physiological properties and have different distributions in the brain, pharmacological agents that target NMDARs with specific GluN2 subunits have significant potential for therapeutic applications. In our previous research, we have identified a family of novel allosteric modulators that differentially potentiate and/or inhibit NMDARs of differing GluN2 subunit composition. To further elucidate their molecular mechanisms, in the present study, we have identified four potential binding sites for novel allosteric modulators by performing molecular modeling, docking, and in silico mutations. The molecular determinants of the modulator binding sites (MBS), analysis of particular MBS electrostatics, and the specific loss or gain of binding after mutations have revealed modulators that have strong potential affinities for specific MBS on given subunits and the role of key amino acids in either promoting or obstructing modulator binding. These findings will help design higher affinity GluN2 subunit-selective pharmaceuticals, which are currently unavailable to treat psychiatric and neurological disorders.     

A mechanistic approach to explore novel HDAC1 inhibitor using pharmacophore modeling, 3D- QSAR analysis,molecular docking,density functional and molecular dynamics simulation study

Deregulated epigenetic activity of Histone deacetylase 1 (HDAC1) in tumor development and carcinogenesis pronounces it as promising therapeutic target for cancer treatment. HDAC1 has recently captured the attention of researchers owing to its decisive role in multiple types of cancer. In the present study a multistep framework combining ligand based 3D-QSAR, molecular docking and Molecular Dynamics (MD) simulation studies were performed to explore potential compound with good HDAC1 binding affinity. Four different pharmacophore hypotheses Hypo1 (AADR), Hypo2 (AAAH), Hypo3 (AAAR) and Hypo4 (ADDR) were obtained. The hypothesis Hypo1 (AADR) with two hydrogen bond acceptors (A), one hydrogen bond donor (D) and one aromatics ring (R) was selected to build 3D-QSAR model on the basis of statistical parameter. The pharmacophore hypothesis produced a statistically significant QSAR model, with co-efficient of correlation r² = 0.82 and cross validation correlation co-efficient q² = 0.70. External validation result displays high predictive power with r² (o) value of 0.88 and r² (m) value of 0.58 to carry out further in silico studies. Virtual screening result shows ZINC70450932 as the most promising lead where HDAC1 interacts with residues Asp99, His178, Tyr204, Phe205 and Leu271 forming seven hydrogen bonds. A high docking score (−11.17 kcal/mol) and lower docking energy −37.84 kcal/mol) displays the binding efficiency of the ligand. Binding free energy calculation was done using MM/GBSA to access affinity of ligands towards protein. Density Functional Theory was employed to explore electronic features of the ligands describing intramolcular charge transfer reaction. Molecular dynamics simulation studies at 50 ns display metal ion (Zn)-ligand interaction which is vital to inhibit the enzymatic activity of the protein.     

Q: a molecular dynamics program for free energy calculations and empirical valence bond simulations in biomolecular systems

A new molecular dynamics program for free energy calculations in biomolecular systems is presented. It is principally designed for free energy perturbation simulations, empirical valence bond calculations, and binding affinity estimation by linear interaction energy methods. Evaluation of ligand-binding selectivity and free energy profiles for nucleophile activation in two protein tyrosine phosphatases as well as absolute binding affinity estimation for a lysine-binding protein are given as examples.     

A comparative study based on docking and molecular dynamics simulations over HDAC-tubulin dual inhibitors

Nowadays the ability to prediction of complex behavior rationally based on the computational approaches has been a successful technique in drug discovery. In the present study interactions of a new series of hybrids, which were made by linking colchicine as a tubulin inhibitor and suberoylanilide hydroxamic acid (SAHA) as a HDAC inhibitor, with HDAC8 and HDAC1 were investigated and compared. This research has been facilitated by the availability of experimental information besides employing docking methodology as well as classical molecular dynamics simulations and binding free energy calculation were performed. The obtained findings indicate different modes of interactions and inhibition strengths of the studied inhibitors for HDAC8 and HDAC1. HDAC8 binding free energies (−34.35 to −26.27 kcal/mol) revealed higher binding affinity to HDAC8 compared to HDAC1 (−33.17 to −7.99 kcal/mol). The binding energy contribution of each residue with the hybrid compounds 4a-4e within the active site of HDAC1 and HDAC8 was analyzed and the results confirmed the rule of key amino acids in interaction with the hybrid compounds.     

An effective HIV-1 integrase inhibitor screening platform: Rationality validation of drug screening,conformational mobility and molecular recognition analysis for PFV integrase complex with viral DNA

As an important target for the development of novel anti-AIDS drugs, HIV-1 integrase (IN) has been widely concerned. However, the lack of a complete accurate crystal structure of HIV-1 IN greatly blocks the discovery of novel inhibitors. In this work, an effective HIV-1 IN inhibitor screening platform, namely PFV IN, was filtered from all species of INs. Next, the 40.8% similarity with HIV-1 IN, as well as the high efficiency of virtual screening and the good agreement between calculated binding free energies and experimental ones all proved PFV IN is a promising screening platform for HIV-1 IN inhibitors. Then, the molecular recognition mechanism of PFV IN by its substrate viral DNA and six naphthyridine derivatives (NRDs) inhibitors was investigated through molecular docking, molecular dynamics simulations and water-mediated interactions analyses. The functional partition of NRDs IN inhibitors could be divided into hydrophobic and hydrophilic ones, and the Mg²⁺ ions, water molecules and conserved DDE motif residues all interacted with the hydrophilic partition, while the bases in viral DNA and residues like Tyr212, Pro214 interacted with the hydrophobic one. Finally, the free energy landscape (FEL) and cluster analyses were performed to explore the molecular motion of PFV IN-DNA system. It is found that the association with NRDs inhibitors would obviously decrease the motion amplitude of PFV IN-DNA, which may be one of the most potential mechanisms of IN inhibitors. This work will provide a theoretical basis for the inhibitor design based on the structure of HIV-1 IN.     

Computational determination of binding structures and free energies of glucose 6-phosphate dehydrogenase with novel steroid inhibitors

Glucose 6-phosphate dehydrogenase (G6PD), the first and the rate-limiting enzyme in the pentose phosphate pathway (PPP), catalyzes the oxidation of G6P to 6-phosphogluconolactone and the reduction of NADP⁺ to NADPH. Its key role in cancer promotes the development of a potent and selective inhibitor that might increase cancer cell death when combined with radiotherapy. In the present study, we investigated the detailed binding modes and binding free energies for G6PD interacting with a promising series of recently developed inhibitors, i.e., the steroid derivatives, by performing molecular docking, molecular dynamics (MD) simulations, and binding free energy calculations. The docking indicates that the inhibitors occupy the binding sites of both G6P and NADP⁺. The calculated binding free energies on the basis of the MD-simulated enzyme–inhibitor complexes are in good agreement with the experimental activity data for all of the examined inhibitors. The valuable insights into the detailed enzyme–inhibitor binding including the important intermolecular interactions, e.g., the hydrogen bond interaction and the hydrophobic interaction, have been provided. The computational results provide new insights into future rational design of more potent inhibitors of G6PD as a treatment for cancer.     

Structural insights into the camel milk lactoperoxidase: Homology modeling and molecular dynamics simulation studies

Lactoperoxodase (LPO) is a heme peroxidase enzyme present in mammalian milk. It is an antimicrobial protein with wide range of industrial applications. Although the three dimensional structure of LPO from various mammalian species has been reported, but its structure from camel source is still unknown. So far, the crystallization attempts have not been successful in determining camel LPO (cLPO) structure. Herein, we developed the three dimensional structure of cLPO by homology modeling approach using prime module available in Schrodinger suite. The developed model in complex with ligand hypothiocyanate (OSCN⁻) was further validated by Ramachandran plot followed by molecular dynamics (MD) simulation studies using Desmond module of Schrodinger. cLPO model exhibited overall structural similarity with template crystal structure, however, it displayed different interaction pattern of amino acid residues with ligand OSCN⁻ in comparison to template crystal structure. Moreover, the ligand binding site environment in cLPO is more polar, less hydrophobic, and harbours more number of charged residues than template crystal structure. The substrate binding pocket environment of cLPO shows a considerable difference from template crystal structure. This subsequently resulted in dissimilar behaviour of ligand during the course of MD simulation studies.     

Combined in silico approaches for the identification of novel inhibitors of human islet amyloid polypeptide (hIAPP) fibrillation

Human islet amyloid polypeptide (hIAPP) is a natively unfolded polypeptide hormone of glucose metabolism, which is co-secreted with insulin by the β-cells of the pancreas. In patients with type 2 diabetes, IAPP forms amyloid fibrils because of diabetes-associated β-cells dysfunction and increasing fibrillation, in turn, lead to failure of secretory function of β-cells. This provides a target for the discovery of small organic molecules against protein aggregation diseases. However, the binding mechanism of these molecules with monomers, oligomers and fibrils to inhibit fibrillation is still an open question. In this work, ligand and structure-based in silico approaches were used to identify novel fibrillation inhibitors and/or fibril binding compounds. The best pharmacophore model was used as a 3D search query for virtual screening of a compound database to identify novel molecules having the potential to be therapeutic agents against protein aggregation diseases. Docking and molecular dynamics simulation studies were used to explore the interaction pattern and mechanism of the identified novel small molecules with predicted hIAPP structure, its aggregation prone conformation and fibril forming segments. We show that catechins with galloyl group and molecules having two to three planar apolar rings bind to hIAPP structures and fibril forming segments with greater affinity. The differences in binding affinities of different compounds against several fibril forming segments of the peptide suggest that a mixture of active compounds may be required for treatment of aggregation diseases.     

Computational docking,molecular dynamics simulation and subsite structure analysis of a maltogenic amylase from Bacillus lehensis G1 provide insights into substrate and product specificity

Maltogenic amylase (MAG1) from Bacillus lehensis G1 displayed the highest hydrolysis activity on β-cyclodextrin (β-CD) to produce maltose as a main product and exhibited high transglycosylation activity on malto-oligosaccharides with polymerization degree of three and above. These substrate and product specificities of MAG1 were elucidated from structural point of view in this study. A three-dimensional structure of MAG1 was constructed using homology modeling. Docking of β-CD and malto-oligosaccharides was then performed in the MAG1 active site. An aromatic platform in the active site was identified which is responsible in substrate recognition especially in determining the enzyme’s preference toward β-CD. Molecular dynamics (MD) simulation showed MAG1 structure is most stable when docked with β-CD and least stable when docked with maltose. The docking analysis and MD simulation showed that the main subsites for substrate stabilization in the active site are −2, −1, +1 and +2. A bulky residue, Trp359 at the +2 subsite was identified to cause steric interference to the bound linear malto-oligosaccharides thus prevented it to occupy subsite +3, which can only be reached by a highly bent glucose molecule such as β-CD. The resulted modes of binding from docking simulation show a good correlation with the experimentally determined hydrolysis pattern. The subsite structure generated from this study led to a possible mode of action that revealed how maltose was mainly produced during hydrolysis. Furthermore, maltose only occupies subsite +1 and +2, therefore could not be hydrolyzed or transglycosylated by the enzyme. This important knowledge has paved the way for a novel structure-based molecular design for modulation of its catalytic activities.     

Structural analysis of a novel N-carbamoyl-d-amino acid amidohydrolase from a Brazilian Bradyrhizobium japonicum strain: In silico insights by molecular modelling,docking and molecular dynamics

In this work we performed several in silico analyses to describe the relevant structural aspects of an enzyme N-Carbamoyl-d-amino acid amidohydrolase (d-NCAase) encoded on the genome of the Brazilian strain CPAC 15 (=SEMIA 5079) of Bradyrhizobium japonicum, a nonpathogenic species belonging to the order Rhizobiales. d-NCAase has wide applications particularly in the pharmaceutical industry, since it catalyzes the production of d-amino acids such as D-p-hydroxyphenylglycine (D-HPG), an intermediate in the synthesis of β-lactam antibiotics. We applied a homology modelling approach and 50 ns of molecular dynamics simulations to predict the structure and the intersubunit interactions of this novel d-NCAase. Also, in order to evaluate the substrate binding site, the model was subjected to 50 ns of molecular dynamics simulations in the presence of N-Carbamoyl-d-p-hydroxyphenylglycine (Cp-HPG) (a d-NCAase canonical substrate) and water-protein/water-substrate interactions analyses were performed. Overall, the structural analysis and the molecular dynamics simulations suggest that d-NCAase of B. japonicum CPAC-15 has a homodimeric structure in solution. Here, we also examined the substrate specificity of the catalytic site of our model and the interactions with water molecules into the active binding site were comprehensively discussed. Also, these simulations showed that the amino acids Lys123, His125, Pro127, Cys172, Asp174 and Arg176 are responsible for recognition of ligand in the active binding site through several chemical associations, such as hydrogen bonds and hydrophobic interactions. Our results show a favourable environment for a reaction of hydrolysis that transforms N-Carbamoyl-d-p-hydroxyphenylglycine (Cp-HPG) into the active compound D-p-hydroxyphenylglycine (D-HPG). This work envisage the use of d-NCAase from the Brazilian Bradyrhizobium japonicum strain CPAC-15 (=SEMIA 5079) for the industrial production of D-HPG, an important intermediate for semi-synthesis of β-lactam antibiotics such as penicillins, cephalosporins and amoxicillin.     

Virtual screening of eighteen million compounds against dengue virus: Combined molecular docking and molecular dynamics simulations study

Dengue virus is a major issue of tropical and sub-tropical regions. Dengue virus has been the cause behind the major alarming epidemics in the history with mass causalities from the decades. Unavailability of on-shelf drugs for the prevention of further proliferation of virus inside the human body results in immense number of deaths each year. This issue necessitates the design of novel anti-dengue drug. The protease enzyme pathway is the critical target for drug design due to its significance in the replication, survival and other cellular activities of dengue virus. Therefore, approximately eighteen million compounds from the ZINC database have been virtually screened against nonstructural protein 3 (NS3). The incremental construction algorithm of Glide docking program has been used with its features high throughput virtual screening (HTVS), standard precision (SP), extra precision (XP) and in combination of Prime module, induced fit docking (IFD) approach has also been applied. Five top-ranked compounds were then selected from the IFD results with better predicted binding energies with the catalytic triad residues (His51, Asp75, and Ser135) that may act as potential inhibitors for the underlying target protease enzyme. The top-ranked compounds ZINC95518765, ZINC44921800, ZINC71917414, ZINC39500661, ZINC36681949 have shown the predicted binding energies of −7.55, −7.36, −8.04, −8.41, −9.18 kcal/mol, respectively, forming binding interactions with three catalytically important amino acids. Top-docking poses of compounds are then used in molecular dynamics (MD) simulations. In computational studies, our proposed compounds confirm promising results against all the four serotypes of dengue virus, strengthening the opportunity of these compounds to work as potential on-shelf drugs against dengue virus. Further experimentation on the proposed compounds can result in development of strong inhibitors.     

Multicomplex-based pharmacophore modeling coupled with molecular dynamics simulations: An efficient strategy for the identification of novel inhibitors of PfDHODH

Enormous efforts have been made in the past to identify novel scaffolds against the potential therapeutic target, Plasmodium falciparum dihydroorotate dehydrogenase (PfDHODH). Fourteen different organic molecules have been crystallized to understand the structural basis of the inhibition. However, the pharmacophoric studies carried out so far, have not exploited all the structural information simultaneously to identify the novel inhibitors. Therefore, an attempt was made to construct the pharmacophore hypotheses from the available PfDHODH structural proteome. Among the generated hypotheses, a representative hypothesis was employed as a primary filter to list the molecules with complimentary features accountable for inhibition. Moreover, the auxiliary evaluations of the filtered molecules were accomplished via docking and drug-likeness studies. Subsequently, the stability of the protein-ligand complex was evaluated by using molecular dynamics simulations (MDs). The molecular details of binding interactions of the potential hits were compared with the highly active experimental structure (5FI8) to seek the more potent candidates that can be targeted against PfDHODH. Overall, the combination of screening and stability procedures resulted in the identification of three potent candidates. The drug-likeness of these molecules lie within the acceptable range and consequently increased the opportunities for their development to new anti-malarials.     

Effects of different solvents on the conformations of apoptotic cytochrome c: Structural insights from molecular dynamics simulation

Cytochrome c (cyt-c) upon binding with cardiolipin acquires peroxidase activity and is strictly connected to the pathogenesis of many human diseases including neurodegenerative and cardiovascular diseases. Interaction of cyt-c with cardiolipin mimics partial unfolding/conformational changes of cyt-c in different solvent environments. Dynamic pictures of these conformational changes of cyt-c are crucial in understanding their physiological roles in mitochondrial functions. Therefore, atomistic molecular dynamics (MD) simulations have been carried out to investigate the effect of different solvents (water, urea/water, MeOH and DMSO) on the structure and conformations of apoptotic cyt-c (Fe³⁺). Our study demonstrates that the structural changes in the protein are solvent dependent. The structural differences are observed majorly on the β-sheets and α-helical conformations and the degree of their perturbation are specific to the solvent. Although a complete loss of β-sheets (0%) is observed in MeOH and DMSO, by contrast, well preserved β-sheets (3.84%) are observed in water and urea/water. A significant decrease in the α-helical contents is observed in MeOH (41.34%) and water (42.46%), a negligible alteration in DMSO (44.25%) and well preserved α-helical (45.19%) contents in urea/water. The distances between the residues critical for electron transfer are decreased considerably for DMSO. Further, the reduction in residue flexibility and the conformational space indicate that the collective motions of cyt-c are reduced when compared to other cosolvents. Essential dynamics analysis implies that the overall motions of cyt-c in water, MeOH and urea/water are involved in three to four eigenvectors and in first eigenvector in DMSO. Overall, we believe that MD simulations of cyt-c in different solvents can provide a detailed microscopic understanding of the physiological roles, electron transport and peroxidase function in the early events of apoptosis which are hard to probe experiments.     

Exploring the molecular basis of the enantioselective binding of penicillin G acylase towards a series of 2-aryloxyalkanoic acids: a docking and molecular dynamics study

In the present paper, molecular modeling studies were undertaken in order to shed light on the molecular basis of the observed enantioselectivity of penicillin G acylase (PGA), a well known enzyme for its industrial applications, towards 16 racemic 2-aryloxyalkanoic acids, which have been reported to affect several biological systems. With this intention docking calculations and MD simulations were performed. Docking results indicated that the (S)-enantiomers establish several electrostatic interactions with SerB1, SerB386 and ArgB263 of PGA. Conversely, the absence of specific polar interactions between the (R)-enantiomers and ArgB263 seems to be the main reason for the different binding affinities observed between the two enantiomers. Results of molecular dynamics simulations demonstrated that polar interactions are responsible for both the ligand affinity and PGA enantiospecificity. Modeling calculations provided possible explanations for the observed enantioselectivity of the enzyme that rationalize available experimental data and could be the basis for future protein engineering efforts.     

Structural dynamics of wild type and mutated forms of human L1 endonuclease and insights into its sequence specific nucleic acid binding mechanism: A molecular dynamics study

Biomolecular recognition of proteins and nucleic acids is mainly mediated by their structural features and the molecular dynamics simulations approach has been used to explore this recognition processes at the atomic level. L1-Endonuclease, an enzyme involved in L1 retrotransposition, cleaves the TA junction DNA (5′-TTTT/AA-3′) and expresses high specificity for target site recognition. The present study highlights the structural features of L1-endonuclease as well as DNA responsible for such specific recognition. Especially, the importance of βB6-B5 hairpin loop in DNA recognition has been elucidated by analyzing the dynamics of Thr192 mutated L1-endonuclease. In addition, simulations of the endonuclease complexed with DNA substrates (sequences having TA and CG junctions) revealed the specificity of L1 endonuclease towards TA junction. Molecular dynamics simulations revealed that the βB6-B5 hairpin loop protrudes well into the minor groove of DNA having TA junction and induces DNA bending such that the width of minor groove is increased. Such endonuclease induced bending of TA junction DNA sequence positions the scissile phosphodiester bond of DNA for cleavage. The innate property of minor groove widening in TA junction than in CG junction is utilized by the βB6-βB5 hairpin loop of endonuclease while recognizing the DNA sequences. The present study also highlights the role of Mg²⁺ cation in catalysis and attempts to explore the possible target site DNA cleavage mechanism.