Recognition of corn defense chitinases by fungal polyglycine hydrolases

Polyglycine hydrolases (PGH)s are secreted fungal endoproteases that cleave peptide bonds in the polyglycine interdomain linker of ChitA chitinase, an antifungal protein from domesticated corn (Zea mays ssp. mays). These target‐specific endoproteases are unusual because they do not cut a specific peptide bond but select one of many Gly‐Gly bonds within the polyglycine region. Some Gly‐Gly bonds are cleaved frequently while others are never cleaved. Moreover, we have previously shown that PGHs from different fungal pathogens prefer to cleave different Gly‐Gly peptide bonds. It is not understood how PGHs selectively cleave the ChitA linker, especially because its polyglycine structure lacks peptide sidechains. To gain insights into this process we synthesized several peptide analogs of ChitA to evaluate them as potential substrates and inhibitors of Es‐cmp, a PGH from the plant pathogenic fungus Epicoccum sorghi. Our results showed that part of the PGH recognition site for substrate chitinases is adjacent to the polyglycine linker on the carboxy side. More specifically, four amino acid residues were implicated, each spaced four residues apart on an alpha helix. Moreover, analogous peptides with selective Gly‐>sarcosine (N‐methylglycine) mutations or a specific Ser‐>Thr mutation retained inhibitor activity but were no longer cleaved by PGH. Additonally, our findings suggest that peptide analogs of ChitA that inhibit PGH activity could be used to strengthen plant defenses.     

Association of PepT1 with lipid rafts differently modulates its transport activity in polarized and nonpolarized cells

The transporter PepT1, apically expressed in intestinal epithelial cells, is responsible for the uptake of di/tripeptides. PepT1 is also expressed in nonpolarized immune cells. Here we investigated the localization of PepT1 in lipid rafts in small intestinal brush border membranes (BBMs) and polarized and nonpolarized cells, as well as functional consequences of the association of PepT1 with lipid rafts. Immunoblot analysis showed the presence of PepT1 in low-density fractions isolated from mouse intestinal BBMs, polarized intestinal Caco2-BBE cells, and nonpolarized Jurkat cells by solubilization in ice-cold 0.5% Triton X-100 and sucrose gradient fractionation. PepT1 colocalized with lipid raft markers GM1 and N-aminopeptidase in intestinal BBMs and Caco2-BBE cell membranes. Disruption of lipid rafts with methyl-beta-cyclodextrin (MbetaCD) shifted PepT1 from low- to high-density fractions. Remarkably, we found that MbetaCD treatment increased PepT1 transport activity in polarized intestinal epithelia but decreased that in intestinal BBM vesicles and nonpolarized immune cells. Mutational analysis showed that phenylalanine 293, phenylalanine 297, and threonine 281 in transmembrane segment 7 of the human di/tripeptide transporter, hPepT1, are important for the targeting to lipid rafts and transport activity of hPepT1. In conclusion, the association of PepT1 with lipid rafts differently modulates its transport activity in polarized and nonpolarized cells.     

Silencing of Mammalian Sar1 Isoforms Reveals COPII‐Independent Protein Sorting and Transport

The Sar1 GTPase coordinates the assembly of coat protein complex‐II (COPII) at specific sites of the endoplasmic reticulum (ER). COPII is required for ER‐to‐Golgi transport, as it provides a structural and functional framework to ship out protein cargoes produced in the ER. To investigate the requirement of COPII‐mediated transport in mammalian cells, we used small interfering RNA (siRNA)‐mediated depletion of Sar1A and Sar1B. We report that depletion of these two mammalian forms of Sar1 disrupts COPII assembly and the cells fail to organize transitional elements that coordinate classical ER‐to‐Golgi protein transfer. Under these conditions, minimal Golgi stacks are seen in proximity to juxtanuclear ER membranes that contain elements of the intermediate compartment, and from which these stacks coordinate biosynthetic transport of protein cargo, such as the vesicular stomatitis virus G protein and albumin. Here, transport of procollagen‐I is inhibited. These data provide proof‐of‐principle for the contribution of alternative mechanisms that support biosynthetic trafficking in mammalian cells, providing evidence of a functional boundary associated with a bypass of COPII .     

Transepithelial transport of the bioactive tripeptide, Val-Pro-Pro, in human intestinal Caco-2 cell monolayers

Some of the food-derived tripeptides with angiotensin converting enzyme (ACE)-inhibitory activity have been reported to be hypotensive after being orally administered. The mechanism for the intestinal transport of these tripeptides was studied by using monolayer-cultured human intestinal Caco-2 cells which express many enterocyte-like functions including the peptide transporter (PepT1)-mediated transport system. Val-Pro-Pro, an ACE-inhibitory peptide from fermented milk, was used as a model tripeptide. A significant amount of intact Val-Pro-Pro was transported across the Caco-2 cell monolayer. This transport was hardly inhibited by a competitive substrate for PepT1. Since no intact Val-Pro-Pro was detected in the cells, Val-Pro-Pro apically taken by Caco-2 cells via PepT1 was likely to have been quickly hydrolyzed by intracellular peptidases, producing free Val and Pro. These findings suggest that PepT1-mediated transport was not involved in the transepithelial transport of intact Val-Pro-Pro. Paracellular diffusion is suggested to have been the main mechanism for the transport of intact Val-Pro-Pro across the Caco-2 cell monolayer.     

Molecular simulation of N-acetylneuraminic acid analogs and molecular dynamics studies of cholera toxin-Neu5Gc complex

Cholera toxin (CT) is an AB₅ protein complex secreted by the pathogen Vibrio cholera, which is responsible for cholera infection. N-acetylneuraminic acid (NeuNAc) is a derivative of neuraminic acid with nine-carbon backbone. NeuNAc is distributed on the cell surface mainly located in the terminal components of glycoconjugates, and also plays an important role in cell–cell interaction. In our current study, molecular docking and molecular dynamic (MD) simulations were implemented to identify the potent NeuNAc analogs with high-inhibitory activity against CT protein. Thirty-four NeuNAc analogs, modified in different positions C-1/C-2/C-4/C-5/C-7/C-8/C-9, were modeled and docked against the active site of CT protein. Among the 34 NeuNAc analogs, the analog Neu5Gc shows the least extra precision glide score of ?9.52 and glide energy of ?44.71?kcal/mol. NeuNAc analogs block the CT active site residues HIS:13, ASN:90, LYS:91, GLN:56, GLN:61, and TRP:88 through intermolecular hydrogen bonding. The MD simulation for CT-Neu5Gc docking complex was performed using Desmond. MD simulation of CT-Neu5Gc complex reveals the stable nature of docking interaction.     

Sarcosine1,glycine8 angiotensin II is an AT1 angiotensin II receptor subtype selective antagonist

Studies predating the discovery of the two major subtypes of angiotensin II (Ang II) receptors, AT1 and AT2, revealed anomalous characteristics of sarcosine1,glycine8 Ang II (Sar1,Gly8 Ang II). It competed poorly for 125I-Ang II binding in bovine brain but potently antagonized dipsogenic responses to intracerebroventricularly administered Ang II. Subsequent recognition that bovine brain contains AT(2) receptors, while dipsogenic responses to Ang II are mediated by AT1 receptors, suggests that Sar1,Gly(8) Ang II is AT1 selective. Sar1,Gly8 Ang II competed for 125I-sarcosine1,isoleucine8 Ang II binding to AT1 receptors in pituitary, liver and adrenal (the latter with the AT2 selective antagonist PD 123,319) with Ki's of 0.66, 1.40 and 1.36 nM, respectively. In contrast, the Ki of Sar1,Gly8 Ang II for AT2 receptors in rat adrenal (with the selective AT1 antagonist losartan) was 52 nM. 125I-Sar1,Gly8 Ang II (0.5-3 nM) bound to AT1 receptors in pituitary, liver, heart, adrenal, and hypothalamic membranes with high affinity (Kd=0.43, 1.6, 2.3, 0.96 and 1.8 nM, respectively), but showed no saturable binding to the adrenal AT2 receptor. 125I-Sar1,Gly8 Ang II selectively labeled AT1 receptors in sections of adrenal using receptor autoradiography. Thus, binding studies reveal Sar1,Gly8 Ang II to be the first angiotensin peptide analog to show AT1 receptor selectivity. 125I-Sar1,Gly8 Ang II offers a new means to selectively radiolabel AT1 receptors and may help to characterize ligand docking sites and agonist switches for AT1 versus AT2 receptors.     

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.     

The transmembrane tyrosines Y56, Y91 and Y167 play important roles in determining the affinity and transport rate of the rabbit proton-coupled peptide transporter PepT1

The mammalian proton-coupled peptide transporter PepT1 is widely accepted as the major route of uptake for dietary nitrogen, as well as being responsible for the oral absorption of a number of classes of drugs, including β-lactam antibiotics and angiotensin-converting enzyme (ACE) inhibitors. Using site-directed mutagenesis and zero-trans transport assays, we investigated the role of conserved tyrosines in the transmembrane domains (TMDs) of rabbit PepT1 as predicted by hydropathy plots.All the individual TMD tyrosines were substituted with phenylalanine and shown to retain the ability to traffic to the plasma membrane of Xenopus laevis oocytes. These single substitutions of TMD tyrosines by phenylalanine residues did not affect the proton dependence of peptide uptake, with all retaining wild-type PepT1-like pH dependence. Individual mutations of four of the nine TMD residue tyrosines (Y64, Y287, Y345 and Y587) were without measurable effect on PepT1 function, whereas the other five (Y12, Y56, Y91, Y167 and Y345) were shown to result in altered transport function compared to the wild-type PepT1.Intriguingly, the affinity of Y56F-PepT1 was found to be dramatically increased (approximately 100-fold) in comparison to that of the wild-type rabbit PepT1. Y91 mutations also affected the substrate affinity of the transporter, which increased in line with the hydrophilicity of the substituted amino acid (F > Y > Q > R). Y167 was demonstrated to play a pivotal role in rabbit PepT1 function since Y167F, Y167R and Y167Q demonstrated very little transport function. These results are discussed with regard to a proposed mechanism for PepT1 substrate binding.     

Transport of Charged Dipeptides by the Intestinal H+/Peptide Symporter PepT1 Expressed in Xenopus laevis Oocytes

The cloned intestinal peptide transporter is capable of electrogenic H⁺-coupled cotransport of neutral di- and tripeptides and selected peptide mimetics. Since the mechanism by which PepT1 transports substrates that carry a net negative or positive charge at neutral pH is poorly understood, we determined in Xenopus oocytes expressing PepT1 the characteristics of transport of differently charged glycylpeptides. Transport function of PepT1 was assessed by flux studies employing a radiolabeled dipeptide and by the two-electrode voltage-clamp-technique. Our studies show, that the transporter is capable of translocating all substrates by an electrogenic process that follows Michaelis Menten kinetics. Whereas the apparent K_0.5 value of a zwitterionic substrate is only moderately affected by alterations in pH or membrane potential, K_0.5 values of charged substrates are strongly dependent on both, pH and membrane potential. Whereas the affinity of the anionic dipeptide increased dramatically by lowering the pH, a cationic substrate shows only a weak affinity for PepT1 at all pH values (5.5–8.0). The driving force for uptake is provided mainly by the inside negative transmembrane electrical potential. In addition, affinity for proton interaction with PepT1 was found to depend on membrane potential and proton binding subsequently affects the substrate affinity. Furthermore, our studies suggest, that uptake of the zwitterionic form of a charged substrate contributes to overall transport and that consequently the stoichiometry of the flux-coupling ratios for peptide: H⁺/H₃O⁺ cotransport may vary depending on pH. Received: 19 August 1996/Revised: 10 October 1996     

Crystal structure analysis,covalent docking,and molecular dynamics calculations reveal a conformational switch in PhaZ7 PHB depolymerase

An open and a closed conformation of a surface loop in PhaZ7 extracellular poly(3‐hydroxybutyrate) depolymerase were identified in two high‐resolution crystal structures of a PhaZ7 Y105E mutant. Molecular dynamics (MD) simulations revealed high root mean square fluctuations (RMSF) of the 281–295 loop, in particular at residue Asp289 (RMSF 7.62 Å). Covalent docking between a 3‐hydroxybutyric acid trimer and the catalytic residue Ser136 showed that the binding energy of the substrate is significantly more favorable in the open loop conformation compared to that in the closed loop conformation. MD simulations with the substrate covalently bound depicted 1 Å RMSF higher values for the residues 281–295 in comparison to the apo (substrate‐free) form. In addition, the presence of the substrate in the active site enhanced the ability of the loop to adopt a closed form. Taken together, the analysis suggests that the flexible loop 281–295 of PhaZ7 depolymerase can act as a lid domain to control substrate access to the active site of the enzyme. Proteins 2017; 85:1351–1361. © 2017 Wiley Periodicals, Inc.     

Novel urushiol derivatives as HDAC8 inhibitors: rational design,virtual screening,molecular docking and molecular dynamics studies

Three series of novel urushiol derivatives were designed by introducing a hydroxamic acid moiety into the tail of an alkyl side chain and substituents with differing electronic properties or steric bulk onto the benzene ring and alkyl side chain. The compounds’ binding affinity toward HDAC8 was screened by Glide docking. The highest-scoring compounds were processed further with molecular docking, MD simulations, and binding free energy studies to analyze the binding modes and mechanisms. Ten compounds had Glide scores of ?8.2 to ?10.2, which revealed that introducing hydroxy, carbonyl, amino, or methyl ether groups into the alkyl side chain or addition of –F, –Cl, sulfonamide, benzamido, amino, or hydroxy substituents on the benzene ring could significantly increase binding affinity. Molecular docking studies revealed that zinc ion coordination, hydrogen bonding, and hydrophobic interactions contributed to the high calculated binding affinities of these compounds toward HDAC8. MD simulations and binding free energy studies showed that all complexes possessed good stability, as characterized by low RMSDs, low RMSFs of residues, moderate hydrogen bonding and zinc ion coordination and low values of binding free energies. Hie147, Tyr121, Phe175, Hip110, Phe119, Tyr273, Lys21, Gly118, Gln230, Leu122, Gly269, and Gly107 contributed favorably to the binding; and Van der Waals and electrostatic interactions provided major contributions to the stability of these complexes. These results show the potential of urushiol derivatives as HDAC8 binding lead compounds, which have great therapeutic potential in the treatment of various malignancies, neurological disorders, and human parasitic diseases.     

Holo‐BtuF stabilizes the open conformation of the vitamin B12 ABC transporter BtuCD

While there is evidence that other ABC transporters can tell between empty and loaded substrate binding protein, reconstitution experiments suggest otherwise for the Escherichia coli vitamin B12 importer BtuCD‐F. Here, we address the question of BtuCD‐F substrate sensitivity in a combined protein–protein docking and molecular dynamics simulation approach. Starting from the BtuCD and holo‐BtuF crystal structures, we model two holo‐BtuCD‐F docking complexes differing by a 180° orientation of BtuF. One of these is similar to the apo‐BtuCD‐F crystal structure. Both docking complexes were embedded in a lipid/water environment to investigate their dynamics and BtuCD's conformational response to the presence and absence of BtuF, vitamin B12, and Mg‐ATP in a series of 28 independent MD simulations. We find holo‐BtuF stabilizing the open conformation of BtuCD, whereas the transporter begins to close again when BtuF or vitamin B12 is removed—suggesting BtuCD‐F is capable of substrate sensitivity. We identified BtuC transmembrane helices 3 and 5, the L‐loops and the adjacent helices comprised of BtuC residues 170–180 as hotspots of conformational change. We propose the latter to act as substrate sensors. BtuF‐Trp44 appears to act as a lid on the vitamin B12 binding cleft in BtuF X‐ray structures and protrudes into the BtuCD transport channel in one of our simulations, which might represent an initial step in vitamin B12 uptake. On an average, we observe subunit motions where the nucleotide binding domains approach each other while the transmembrane domains display an opening trend toward the periplasm. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Comparative study of stability and transport of molecules through cyclic peptide nanotube and aquaporin: a molecular dynamics simulation approach

Abstract

The structural stability and transport properties of the cyclic peptide nanotube (CPN) 8?×?[Cys–Gly–Met–Gly]₂ in different phospholipid bilayers such as POPA (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidic acid), POPE (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine), POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine), POPG (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol) and POPS (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoserine) with water have been investigated using molecular dynamics (MD) simulation. The hydrogen bonds and non-bonded interaction energies were calculated to study the stability in different bilayers. One µs MD simulation in POPA lipid membrane reveals the stability of the cyclic peptide nanotube, and the simulations at various temperatures manifest the higher stability of 8?×?[Cys–Gly–Met–Gly]₂. We demonstrated that the presence of sulphur-containing amino acids in CPN enhances the stability through disulphide bonds between the adjacent rings. Further, the water permeation coefficient of the CPN is calculated and compared with human aquaporin-2 (AQP2) channel protein. It is found that the coefficients are highly comparable to the AQP2 channel though the mechanism of water transport is not similar to AQP 2; the flow of water in the CPN is taking place as a two-line 1–2–1–2 file fashion. In addition to that, the transport behavior of Na⁺ and K⁺ ions, single water molecule, urea and anti-cancer drug fluorouracil were investigated using pulling simulation and potential of mean force calculation. The above transport behavior shows that Na⁺ is trapped in CPN for a longer time than other molecules. Also, the interactions of the ions and molecules in Cα and mid-Cα plane were studied to understand the transport behavior of the CPN. Abbreviations AQP2 Aquaporin-2

CPN Cyclic peptide nanotube

MD Molecular dynamics

POPA 1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphatidic acid

POPE 1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine

POPG 1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol

POPS 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoserine

Communicated by Ramaswamy H. Sarma     

New insight into the architecture of oxy‐anion pocket in unliganded conformation of GAT domains: A MD‐simulation study

Human Guanine Monophosphate Synthetase (hGMPS) converts XMP to GMP, and acts as a bifunctional enzyme with N‐terminal “glutaminase” (GAT) and C‐terminal “synthetase” domain. The enzyme is identified as a potential target for anti‐cancer and immunosuppressive therapies. GAT domain of enzyme plays central role in metabolism, and contains conserved catalytic residues Cys104, His190, and Glu192. MD simulation studies on GAT domain suggest that position of oxyanion in unliganded conformation is occupied by one conserved water molecule (W1), which also stabilizes that pocket. This position is occupied by a negatively charged atom of the substrate or ligand in ligand bound crystal structures. In fact, MD simulation study of Ser75 to Val indicates that W1 conserved water molecule is stabilized by Ser75, while Thr152, and His190 also act as anchor residues to maintain appropriate architecture of oxyanion pocket through water mediated H‐bond interactions. Possibly, four conserved water molecules stabilize oxyanion hole in unliganded state, but they vacate these positions when the enzyme (hGMPS)‐substrate complex is formed. Thus this study not only reveals functionally important role of conserved water molecules in GAT domain, but also highlights essential role of other non‐catalytic residues such as Ser75 and Thr152 in this enzymatic domain. The results from this computational study could be of interest to experimental community and provide a testable hypothesis for experimental validation. Conserved sites of water molecules near and at oxyanion hole highlight structural importance of water molecules and suggest a rethink of the conventional definition of chemical geometry of inhibitor binding site. Proteins 2016; 84:360–373. © 2016 Wiley Periodicals, Inc.     

An explorative study on Staphylococcus aureus MurE inhibitor: induced fit docking,binding free energy calculation,and molecular dynamics

Staphylococcus aureus MurE enzyme catalyzes the addition of l-lysine as third residue of the peptidoglycan peptide moiety. Due to the high substrate specificity and its ubiquitous nature among bacteria, MurE enzyme is considered as one of the potential target for the development of new therapeutic agents. In the present work, induced fit docking (IFD), binding free energy calculation, and molecular dynamics (MD) simulation were carried out to elucidate the inhibition potential of 2-thioxothiazolidin-4-one based inhibitor 1 against S. aureus MurE enzyme. The inhibitor 1 formed majority of hydrogen bonds with the central domain residues Asn151, Thr152, Ser180, Arg187, and Lys219. Binding free-energy calculation by MM-GBSA approach showed that van der Waals (ΔG_vdW, ?57.30?kcal/mol) and electrostatic solvation (ΔG_solv, ?36.86?kcal/mol) energy terms are major contributors for the inhibitor binding. Further, 30-ns MD simulation was performed to validate the stability of ligand–protein complex and also to get structural insight into mode of binding. Based on the IFD and MD simulation results, we designed four new compounds D1–D4 with promising binding affinity for the S. aureus MurE enzyme. The designed compounds were subjected to the extra-precision docking and binding free energy was calculated for complexes. Further, a 30-ns MD simulation was performed for D1/4C13 complex.     

Structure‐guided screening of chemical database to identify NS3‐NS2B inhibitors for effective therapeutic application in dengue infection

Dengue infection is the most common arthropod‐borne disease caused by dengue viruses, predominantly affecting millions of human beings annually. To find out promising chemical entities for therapeutic application in Dengue, in the current research, a multi‐step virtual screening effort was conceived to screen out the entire “screening library” of the Asinex database. Initially, through “Lipinski rule of five” filtration criterion almost 0.6 million compounds were collected and docked with NS3‐NS2B protein. Thereby, the chemical space was reduced to about 3500 compounds through the analysis of binding affinity obtained from molecular docking study in AutoDock Vina. Further, the “Virtual Screening Workflow” (VSW) utility of Schrödinger suite was used, which follows a stepwise multiple docking programs such as ‐ high‐throughput virtual screening (HTVS), standard precision (SP), and extra precision (XP) docking, and in postprocessing analysis the MM‐GBSA based free binding energy calculation. Finally, five potent molecules were proposed as potential inhibitors for the dengue NS3‐NS2B protein based on the investigation of molecular interactions map and protein‐ligand fingerprint analyses. Different pharmacokinetics and drug‐likeness parameters were also checked, which favour the potentiality of selected molecules for being drug‐like candidates. The molecular dynamics (MD) simulation analyses of protein‐ligand complexes were explained that NS3‐NS2B bound with proposed molecules quite stable in dynamic states as observed from the root means square deviation (RMSD) and root means square fluctuation (RMSF) parameters. The binding free energy was calculated using MM‐GBSA method from the MD simulation trajectories revealed that all proposed molecules possess such a strong binding affinity towards the dengue NS3‐NS2B protein. Therefore, proposed molecules may be potential chemical components for effective inhibition of dengue NS3‐NS2B protein subjected to experimental validation.     

Insights into regioselective metabolism of mefenamic acid by cytochrome P450 BM3 mutants through crystallography,docking, molecular dynamics,and free energy calculations

Cytochrome P450 BM3 (CYP102A1) mutant M11 is able to metabolize a wide range of drugs and drug‐like compounds. Among these, M11 was recently found to be able to catalyze formation of human metabolites of mefenamic acid and other nonsteroidal anti‐inflammatory drugs (NSAIDs). Interestingly, single active‐site mutations such as V87I were reported to invert regioselectivity in NSAID hydroxylation. In this work, we combine crystallography and molecular simulation to study the effect of single mutations on binding and regioselective metabolism of mefenamic acid by M11 mutants. The heme domain of the protein mutant M11 was expressed, purified, and crystallized, and its X‐ray structure was used as template for modeling. A multistep approach was used that combines molecular docking, molecular dynamics (MD) simulation, and binding free‐energy calculations to address protein flexibility. In this way, preferred binding modes that are consistent with oxidation at the experimentally observed sites of metabolism (SOMs) were identified. Whereas docking could not be used to retrospectively predict experimental trends in regioselectivity, we were able to rank binding modes in line with the preferred SOMs of mefenamic acid by M11 and its mutants by including protein flexibility and dynamics in free‐energy computation. In addition, we could obtain structural insights into the change in regioselectivity of mefenamic acid hydroxylation due to single active‐site mutations. Our findings confirm that use of MD and binding free‐energy calculation is useful for studying biocatalysis in those cases in which enzyme binding is a critical event in determining the selective metabolism of a substrate. Proteins 2016; 84:383–396. © 2016 Wiley Periodicals, Inc.     

Diabetes mellitus caused by mutations in human insulin: analysis of impaired receptor binding of insulins Wakayama,Los Angeles and Chicago using pharmacoinformatics

Several naturally occuring mutations in the human insulin gene are associated with diabetes mellitus. The three known mutant molecules, Wakayama, Los Angeles and Chicago were evaluated using molecular docking and molecular dynamics (MD) to analyse mechanisms of deprived binding affinity for insulin receptor (IR). Insulin Wakayama, is a variant in which valine at position A3 is substituted by leucine, while in insulin Los Angeles and Chicago, phenylalanine at positions B24 and B25 is replaced by serine and leucine, respectively. These mutations show radical changes in binding affinity for IR. The ZDOCK server was used for molecular docking, while AMBER 14 was used for the MD study. The published crystal structure of IR bound to natural insulin was also used for MD. The binding interactions and MD trajectories clearly explained the critical factors for deprived binding to the IR. The surface area around position A3 was increased when valine was substituted by leucine, while at positions B24 and B25 aromatic amino acid phenylalanine replaced by non-aromatic serine and leucine might be responsible for fewer binding interactions at the binding site of IR that leads to instability of the complex. In the MD simulation, the normal mode analysis, rmsd trajectories and prediction of fluctuation indicated instability of complexes with mutant insulin in order of insulin native insulin < insulin Chicago < insulin Los Angeles < insulin Wakayama molecules which corresponds to the biological evidence of the differing affinities of the mutant insulins for the IR.