Flap flexibility amongst plasmepsins I,II, III,IV, and V: Sequence,structural, and molecular dynamics analyses

Herein, for the first time, we comparatively report the opening and closing of apo plasmepsin I – V. Plasmepsins belong the aspartic protease family of enzymes, and are expressed during the various stages of the P. falciparum lifecycle, the species responsible for the most lethal and virulent malaria to infect humans. Plasmepsin I, II, IV and HAP degrade hemoglobin from infected red blood cells, whereas plasmepsin V transport proteins crucial to the survival of the malaria parasite across the endoplasmic reticulum. Flap‐structures covering the active site of aspartic proteases (such as HIV protease) are crucial to the conformational flexibility and dynamics of the protein, and ultimately control the binding landscape. The flap‐structure in plasmepsins is made up of a flip tip in the N‐terminal lying perpendicular to the active site, adjacent to the flexible loop region in the C‐terminal. Using molecular dynamics, we propose three parameters to better describe the opening and closing of the flap‐structure in apo plasmepsins. Namely, the distance, d₁, between the flap tip and the flexible region; the dihedral angle, ?, to account for the twisting motion; and the TriCα angle, θ₁. Simulations have shown that as the flap‐structure twists, the flap and flexible region move apart opening the active site, or move toward each other closing the active site. The data from our study indicate that of all the plasmepsins investigated in the present study, Plm IV and V display the highest conformational flexibility and are more dynamic structures versus Plm I, II, and HAP. Proteins 2015; 83:1693–1705. © 2015 Wiley Periodicals, Inc.     

Molecular recognition and binding of beta‐lactamase II from Bacillus cereus with penicillin V and sulbactam by spectroscopic analysis in combination with docking simulation

The molecular recognition and binding interaction of beta‐lactamase II from Bacillus cereus (Bc II) with penicillin V (PV) and sulbactam (Sul) at 277 K were studied by spectroscopic analysis and molecular docking. The results showed that a non‐fluorescence static complex was separately formed between Bc II and two ligands, the molecular ratio of Bc II to PV or Sul was both 1:1 in the binding and the binding constants were 2.00 × 10⁶ and 3.98 × 10⁵ (L/mol), respectively. The negative free energy changes and apparent activation energies indicated that both the binding processes were spontaneous. Molecular docking showed that in the binding process, the whole Sul molecule entered into the binding pocket of Bc II while only part of the whole PV molecule entered into the pocket due to a long side chain, and electrostatic interactions were the major contribution to the binding processes. In addition, a weak conformational change of Bc II was also observed in the molecular recognition and binding process of Bc II with PV or Sul. This study may provide some valuable information for exploring the recognition and binding of proteins with ligands in the binding process and for the design of novel super‐antibiotics.     

Predicting Ligand Binding Affinity with Alchemical Free Energy Methods in a Polar Model Binding Site

We present a combined experimental and modeling study of organic ligand molecules binding to a slightly polar engineered cavity site in T4 lysozyme (L99A/M102Q). For modeling, we computed alchemical absolute binding free energies. These were blind tests performed prospectively on 13 diverse, previously untested candidate ligand molecules. We predicted that eight compounds would bind to the cavity and five would not; 11 of 13 predictions were correct at this level. The RMS error to the measurable absolute binding energies was 1.8 kcal/mol. In addition, we computed “relative” binding free energies for six phenol derivatives starting from two known ligands: phenol and catechol. The average RMS error in the relative free energy prediction was 2.5 kcal/mol (phenol) and 1.1 kcal/mol (catechol). To understand these results at atomic resolution, we obtained x-ray co-complex structures for nine of the diverse ligands and for all six phenol analogs. The average RMSD of the predicted pose to the experiment was 2.0 Å (diverse set), 1.8 Å (phenol-derived predictions), and 1.2 Å (catechol-derived predictions). We found that predicting accurate affinities and rank-orderings required near-native starting orientations of the ligand in the binding site. Unanticipated binding modes, multiple ligand binding, and protein conformational change all proved challenging for the free energy methods. We believe that these results can help guide future improvements in physics-based absolute binding free energy methods.     

Zinc finger protein binding to DNA: an energy perspective using molecular dynamics simulation and free energy calculations on mutants of both zinc finger domains and their specific DNA bases

Energy calculations based on MM-GBSA were employed to study various zinc finger protein (ZF) motifs binding to DNA. Mutants of both the DNA bound to their specific amino acids were studied. Calculated energies gave evidence for a relationship between binding energy and affinity of ZF motifs to their sites on DNA. ΔG values were ?15.82(12), ?3.66(12), and ?12.14(11.6) kcal/mol for finger one, finger two, and finger three, respectively. The mutations in the DNA bases reduced the value of the negative energies of binding (maximum value for ΔΔG = 42Kcal/mol for F1 when GCG mutated to GGG, and ΔΔG = 22 kcal/mol for F2, the loss in total energy of binding originated in the loss in electrostatic energies upon mutation (r = .98). The mutations in key amino acids in the ZF motif in positions-1, 2, 3, and 6 showed reduced binding energies to DNA with correlation coefficients between total free energy and electrostatic was .99 and with Van der Waal was .93. Results agree with experimentally found selectivity which showed that Arginine in position-1 is specific to G, while Aspartic acid (D) in position 2 plays a complicated role in binding. There is a correlation between the MD calculated free energies of binding and those obtained experimentally for prepared ZF motifs bound to triplet bases in other reports (), our results may help in the design of ZF motifs based on the established recognition codes based on energies and contributing energies to the total energy.     

Theoretical investigation on the glycan‐binding specificity of Agrocybe cylindracea galectin using molecular modeling and molecular dynamics simulation studies

Galectins are β‐galactoside binding proteins which have the ability to serve as potent antitumor, cancer biomarker, and induce tumor cell apoptosis. Agrocybe cylindracea galectin (ACG) is a fungal galectin which specifically recognizes α(2,3)‐linked sialyllactose at the cell surface that plays extensive roles in the biological recognition processes. To investigate the change in glycan‐binding specificity upon mutations, single point and double point site‐directed in silico mutations are performed at the binding pocket of ACG. Molecular dynamics (MD) simulation studies are carried out for the wild‐type (ACG) and single point (ACG1) and double point (ACG2) mutated ACGs to investigate the dynamics of substituted mutants and their interactions with the receptor sialyllactose. Plausible binding modes are proposed for galectin–sialylglycan complexes based on the analysis of hydrogen bonding interactions, total pair‐wise interaction energy between the interacting binding site residues and sialyllactose and binding free energy of the complexes using molecular mechanics–Poisson–Boltzmann surface area. Our result shows that high contribution to the binding in different modes is due to the direct and water‐mediated hydrogen bonds. The binding specificity of double point mutant Y59R/N140Q of ACG2 is found to be high, and it has 26 direct and water‐mediated hydrogen bonds with a relatively low‐binding free energy of −47.52 ± 5.2 kcal/mol. We also observe that the substituted mutant Arg59 is crucial for glycan‐binding and for the preference of α(2,3)‐linked sialyllactose at the binding pocket of ACG2 galectin. When compared with the wild‐type and single point mutant, the double point mutant exhibits enhanced affinity towards α(2,3)‐linked sialyllactose, which can be effectively used as a model for biological cell marker in cancer therapeutics. Copyright © 2015 John Wiley & Sons, Ltd.     

Exploration and validation of diphosphate‐based Plasmodium LytB inhibitors using computational approaches

LytB or IspH is an indispensable enzyme and a suitable drug target of Plasmodium falciparum that participate in isoprenoid biosynthesis of nonmevalonate pathway (MEP). Recently, we have investigated the structural dynamics of Plasmodium LytB and proposed some novel diphosphate‐based inhibitors using molecular modeling and docking studies. Here, we have tried to characterize those previously screened molecules by quantitative structure activity relationships and pharmacophore‐based analyses, as well as validated the dynamics of their interactions with LytB protein. Five total compounds having PubChem CID 516 , 125696 , 165275 , 448012 , and 9921431 were predicted with significant inhibitory activity by quantitative structure activity relationships and pharmacophore models. Again, the molecular dynamics simulation results showed that these five compounds are able to form stable complexes with the receptor through many direct and water mediated interactions. The binding free energies calculated by Poisson‐Boltzmann surface area method resulted within the range between ?99.77 and ? 43.74 kcal/mol, which favoured their profound inhibitory affinity. Residues of LytB like His41, His74, Ser222, Ser223, and Asn224 in LytB were the main protagonists in contributing the majority of interaction energies to the ligands. Finally, the ADMET, toxicity, and drug‐likeness scores also affirm these compounds to be considered for further development of new antimalarial inhibitor in the future.     

Molecular insights into Aβ42 protofibril destabilization with a fluorinated compound D744: A molecular dynamics simulation study

The aggregation of amyloid β‐peptide (Aβ₄₂) into toxic oligomers, fibrils, has been identified as a key process in Alzheimer's disease (AD) progression. The role of halogen‐substituted compounds have been highlighted in the disassembly of Aβ protofibril. However, the underlying inhibitory mechanism of Aβ₄₂ protofibril destabilization remains elusive. In this regard, a combined molecular docking and molecular dynamics (MD) simulations were performed to elucidate the inhibitory mechanism of a fluorinated compound, D744 , which has been reported previously for potential in vitro and in vivo inhibitory activity against Aβ₄₂ aggregation and reduction in the Aβ‐induced cytotoxicity. The molecular docking analysis highlights that D744 binds and interacts with chain A of the protofibril structure with hydrophobic contacts and orthogonal multipolar interaction. MD simulations reveal destabilization of the protofibril structure in the presence of D744 due to the decrease in β‐sheet content and a concomitant increase of coil and bend structures, increase in the interchain D23‐K28 salt bridge distance, decrease in the number of backbone hydrogen bonds, increase in the average distance between Cα atoms, and decrease in the binding affinity between chains A and B of the protofibril structure. The binding free‐energy analysis between D744 and the protofibril structure with Molecular Mechanics Poisson‐Boltzmann Surface Area (MM‐PBSA) reveal that residues Leu17, Val18, Phe19, Phe20, Ala21, Glu22, Asp23, Leu34, Val36, Gly37, and Gly38 of chain A of the protofibril structure contribute maximum towards binding free energy (ΔG _binding  = −44.87 kcal/mol). The insights into the underlying inhibitory mechanism of small molecules that show potential in vitro anti‐aggregation activity against Aβ₄₂ will be beneficial for the current and future AD therapeutic studies.     

Theoretical study of the ligand-CYP2B4 complexes: effect of structure on binding free energies and heme spin state

The molecular origins of temperature-dependent ligand-binding affinities and ligand-induced heme spin state conversion have been investigated using free energy analysis and DFT calculations for substrates and inhibitors of cytochrome P450 2B4 (CYP2B4), employing models of CYP2B4 based on CYP2C5(3LVdH)/CYP2C9 crystal structures, and the results compared with experiment. DFT calculations indicate that large heme-ligand interactions (ca. -15 kcal/mol) are required for inducing a high to low spin heme transition, which is correlated with large molecular electrostatic potentials (approximately -45 kcal/mol) at the ligand heteroatom. While type II ligands often contain oxygen and nitrogen heteroatoms that ligate heme iron, DFT results indicate that BP and MF heme complexes, with weak substrate-heme interactions (ca. -2 kcal/mol), and modest MEPS minima (>-35 kcal/mol) are high spin. In contrast, heme complexes of the CYP2B4 inhibitor, 4PI, the product of benzphetamine metabolism, DMBP, and water are low spin, have substantial heme-ligand interaction energies (<-15 kcal/mol) and deep MEPS minima (<-45 kcal/mol) near their heteroatoms. MMPBSA analysis of MD trajectories were made to estimate binding free energies of these ligands at the heme binding site of CYP2B4. In order to initially assess the realism of this approach, the binding free energy of 4PI inhibitor was computed and found to be a reasonable agreement with experiment: -7.7 kcal/mol [-7.2 kcal/mol (experiment)]. BP was determined to be a good substrate [-6.3 kcal/mol (with heme-ligand water), -7.3 kcal/mol (without ligand water)/-5.8 kcal/mol (experiment)], whereas the binding of MF was negligible, with only marginal binding binding free energy of -1.7 kcal/mol with 2-MF bound [-3.8 kcal/mol (experiment)], both with and without retained heme-ligand water. Analysis of the free energy components reveal that hydrophobic/nonpolar contributions account for approximately 90% of the total binding free energy of these substrates and are the source of their differential and temperature-dependent CYP2B4 binding. The results indicate the underlying origins of the experimentally observed differential binding affinities of BP and MF, and indicate the plausibility of the use of models derived from moderate sequence identity templates in conjunction with approximate free energy methods in the estimation of ligand-P450 binding affinities.     

Crystal structure and molecular modeling study of N-carbamoylsarcosine amidase Ta0454 from Thermoplasma acidophilum

A crystal structure of the putative N-carbamoylsarcosine amidase (CSHase) Ta0454 from Thermoplasma acidophilum was solved by single-wavelength anomalous diffraction and refined at a resolution of 2.35 Å. CSHases are involved in the degradation of creatinine. Ta0454 shares a similar fold and a highly conserved C-D-K catalytic triad (Cys123, Asp9, and Lys90) with the structures of three cysteine hydrolases (PDB codes 1NBA, 1IM5, and 2H0R). Molecular dynamics (MD) simulations of Ta0454/N-carbamoylsarcosine and Ta0454/pyrazinamide complexes were performed to determine the structural basis of the substrate binding pattern for each ligand. Based on the MD-simulated trajectories, the MM/PBSA method predicts binding free energies of ?24.5 and ?17.1 kcal/mol for the two systems, respectively. The predicted binding free energies suggest that Ta0454 is selective for N-carbamoylsarcosine over pyrazinamide, and zinc ions play an important role in the favorable substrate bound states.     

Analysis of cyclosporin A and a set of analogs as inhibitors of a T. cruzi cyclophilin by docking and molecular dynamics

Cyclophilins (CyPs) are enzymes involved in protein folding. In Trypanosoma cruzi (T. cruzi), the most abundantly expressed CyP is the isoform TcCyP19. It has been shown that TcCyP19 is inhibited by the immunosuppressive drug cyclosporin A (CsA) and analogs, which also proved to have potent trypanosomicidal activity in vitro. In this work, we continue and expand a previous study on the molecular interactions of CsA, and a set of analogs modeled in complexes with TcCyP19. The modeled complexes were used to evaluate binding free energies by molecular dynamics (MD), applying the Linear Interaction Energy (LIE) method. In addition, putative binding sites were identified by molecular docking. In our analysis, the binding free energy calculations did not correlate with experimental data. The heterogeneity of the non-bonded energies and the variation in the pattern of hydrogen bonds suggest that the systems may not be suitable for the application of the LIE method. Further, the docking calculations identified two other putative binding sites with comparable scoring energies to the active site, a fact that may also explain the lack of correlation found. Kinetic experiments are needed to confirm or reject the multiple binding sites hypothesis. In the meantime, MD simulations at the alternative sites, employing other methods to compute binding free energies, might be successful at finding good correlations with the experimental data.     

Molecular basis of endosomal-membrane association for the dengue virus envelope protein

Dengue virus is coated by an icosahedral shell of 90 envelope protein dimers that convert to trimers at low pH and promote fusion of its membrane with the membrane of the host endosome. We provide the first estimates for the free energy barrier and minimum for two key steps in this process: host membrane bending and protein–membrane binding. Both are studied using complementary membrane elastic, continuum electrostatics and all-atom molecular dynamics simulations. The predicted host membrane bending required to form an initial fusion stalk presents a 22–30 kcal/mol free energy barrier according to a constrained membrane elastic model. Combined continuum and molecular dynamics results predict a 15 kcal/mol free energy decrease on binding of each trimer of dengue envelope protein to a membrane with 30% anionic phosphatidylglycerol lipid. The bending cost depends on the preferred curvature of the lipids composing the host membrane leaflets, while the free energy gained for protein binding depends on the surface charge density of the host membrane. The fusion loop of the envelope protein inserts exactly at the level of the interface between the membrane's hydrophobic and head-group regions. The methods used in this work provide a means for further characterization of the structures and free energies of protein-assisted membrane fusion.     

Incorporating replacement free energy of binding‐site waters in molecular docking

Binding‐site water molecules play a crucial role in protein‐ligand recognition, either being displaced upon ligand binding or forming water bridges to stabilize the complex. However, rigorously treating explicit binding‐site waters is challenging in molecular docking, which requires to fully sample ensembles of waters and to consider the free energy cost of replacing waters. Here, we describe a method to incorporate structural and energetic properties of binding‐site waters into molecular docking. We first developed a solvent property analysis (SPA) program to compute the replacement free energies of binding‐site water molecules by post‐processing molecular dynamics trajectories obtained from ligand‐free protein structure simulation in explicit water. Next, we implemented a distance‐dependent scoring term into DOCK scoring function to take account of the water replacement free energy cost upon ligand binding. We assessed this approach in protein targets containing important binding‐site waters, and we demonstrated that our approach is reliable in reproducing the crystal binding geometries of protein‐ligand‐water complexes, as well as moderately improving the ligand docking enrichment performance. In addition, SPA program (free available to academic users upon request) may be applied in identifying hot‐spot binding‐site residues and structure‐based lead optimization. Proteins 2014; 82:1765–1776. © 2014 Wiley Periodicals, Inc.     

A dual target of Plasmepsin IX and X: Unveiling the atomistic superiority of a core chemical scaffold in malaria therapy

Plasmepsin IX and X, members of the prominent aspartic family of proteases whose function were hitherto unknown have only recently been established as key mediators of erythrocyte invasion and egress of the virulent malarial parasite. Inhibitor 49c, a potent antimalarial peptidomimetic inhibitor initially developed to target Plasmepsin II has lately been proven to exhibit potent inhibitory activity against Plasmepsin IX and X. However, the molecular and structural dynamics supporting its inhibitory activity remain inconclusive. Hindering the motion of the flap and hinge region of an aspartic protease remains essential for disabling the catalytic activity of the enzyme. Integrating molecular dynamic simulations coupled with other advanced biocomputational tools, we reveal the enhanced structural mechanistic competence of 49c in complex with Plasmepsin IX and X relative to Pepstatin. Pepstatin, a known aspartic protease inhibitor which actively hinders the opening and closing of the flap tip and flexible loop and consequently limits access to the catalytic aspartic residues, however, its administration has been related to elevated levels of toxicity. Thermodynamic calculations reveal a higher relative binding free energy associated with Plasmepsin IX and X in complex with 49c as opposed to Pepstatin. A relatively compact and structurally rigid 49c bound complexes sequel into the restriction of the flap and hinge residues by restraining cohesive movement, consequently hindering their “twisting motion” from transpiring. Findings unveil an atomistic perspective into the structural superiority of 49c in complex with Plasmepsin IX and X.     

Synthesis and multi‐spectroscopic study on DNA‐binding,cleavage and biological properties of M(II) complexes based on N2O2 donor Schiff base ligand

A novel Schiff base, (S,Z)‐4‐(methylthio)‐2‐((3‐oxo‐2,3‐dihydro‐1H‐inden‐1‐ylidene)amino)butanoic acid (L) and four M(II) complexes (where M = Co, Cu, Ni and Zn) were synthesized and characterized. The DNA‐binding characteristics of the complexes were investigated using various spectroscopic methods and viscosity measurements. Analysis of the results suggests that all the complexes bind to calf thymus DNA via intercalation. Among the four, Cu(II) complex was found to promote the photocleavage of plasmid DNA pBR322 under irradiation at 365 nm. These complexes also exhibit good antioxidant activities against 2,2‐diphenyl‐1‐picrylhydrazyl radical. In vitro antibacterial and antifungal assay indicates that these complexes are good antimicrobial agents.     

Isoform-selective inhibition of chrysin towards human cytochrome P450 1A2. Kinetics analysis,molecular docking,and molecular dynamics simulations

Our kinetics studies demonstrated that the nature product chrysin exhibited a high inhibitory affinity of 54 nM towards human cytochrome P450 1A2 and was comparable to α-naphthoflavone (49 nM), whereas it represented a moderate affinity of 5225 nM against human cytochrome P450 2C9. However, it remains unclear how this inhibitor selectively binds 1A2. To better understand the isoform selectivity of chrysin, molecular docking and molecular dynamics simulations were performed. Chrysin formed a strong H-bond with Asp313 of 1A2. The stacking interactions with Phe226 also contributed to its tight binding to 1A2. The larger and much more open active site architectures of 2C9 may explain the weaker inhibitory affinity of chrysin towards 2C9. The predicted binding free energies suggest that chrysin preferred 1A2 (ΔG_{bind, pred} = ?23.11 kcal/mol) to 2C9 (?20.41 kcal/mol). Additionally, the present work revealed that 7-hydroxy-flavone bound to 1A2 in a similar pattern as chrysin and represented a slightly less negative predicted binding free energy, which was further validated by our kinetics analysis (IC₅₀ = 240 nM). Results of the study can provide insight for designing novel isoform-selective 1A2 inhibitors.     

Interaction of polyethyleneimine‐anchored copper(II) complexes with tRNA studied by spectroscopy methods and biological activities

Ultraviolet–visible, emission and circular dichroism spectroscopic methods were used in transfer RNA (tRNA) interaction studies performed for polyethyleneimine–copper(II) complexes [Cu(phen)(l ‐Tyr)BPEI]ClO₄ (where phen =1,10‐phenanthroline, l ‐Tyr = l ‐tyrosine and BPEI = branched polyethyleneimine) with various degrees of coordination (x = 0.059, 0.149, 0.182) in the polymer chain. The results indicated that polyethyleneimine–copper(II) complexes bind with tRNA mostly through surface binding, although other binding modes, such as hydrogen bonding and van der Waals interactions, might also be present. Dye‐exclusion, sulforhodamine B and 3‐(4,5‐dimethylthiazol‐2‐yl)‐2,5‐diphenyltetrazolium bromide assays of a polyethyleneimine–copper(II) complex with a higher degree of coordination against different cancer cell lines proved that the complex exhibited cytotoxic specificity and a significant cancer cell inhibition rate. Antimicrobial screening showed activity against some human pathogens.     

Atomic Level Rendering of DNA-Drug Encounter

Computer simulations have been demonstrated to be important for unraveling atomic mechanisms in biological systems. In this study, we show how combining unbiased molecular dynamic simulations with appropriate analysis tools can successfully describe metal-based drug interactions with DNA. To elucidate the noncovalent affinity of cisplatin’s family to DNA, we performed extensive all-atom molecular dynamics simulations (3.7 μs total simulation length). The results show that the parent drug, cisplatin, has less affinity to form noncovalent adducts in the major groove than its aquo complexes. Furthermore, the relative position in which the drugs enter the major groove is dependent on the compound’s net charge. Based on the simulations, we estimated noncovalent binding free energies through the use of Markov state models. In addition, and to overcome the lack of experimental information, we employed two additional methods: Molecular Mechanics Poisson-Boltzmann Surface Area (MMPB-SA) and steered molecular dynamics with the Jarzynski estimator, with an overall good agreement between the three methods. All complexes show interaction energies below 3 kcal/mol with DNA but the charged hydrolysis products have slightly more favorable binding free energies than the parent drug. Moreover, this study sets the precedent for future unbiased DNA-ligand simulations of more complex binders.     

Nucleotide recognition by the initiation factor aIF5B: Free energy simulations of a neoclassical GTPase

The GTPase aIF5B is a universally conserved initiation factor that assists ribosome assembly. Crystal structures of its nucleotide complexes, X‐ray(GTP) and X‐ray(GDP), are similar in the nucleotide vicinity, but differ in the orientation of a distant domain IV. This has led to two, contradictory, mechanistic models. One postulates that X‐ray(GTP) and X‐ray(GDP) are, respectively, the active, “ON” and the inactive, “OFF” states; the other postulates that both structures are OFF, whereas the ON state is still uncharacterized. We study GTP/GDP binding using molecular dynamics and a continuum electrostatic free energy method. We predict that X‐ray(GTP) has a ≈ 3 kcal/mol preference to bind GDP, apparently contradicting its assignment as ON. However, the preference arises mainly from a single, nearby residue from the switch 2 motif: Glu81, which becomes protonated upon GTP binding, with a free energy cost of about 4 kcal/mol. We then propose a different model, where Glu81 protonation/deprotonation defines the ON/OFF states. With this model, the X‐ray(GTP):GTP complex, with its protonated Glu81, is ON, whereas X‐ray(GTP):GDP is OFF. The model postulates that distant conformational changes such as domain IV rotation are “uncoupled” from GTP/GDP exchange and do not affect the relative GTP/GDP binding affinities. We analyze the model using a general thermodynamic framework for GTPases. It yields rather precise predictions for the nucleotide specificities of each state, and the state specificities of each nucleotide, which are roughly comparable to the homologues IF2 and aIF2, despite the lack of any conformational switching in the model. © 2012 Wiley Periodicals, Inc.