Molecular dynamics simulation studies and in vitro site directed mutagenesis of avian beta-defensin Apl_AvBD2
© Soman et al. 2010
Published: 18 January 2010
Skip to main content
© Soman et al. 2010
Published: 18 January 2010
Defensins comprise a group of antimicrobial peptides, widely recognized as important elements of the innate immune system in both animals and plants. Cationicity, rather than the secondary structure, is believed to be the major factor defining the antimicrobial activity of defensins. To test this hypothesis and to improve the activity of the newly identified avian β-defensin Apl_AvBD2 by enhancing the cationicity, we performed in silico site directed mutagenesis, keeping the predicted secondary structure intact. Molecular dynamics (MD) simulation studies were done to predict the activity. Mutant proteins were made by in vitro site directed mutagenesis and recombinant protein expression, and tested for antimicrobial activity to confirm the results obtained in MD simulation analysis.
MD simulation revealed subtle, but critical, structural variations between the wild type Apl_AvBD2 and the more cationic in silico mutants, which were not detected in the initial structural prediction by homology modelling. The C-terminal cationic 'claw' region, important in antimicrobial activity, which was intact in the wild type, showed changes in shape and orientation in all the mutant peptides. Mutant peptides also showed increased solvent accessible surface area and more number of hydrogen bonds with the surrounding water molecules. In functional studies, the Escherichia coli expressed, purified recombinant mutant proteins showed total loss of antimicrobial activity compared to the wild type protein.
The study revealed that cationicity alone is not the determining factor in the microbicidal activity of antimicrobial peptides. Factors affecting the molecular dynamics such as hydrophobicity, electrostatic interactions and the potential for oligomerization may also play fundamental roles. It points to the usefulness of MD simulation studies in successful engineering of antimicrobial peptides for improved activity and other desirable functions.
Defensins are recognized as important elements of the innate immune system in almost all living beings . The most distinct molecular feature of defensins is their high pI value, ranging from +6 to +12 as monomers, manifested by abundant arginine and lysine residues in their sequences . They kill microorganisms through permeabilization of the microbial membrane composed of negatively charged components such as phospholipids, teichoic acids and lipopolysaccharides . It is believed that electrostatic interactions dictate not only the uptake of cationic defensins across the bacterial cell wall but also their ability to permeabilize the cytoplasmic membrane and to induce leakage of cellular contents .
Many previous studies hypothesized that the major factor defining the antimicrobial activity of defensin is its cationicity or its isoelectric point (pI) value, than its secondary structure . The present study attempted to verify this hypothesis and generate derivatives of a beta-defensin with more microbicidal activity by engineering the protein. Apl_AvBD2 is a β-defensin homologue identified from domestic duck, and was found to exhibit antibacterial and immunomodulatory properties [6, 7]. In the present study, we made in silico mutants of Apl_AvBD2 with higher pI values than the wild type, and analyzed them by Molecular Dynamic (MD) simulation analysis to find the structure-function relationship. Subsequently, mutated recombinant proteins were made in vitro and were evaluated for antibacterial activity to confirm the observations from the computational studies. Our results indicated that subtle structural differences in critical areas of the molecule can drastically alter the antibacterial potential of β-defensin molecules.
The predicted amino acid sequence of Apl_AvBD2 [GenBank: AY641439] was subjected to a homology search using BLAST ( and PSI-BLAST  against NCBI PDB database. The top hits were aligned against the query sequence in a multiple sequence alignment using Clustal W . We chose the top scoring sequence Apa_AvBD2 (Spheniscin 2: showed 35% sequence identity), originating from King Penguin as the template for further study and the PDB co-ordinate ((PDB id: 1ut3) was retrieved from Protein Data Bank. Homology modelling was carried out using MODELLER package  based on the sequence alignment generated between template and target sequences. The atomic coordinates were obtained from the template structures to model Apl_AvBD2. Care was taken to make the coordination geometry of side chain atoms most favorable. Conformations of a few residues were therefore adjusted using loop refinement programme within MODELLER package. Energy minimization of the top scored model was carried out with GROMACS 3.3.1 (The Groningen Machine for Chemical Simulations)  using OPLSAA force field. The minimization was set to run for 5000 steps or until convergence to machine precision. PROSA2003  program was used for validation of the model, by analyzing residue interaction energy and z-score. These procedures were iterated several times until a good quality model was obtained.
Molecular dynamics (MD) simulations were performed using the GROMACS version 3.3.1 and OPLSAA force field. The initial structures were solvated with three-point transferable intermolecular potential (TIP3P) water molecules  and appropriate number of Chloride ions in a rectangular box to neutralize the system; the box dimensions ensured that any protein atom was at least 8 Å away from the wall of the box. After energy minimization, MD simulations were performed for 100 ps at constant temperature (300 K) and pressure (1 atm) with periodic boundary conditions, particle-mesh Ewald summation, and a 1-fs time step to heat and equilibrate the system. This was followed by production runs of 10 ns duration for each simulation. Structures were saved every 10 ps for analysis. The relative binding energies computed using the tool g_energy module of 'GROMACS 3.3.1, employing molecular mechanics and a continuum solvent model. The output files (.xvg) from the GROMACS 3.3.1 was analysed in XMGRACE  software.
The parameters analysed were: area per atom, area per residue, energy variations (kinetic energy, potential energy and total energy), van der Waals interactions, intra-molecular hydrogen bonds, inter-molecular hydrogen bonds with surrounding water molecules, radius of gyration, root mean square deviations (RMSD) of each amino acid residue, root mean square fluctuations (RMSF) of Cα atoms of each amino acid residue, solvent accessible surface area, and hydrophobic and hydrophilic interactions. Pressure, volume, temperature, minimum distance to periodic image and maximum internal distance of the MD simulation system were examined. The obtained parameters for wild type and mutant peptides were compared.
Mutagenesis reactions were carried out using the modified mutagenic primers as per the protocol by Quick Change Mutagenesis Kit (Stratagene, La Jolla, CA). Apl_AvBD2-pET-32 gene-vector construct was used for this purpose. The cloning strategy and recombinant expression protocols are described elsewhere . Selected clones were sequenced directly using an ABI 3730 Genetic Analyzer automated DNA sequencer (PE Applied Biosystems, Foster City, CA) for confirmation of mutation. Recombinant proteins of single amino acid mutants and serial progressive mutants of Apl_AvBD2, which were selected for the MD analysis (Figure 1), were made in vitro. These would make variants of the protein with pI value ranging from 8.69 to 9.97. These proteins were expressed in BL21DE3 pLysS bacterial cells, purified and used for antibacterial assay .
Apl_AvBD2 was modelled using the known three-dimensional structure of the Apa_AvBD2 as the template for homology modelling. Superimposition of the Cα atoms of the template and target protein gave a calculated RMSD of 0.39 Å. Analysis of the secondary structure revealed the presence of mainly two β-strands in Apl_AvBD2. PROSA2003 Z-scores pointed to the compatibility between the model and the template. The Z-scores of both the structures were more or less similar. PROSA2003 quality check showed that the model of the Apl_AvBD2 was of good quality, wherein the interaction energy of each residue with the residual of the protein was negative. The Apl_AvBD2 model had a Z-score of -5.55 compared to that of the template (-5.83). Inspection of PROSA2003 plots revealed no region of the model with positive PROSA2003 energies (data not shown).
In the MD simulation experiments, we selected dimeric state of the Apl_AvBD2 and its mutants for analysis. This was because several previous studies have shown that the β-defensins function in the form of dimers, which are their most stable form . We also observed that in MD simulation analysis using the monomer form of wild type Apl_AvBD2, the system was not stable after 5 ns, confirming the monomer's instability (data not shown).
From the results obtained in this study, it can be concluded that the increase in cationicity alone may not enhance the antibacterial activity of defensins. The antimicrobial activity of these peptides requires a balance between its cationicity and hydrophobicity. The substitution of hydrophobic residues with more cationic hydrophilic residues leads to complete loss of activity. Moreover, even a single amino acid change can cause deleterious effect in the antibacterial activity. Most of the naturally occurring antimicrobial peptides have undergone evolutionary selection to perform best in the host. Manipulations to enhance its effects may not always yield viable results and has to be done prudently. MD simulation experiments would be useful while manipulating antimicrobial peptides to improve its activity by mutagenesis or by chemical synthesis.
The authors acknowledge the Department of Science and Technology, Government of India for financial assistance to this project. Soja Saghar Soman was supported with Senior Research Fellowship from Council of Scientific and Industrial Research (CSIR), Government of India. Bioinformatics facility at RGCB is supported by Department of Biotechnology, Government of India.
This article has been published as part of BMC Bioinformatics Volume 11 Supplement 1, 2010: Selected articles from the Eighth Asia-Pacific Bioinformatics Conference (APBC 2010). The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2105/11?issue=S1.
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.