An in silico study of the citrus dioxygenases CCD4 family substrates.

The coloration of Citrus fruits is related with the concentration of carotenoids, isoprenoid pigments of 40 carbon atoms (C40). Rodrigo et al. and Ma et al. reported a CCD4-type citrus dioxygenase responsible for the generation of C30 apocarotenoids providing a reddish-orange pigmentation to the peel of many mandarins and oranges. Among them, CCD4b was the first case described of a dioxygenase that cleaves carotenoids C40 in the double bond 7', 8' or 7, 8, generating ß-citraurin or 8-ß-apocarotenal. Here we report the three-dimensional structures of CCD4a and CCD4b, modeled by sequence homology (2BIW) and validated by molecular dynamics (MD). Docking calculations were performed in CCD4a and CCD4b structures with thousands of rotated initial carotenoid conformations and all the possible poses in the active site were found. The interaction energy was measured by means of ASE scoring, Amber99 refinement and London ΔG rescoring. For the case of CCD4a model, the results showed London ΔG score of -19, -17 and -15 kcal/mol for zeaxanthin, ß-cryptoxanthin and ß-carotene, respectively. The same sequence in the estimated interaction strength for the three ligands was obtained using MD. The interaction energy of CCD4b indicated that, in agreement with experimental data, zeaxanthin and ß-cryptoxanthin could be cleaved by the enzyme, ß- and α-carotene have chances to be oxidized and lycopene has not good interaction energy to be predicted as substrate. These findings will be discussed considering the potential in vivo substrates and products, and the physiological role in Citrus fruits. Communicated by Ramaswamy H. Sarma.

In silico studies of Echinococcus granulosus FABPs.

Fatty acid (FA) binding proteins are small intracellular proteins whose members exhibit great diversity and low similarity at the primary structure level, but a highly conserved three-dimensional structure. Characterised by a high-affinity non-covalent binding of hydrophobic ligands, these proteins have a molecular mass of 14-15 kDa with a characteristic ß-barrel structure. Members of this family have been identified along the zoological scale, with Platyhelminthes being the more primitive organisms where they have been reported. Two FA binding proteins (FABPs), EgFABP1 and EgFABP2, with 88% similarity have been identified in Echinococcus granulosus. In an effort to understand why two such similar proteins are expressed by this organism, we performed an in silico analysis of the binding capabilities of both proteins. The crystallographic structure of EgFABP1 was utilised as a template to model EgFABP2, and both were docked against palmitate, oleate, linoleate and arachidonate. The docked structures were submitted to 4 ns molecular dynamics simulations, and their protein-ligand interaction energies were measured. The collected data demonstrated that linoleate and arachidonate had the higher interaction energies when bound to EgFABP1 and that palmitate and linoleate had the higher interaction energies when bound to EgFABP2. External and internal binding surfaces were analysed, showing differences at both levels. Internal surface compositions suggested that both proteins could have preferences for certain FAs. Comparisons of the holo and apo forms of each protein indicated that the ligand imposed subtle, but specific modifications that could trigger surface signals. The differences found between the proteins under study suggest that they could have functional uniqueness in the parasite's metabolism.

In silico characterization of cytisinoids docked into an acetylcholine binding protein.

Homology models of nicotinic acetylcholine receptors (nAChRs) suggest that subtype specificity is due to non-conserved residues in the complementary subunit of the ligand-binding pocket. Cytisine and its derivatives generally show a strong preference for heteromeric alpha4beta2* nAChRs over the homomeric alpha7 subtype, and the structural modifications studied do not cause large changes in their nAChR subtype selectivity. In the present work we docked cytisine, N-methylcytisine, and several pyridone ring-substituted cytisinoids into the crystallographic structure of the Lymnaea stagnalis acetylcholine binding protein (AChBP) co-crystallized with nicotine (1UW6). The graphical analysis of the best poses showed that cytisinoids have weak interactions with the side chains of the non-conserved amino acids in the complementary subunit justifying the use of PDB 1UWB as a surrogate for nAChR. Furthermore, we found a high correlation (R(2)=0.96) between the experimental pIC(50) values at alpha4beta2* nAChR and docking energy (S) of the best cytisinoid poses within the AChBP. Due to the quality of the correlation we suggest that this equation might be used as a predictive model to propose new cytisine-derived nAChRs ligands. Our docking results also suggest that further structural modifications of these cytisinoids will not greatly alter their alpha4beta2*/alpha7 selectivity.