Molecular modeling of antibodies for the treatment of TNFα-related immunological diseases.

Therapeutic monoclonal antibodies (mAbs) have high efficacy in treating TNF α-related immunological diseases. Other than neutralizing TNF α, these IgG1 antibodies exert Fc receptor-mediated effector functions such as the complement-dependent cytotoxicity (CDC) and antibody-dependent cell cytotoxicity (ADCC). The crystallizable fragment (Fc) of these IgG1 contains a single glycosylation site at Asn 297/300 that is essential for the CDC and ADCC. Glycosylated antibodies lacking core fucosylation showed an improved ADCC. However, no structural data are available concerning the ligand-binding interaction of these mAbs used in TNF α-related diseases and the role of the fucosylation. We therefore used comparative modeling for generating complete 3D mAb models that include the antigen-binding fragment (Fab) portions of infliximab, complexed with TNF α (4G3Y.pdb), the Fc region of the human IGHG1 fucosylated (3SGJ) and afucosylated (3SGK) complexed with the Fc receptor subtype Fcγ RIIIA, and the Fc region of a murine immunoglobulin (1IGT). After few thousand steps of energy minimization on the resulting 3D mAb models, minimized final models were used to quantify interactions occurring between Fcγ RIIIA and the fucosylated/afucosylated Fc fragments. While fucosylation does not affect Fab-TNF α interactions, we found that in the absence of fucosylation the Fc-mAb domain and Fcγ RIIIA are closer and new strong interactions are established between G129 of the receptor and S301 of the Chimera 2 Fc mAb; new polar interactions are also established between the Chimera 2 Fc residues Y299, N300, and S301 and the Fcγ RIIIA residues K128, G129, R130, and R155. These data help to explain the reduced ADCC observed in the fucosylated mAbs suggesting the specific AA residues involved in binding interactions.

Distance-dependent hydrophobic-hydrophobic contacts in protein folding simulations.

Successful prediction of protein folding from an amino acid sequence is a challenge in computational biology. In order to reveal the geometric constraints that drive protein folding, highlight those constraints kept or missed by distinct lattices and for establishing which class of intra- and inter-secondary structure element interactions is the most relevant for the correct folding of proteins, we have calculated inter-alpha carbon distances in a set of 42 crystal structures consisting of mainly helix, sheet or mixed conformations. The inter-alpha carbon distances were also calculated in several lattice "hydrophobic-polar" models built from the same protein set. We found that helix structures are more prone to form "hydrophobic-hydrophobic" contacts than beta-sheet structures. At a distance lower than or equal to 3.8 Å (very short-range interactions), "hydrophobic-hydrophobic" contacts are almost absent in the native structures, while they are frequent in all the analyzed lattice models. At distances in-between 3.8 and 9.5 Å (short-/medium-range interactions), the best performing lattice for reproducing mainly helix structures is the body-centered-cubic lattice. If protein structures contain sheet portions, lattice performances get worse, with few exceptions observed for double-tetrahedral and body-centered-cubic lattices. Finally, we can observe that ab initio protein folding algorithms, i.e. those based on the employment of lattices and Monte Carlo simulated annealings, can be improved simply and effectively by preventing the generation of "hydrophobic-hydrophobic" contacts shorter than 3.8 Å, by monitoring the "hydrophobic-hydrophobic/polar-polar" contact ratio in short-/medium distance ranges and by using preferentially a body-centered-cubic lattice.

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.

The human extended mitochondrial metabolic network: new hubs from lipids.

Even if systems thinking is not new in biology, rationalizing the explosively growing amount of knowledge has been the compelling reason for the sudden rise and spreading of systems biology. Based on 'omics' data, several genome-scale metabolic networks have been reconstructed and validated. One of the most striking aspects of complex metabolic networks is the pervasive power-law appearance of metabolite connectivity. However, the combinatorial diversity of some classes of compounds, such as lipids, has been scarcely considered so far. In this work, a lipid-extended human mitochondrial metabolic network has been built and analyzed. It is shown that, considering combinatorial diversity of lipids and multipurpose enzymes, an intimate connection between membrane lipids and oxidative phosphorilation appears. This finding leads to some biomedical considerations on diseases involving mitochondrial enzymes. Moreover, the lipid-extended network still shows power-law features. Power-law distributions are intrinsic to metabolic network organization and evolution. Hubs in the lipid-extended mitochondrial network strongly suggest that the "RNA world" and the "lipid world" hypothesis are both correct.

Molecular dynamics in cytochrome c oxidase Mössbauer spectra deconvolution.

In this work low temperature molecular dynamics simulations of cytochrome c oxidase are used to predict an experimentally observable, namely Mössbauer spectra width. Predicted lineshapes are used to model Lorentzian doublets, with which published cytochrome c oxidase Mössbauer spectra were simulated. Molecular dynamics imposed constraints to spectral lineshapes permit to obtain useful information, like the presence of multiple chemical species in the binuclear center of cytochrome c oxidase. Moreover, a benchmark of quality for molecular dynamic simulations can be obtained. Despite the overwhelming importance of dynamics in electron-proton transfer systems, limited work has been devoted to unravel how much realistic are molecular dynamics simulations results. In this work, molecular dynamics based predictions are found to be in good agreement with published experimental spectra, showing that we can confidently rely on actual simulations. Molecular dynamics based deconvolution of Mössbauer spectra will lead to a renewed interest for application of this approach in bioenergetics.