The quaternary structure of the Saccharomyces cerevisiae succinate dehydrogenase. Homology modeling, cofactor docking, and molecular dynamics simulation studies

Succinate dehydrogenases and fumarate reductases are complex mitochondrial or bacterial respiratory chain proteins with remarkably similar structures and functions. Succinate dehydrogenase oxidizes succinate and reduces ubiquinone using a flavin adenine dinucleotide cofactor and iron-sulfur clusters to transport electrons. A model of the quaternary structure of the tetrameric Saccharomyces cerevisiae succinate dehydrogenase was constructed based on the crystal structures of the Escherichia coli succinate dehydrogenase, the E. coli fumarate reductase, and the Wolinella succinogenes fumarate reductase. One FAD and three iron-sulfur clusters were docked into the Sdh1p and Sdh2p catalytic dimer. One b-type heme and two ubiquinone or inhibitor analog molecules were docked into the Sdh3p and Sdh4p membrane dimer. The model is consistent with numerous experimental observations. The calculated free energies of inhibitor binding are in excellent agreement with the experimentally determined inhibitory constants. Functionally important residues identified by mutagenesis of the SDH3 and SDH4 genes are located near the two proposed quinone-binding sites, which are separated by the heme. The proximal quinone-binding site, located nearest the catalytic dimer, has a considerably more polar environment than the distal site. Alternative low energy conformations of the membrane subunits were explored in a molecular dynamics simulation of the dimer embedded in a phospholipid bilayer. The simulation offers insight into why Sdh4p Cys-78 may be serving as the second axial ligand for the heme instead of a histidine residue. We discuss the possible roles of heme and of the two quinone-binding sites in electron transport.     

Molecular modeling studies of atorvastatin analogues as HMGR inhibitors using 3D-QSAR,molecular docking and molecular dynamics simulations

3-Hydroxy-3-methylglutaryl coenzyme-A reductase (HMGR) is generally regarded as targets for the treatment of hypercholesterolemia. HMGR inhibitors (more commonly known as statins) are discovered as plasma cholesterol lowering molecules. In this work, 120 atorvastatin analogues were studied using a combination of molecular modeling techniques including three-dimensional quantitative structure–activity relationship (3D-QSAR), molecular docking and molecular dynamics (MD) simulation. The results show that the best CoMFA (comparative molecular field analysis) model has q² = 0.558 and r² = 0.977, and the best CoMSIA (comparative molecular similarity indices analysis) model has q² = 0.582 and r² = 0.919. Molecular docking and MD simulation explored the binding relationship of the ligand and the receptor protein. The calculation results indicated that the hydrophobic and electrostatic fields play key roles in QSAR model. After MD simulation, we found four vital residues (Lys735, Arg590, Asp690 and Asn686) and three hydrophobic regions in HMGR binding site. The calculation results show that atorvastatin analogues obtained by introduction of F atoms or gem-difluoro groups could obviously improve the inhibitory activity. The new HMGR inhibitor analogues design in this Letter had been submitted which is being currently synthesized by our laboratories.     

Deciphering the GPER/GPR30-agonist and antagonists interactions using molecular modeling studies,molecular dynamics,and docking simulations

The G-protein coupled estrogen receptor 1 GPER/GPR30 is a transmembrane seven-helix (7TM) receptor involved in the growth and proliferation of breast cancer. Due to the absence of a crystal structure of GPER/GPR30, in this work, molecular modeling studies have been carried out to build a three-dimensional structure, which was subsequently refined by molecular dynamics (MD) simulations (up to 120 ns). Furthermore, we explored GPER/GPR30’s molecular recognition properties by using reported agonist ligands (G1, estradiol (E2), tamoxifen, and fulvestrant) and the antagonist ligands (G15 and G36) in subsequent docking studies. Our results identified the E2 binding site on GPER/GPR30, as well as other receptor cavities for accepting large volume ligands, through GPER/GPR30 π–π, hydrophobic, and hydrogen bond interactions. Snapshots of the MD trajectory at 14 and 70 ns showed almost identical binding motifs for G1 and G15. It was also observed that C107 interacts with the acetyl oxygen of G1 (at 14 ns) and that at 70 ns the residue E275 interacts with the acetyl group and with the oxygen from the other agonist whereas the isopropyl group of G36 is oriented toward Met141, suggesting that both C107 and E275 could be involved in the protein activation. This contribution suggest that GPER1 has great structural changes which explain its great capacity to accept diverse ligands, and also, the same ligand could be recognized in different binding pose according to GPER structural conformations.     

Homology modeling,molecular dynamic simulation and docking studies of cyclin dependent kinase 1

In order to develop promising cyclin dependent kinase 1 inhibitors, homology modeling, docking and molecular dynamic simulation techniques were applied to get insight into the functional and structural properties of cyclin dependent kinase 1 (CDK1). Since there is no reported CDK1 crystal structural data, the three dimensional structure of CDK1 was constructed based on homology modeling. An extensive dynamic simulation was also performed on a Flavopiridol-CDK1 complex for probing the binding pattern of Flavopiridol in the active site of CDK1. The binding modes of other inhibitors to CDK1 were also proposed by molecular docking. The structural requirement for developing more potent CDK1 inhibitors was obtained by the above-mentioned molecular simulations and pharmacophore modeling.     

Homology modeling and molecular dynamics simulations of transmembrane domain structure of human neuronal nicotinic acetylcholine receptor

A three-dimensional model of the transmembrane domain of a neuronal-type nicotinic acetylcholine receptor (nAChR), (alpha4)2(beta2)3, was constructed from a homology structure of the muscle-type nAChR recently determined by cryo-electron microscopy. The neuronal channel model was embedded in a fully hydrated DMPC lipid bilayer, and molecular-dynamics simulations were performed for 5 ns. A comparative analysis of the neuronal- versus muscle-type nAChR models revealed many conserved pore-lining residues, but an important difference was found near the periplasmic mouth of the pore. A flickering salt-bridge of alpha4-E266 with its adjacent beta2-K260 was observed in the neuronal-type channel during the course of the molecular-dynamics simulations. The narrowest region, with a pore radius of approximately 2 A formed by the salt-bridges, does not seem to be the restriction site for a continuous water passage. Instead, two hydrophobic rings, formed by alpha4-V259, alpha4-L263, and the homologous residues in the beta2-subunits, act as the gates for water flow, even though the region has a slightly larger pore radius. The model offers new insight into the water transport across the (alpha4)2(beta2)3 nAChR channel, and may lead to a better understanding of the structures, dynamics, and functions of this family of ion channels.     

Homology modelling and molecular dynamics simulations: comparative studies of human aquaporin-1

The structures of the mammalian water transport protein Aqp1 and of its bacterial homologue GlpF enables us to test whether homology models can be used to explore relationships between structure, dynamics and function in mammalian transport proteins. Molecular dynamics simulations (totalling almost 40 ns) were performed starting from: the X-ray structure of Aqp1; a homology model of Aqp1 based on the GlpF structure; and intermediate resolution structures of Aqp1 derived from electron microscopy. Comparisons of protein RMSDs vs. time suggest that the homology models are of comparable conformational stability to the X-ray structure, whereas the intermediate resolution structures exhibit significant conformation drift. For simulations based on the X-ray structure and on homology models, the flexibility profile vs. residue number correlates well with the crystallographic B-values for each residue. In the simulations based on intermediate resolution structures, mobility of the highly conserved NPA loops is substantially higher than in the simulations based on the X-ray structure or the homology models. Pore radius profiles remained relatively constant in the X-ray and homology model simulations but showed substantial fluctuations (reflecting the higher NPA loop mobility) in the intermediate resolution simulations. The orientation of the dipoles of water molecules within the pore is of key importance in maintaining low proton permeability through Aqp1. This property seems to be quite robust to the starting model used in the simulation. These simulations suggest that homology models based on bacterial homologues may be used to derive functionally relevant information on the structural dynamics of mammalian transport proteins.     

Homology modeling of dopamine d(2) and d(3) receptors: molecular dynamics refinement and docking evaluation

Dopamine (DA) receptors, a class of G-protein coupled receptors (GPCRs), have been targeted for drug development for the treatment of neurological, psychiatric and ocular disorders. The lack of structural information about GPCRs and their ligand complexes has prompted the development of homology models of these proteins aimed at structure-based drug design. Crystal structure of human dopamine D(3) (hD(3)) receptor has been recently solved. Based on the hD(3) receptor crystal structure we generated dopamine D(2) and D(3) receptor models and refined them with molecular dynamics (MD) protocol. Refined structures, obtained from the MD simulations in membrane environment, were subsequently used in molecular docking studies in order to investigate potential sites of interaction. The structure of hD(3) and hD(2L) receptors was differentiated by means of MD simulations and D(3) selective ligands were discriminated, in terms of binding energy, by docking calculation. Robust correlation of computed and experimental K(i) was obtained for hD(3) and hD(2L) receptor ligands. In conclusion, the present computational approach seems suitable to build and refine structure models of homologous dopamine receptors that may be of value for structure-based drug discovery of selective dopaminergic ligands.     

Homology modeling,substrate docking,and molecular simulation studies of mycobacteriophage Che12 lysin A

Mycobacteriophages produce lysins that break down the host cell wall at the end of lytic cycle to release their progenies. The ability to lyse mycobacterial cells makes the lysins significant. Mycobacteriophage Che12 is the first reported temperate phage capable of infecting and lysogenising Mycobacterium tuberculosis. Gp11 of Che12 was found to have Chitinase domain that serves as endolysin (lysin A) for Che12. Structure of gp11 was modeled and evaluated using Ramachandran plot in which 98 % of the residues are in the favored and allowed regions. Che12 lysin A was predicted to act on NAG-NAM-NAG molecules in the peptidoglycan of cell wall. The tautomers of NAG-NAM-NAG molecule were generated and docked with lysin A. The stability and binding affinity of lysin A – NAG-NAM-NAG tautomers were studied using molecular dynamics simulations.     

Homology modeling, molecular dynamics, e-pharmacophore mapping and docking study of Chikungunya virus nsP2 protease

To date, no suitable vaccine or specific antiviral drug is available to treat Chikungunya viral (CHIKV) fever. Hence, it is essential to identify drug candidates that could potentially impede CHIKV infection. Here, we present the development of a homology model of nsP2 protein based on the crystal structure of the nsP2 protein of Venezuelan equine encephalitis virus (VEEV). The protein modeled was optimized using molecular dynamics simulation; the junction peptides of a nonstructural protein complex were then docked in order to investigate the possible protein–protein interactions between nsP2 and the proteins cleaved by nsP2. The modeling studies conducted shed light on the binding modes, and the critical interactions with the peptides provide insight into the chemical features needed to inhibit the CHIK virus infection. Energy-optimized pharmacophore mapping was performed using the junction peptides. Based on the results, we propose the pharmacophore features that must be present in an inhibitor of nsP2 protease. The resulting pharmacophore model contained an aromatic ring, a hydrophobic and three hydrogen-bond donor sites. Using these pharmacophore features, we screened a large public library of compounds (Asinex, Maybridge, TOSLab, Binding Database) to find a potential ligand that could inhibit the nsP2 protein. The compounds that yielded a fitness score of more than 1.0 were further subjected to Glide HTVS and Glide XP. Here, we report the best four compounds based on their docking scores; these compounds have IDs of 27943, 21362, ASN 01107557 and ASN 01541696. We propose that these compounds could bind to the active site of nsP2 protease and inhibit this enzyme. Furthermore, the backbone structural scaffolds of these four lead compounds could serve as building blocks when designing drug-like molecules for the treatment of Chikungunya viral fever.     

An approach to prevent aggregation during the purification and crystallization of wild type acyl coenzyme A: isopenicillin N acyltransferase from Penicillium chrysogenum

Acyl coenzyme A: isopenicillin N acyltransferase (AT) from Penicillium chrysogenum is an enzyme of interest for the biosynthesis of beta-lactam antibiotics. Severe aggregation problems with wild type AT have, however, prevented significant progress in the structure-function analysis of this enzyme for a decade. In this study, we show an approach to solve this aggregation problem by using dynamic light scattering (DLS) analysis to probe the aggregation state of the protein in the presence of various additives. After a one-step purification of recombinant wild type AT with a C-terminal His-tag using Ni2+ affinity chelate chromatography, addition of a combination of 5 mM DTT, 250 mM NaCl, and 5 mM EDTA to the purified AT effectively prevented aggregation. In the presence of these additives, the DLS profile of AT shows a narrow size distribution indicative of a homogeneous protein solution and the absence of aggregation. The purity and mono-dispersity of wild type AT was sufficient for the growth of high quality crystals diffracting to 1.64 A resolution.     

Insights into the structure and inhibition of Giardia intestinalis arginine deiminase: homology modeling,docking, and molecular dynamics studies

Giardia intestinalis arginine deiminase (GiADI) is an important metabolic enzyme involved in the energy production and defense of this protozoan parasite. The lack of this enzyme in the human host makes GiADI an attractive target for drug design against G. intestinalis. One approach in the design of inhibitors of GiADI could be computer-assisted studies of its crystal structure, such as docking; however, the required crystallographic structure of the enzyme still remains unresolved. Because of its relevance, in this work, we present a three-dimensional structure of GiADI obtained from its amino acid sequence using the homology modeling approximation. Furthermore, we present an approximation of the most stable dimeric structure of GiADI identified through molecular dynamics simulation studies. An in silico analysis of druggability using the structure of GiADI was carried out in order to know if it is a good target for design and optimization of selective inhibitors. Potential GiADI inhibitors were identified by docking of a set of 3196 commercial and 19 in-house benzimidazole derivatives, and molecular dynamics simulation studies were used to evaluate the stability of the ligand–enzyme complexes.     

Anesthetic interaction with ketosteroid isomerase: insights from molecular dynamics simulations

The nature and the sites of interactions between anesthetic halothane and homodimeric Delta5-3-ketosteroid isomerase (KSI) are characterized by flexible ligand docking and confirmed by 1H-15N NMR. The dynamics consequence of halothane interaction and the implication of the dynamic changes to KSI function are studied by multiple 5-ns molecular dynamics simulations in the presence and absence of halothane. Both docking and MD simulations show that halothane prefer the amphiphilic dimeric interface to the hydrophobic active site of KSI. Halothane occupancy at the dimer interface disrupted the intersubunit hydrogen bonding formed either directly through side chains of polar residues or indirectly through the mediation of the interfacial water molecules. Moreover, in the presence of halothane, the exchange rate of the bound waters with bulk water was increased. Halothane perturbation to the dimer interface affected the overall flexibility of the active site. This action is likely to contribute to the halothane-induced reduction of the KSI activity. The allosteric halothane modulation of the dynamics-function relationship of KSI without direct competition at the enzymatic active sites may be generalized to offer a unifying explanation of anesthetic action on a diverse range of multidomain neuronal proteins that are potentially relevant to clinical general anesthesia.     

Homology modeling, molecular docking and electrostatic potential analysis of MurF ligase from Klebsiella pneumonia

In spite of availability of moderately protective vaccine and antibiotics, new antibacterial agents are urgently needed to decrease the global incidence of Klebsiella pneumonia infections. MurF ligase, a key enzyme, which participates in the bacterial cell wall assembly, is indispensable to existence of K. pneumonia. MurF ligase lack mammalian vis-à-vis and have high specificity, uniqueness, and occurrence only in eubacteria, epitomizing them as promising therapeutic targets for intervention. In this study, we present a unified approach involving homology modeling and molecular docking studies on MurF ligase enzyme. As part of this study, a homology model of K. pneumonia (MurF ligase) enzyme was predicted for the first time in order to carry out structurebased drug design. The accuracy of the model was further validated using different computational approaches. The comparative molecular docking study on this enzyme was undertaken using different phyto-ligands from Desmodium sp. and a known antibiotic Ciprofloxacin. The docking analysis indicated the importance of hotspots (HIS 281 and ASN 282) within the MurF binding pocket. The Lipinski's rule of five was analyzed for all ligands considered for this study by calculating the ADME/Tox, drug likeliness using Qikprop simulation. Only ten ligands were found to comply with the Lipinski rule of five. Based on the molecular docking results and Lipinki values 6-Methyltetrapterol A was confirmed as a promising lead compound. The present study should therefore play a guiding role in the experimental design and development of 6-Methyltetrapterol A as a bactericidal agent.     

Homology modelling of human DHCR24 (seladin-1) and analysis of its binding properties through molecular docking and dynamics simulations

Recent biochemical and clinical evidences unveiled that DHCR24 enzyme (3-beta-hydoxysterol-Delta(24)-reductase, also named seladin-1), which catalyzes the last step of the cholesterol biosynthesis, is implicated in relevant neuroprotective processes by modulating the level of cholesterol in membrane. The present study was undertaken with a view to model the DHCR24 enzyme and its catalytic site, analyzing the substrate recognition at an atomic level. A homology model of the enzyme was obtained based on plant Cytokinin Dehydrogenase, and its active site was found to bind the desmosterol plus a set of post-squalenic intermediates of the cholesterol biosynthesis in a binding mode conducive to catalysis, even if the docking results suggested that the enzyme has a clear preference for the last intermediates of such biosynthetic pathway. Since DHCR24 possesses a putative transmembrane segment, the enzyme was, then, inserted in a suitable membrane model and the membrane-anchored structure in complex with desmosterol and cholesterol underwent 10ns MD simulations. Such simulations evidenced a clearly different behavior between substrate and product since the product only completely leaves the catalytic cavity whereas desmosterol firmly conserves its pivotal interactions during all simulation time. This is one of the first reports documenting the enzymatic product egress using simple MD simulations in which all atoms are free to move.     

Full length Vpu from HIV-1: combining molecular dynamics simulations with NMR spectroscopy

Based on structures made available by solution NMR, molecular models of the protein Vpu from HIV-1 were built and refined by 6 ns MD simulations in a fully hydrated lipid bilayer. Vpu is an 81 amino acid type I integral membrane protein encoded by the human immunodeficiency virus type-1 (HIV-1) and closely related simian immunodeficiency viruses (SIVs). Its role is to amplify viral release. Upon phosphorylation, the cytoplasmic domain adopts a more compact shape with helices 2 and 3 becoming almost parallel to each other. A loss of helicity for several residues belonging to the helices adjacent to both ends of the loop region containing serines 53 and 57 is observed. A fourth helix, present in one of the NMR-based structures of the cytoplasmic domain and located near the C-terminus, is lost upon phosphorylation.     

Transmembrane helix structure, dynamics, and interactions: multi-nanosecond molecular dynamics simulations.

To probe the fundamentals of membrane/protein interactions, all-atom multi-nanosecond molecular dynamics simulations were conducted on a single transmembrane poly(32)alanine helix in a fully solvated dimyristoyphosphatidylcholine (DMPC) bilayer. The central 12 residues, which interact only with the lipid hydrocarbon chains, maintained a very stable helical structure. Helical regions extended beyond these central 12 residues, but interactions with the lipid fatty-acyl ester linkages, the lipid headgroups, and water molecules made the helix less stable in this region. The C and N termini, exposed largely to water, existed as random coils. As a whole, the helix tilted substantially, from perpendicular to the bilayer plane (0 degree) to a 30 degrees tilt. The helix experienced a bend at its middle, and the two halves of the helix at times assumed substantially different tilts. Frequent hydrogen bonding, of up to 0.7 ns in duration, occurred between peptide and lipid molecules. This resulted in correlated translational diffusion between the helix and a few lipid molecules. Because of the large variation in lipid conformation, the lipid environment of the peptide was not well defined in terms of "annular" lipids and on average consisted of 18 lipid molecules. When compared with a "neat" bilayer without peptide, no significant difference was seen in the bilayer thickness, lipid conformations or diffusion, or headgroup orientation. However, the lipid hydrocarbon chain order parameters showed a significant decrease in order, especially in those methylene groups closest to the headgroup.     

Molecular basis of Bcl-X(L)-p53 interaction: insights from molecular dynamics simulations

Bcl-X(L), an antiapoptotic Bcl-2 family protein, plays a central role in the regulation of the apoptotic pathway. Heterodimerization of the antiapoptotic Bcl-2 family proteins with the proapoptotic family members such as Bad, Bak, Bim and Bid is a crucial step in the apoptotic regulation. In addition to these conventional binding partners, recent evidences reveal that the Bcl-2 family proteins also interact with noncanonical binding partners such as p53. Our previous NMR studies showed that Bcl-X(L): BH3 peptide and Bcl-X(L): SN15 peptide (a peptide derived from residues S15-N29 of p53) complex structures share similar modes of bindings. To further elucidate the molecular basis of the interactions, here we have employed molecular dynamics simulations coupled with MM/PBSA approach. Bcl-X(L) and other Bcl-2 family proteins have 4 hydrophobic pockets (p1-p4), which are occupied by four systematically spaced hydrophobic residues (h1-h4) of the proapoptotic Bad and Bak BH3 peptides. We observed that three conserved hydrophobic residues (F19, W23 and L26) of p53 (SN15) peptide anchor into three hydrophobic pockets (p2-p4) of Bcl-X(L) in a similar manner as BH3 peptide. Our results provide insights into the novel molecular recognition by Bcl-X(L) with p53.     

Homology modeling, simulation and molecular docking studies of catechol-2, 3-Dioxygenase from Burkholderia cepacia: Involved in degradation of Petroleum hydrocarbons

Catechol 2, 3-dioxygenase is present in several types of bacteria and undergoes degradation of environmental pollutants through an important key biochemical pathways. Specifically, this enzyme cleaves aromatic rings of several environmental pollutants such as toluene, xylene, naphthalene and biphenyl derivatives. Hence, the importance of Catechol 2, 3-dioxygenase and its role in the degradation of environmental pollutants made us to predict the three-dimensional structure of Catechol 2, 3-dioxygenase from Burkholderia cepacia. The 10ns molecular dynamics simulation was carried out to check the stability of the modeled Catechol 2, 3- dioxygenase. The results show that the model was energetically stable, and it attains their equilibrium within 2000 ps of production MD run. The docking of various petroleum hydrocarbons into the Catechol 2,3-dioxygenase reveals that the benzene, O-xylene, Toluene, Fluorene, Naphthalene, Carbazol, Pyrene, Dibenzothiophene, Anthracene, Phenanthrene, Biphenyl makes strong hydrogen bond and Van der waals interaction with the active site residues of H150, L152, W198, H206, H220, H252, I254, T255, Y261, E271, L276 and F309. Free energy of binding and estimated inhibition constant of these compounds demonstrates that they are energetically stable in their binding cavity. Chrysene shows positive energy of binding in the active site atom of Fe. Except Pyrene all the substrates made close contact with Fe atom by the distance ranges from 1.67 to 2.43 Å. In addition to that, the above mentioned substrate except pyrene all other made π-π stacking interaction with H252 by the distance ranges from 3.40 to 3.90 Å. All these docking results reveal that, except Chrysene all other substrate has good free energy of binding to hold enough in the active site and makes strong VdW interaction with Catechol-2,3-dioxygenase. These results suggest that, the enzyme is capable of catalyzing the above-mentioned substrate.     

Multicomplex-based pharmacophore modeling in conjunction with multi-target docking and molecular dynamics simulations for the identification of PfDHFR inhibitors

Abstract

Plasmodium falciparum dihydrofolate reductase enzyme (PfDHFR) is counted as one of the attractive and validated antimalarial drug targets. However, the point mutations in the active site of wild-type PfDHFR have developed resistance against the well-known antifolates. Therefore, there is a dire need for the development of inhibitors that can inhibit both wild-type and mutant-type DHFR enzyme. In the present contribution, we have constructed the common feature pharmacophore models from the available PfDHFR. A representative hypothesis was prioritized and then employed for the screening of natural product library to search for the molecules with complementary features responsible for the inhibition. The screened candidates were processed via drug-likeness filters and molecular docking studies. The docking was carried out on the wild-type PfDHFR (3QGT); double-mutant PfDHFR (3UM5 and 1J3J) and quadruple-mutant PfDHFR (1J3K) enzymes. A total of eight common hits were obtained from the docking calculations that could be the potential inhibitors for both wild and mutant type DHFR enzymes. Eventually, the stability of these candidates with the selected proteins was evaluated via molecular dynamics simulations. Except for SPECS14, all the prioritized candidates were found to be stable throughout the simulation run. Overall, the strategy employed in the present work resulted in the retrieval of seven candidates that may show inhibitory activity against PfDHFR and could be further exploited as a scaffold to develop novel antimalarials.

Communicated by Ramaswamy H. Sarma