Revised CHARMM force field parameters for iron-containing cofactors of photosystem II

Structural basis of human NOX5 activation

NADPH oxidase 5 (NOX5) catalyzes the production of superoxide free radicals and regulates physiological processes from sperm motility to cardiac rhythm. Overexpression of NOX5 leads to cancers, diabetes, and cardiovascular diseases. NOX5 is activated by intracellular calcium signaling, but the underlying molecular mechanism of which - in particular, how calcium triggers electron transfer from NADPH to FAD - is still unclear. Here we capture motions of full-length human NOX5 upon calcium binding using single-particle cryogenic electron microscopy (cryo-EM). By combining biochemistry, mutagenesis analyses, and molecular dynamics (MD) simulations, we decode the molecular basis of NOX5 activation and electron transfer. We find that calcium binding to the EF-hand domain increases NADPH dynamics, permitting electron transfer between NADPH and FAD and superoxide production. Our structural findings also uncover a zinc-binding motif that is important for NOX5 stability and enzymatic activity, revealing modulation mechanisms of reactive oxygen species (ROS) production.

Perspective on Integrative Simulations of Bioenergetic Domains

Bioenergetic processes in cells, such as photosynthesis or respiration, integrate many time and length scales, which makes the simulation of energy conversion with a mere single level of theory impossible. Just like the myriad of experimental techniques required to examine each level of organization, an array of overlapping computational techniques is necessary to model energy conversion. Here, a perspective is presented on recent efforts for modeling bioenergetic phenomena with a focus on molecular dynamics simulations and its variants as a primary method. An overview of the various classical, quantum mechanical, enhanced sampling, coarse-grained, Brownian dynamics, and Monte Carlo methods is presented. Example applications discussed include multiscale simulations of membrane-wide electron transport, rate kinetics of ATP turnover from electrochemical gradients, and finally, integrative modeling of the chromatophore, a photosynthetic pseudo-organelle.

GalaxyDock2-HEME: Protein-ligand docking for heme proteins

Prediction of protein-ligand binding poses is an essential component for understanding protein-ligand interactions and computer-aided drug design. Various proteins involve prosthetic groups such as heme for their functions, and adequate consideration of the prosthetic groups is vital for protein-ligand docking. Here, we extend the GalaxyDock2 protein-ligand docking algorithm to handle ligand docking to heme proteins. Docking to heme proteins involves increased complexity because the interaction of heme iron and ligand has covalent nature. GalaxyDock2-HEME, a new protein-ligand docking program for heme proteins, has been developed based on GalaxyDock2 by adding an orientation-dependent scoring term to describe heme iron-ligand coordination interaction. This new docking program performs better than other noncommercial docking programs such as EADock with MMBP, AutoDock Vina, PLANTS, LeDock, and GalaxyDock2 on a heme protein-ligand docking benchmark set in which ligands are known to bind iron. In addition, docking results on two other sets of heme protein-ligand complexes in which ligands do not bind iron show that GalaxyDock2-HEME does not have a high bias toward iron binding compared to other docking programs. This implies that the new docking program can distinguish iron binders from noniron binders for heme proteins.

In silico investigation of cytochrome bc1 molecular inhibition mechanism against Trypanosoma cruzi

Chagas' disease is a neglected tropical disease caused by the kinetoplastid protozoan Trypanosoma cruzi. The only therapies are the nitroheterocyclic chemicals nifurtimox and benznidazole that cause various adverse effects. The need to create safe and effective medications to improve medical care remains critical. The lack of verified T. cruzi therapeutic targets hinders medication research for Chagas' disease. In this respect, cytochrome bc1 has been identified as a promising therapeutic target candidate for antibacterial medicines of medical and agricultural interest. Cytochrome bc1 belongs to the mitochondrial electron transport chain and transfers electrons from ubiquinol to cytochrome c1 by the action of two catalytic sites named Qi and Qo. The two binding sites are highly selective, and specific inhibitors exist for each site. Recent studies identified the Qi site of the cytochrome bc1 as a promising drug target against T. cruzi. However, a lack of knowledge of the drug mechanism of action unfortunately hinders the development of new therapies. In this context, knowing the cause of binding site selectivity and the mechanism of action of inhibitors and substrates is crucial for drug discovery and optimization processes. In this paper, we provide a detailed computational investigation of the Qi site of T. cruzi cytochrome b to shed light on the molecular mechanism of action of known inhibitors and substrates. Our study emphasizes the action of inhibitors at the Qi site on a highly unstructured portion of cytochrome b that could be related to the biological function of the electron transport chain complex.

A versatile artificial metalloenzyme scaffold enabling direct bioelectrocatalysis in solution

Artificial metalloenzymes (ArMs) are commonly designed with protein scaffolds containing buried coordination pockets to achieve substrate specificity and product selectivity for homogeneous reactions. However, their reactivities toward heterogeneous transformations are limited because interfacial electron transfers are hampered by the backbone shells. Here, we introduce bacterial small laccase (SLAC) as a new protein scaffold for constructing ArMs to directly catalyze electrochemical transformations. We use molecular dynamics simulation, x-ray crystallography, spectroscopy, and computation to illustrate the scaffold-directed assembly of an oxo-bridged dicobalt motif on protein surface. The resulting ArM in aqueous phase catalyzes electrochemical water oxidation without mediators or electrode modifications. Mechanistic investigation reveals the role of SLAC scaffold in defining the four-electron transfer pathway from water to oxygen. Furthermore, we demonstrate that SLAC-based ArMs implemented with Ni²⁺, Mn²⁺, Ru³⁺, Pd²⁺, or Ir³⁺ also enable direct bioelectrocatalysis of water electrolysis. Our study provides a versatile and generalizable route to complement heterogeneous repertoire of ArMs for expanded applications.

Computer-Aided Rational Construction of Mediated Bioelectrocatalysis with π-Conjugated (Hetero)cyclic Molecules: Toward Promoted Distant Electron Tunneling and Improved Biosensing

Highly π-conjugated (hetero)cyclic molecules having delocalized orbitals and tunable charge mobilities are attractive redox relays for mediated bioelectrocatalysis in manifold applications. As rigid molecules, their dynamics within the soft but confined intraprotein space becomes the crucial determinant of the enzyme-mediator electron-tunneling efficiency. However, it is rarely investigated in designing the mediated interface with a particular biocatalyst (e.g., oxidoreductase), which remains an empirical but try-and-error process. Herein, we present the computer-aided exploration of interactions between a flavin-containing reductive synthase and structurally diverse π-extended (hetero)cyclic mediators to realize reversed bioelectrocatalytic oxidation at low overpotentials. Compared to ring-fused systems with unbroken molecular planarity, heteroatom-bridged cyclic, and rotatable conjugated structures (e.g., indophenols) can experience unusually large dynamic torsion under biased forces of hydrogen bonding with enzyme residues. This behavior led to fast intraprotein reorientation (<50 ps) that shortened the electron-tunneling distance from 12 to 9 Å. Meanwhile, the lowest unoccupied molecular orbital level upon molecular torsion was decreased by 0.5 eV to further promote electron abstraction from the reduced flavin cofactor. An efficient distant electron tunneling also obviated mediator transport through the substrate channel, thus avoiding competitive inhibition on enzyme kinetics to broaden the operating concentration range. The resulting bioelectrocatalytic interface enables low-potential biosensing of glutamate with improved selectivity. Our finding provides new structural insights into constructing efficient long-range heterogeneous charge transport with biomacromolecular catalysts.

Exploring the folding energy landscapes of heme proteins using a hybrid AWSEM-heme model

Heme is an active center in many proteins. Here we explore computationally the role of heme in protein folding and protein structure. We model heme proteins using a hybrid model employing the AWSEM Hamiltonian, a coarse-grained forcefield for the protein chain along with AMBER, an all-atom forcefield for the heme. We carefully designed transferable force fields that model the interactions between the protein and the heme. The types of protein-ligand interactions in the hybrid model include thioester covalent bonds, coordinated covalent bonds, hydrogen bonds, and electrostatics. We explore the influence of different types of hemes (heme b and heme c) on folding and structure prediction. Including both types of heme improves the quality of protein structure predictions. The free energy landscape shows that both types of heme can act as nucleation sites for protein folding and stabilize the protein folded state. In binding the heme, coordinated covalent bonds and thioester covalent bonds for heme c drive the heme toward the native pocket. The electrostatics also facilitates the search for the binding site.

Decoding the Molecular Effects of Atovaquone Linked Resistant Mutations on Plasmodium falciparum Cytb-ISP Complex in the Phospholipid Bilayer Membrane

Atovaquone (ATQ) is a drug used to prevent and treat malaria that functions by targeting the Plasmodium falciparum cytochrome b (PfCytb) protein. PfCytb catalyzes the transmembrane electron transfer (ET) pathway which maintains the mitochondrial membrane potential. The ubiquinol substrate binding site of the protein has heme bL, heme bH and iron-sulphur [2FE-2S] cluster cofactors that act as redox centers to aid in ET. Recent studies investigating ATQ resistance mechanisms have shown that point mutations of PfCytb confer resistance. Thus, understanding the resistance mechanisms at the molecular level via computational approaches incorporating phospholipid bilayer would help in the design of new efficacious drugs that are also capable of bypassing parasite resistance. With this knowledge gap, this article seeks to explore the effect of three drug resistant mutations Y268C, Y268N and Y268S on the PfCytb structure and function in the presence and absence of ATQ. To draw reliable conclusions, 350 ns all-atom membrane (POPC:POPE phospholipid bilayer) molecular dynamics (MD) simulations with derived metal parameters for the holo and ATQ-bound -proteins were performed. Thereafter, simulation outputs were analyzed using dynamic residue network (DRN) analysis. Across the triplicate MD runs, hydrophobic interactions, reported to be crucial in protein function were assessed. In both, the presence and absence of ATQ and a loss of key active site residue interactions were observed as a result of mutations. These active site residues included: Met 133, Trp136, Val140, Thr142, Ile258, Val259, Pro260 and Phe264. These changes to residue interactions are likely to destabilize the overall intra-protein residue communication network where the proteins' function could be implicated. Protein dynamics of the ATQ-bound mutant complexes showed that they assumed a different pose to the wild-type, resulting in diminished residue interactions in the mutant proteins. In summary, this study presents insights on the possible effect of the mutations on ATQ drug activity causing resistance and describes accurate MD simulations in the presence of the lipid bilayer prior to conducting inhibitory drug discovery for the PfCytb-iron sulphur protein (Cytb-ISP) complex.

Molecular dynamics simulations in photosynthesis

Photosynthesis is regulated by a dynamic interplay between proteins, enzymes, pigments, lipids, and cofactors that takes place on a large spatio-temporal scale. Molecular dynamics (MD) simulations provide a powerful toolkit to investigate dynamical processes in (bio)molecular ensembles from the (sub)picosecond to the (sub)millisecond regime and from the Å to hundreds of nm length scale. Therefore, MD is well suited to address a variety of questions arising in the field of photosynthesis research. In this review, we provide an introduction to the basic concepts of MD simulations, at atomistic and coarse-grained level of resolution. Furthermore, we discuss applications of MD simulations to model photosynthetic systems of different sizes and complexity and their connection to experimental observables. Finally, we provide a brief glance on which methods provide opportunities to capture phenomena beyond the applicability of classical MD.

A sixty-year tryst with photosynthesis and related processes: an informal personal perspective

After briefly describing my early collaborative work at the University of Allahabad, that had laid the foundation of my research life, I present here some of our research on photosynthesis at the University of Illinois at Urbana-Champaign, randomly selected from light absorption to NADP⁺ reduction in plants, algae, and cyanobacteria. These include the fact that (i) both the light reactions I and II are powered by light absorbed by chlorophyll (Chl) a of different spectral forms; (ii) light emission (fluorescence, delayed fluorescence, and thermoluminescence) by plants, algae, and cyanobacteria provides detailed information on these reactions and beyond; (iii) primary photochemistry in both the photosystems I (PS I) and II (PS II) occurs within a few picoseconds; and (iv) most importantly, bicarbonate plays a unique role on the electron acceptor side of PS II, specifically at the two-electron gate of PS II. Currently, the ongoing research around the world is, and should be, directed towards making photosynthesis better able to deal with the global issues (such as increasing population, dwindling resources, and rising temperature) particularly through genetic modification. However, basic research is necessary to continue to provide us with an understanding of the molecular mechanism of the process and to guide us in reaching our goals of increasing food production and other chemicals we need for our lives.

Dynamic multiple-scattering treatment of X-ray absorption: Parameterization of a new molecular dynamics force field for myoglobin

We present a detailed analysis of the X-ray absorption near-edge structure (XANES) data on the Fe K-edge of CO Myoglobin based on a combined procedure of Molecular Dynamics (MD) calculations and MXAN (Minuit XANes) data analysis that we call D-MXAN. The ability of performing quantitative XANES data analysis allows us to refine classical force field MD parameters, thus obtaining a reliable tool for the atomic investigation of this important model system for biological macromolecules. The iterative procedure here applied corrects the greatest part of the structural discrepancy between classical MD sampling and experimental determinations. Our procedure, moreover, is able to discriminate between different heme conformational basins visited during the MD simulation, thus demonstrating the necessity of a sampling on the order of tens of nanoseconds, even for an application such X-ray absorption spectroscopy data analysis.

AMOEBA Polarizable Force Field Parameters of the Heme Cofactor in Its Ferrous and Ferric Forms

We report the first parameters of the heme redox cofactors for the polarizable AMOEBA force field in both the ferric and ferrous forms. We consider two types of complexes, one with two histidine side chains as axial ligands and one with a histidine and a methionine side chain as ligands. We have derived permanent multipoles from second-order Møller-Plesset perturbation theory (MP2). The sets of parameters have been validated in a first step by comparison of AMOEBA interaction energies of heme and a collection of biologically relevant molecules with MP2 and Density Functional Theory (DFT) calculations. In a second validation step, we consider interaction energies with large aggregates comprising around 80 H₂O molecules. These calculations are repeated for 30 structures extracted from semiempirical PM7 DM simulations. Very encouraging agreement is found between DFT and the AMOEBA force field, which results from an accurate treatment of electrostatic interactions. We finally report long (10 ns) MD simulations of cytochromes in two redox states with AMOEBA testing both the 2003 and 2014 AMOEBA water models. These simulations have been carried out with the TINKER-HP (High Performance) program. In conclusion, owing to their ubiquity in biology, we think the present work opens a wide array of applications of the polarizable AMOEBA force field on hemeproteins.

Development of CHARMM-Compatible Force-Field Parameters for Cobalamin and Related Cofactors from Quantum Mechanical Calculations

Corrinoid cofactors such as cobalamin are used by many enzymes and are essential for most living organisms. Therefore, there is broad interest in investigating cobalamin-protein interactions with molecular dynamics simulations. Previously developed parameters for cobalamins are based mainly on crystal structure data. Here, we report CHARMM-compatible force field parameters for several corrinoids developed from quantum mechanical calculations. We provide parameters for corrinoids in three oxidation states, Co³⁺, Co²⁺, and Co¹⁺, and with various axial ligands. Lennard-Jones parameters for the cobalt center in the Co(II) and Co(I) states were optimized using a helium atom probe, and partial atomic charges were obtained with a combination of natural population analysis (NPA) and restrained electrostatic potential (RESP) fitting approaches. The Force Field Toolkit was used to optimize all bonded terms. The resulting parameters, determined solely from calculations of cobalamin alone or in water, were then validated by assessing their agreement with density functional theory geometries and by analyzing molecular dynamics simulation trajectories of several corrinoid proteins for which X-ray crystal structures are available. In each case, we obtained excellent agreement with the reference data. In comparison to previous CHARMM-compatible parameters for cobalamin, we observe a better agreement for the fold angle and lower RMSD in the cobalamin binding site. The approach described here is readily adaptable for developing CHARMM-compatible force-field parameters for other corrinoids or large biomolecules.