Neration of the nearUV CD spectrum by means of interactions between the transitions of the aromatic chromophores; evaluating the impact of the protein BIBS39 biological activity conformational flexibility on the quality of the calculated spectra; exploring the sensitivity of chromophore interactions identified in the near-UV to the effect of the protein conformational dynamics; computing the effects of tryptophan mutations on the CD spectra in correlation with the experimental ones; evaluating the applicability of restricted structural model including only the tryptophan and tyrosine chromophores at both semiempirical level (using the matrix method) and Time-Dependent Density Functional Theory (TDDFT);ii) iii)iv) v)This study is focused mainly on the aromatic contributions (Lb and La transitions) in the near-UV CD. Indeed the higher energy aromatic transitions (Bb and Ba) might contribute sensitively to the far-UV [3,10] where they mix with a huge number of peptide transitions. The analysis of the interactions would be therefore complicated and is not present here.MethodsThree levels of modelling methods were carried out in the study of HCAII CD spectral features: i) Atomistic Molecular Dynamics (MD) simulations [13,14]; ii) Approximate Quantum Mechanical CD calculations using the Matrix Salmon calcitonin biological activity method [15] and iii) Time Dependent Density Functional Theory (TDDFT) calculations [16]. Tryptophan mutant structures were prepared by in silico mutagenesis from the crystal structure of the wild-type of HCAII taken from Protein Data Bank (Berman and others 2000) (PDB ID code 2cba) (Hakansson and others 1992), and structural snapshots of the wild-type protein and tryptophan mutant forms were taken from MD simulations. The CD calculations with the matrix method were performed incorporating all peptides and side chain chromophores. The matrix method calculations were performed using the Dichrocalc web interface [17]. This method [15] in its origin-independent form [18] considers the protein as a system of M independent chromophoric groups. The wave function of the entire molecule is represented as a linear superposition of basis functions. Every basis function is a product of all monomer wave functions where only one group is in an excited state. This way the matrix method incorporates all mechanisms of generation of the rotational strengths (m-m, m-m and the static field effect). The interactions between the chromophores are considered to be purely electrostatic and therefore the permanent and transition electron densities (represented 1326631 by monopoles) are implemented from electronic structure calculations on model systems. Finally, the Hamiltonian matrix is diagonalized by unitary transformation in order to represent the excited states in the interacting system. More details about the matrix method can be found in [5,19,20]. The monopoles for the side chain chromophores (including the aromatic ones) are taken from ab initio calculations [21] and the monopoles for the peptide chromophores are taken from ab intio [22] and semi-empirical calculations [23]. TDDFT calculations were done with Gaussian09 code [24] and to the best of our knowledge represent one of the largest biomolecular TDDFT calculations. The system included only 3methylindole parts from the side chains of the tryptophans and the phenol parts from the side chains for the tyrosines kept at theirFigure 1. Structure of HCAII. The tryptophan chromophores are shown in blue licorice. doi:10.1371/journal.pone.0056874.gConformat.Neration of the nearUV CD spectrum by means of interactions between the transitions of the aromatic chromophores; evaluating the impact of the protein conformational flexibility on the quality of the calculated spectra; exploring the sensitivity of chromophore interactions identified in the near-UV to the effect of the protein conformational dynamics; computing the effects of tryptophan mutations on the CD spectra in correlation with the experimental ones; evaluating the applicability of restricted structural model including only the tryptophan and tyrosine chromophores at both semiempirical level (using the matrix method) and Time-Dependent Density Functional Theory (TDDFT);ii) iii)iv) v)This study is focused mainly on the aromatic contributions (Lb and La transitions) in the near-UV CD. Indeed the higher energy aromatic transitions (Bb and Ba) might contribute sensitively to the far-UV [3,10] where they mix with a huge number of peptide transitions. The analysis of the interactions would be therefore complicated and is not present here.MethodsThree levels of modelling methods were carried out in the study of HCAII CD spectral features: i) Atomistic Molecular Dynamics (MD) simulations [13,14]; ii) Approximate Quantum Mechanical CD calculations using the Matrix Method [15] and iii) Time Dependent Density Functional Theory (TDDFT) calculations [16]. Tryptophan mutant structures were prepared by in silico mutagenesis from the crystal structure of the wild-type of HCAII taken from Protein Data Bank (Berman and others 2000) (PDB ID code 2cba) (Hakansson and others 1992), and structural snapshots of the wild-type protein and tryptophan mutant forms were taken from MD simulations. The CD calculations with the matrix method were performed incorporating all peptides and side chain chromophores. The matrix method calculations were performed using the Dichrocalc web interface [17]. This method [15] in its origin-independent form [18] considers the protein as a system of M independent chromophoric groups. The wave function of the entire molecule is represented as a linear superposition of basis functions. Every basis function is a product of all monomer wave functions where only one group is in an excited state. This way the matrix method incorporates all mechanisms of generation of the rotational strengths (m-m, m-m and the static field effect). The interactions between the chromophores are considered to be purely electrostatic and therefore the permanent and transition electron densities (represented 1326631 by monopoles) are implemented from electronic structure calculations on model systems. Finally, the Hamiltonian matrix is diagonalized by unitary transformation in order to represent the excited states in the interacting system. More details about the matrix method can be found in [5,19,20]. The monopoles for the side chain chromophores (including the aromatic ones) are taken from ab initio calculations [21] and the monopoles for the peptide chromophores are taken from ab intio [22] and semi-empirical calculations [23]. TDDFT calculations were done with Gaussian09 code [24] and to the best of our knowledge represent one of the largest biomolecular TDDFT calculations. The system included only 3methylindole parts from the side chains of the tryptophans and the phenol parts from the side chains for the tyrosines kept at theirFigure 1. Structure of HCAII. The tryptophan chromophores are shown in blue licorice. doi:10.1371/journal.pone.0056874.gConformat.