Molecular Modeling of the Ternary Complex of Rusticyanin-Cytochrome c4-Cytochrome Oxidase: An Insight to Possible H-Bond Mediated Recognition and Electron Transfer Reaction inT.ferrooxidans

Evaluation of Cyc1 protein stability in Acidithiobacillus ferrooxidans bacterium after E121D mutation by molecular dynamics simulation to improve electron transfer

Cyc₁ (Cytochrome c₅₅₂) is a protein in the electron transport chain of the Acidithiobacillus ferrooxidans (Af) bacteria which obtain their energy from oxidation Fe²⁺ to Fe³⁺. The electrons are directed through Cyc₂, RCY (rusticyanin), Cyc₁ and Cox aa₃ proteins to O₂. Cyc₁ protein consists of two chains, A and B. In the present study, a novel mutation (E121D) in the A chain of Cyc₁ protein was selected due to electron receiving from Histidine 143 of RCY. Then, the changes performed in the E121D mutant were evaluated by MD simulations analyzes. Cyc₁ and RCY proteins were docked by a Patchdock server. By E121D mutation, the connection between Zn 1388 of chain B and aspartate 121 of chain A weaken. Asp 121 gets farther from Zn 1388. Therefore, the aspartate gets closer to Cu 1156 of the RCY leading to the higher stability of the RCY/Cyc₁ complex. Further, an acidic residue (Glu121) becomes a more acidic residue (Asp121) and improves the electron transfer to Cyc₁ protein. The results of RMSF analysis showed further ligand flexibility in mutation. This leads to fluctuation of the active site and increases redox potential at the mutation point and the speed of electron transfer. This study also predicts that in all respiratory chain proteins, electrons probably enter the first active site via glutamate and exit histidine in the second active site of each respiratory chain protein.

Computational structure prediction provides a plausible mechanism for electron transfer by the outer membrane protein Cyc2 from Acidithiobacillus ferrooxidans

Cyc2 is the key protein in the outer membrane of Acidithiobacillus ferrooxidans that mediates electron transfer between extracellular inorganic iron and the intracellular central metabolism. This cytochrome c is specific for iron and interacts with periplasmic proteins to complete a reversible electron transport chain. A structure of Cyc2 has not yet been characterized experimentally. Here we describe a structural model of Cyc2, and associated proteins, to highlight a plausible mechanism for the ferrous iron electron transfer chain. A comparative modeling protocol specific for trans membrane beta barrel (TMBB) proteins in acidophilic conditions (pH ~ 2) was applied to the primary sequence of Cyc2. The proposed structure has three main regimes: Extracellular loops exposed to low-pH conditions, a TMBB, and an N-terminal cytochrome-like region within the periplasmic space. The Cyc2 model was further refined by identifying likely iron and heme docking sites. This represents the first computational model of Cyc2 that accounts for the membrane microenvironment and the acidity in the extracellular matrix. This approach can be used to model other TMBBs which can be critical for chemolithotrophic microbial growth.

Investigation of Cyc1 protein structure stability after H53I mutation using computational approaches to improve redox potential

Acidithiobacillus ferrooxidans (Af) is an acidophilic bacterium that grows in rigid surroundings and gets its own energy from the oxidation of Fe²⁺ to Fe³⁺. These bacteria are involved in the bioleaching process. Cyc₁ is a periplasmic protein with a crucial role in electron transportation in the respiratory chain. His53 of the Cyc₁ protein, involved in electron transfer to CoxB, was selected for mutation and bioinformatics studies. His53 was substituted by Ile using PyMol software. Molecular dynamics simulations were performed for wild and mutant types of Cyc₁ protein. The conformational changes of mutated protein were studied by analyzing RMSD, RMSF, SASA, Rg, H Bond, and DSSP. The results of the RMSF analysis indicated an increase in the flexibility of the ligand in the mutant. Finally, active site instability leads to an increase in the value of E₀ at the mutation point and improving electron transfer. On the other, His53 in Cyc₁ is interconnected to Glu126 in CoxB through the water molecule (W76) and hydrogen bonding. In the H53I mutation, there was a decrease in the distance between H2O 2030, 2033, and isoleucine 53, and subsequently, the distance to the water molecule 76 between the two proteins was reduced and strengthens the hydrogen bond between Cyc₁ and CoxB, finally improves electron transfer and the bioleaching process.

Recognition dynamics of trinuclear copper cluster and associated histidine residues through conserved or semi-conserved water molecules in human Ceruloplasmin: The involvement of aspartic and glutamic acid gates

Human Ceruloplasmin belongs to the family of multi-copper oxidases and it is involved in different physiological processes, copper ion transport, iron metabolism, iron homeostasis, and biogenic amine metabolism. MD-simulation studies have indicated the higher hydrophilic susceptibility of the trinuclear copper cluster in native CP compared to its oxygen bound form. The copper (T2/T3) atom Cu3047 of the cluster, which is close to T1 copper center Cu3052 (~13 Å) has a higher affinity for water molecules compared to other copper centers. The water molecules of W3, W4, W5, W9, and W12 conserved water sites are coordinated to Cu3047, where W3, W9, and W12 centers are found to play some crucial role in the stabilization of native trinuclear copper cluster. The hydrogen bonding interaction of Asp169, Glu112, Asp995, and Glu1032 residues with the copper-bound conserved water molecules (W3, W4, W5, W10, and W11) in native CP is observed to be unique. The conformational flexibility of Asp169 and Glu112 and their association with the copper-bound water molecules, but the absence of such interaction in O₂-bound simulated structure of the enzyme is indicating some plausible rational on the role of these acidic residues in the gating of O₂ molecule in the native trinuclear Cu cluster of CP. The simulation results may shade some new light on the biochemistry/chemistry of CP, specially on the hydration dynamics of the trinuclear copper cluster.

Homology modeling and docking studies of FabH (β-ketoacyl-ACP synthase III) enzyme involved in type II fatty acid biosynthesis of Chlorella variabilis: a potential algal feedstock for biofuel production

The concept of using microalgae as an alternative renewable source of biofuel has gained much importance in recent years. However, its commercial feasibility is still an area of concern for researchers. Unraveling the fatty acid metabolic pathway and understanding structural features of various key enzymes regulating the process will provide valuable insights to target microalgae for augmented oil content. FabH (β-ketoacyl-acyl carrier protein synthase; KAS III) is a condensing enzyme catalyzing the initial elongation step of type II fatty acid biosynthetic process and acyl carrier protein (ACP) facilitates the shuttling of the fatty acyl intermediates to the active site of the respective enzymes in the pathway. In the present study, a reliable three-dimensional structure of FabH from Chlorella variabilis, an oleaginous green microalga was modeled and subsequently the key residues involved in substrate binding were determined by employing protein-protein docking and molecular dynamics (MD) simulation protocols. The FabH-ACP complex having the lowest docking energy score showed the binding of ACP to the electropositive FabH surface with strong hydrogen bond interactions. The MD simulation results indicated that the substrate-complexed FabH adopted a more stable conformation than the free enzyme. Further, the FabH structure retained its stability throughout the simulation although noticeable displacements were observed in the loop regions. Molecular simulation studies suggested the importance of crucial hydrogen bonding of the conserved Arg(91) of FabH with Glu(53) and Asp(56) of ACP for exhibiting high affinity between the enzyme and substrate. The molecular modeling results are consistent with available experimental results on the flexibility of FabH and the present study provides first in silico insights into the structural and dynamical aspect of catalytic mechanism of FabH, which could be used for further site-specific mutagenic experiments to develop engineered high oil-yielding microalgal strains for biofuel production.

Homology modeling of Cry10Aa toxin from B. thuringiensis israelensis and B. thuringiensis subsp. LDC-9

A three dimensional model was developed for Cry10Aa protein sequence of B. thuringiensis LDC-9 and B. thuringiensis israelensis that has not been solved empirically by X-ray crystallography or NMR. Homology modeling was employed for the structure prediction using Cry2Aa as template protein, a high-resolution X-ray crystallography structure. The model predicted for the B. thuringiensis LDC-9 Cry10Aa protein reveals a partial N-terminal domain only due to its partial sequence of 104 amino acids. B. thuringiensis israelensis Cry10Aa model contains three domains such as domain I, a bundle of eight alpha helices with the central relatively hydrophobic helix surrounded by amphipathic helices while domain II and III contain mostly beta-sheets. Significant structural differences within domain II in this model among all Cry protein structures indicates that it is involved in recognition and binding to cell surfaces. Comparison of B. thuringiensis israelensis predicted structure with available experimentally determined Cry structures reveals identical folds. The distribution of electrostatic potential on the surface of the molecules in the model is non-uniform and identifies one side of the alpha-helical domain as negatively charged indicating orientation of toxic molecules toward the cell membrane during the initial binding with a cell surface receptor. The collective knowledge of Cry toxin structures will lead to a more critical understanding of the structural basis for receptor binding and pore formation, as well as allowing the scope of diversity to be better appreciated. This model will serve as a starting point for the design of mutagenesis experiments aimed to improve the toxicity and to provide a new tool for the elucidation of the mechanism of action of these mosquitocidal proteins.

Anti-infectious bronchitis virus (IBV) activity of 1,8-cineole: effect on nucleocapsid (N) protein

In the present study, anti-IBV (infectious bronchitis virus) activity of 1,8-cineole was studied by MTT assay, as well as docking and molecular dynamic (MD) simulations. The CC50 of 1,8-cineole was above 10 mM. And the maximum noncytotoxic concentration (TD0) of 1,8-cineole was determined to be 3.90 ± 0.22 mM, which was much higher than that of ribavirin (0.78 ± 0.15 mM). 1,8-cineole could inhibit IBV with an IC(50) of 0.61 mM. MTT assay showed that the inhibition of IBV by 1, 8-cineole appears to occur moderately before entering the cell but much strongly after penetration of the virus into the cell. In silico simulations indicated that the binding site of 1,8-cineole was located at the N terminus of phosphorylated nucleocapsid (N) protein, with interaction energy equaling -40.33 kcal mol(-1). The residues TyrA92, ProA134, PheA137, AspA138 and TyrA140 had important roles during the binding process and are fully or partially conserved in various IBV strains. Based on spatial and energetic criteria, 1,8-cineole interfered with the binding between RNA and IBV N-protein. Results presented here may suggest that 1,8-cineole possesses anti-IBV properties, and therefore is a potential source of anti-IBV ingredients for the pharmaceutical industry.

Cupredoxins--a study of how proteins may evolve to use metals for bioenergetic processes

Cupredoxins are small proteins that contain type I copper centers, which are ubiquitous in nature. They function as electron transfer shuttles between proteins. This review of the structure and properties of native cupredoxins, and those modified by site-directed mutagenesis, illustrates how these proteins may have evolved to specifically bind copper, develop recognition sites for specific redox partners, tune redox potential for a particular function, and allow for efficient electron transfer through the protein matrix. This is relevant to the general understanding of the roles of metals in energy metabolism, respiration and photosynthesis.

The diheme cytochrome c(4) from Vibrio cholerae is a natural electron donor to the respiratory cbb(3) oxygen reductase

The respiratory chain of Vibrio cholerae contains three bd-type quinol oxygen reductases as well as one cbb(3) oxygen reductase. The cbb(3) oxygen reductase has been previously isolated and characterized; however, the natural mobile electron donor(s) that shuttles electrons between the bc(1) complex and the cbb(3) oxygen reductase is not known. The most likely candidates are the diheme cytochrome c(4) and monoheme cytochrome c(5), which have been previously shown to be present in the periplasm of aerobically grown cultures of V. cholerae. Both cytochromes c(4) and c(5) from V. cholerae have been cloned and expressed heterologously in Escherichia coli. It is shown that reduced cytochrome c(4) is a substrate for the purified cbb(3) oxygen reductase and can support steady state oxygen reductase activity of at least 300 e(-1)/s. In contrast, reduced cytochrome c(5) is not a good substrate for the cbb(3) oxygen reductase. Surprisingly, the dependence of the oxygen reductase activity on the concentration of cytochrome c(4) does not exhibit saturation. Global spectroscopic analysis of the time course of the oxidation of cytochrome c(4) indicates that the apparent lack of saturation is due to the strong dependence of K(M) and V(max) on the concentration of oxidized cytochrome c(4). Whether this is an artifact of the in vitro assay or has physiological significance remains unknown. Cyclic voltammetry was used to determine that the midpoint potentials of the two hemes in cytochrome c(4) are 240 and 340 mV (vs standard hydrogen electrode), similar to the electrochemical properties of other c(4)-type cytochromes. Genomic analysis shows a strong correlation between the presence of a c(4)-type cytochrome and a cbb(3) oxygen reductase within the beta- and gamma-proteobacterial clades, suggesting that cytochrome c(4) is the likely natural electron donor to the cbb(3) oxygen reductases within these organisms. These would include the beta-proteobacteria Neisseria meningitidis and Neisseria gonnorhoeae, in which the cbb(3) oxygen reductases are the only terminal oxidases in their respiratory chains, and the gamma-proteobacterium Pseudomonas stutzeri.

Structural model of the pre-pore ring-like structure of Panton-Valentine leukocidin: providing dimensionality to biophysical and mutational data

Panton-Valentine leukocidin (PVL) is a bipartite toxin that plays an important role in the pathogenesis of methicillin-resistant Staphylococcus aureus. Recent clinical data suggest a correlation between PVL and severe cases of S. aureus pneumonia. A clear understanding of the structure and function of PVL is critical to the development of novel, effective treatments. Here, we report an all-atom model of the macromolecular structure of Panton-Valentine leukocidin in its octameric, pre-pore conformation that confirms and extends our understanding of the toxin's mechanism of action.