Molecular dynamics in cytochrome c oxidase Mössbauer spectra deconvolution

The Dynamics of OXA-23 β-Lactamase from Acinetobacter baumannii

Antibiotic resistance is a pressing topic, which also affects β-lactam antibiotic molecules. Until a few years ago, it was considered no more than an interesting species from an academic point of view, Acinetobacter baumanii is today one of the most serious threats to public health, so much so that it has been declared one of the species for which the search for new antibiotics, or new ways to avoid its resistance, is an absolute priority according to WHO. Although there are several molecular mechanisms that are responsible for the extreme resistance of A. baumanii to antibiotics, a class D β-lactamase is the main cause for the clinical concern of this bacterial species. In this work, we analyzed the A. baumanii OXA-23 protein via molecular dynamics. The results obtained show that this protein is able to assume different conformations, especially in some regions around the active site. Part of the OXA-23 protein has considerable conformational motility, while the rest is less mobile. The importance of these observations for understanding the functioning mechanism of the enzyme as well as for designing new effective molecules for the treatment of A. baumanii is discussed.

AI-Aided Search for New HIV-1 Protease Ligands

The availability of drugs capable of blocking the replication of microorganisms has been one of the greatest triumphs in the history of medicine, but the emergence of an ever-increasing number of resistant strains poses a serious problem for the treatment of infectious diseases. The search for new potential ligands for proteins involved in the life cycle of pathogens is, therefore, an extremely important research field today. In this work, we have considered the HIV-1 protease, one of the main targets for AIDS therapy. Several drugs are used today in clinical practice whose mechanism of action is based on the inhibition of this enzyme, but after years of use, even these molecules are beginning to be interested by resistance phenomena. We used a simple artificial intelligence system for the initial screening of a data set of potential ligands. These results were validated by docking and molecular dynamics, leading to the identification of a potential new ligand of the enzyme which does not belong to any known class of HIV-1 protease inhibitors. The computational protocol used in this work is simple and does not require large computational power. Furthermore, the availability of a large number of structural information on viral proteins and the presence of numerous experimental data on their ligands, with which it is possible to compare the results obtained with computational methods, make this research field the ideal terrain for the application of these new computational techniques.

Targeting mitochondrial impairment for the treatment of cardiovascular diseases: From hypertension to ischemia-reperfusion injury, searching for new pharmacological targets

Mitochondria and mitochondrial proteins represent a group of promising pharmacological target candidates in the search of new molecular targets and drugs to counteract the onset of hypertension and more in general cardiovascular diseases (CVDs). Indeed, several mitochondrial pathways result impaired in CVDs, showing ATP depletion and ROS production as common traits of cardiac tissue degeneration. Thus, targeting mitochondrial dysfunction in cardiomyocytes can represent a successful strategy to prevent heart failure. In this context, the identification of new pharmacological targets among mitochondrial proteins paves the way for the design of new selective drugs. Thanks to the advances in omics approaches, to a greater availability of mitochondrial crystallized protein structures and to the development of new computational approaches for protein 3D-modelling and drug design, it is now possible to investigate in detail impaired mitochondrial pathways in CVDs. Furthermore, it is possible to design new powerful drugs able to hit the selected pharmacological targets in a highly selective way to rescue mitochondrial dysfunction and prevent cardiac tissue degeneration. The role of mitochondrial dysfunction in the onset of CVDs appears increasingly evident, as reflected by the impairment of proteins involved in lipid peroxidation, mitochondrial dynamics, respiratory chain complexes, and membrane polarization maintenance in CVD patients. Conversely, little is known about proteins responsible for the cross-talk between mitochondria and cytoplasm in cardiomyocytes. Mitochondrial transporters of the SLC25A family, in particular, are responsible for the translocation of nucleotides (e.g., ATP), amino acids (e.g., aspartate, glutamate, ornithine), organic acids (e.g. malate and 2-oxoglutarate), and other cofactors (e.g., inorganic phosphate, NAD⁺, FAD, carnitine, CoA derivatives) between the mitochondrial and cytosolic compartments. Thus, mitochondrial transporters play a key role in the mitochondria-cytosol cross-talk by leading metabolic pathways such as the malate/aspartate shuttle, the carnitine shuttle, the ATP export from mitochondria, and the regulation of permeability transition pore opening. Since all these pathways are crucial for maintaining healthy cardiomyocytes, mitochondrial carriers emerge as an interesting class of new possible pharmacological targets for CVD treatments.

57 Fe Mössbauer Spectroscopy Characterization of Electrocatalysts

This work addresses the importance of Mössbauer spectroscopy for the characterization of iron-containing electrocatalysts. The most important aspects of electrocatalysis and Mössbauer spectroscopy are summarized. Next, Fe-N-C catalysts and important conclusions made by this technique on preparation, active site identification and degradation are summarized. Furthermore, recent highlights derived for other iron-containing electrocatalysts are summarized.

Cytochrome c oxidase structures suggest a four-state stochastic pump mechanism

A simple stochastic model for a cytochrome c oxidase proton pump.

A random version of principal component analysis in data clustering

Principal component analysis (PCA) is a widespread technique for data analysis that relies on the covariance/correlation matrix of the analyzed data. However, to properly work with high-dimensional data sets, PCA poses severe mathematical constraints on the minimum number of different replicates, or samples, that must be included in the analysis. Generally, improper sampling is due to a small number of data respect to the number of the degrees of freedom that characterize the ensemble. In the field of life sciences it is often important to have an algorithm that can accept poorly dimensioned data sets, including degenerated ones. Here a new random projection algorithm is proposed, in which a random symmetric matrix surrogates the covariance/correlation matrix of PCA, while maintaining the data clustering capacity. We demonstrate that what is important for clustering efficiency of PCA is not the exact form of the covariance/correlation matrix, but simply its symmetry.

The reductive phase of Rhodobacter sphaeroides cytochrome c oxidase disentangled by CO ligation

Cytochrome c oxidase (CcO) is a membrane protein of the respiratory chain that catalytically reduces molecular oxygen (O2) to water while translocating protons across the membrane. The enzyme hosts two copper and two heme iron moieties (heme a/heme a3). The atomic details of the sequential steps that go along with this redox-driven proton translocation are a matter of debate. Particularly for the reductive phase of CcO that precedes oxygen binding experimental data are scarce. Here, we use CcO under anaerobic conditions where carbon monoxide (CO) is bound to heme a3 which in tandem with CuB forms the binuclear center (BNC). Fourier-transform infrared (FTIR) absorption spectroscopy is combined with electro-chemistry to probe different redox and protonation states populated by variation of the external electrostatic potential. With this approach, the redox behavior of heme a and the BNC could be separated and the corresponding redox potentials were determined. We also infer the protonation of one of the propionate side chains of heme a3 to correlate with the oxidation of heme a. Experimental changes in the local electric field surrounding CO bound to heme a3 are determined by their vibrational Stark effect and agree well with electrostatic computations. The comparison of experimental and computational results indicates that changes of the heme a3/CuB redox state are coupled to proton transfer towards heme a3. The latter supports the role of the heme a3 propionate D as proton loading site.

Analysis of the conformations of the HIV-1 protease from a large crystallographic data set

The HIV-1 protease performs essential roles in viral maturation by processing specific cleavage sites in the Gag and Gag-Pol precursor polyproteins to release their mature forms. Here the analysis of a large HIV-1 protease data set (containing 552 dimer structures) are reported. These data are related to article entitled “Conformations of the HIV-1 protease: a crystal structure data set analysis” (Palese, 2017) [1].

Conformations of the HIV-1 protease: A crystal structure data set analysis

The HIV protease is an important drug target for HIV/AIDS therapy, and its structure and function have been extensively investigated. This enzyme performs an essential role in viral maturation by processing specific cleavage sites in the Gag and Gag-Pol precursor polyproteins so as to release their mature forms. This 99 amino acid aspartic protease works as a homodimer, with the active site localized in a central cavity capped by two flexible flap regions. The dimer presents closed or open conformations, which are involved in the substrate binding and release. Here the results of the analysis of a HIV-1 protease data set containing 552 dimer structures are reported. Different dimensionality reduction methods have been used in order to get information from this multidimensional database. Most of the structures in the data set belong to two conformational clusters. An interesting observation that comes from the analysis of these data is that some protease sequences are localized preferentially in specific areas of the conformational landscape of this protein.

Protein States as Symmetry Transitions in the Correlation Matrices

Over the last few years, there has been significant progress in the knowledge on protein folding. However, some aspects of protein folding still need further attention. One of these is the exact relationship between the folded and unfolded states and the differences between them. Whereas the folded state is well known, at least from a structural point of view (just think of the thousands of structures in online databases), the unfolded state is more elusive. Also, these are dynamic states of matter, and this aspect cannot be overlooked. Molecular dynamics-derived correlation matrices are an invaluable source of information on the protein dynamics. Here, bulk eigenvalue spectra of the correlation matrices obtained from the Trp-cage dynamics in the folded and unfolded states have been analyzed. The associated modes represent localized vibrations and are significantly affected by the fine details of the structure and interactions. Therefore, these bulk modes can be used as probes of the protein local dynamics in different states. The results of these analyses show that the correlation matrices describing the folded and unfolded dynamics belong to different symmetry classes. This finding provides new support to the phase-transition models of protein folding.

Random Matrix Theory in molecular dynamics analysis

It is well known that, in some situations, principal component analysis (PCA) carried out on molecular dynamics data results in the appearance of cosine-shaped low index projections. Because this is reminiscent of the results obtained by performing PCA on a multidimensional Brownian dynamics, it has been suggested that short-time protein dynamics is essentially nothing more than a noisy signal. Here we use Random Matrix Theory to analyze a series of short-time molecular dynamics experiments which are specifically designed to be simulations with high cosine content. We use as a model system the protein apoCox17, a mitochondrial copper chaperone. Spectral analysis on correlation matrices allows to easily differentiate random correlations, simply deriving from the finite length of the process, from non-random signals reflecting the intrinsic system properties. Our results clearly show that protein dynamics is not really Brownian also in presence of the cosine-shaped low index projections on principal axes.

Prediction of high- and low-affinity quinol-analogue-binding sites in the aa3 and bo3 terminal oxidases from Bacillus subtilis and Escherichia coli1

Haem-copper oxidases are the terminal enzymes in both prokaryotic and eukaryotic respiratory chains. They catalyse the reduction of dioxygen to water and convert redox energy into a transmembrane electrochemical proton gradient during their catalytic activity. Haem-copper oxidases show substantial structure similarity, but spectroscopic and biochemical analyses indicate that these enzymes contain diverse prosthetic groups and use different substrates (i.e. cytochrome c or quinol). Owing to difficulties in membrane protein crystallization, there are no definitive structural data about the quinol oxidase physiological substrate-binding site(s). In the present paper, we propose an atomic structure model for the menaquinol:O2 oxidoreductase of Bacillus subtilis (QOx.aa3). Furthermore, a multistep computational approach is used to predict residues involved in the menaquinol/menaquinone binding within B. subtilis QOx.aa3 as well as those involved in quinol/quinone binding within Escherichia coli QOx.bo3. Two specific sequence motifs, R70GGXDX4RXQX3PX3FX[D/N/E/Q]X2HYNE97 and G159GSPX2GWX2Y169 (B. subtilis numbering), were highlighted within QOx from Bacillales. Specific residues within the first and the second sequence motif participate in the high- and low-affinity substrate-binding sites respectively. Using comparative analysis, two analogous motifs, R71GFXDX4RXQX8[Y/F]XPPHHYDQ101 and G163EFX3GWX2Y173 (E. coli numbering) were proposed to be involved in Enterobacteriales/Rhodobacterales/Rhodospirillales QOx high- and low-affinity quinol-derivative-binding sites. Results and models are discussed in the context of the literature.

Amyloid beta(1-42) in aqueous environments: effects of ionic strength and E22Q (Dutch) mutation

Development of extracellular plaques characteristic of Alzheimer's disease is related to aggregation of amyloid peptides. The Aβ-42 peptide is the most aggregation prone species, and some missense mutant forms increase this aggregation ability. Due to its poor solubility as monomer in aqueous solutions, Aβ-42 conformational transitions in water have been largely investigated by molecular dynamics. Here we report an all-atom molecular dynamics analysis of the Aβ-42 peptide in aqueous environment using as starting conformation a structure obtained in an isotropic, low-polarity medium, representing a plausible model for the membrane-bound species. While previous studies commonly show that Aβ-42 is largely unstructured in aqueous solution, here we report that this peptide can adopt partially folded structures. Importance of ionic strength has been also investigated, showing that at physiological ionic strength condition a loop stabilizing electrostatic interaction involving Lys28 builds up. In addition, besides stable α-helix structures, we observe the appearance of 310 helix, similar to what was reported experimentally for the Aβ-40 species. The effect of E22Q (Dutch) mutation in high ionic strength condition has been explored. We show that this mutation has a dramatic impact on the Aβ-42 structure. Instead of a partially folded, but extended, conformation obtained with the wild type, the E22Q assumes a two-helix collapsed one due to the clustering of hydrophobic residues.