Data for molecular dynamics simulations of B-type cytochrome c oxidase with the Amber force field

Delineating redox cooperativity in water-soluble and membrane multiheme cytochromes through protein design

Nature has evolved diverse electron transport proteins and multiprotein assemblies essential to the generation and transduction of biological energy. However, substantially modifying or adapting these proteins for user-defined applications or to gain fundamental mechanistic insight can be hindered by their inherent complexity. De novo protein design offers an attractive route to stripping away this confounding complexity, enabling us to probe the fundamental workings of these bioenergetic proteins and systems, while providing robust, modular platforms for constructing completely artificial electron-conducting circuitry. Here, we use a set of de novo designed mono-heme and di-heme soluble and membrane proteins to delineate the contributions of electrostatic micro-environments and dielectric properties of the surrounding protein medium on the inter-heme redox cooperativity that we have previously reported. Experimentally, we find that the two heme sites in both the water-soluble and membrane constructs have broadly equivalent redox potentials in isolation, in agreement with Poisson-Boltzmann Continuum Electrostatics calculations. BioDC, a Python program for the estimation of electron transfer energetics and kinetics within multiheme cytochromes, also predicts equivalent heme sites, and reports that burial within the low dielectric environment of the membrane strengthens heme-heme electrostatic coupling. We conclude that redox cooperativity in our diheme cytochromes is largely driven by heme electrostatic coupling and confirm that this effect is greatly strengthened by burial in the membrane. These results demonstrate that while our de novo proteins present minimalist, new-to-nature constructs, they enable the dissection and microscopic examination of processes fundamental to the function of vital, yet complex, bioenergetic assemblies.

Reaction pathways, proton transfer, and proton pumping in ba3 class cytochrome c oxidase: perspectives from DFT quantum chemistry and molecular dynamics

After drawing comparisons between the reaction pathways of cytochrome c oxidase (CcO, Complex 4) and the preceding complex cytochrome bc1 (Complex 3), both being proton pumping complexes along the electron transport chain, we provide an analysis of the reaction pathways in bacterial ba3 class CcO, comparing spectroscopic results and kinetics observations with results from DFT calculations. For an important arc of the catalytic cycle in CcO, we can trace the energy pathways for the chemical protons and show how these pathways drive proton pumping of the vectorial protons. We then explore the proton loading network above the Fe heme a3-CuB catalytic center, showing how protons are loaded in and then released by combining DFT-based reaction energies with molecular dynamics simulations over states of that cycle. We also propose some additional reaction pathways for the chemical and vector protons based on our recent work with spectroscopic support.

An expandable, modular de novo protein platform for precision redox engineering

The electron-conducting circuitry of life represents an as-yet untapped resource of exquisite, nanoscale biomolecular engineering. Here, we report the characterization and structure of a de novo diheme "maquette" protein, 4D2, which we subsequently use to create an expanded, modular platform for heme protein design. A well-folded monoheme variant was created by computational redesign, which was then utilized for the experimental validation of continuum electrostatic redox potential calculations. This demonstrates how fundamental biophysical properties can be predicted and fine-tuned. 4D2 was then extended into a tetraheme helical bundle, representing a 7 nm molecular wire. Despite a molecular weight of only 24 kDa, electron cryomicroscopy illustrated a remarkable level of detail, indicating the positioning of the secondary structure and the heme cofactors. This robust, expressible, highly thermostable and readily designable modular platform presents a valuable resource for redox protein design and the future construction of artificial electron-conducting circuitry.

Current Trends and Changes in Use of Membrane Molecular Dynamics Simulations within Academia and the Pharmaceutical Industry

There has been an almost exponential increase in the use of molecular dynamics simulations in basic research and industry over the last 5 years, with almost a doubling in the number of publications each year. Many of these are focused on neurological membranes, and biological membranes in general, applied to the medical industry. A smaller portion have utilized membrane simulations to answer more basic questions related to the function of specific proteins, chemicals or biological processes. This review covers some newer studies, alongside studies from the last two decades, to determine changes in the field. Some of these are basic, while others are more profound, such as multi-component embedded membrane machinery. It is clear that many facets of the discipline remain the same, while the focus on and uses of the technology are broadening in scope and utilization as a general research tool. Analysis of recent literature provides an overview of the current methodologies, covers some of the recent trends or advances and tries to make predictions of the overall path membrane molecular dynamics will follow in the coming years. In general, the overview presented is geared towards the general scientific community, who may wish to introduce the use of these methodologies in light of these changes, making molecular dynamic simulations more feasible for general scientific or medical research.

Data for molecular dynamics simulations of Escherichia coli cytochrome bd oxidase with the Amber force field

Cytochrome bd-type quinol oxidase is an important metalloenzyme that allows many bacteria to survive in low oxygen conditions. Since bd oxidase is found in many prokaryotes but not in eukaryotes, it has emerged as a promising bacterial drug target. Examples of organisms containing bd oxidases include the Mycobacterium tuberculosis (Mtb) bacterium that causes tuberculosis (TB) in humans, the Vibrio cholerae bacterium that causes cholera, the Pseudomonas aeruginosa bacterium that contributes to antibiotic resistance and sepsis, and the Campylobacter jejuni bacterium that causes food poisoning. Escherichia coli (E. coli) is another organism exhibiting the cytochrome bd oxidase. Since it has the highest sequence identity to Mtb (36%) and we are ultimately interested in finding drug targets for TB, we have built parameters for the E. coli bd oxidase (Protein Data Bank ID number: 6RKO) that are compatible with the all-atom Amber ff14SB force field for molecular dynamics (MD) simulations. Specifically, we built parameters for the three heme cofactors present in all species of bacterial cytochrome bd-type oxidases (heme b558 , heme b595 , and heme d) along with their axial ligands. This data report includes the parameter and library files that can be used with Amber's LEaP program to generate input files for MD simulations using the Amber software package. We also provide the PDB data files of the initial model both by itself and solvated with TIP3P water molecules and counterions.

Subtype-selective mechanisms of negative allosteric modulators binding to group I metabotropic glutamate receptors

Group I metabotropic glutamate receptors (mGlu₁ and mGlu₅) are promising targets for multiple psychiatric and neurodegenerative disorders. Understanding the subtype selectivity of mGlu₁ and mGlu₅ allosteric sites is essential for the rational design of novel modulators with single- or dual-target mechanism of action. In this study, starting from the deposited mGlu₁ and mGlu₅ crystal structures, we utilized computational modeling approaches integrating docking, molecular dynamics simulation, and efficient post-trajectory analysis to reveal the subtype-selective mechanism of mGlu₁ and mGlu₅ to 10 diverse drug scaffolds representing known negative allosteric modulators (NAMs) in the literature. The results of modeling identified six pairs of non-conserved residues and four pairs of conserved ones as critical features to distinguish the selective NAMs binding to the corresponding receptors. In addition, nine pairs of residues are beneficial to the development of novel dual-target NAMs of group I metabotropic glutamate receptors. Furthermore, the binding modes of a reported dual-target NAM (VU0467558) in mGlu₁ and mGlu₅ were predicted to verify the identified residues that play key roles in the receptor selectivity and the dual-target binding. The results of this study can guide rational structure-based design of novel NAMs, and the approach can be generally applicable to characterize the features of selectivity for other G-protein-coupled receptors.

Virtual screening of quinoline derived library for SARS-COV-2 targeting viral entry and replication

The COVID-19 pandemic infection has claimed many lives and added to the social, economic, and psychological distress. The contagious disease has quickly spread to almost 218 countries and territories following the regional outbreak in China. As the number of infected populations increases exponentially, there is a pressing demand for anti-COVID drugs and vaccines. Virtual screening provides possible leads while extensively cutting down the time and resources required for ab-initio drug design. We report structure-based virtual screening of a hundred plus library of quinoline drugs with established antiviral, antimalarial, antibiotic or kinase inhibitor activity. In this study, targets having a role in viral entry, viral assembly, and viral replication have been selected. The targets include: 1) RBD of receptor-binding domain spike protein S 2) M^pro Chymotrypsin main protease 3) P^pro Papain protease 4) RNA binding domain of Nucleocapsid Protein, and 5) RNA Dependent RNA polymerase from SARS-COV-2. An in-depth analysis of the interactions and G-score compared to the controls like hydroxyquinoline and remdesivir has been presented. The salient results are (1) higher scoring of antivirals as potential drugs (2) potential of afatinib by scoring as better inhibitor, and (3) biological explanation of the potency of afatinib. Further MD simulations and MM-PBSA calculations showed that afatinib works best to interfere with the the activity of RNA dependent RNA polymerase of SARS-COV-2, thereby inhibiting replication process of single stranded RNA virus.

Communicated by Ramaswamy H. Sarma

Structure-Based Virtual Screening Identifies Multiple Stable Binding Sites at the RecA Domains of SARS-CoV-2 Helicase Enzyme

With the emergence and global spread of the COVID-19 pandemic, the scientific community worldwide has focused on search for new therapeutic strategies against this disease. One such critical approach is targeting proteins such as helicases that regulate most of the SARS-CoV-2 RNA metabolism. The purpose of the current study was to predict a library of phytochemicals derived from diverse plant families with high binding affinity to SARS-CoV-2 helicase (Nsp13) enzyme. High throughput virtual screening of the Medicinal Plant Database for Drug Design (MPD3) database was performed on SARS-CoV-2 helicase using AutoDock Vina. Nilotinib, with a docking value of -9.6 kcal/mol, was chosen as a reference molecule. A compound (PubChem CID: 110143421, ZINC database ID: ZINC257223845, eMolecules: 43290531) was screened as the best binder (binding energy of -10.2 kcal/mol on average) to the enzyme by using repeated docking runs in the screening process. On inspection, the compound was disclosed to show different binding sites of the triangular pockets collectively formed by Rec1A, Rec2A, and 1B domains and a stalk domain at the base. The molecule is often bound to the ATP binding site (referred to as binding site 2) of the helicase enzyme. The compound was further discovered to fulfill drug-likeness and lead-likeness criteria, have good physicochemical and pharmacokinetics properties, and to be non-toxic. Molecular dynamic simulation analysis of the control/lead compound complexes demonstrated the formation of stable complexes with good intermolecular binding affinity. Lastly, affirmation of the docking simulation studies was accomplished by estimating the binding free energy by MMPB/GBSA technique. Taken together, these findings present further in silco investigation of plant-derived lead compounds to effectively address COVID-19.

Exploring the Binding Mechanism of Metabotropic Glutamate Receptor 5 Negative Allosteric Modulators in Clinical Trials by Molecular Dynamics Simulations

Metabotropic glutamate receptor 5 (mGlu₅) plays a key role in synaptic information storage and memory, which is a well-known target for a variety of psychiatric and neurodegenerative disorders. In recent years, the increasing efforts have been focused on the design of allosteric modulators, and the negative allosteric modulators (NAMs) are the front-runners. Recently, the architecture of the transmembrane (TM) domain of mGlu₅ receptor has been determined by crystallographic experiment. However, it has been not well understood how the pharmacophores of NAMs accommodated into the allosteric binding site. In this study, molecular dynamics (MD) simulations were performed on mGlu₅ receptor bound with NAMs in preclinical or clinical development to shed light on this issue. In order to identify the key residues, the binding free energies as well as per-residue contributions for NAMs binding to mGlu₅ receptor were calculated. Subsequently, the in silico site-directed mutagenesis of the key residues was performed to verify the accuracy of simulation models. As a result, the shared common features of the studied 5 clinically important NAMs (mavoglurant, dipraglurant, basimglurant, STX107, and fenobam) interacting with 11 residues in allosteric site were obtained. This comprehensive study presented a better understanding of mGlu₅ receptor NAMs binding mechanism, which would be further used as a useful framework to assess and discover novel lead scaffolds for NAMs.

Molecular Determinants of the Promiscuity of MexB and MexY Multidrug Transporters of Pseudomonas aeruginosa

Secondary multidrug transporters of the resistance-nodulation-cell division (RND) superfamily contribute crucially to antibiotic resistance in Gram-negative bacteria. Compared to the most studied transporter AcrB of Escherichia coli, little is known about the molecular determinants of distinct polyspecificities of the most important RND transporters MexB and MexY of Pseudomonas aeruginosa. In an effort to add knowledge on this topic, we performed an exhaustive atomic-level comparison of the main putative recognition sites (access and deep binding pockets) in these two Mex transporters. We identified an underlying link between some structural, chemical and dynamical features of the binding pockets and the physicochemical nature of the corresponding substrates recognized by either one or both pumps. In particular, mosaic-like lipophilic and electrostatic surfaces of the binding pockets provide for both proteins several multifunctional sites for diffuse binding of diverse substrates. Specific lipophilicity signatures of the weakly conserved deep pocket suggest a key role of this site as a selectivity filter as in Acr transporters. Finally, the different dynamics of the bottom-loop in MexB and MexY support its possible role in binding of large substrates. Our work represents the first comparative study of the major RND transporters in P. aeruginosa and also the first structure-based study of MexY, for which no experimental structure is available yet.