Next Article in Journal
Structural DNA Nanotechnology: From Design to Applications
Next Article in Special Issue
Spatial Simulations in Systems Biology: From Molecules to Cells
Previous Article in Journal
Structural Analysis of Hypothetical Proteins from Helicobacter pylori: An Approach to Estimate Functions of Unknown or Hypothetical Proteins
Previous Article in Special Issue
Pharmacological Classification and Activity Evaluation of Furan and Thiophene Amide Derivatives Applying Semi-Empirical ab initio Molecular Modeling Methods
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Target Molecular Simulations of RecA Family Protein Filaments

1
Department of Information Management, Chia Nan University of Pharmacy & Science, No. 60, Sec. 1, Erren Rd., Rende Dist., Tainan City 71710, Taiwan
2
National Center for High-Performance Computing, Hsin-Shi, No. 28, Nan-Ke 3rd Rd., Hsin-Shi Dist., Tainan City 74147, Taiwan
3
Department of Physics, National Chung Hsing University, No. 250, Kuo Kuang Rd., Taichung 402, Taiwan
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2012, 13(6), 7138-7148; https://doi.org/10.3390/ijms13067138
Submission received: 11 April 2012 / Revised: 10 May 2012 / Accepted: 29 May 2012 / Published: 11 June 2012
(This article belongs to the Special Issue Advances in Biomolecular Simulation)

Abstract

:
Modeling of the RadA family mechanism is crucial to understanding the DNA SOS repair process. In a 2007 report, the archaeal RadA proteins function as rotary motors (linker region: I71-K88) such as shown in Figure 1. Molecular simulations approaches help to shed further light onto this phenomenon. We find 11 rotary residues (R72, T75-K81, M84, V86 and K87) and five zero rotary residues (I71, K74, E82, R83 and K88) in the simulations. Inclusion of our simulations may help to understand the RadA family mechanism.

1. Introduction

Homologous recombination is a process whereby two DNA duplexes interact to transfer genetic information, then create new genetic linkages and rearrange DNA segments. Genes involved in homologous recombination are important for regulation of gene expression and DNA repair [1,2]. In the recombination process, RecA-like proteins can bring two homologous DNA molecules together and exchange the DNA strands [24]. Examples of RecA-like proteins include archaeal RadA, bacterial RecA, Rad51, meiosis-specific Dmc1 and eukaryotic Rad51. At DNA double strand (dsDNA) break sites, these proteins interact with single-stranded DNA (ssDNA) and form a right-handed helical nucleoprotein filament (presynaptic complex) [5]. The presynaptic complex has both ATPase and DNA strand exchange activities. In the presence of ATP molecules, the DNA strand exchange activity ensures the formation of heteroduplex DNA (hDNA) between ssDNA and its complementary strand in the double-stranded DNA (dsDNA). The molecular mechanism has been proposed such that ATP hydrolysis promotes product release and a rotary mechanism that can solve the DNA topological problem in the strand exchange reaction [68]. In terms of structural functionalism, RecA proteins have similar N-terminal domains (NTD) that may interact with dsDNA, and the C-terminal domain that may hydrolyze ATP molecules and form the new conformation [911]. RecA proteins may form closed rings, as well as left-handed protein filaments and right-handed helical filaments [1217]. The molecular mechanism underlying the conformational flexibility of RecA proteins is still unclear. In a 2007 report [17], Wang et al. proposed that: (1) right-handed RadA proteins will interact with ssDNA and form presynaptic complexes; (2) ATPase of presynaptic complexes decomposes ATP molecules (ATP→ADP). Then the hydrolysis energy can force the DNA strand exchange; (3) the complexes might release the ADP molecules and prepare to dissociate the DNA strands; and (4) after dissociating the DNA strands, left-handed RadA proteins might be formed. The RadA proteins might function as rotary motors (linker region: I71-K88) such as the mechanism shown in Figure 1.
MD techniques can offer a convenient alternative to experimental approaches because they can treat a single macromolecule at an atomic level. TMD (target molecular dynamics) [18] methods can provide calculated reaction paths for most proteins by continuously decreasing the target values. This method can predict reaction paths of ras p21 proteins [19] and chymotrypsin inhibitor 2 proteins [20]. A TIP3P solvent model allows all of the archaeal RadA protein’s rotation trajectories to be sampled.
The cumulative changes in the backbone dihedral angles (CCDA) method can predict important residues in biomolecular dissociation systems [2126]. The backbone dihedral angles of proteins are called phi (ϕ, involving the backbone atoms C′-N-Cα-C′) and psi (ψ, involving the backbone atoms N-Cα-C′-N). Thus, phi controls the C′-C′ distance and psi controls the N-N distance. CCDA is defined as:
C C D A ( α ) = Σ ( α ) j
Here, α is the phi or psi angle, and j is residue number. In biomolecular systems, backbone rotations are provided with higher energy barriers, and backbone dihedral angles are representative of backbone rotations [27,28]. Thus, counting the cumulative changes in backbone dihedral angles can predict important residues of the Abs-Ag complex.
In the present study, we used the Sulfolobus solfataricus (Sso) protein sequence and related 3D structures (closed-ring RadA, MvRadA, right-hand RadA and left-hand RadA) [12,17,29,30] to normalize the 3D archaeal RadA structures. Then the TMD method was used to investigate the conformational mechanism of the archaeal RadA proteins. We also reveal the conformational mechanism, the CCDAs (cumulate changed dihedral angles) of rotational residues (I71-K88), the system’s potential energy, and RMSD of archaeal RadA proteins in the simulation process.

2. Results and Discussion

2.1. Homology Models Construction and Evaluation

From the 2007 report [17], the closed-ring, MvRadA and right-hand RadA proteins can be used as potential templates for homology models. For closed-ring, MvRadA and right-hand RadA proteins, the percentage of sequences identifying with the left-hand RadA protein sequence are 51.6%, 44.0% and 100% (Figure 2). From Prosa method validations, the z-score of the four proteins are −8.91, 0.09, −2.47 and −7.86 for closed-ring, MyRadA, right-hand RadA and left-hand RadA, respectively. These initial structures are a good starting point for building reliable models in the next step.
The four initial proteins were refined by energy minimization and MD simulations. Figure 3 displays the total energy versus frames collected per 1ps during the entire 2 ns MD simulation. Clearly, the total energies of four systems were equilibrated after 0.3 ns. The structures were set as the final models by 2 ns simulations. Final models of the four proteins were assessed by PROSA and PROCHECK. By Prosa method validations, the z-scores of the four proteins are −8.95, 0.01, −2.53 and −7.97, respectively. As to the assessment by PROCHECK, the reliability backbone torsion angles of the four RadA filaments were examined. In the core Ramachandran region, the percentage of dihedral angles is 100.0%, 98.8%, 98.6% and 100.0% for left-hand, closed-ring, MvRadA and right-hand RadA proteins, respectively. The data indicate that these 3D models are reliable for further TMD simulations.

2.2. TMD Conformational Calculations

The potential energy of our left-hand filament is reset to zero, with the energy profiles of the TMD calculations shown in Figure 4. The four potential energy states are 0.00, 254.14, 223.54 and −10.16 kcal/mol. The results in Figure 4 indicate that two energy barriers (0–28 ns and 28–33 ns) occur in the 40 ns simulations processes. Analyzing the energy profile, we find that the major energy barrier occurs during the conformational change between left-hand and pre-right-hand proteins. The major energy barrier height is approximately 280.00 kcal/mol. The minor energy barrier exists as the conformation turns into the right-hand protein. The energy barrier height is approximately 60.00 kcal/mol. The TMD simulation trajectories of the calculations were traced by RMSD and CCDAs methods. Figure 5 shows the 40 ns TMD simulation profiles of RMSD. Due to the simulation’s initial structure of RadA filament (left-hand), the profile of the left-hand RadA filament indicates that the RMSD will increase 27 Å within 10 n, and the value will fluctuate between 25 and 35 Å within 10–30 ns, with the value decreasing to 0 Å at 30–40 ns. The other RMSD profile shows the RMSD values for closed-ring, MvRadA and right-hand RadA filaments.
Figure 6 and Table 1 show the CCDAs of the rotational residues (I71-K88). The results indicate that the CCDAs of the five residues (I71, K74, E82, R83 and K88) are zero in the simulations.
Analyzing the major barrier (0–28 ns), we find that the 11 variations (R72, T75-K81, M84, V86 and K87) are obvious in the CCDAs. At 0–10 ns, the system tends to fold the closed-ring RadA (Figure 5) and the system energy rises to 254.14 kcal/mol. The K81 is the obvious CCDA variation and the event occurs at 6 ns. At 10–20 ns, the system conformation exists between closed-ring and MvRadA proteins. The system energy decreases to 223.54 kcal/mol and three variations (K80, K81 and M84) are obvious in the CCDAs. The K80 variation occurs at 12 ns and the other variations occur at 17 ns. At 20–28 ns, the system energy decreases to −70 kcal/mol and four variations (R72, T75, A76 and E78) emerge at 25 ns in the CCDAs. Analyzing the minor barrier (28–33 ns), we find that the obvious variations of CCDAs are the same as for the major energy barrier analysis.

3. Experimental Section

3.1. Homology Modeling and Refinement of the Models

Homology modeling is a computational approach for three-dimensional protein structure modeling and prediction. Proteins whose structures are still uncharacterized can be built using homology modeling. This method builds an atomic model based on experimentally determined known structures that have sequence homology of more than 40% with the target. Modeling structures with less than 40% template similarity would result in a less reliable model. Homology modeling is also known as comparative modeling. To normalize the 3D archaeal RadA structures, the homology modeling method was applied in the study. All the primary protein sequences and structures can be obtained from the protein data bank [12,17,29,30]. For the four target sequences (closed-ring RadA, MvRadA, right-hand RadA and left-hand RadA), the multiple sequence alignment was calculated by ClustalX2 program [31]. According to the left-hand RadA protein sequence and the template structures (PDB ID: 1PZN and 1T4G), the closed-ring RadA and MvRadA structures were constructed by the Modeller program [32]. The disappearing L2 region [33] of the right-hand RadA (M258-H275) was constructed by the loop optimization method and very slow loop refinements method found in the Modeller program.
The structural refinement process was accomplished in several stages. First, energy minimization of 10,000 steps steepest decent (SD) followed by 30,000 steps conjugated gradient (CG) was carried out while fixing the backbone Cα. Then, we generated the solvent box (TIP3 water molecules and RadA protein) of volume 89.27 × 81.69 × 88.17 Å3. The initial simulation system combined the RadA protein and the solvent box. The water molecules around the RadA protein with 2 Å were deleted. The total number of water molecules was 6104 and the total number of atoms in the system was 53,081. After performing 20,000 step CG energy minimization, MD simulations were performed with 2 ns using an NVT ensemble, periodic boundary conditions box (89.27 × 81.69 × 88.17 Å3), and particle mesh Ewald (PME) method at a temperature of 310 K. The energy minimization and MD simulations mentioned above were accomplished by the Amber program and force fields [34]. The final structures were checked by PROCHECK [35] and PROSA [36].

3.2. TMD Simulations

An additional energy term based on the RMSD of the RadA proteins is relative to a prescribed target structure. The energy term (UTMD) is defined as the formula:
U T M D = 1 2 K N ( R M S D ( t ) - R M S D 0 ( t ) ) 2
where the force constant, K, is 200 kcal·mol−1·Å−2. N is the number of target atoms. RMSD(t) is the RMSD of the simulation structure at time t relative to the prescribed target structure, and RMSD0(t) is the prescribed target RMSD value at time t. As described above, four TMD simulation cases were simulated in the study. The four cases were left-hand/close-ring, close-ring/MvRadA, MvRadA/right-hand and right-hand/left-hand RadA proteins. The target proteins were close-ring, MvRadA, right-hand and left-hand RadA protein, respectively. The TMD simulations were calculated with the NAMD [37] software and 10 ns simulation time.

4. Conclusions

In this article, we propose using molecular simulations techniques to analyze the conformation mechanism of RadA proteins. We report the two energy barriers, 11 high rotational residues and 5 zero rotational residues (I71, K74, E82, R83 and K88) in the TMD simulations. If the initial conformation starts as the left-hand RadA protein, the system should overcome the two energy barriers and the 11 residues (R72, T75-K81, M84, V86 and K87) to provide high conformation degree of freedom. The high rotational residues might increase the flexibility of RadA and make the conformation change easy. The zero rotational residues are the axis of the conformation change mechanism. In the next phase of the project, we will use molecular biology methods to mutate the residues included 11 high rotational residues and 5 zero rotational residues. We hope this finding might help to understand the RadA family mechanism.

Acknowledgment

This work was supported by the National Center for High-Performance Computing and National Science Council (NSC97-2221-E-041-018-MY3 and NSC 100-2113-M-492-001).

References

  1. Seitz, E.M.; Brockman, J.P.; Sandler, S.J.; Clark, A.J.; Kowalczykowski, S.C. RadA protein is an archaeal RecA protein homolog that catalyzes DNA strand exchange. Genes Dev 1998, 12, 1248–1253. [Google Scholar]
  2. Kowalczykowski, S.C.; Dixon, D.A.; Eggleston, A.K.; Lauder, S.D.; Rehrauer, W.M. Biochemistry of homologous recombination in Escherichia coli. Microbiol. Rev 1994, 58, 401–465. [Google Scholar]
  3. Roca, A.I.; Cox, M.M. RecA protein: Structure, function, and role in recombinational DNA repair. Prog. Nucleic Acid Res. Mol. Biol 1997, 56, 129–223. [Google Scholar]
  4. DiRuggiero, J.; Santangelo, N.; Nackerdien, Z.; Ravel, J.; Robb, F.T. Repair of extensive ionizing-radiation DNA damage at 95 degrees C in the hyperthermophilic archaeon Pyrococcus furiosus. J. Bacteriol 1997, 179, 4643–4645. [Google Scholar]
  5. Kil, Y.V.; Glazunov, E.A.; Lanzov, V.A. Characteristic thermodependence of the RadA recombinase from the hyperthermophilic archaeon Desulfurococcus amylolyticus. J. Bacteriol 2005, 187, 2555–2557. [Google Scholar]
  6. Kowalczykowski, S.C.; Eggleston, A.K. Homologous pairing and DNA strand-exchange proteins. Annu. Rev. Biochem 1994, 63, 991–1043. [Google Scholar]
  7. West, S.C. Molecular views of recombination proteins and their control. Nat. Rev. Mol. Cell Biol 2003, 4, 435–445. [Google Scholar]
  8. Cox, M.M. The bacterial reca protein as a motor protein. Annu. Rev. Microbiol 2003, 57, 551–577. [Google Scholar]
  9. Aihara, H.; Ito, Y.; Kurumizaka, H.; Yokoyama, S.; Shibata, T. The N-terminal domain of the human Rad51 protein binds DNA: Structure and a DNA binding surface as revealed by NMR. J. Mol. Biol 1999, 290, 495–504. [Google Scholar]
  10. Kinebuchi, T.; Kagawa, W.; Kurumizaka, H.; Yokoyama, S. Role of the N-terminal domain of the human DMC1 protein in octamer formation and DNA binding. J. Biol. Chem 2005, 280, 28382–28387. [Google Scholar]
  11. Eggler, A.L.; Lusetti, S.L.; Cox, M.M. The C terminus of the Escherichia coli RecA protein modulates the DNA binding competition with single-stranded DNA-binding protein. J. Biol. Chem 2003, 278, 16389–16396. [Google Scholar]
  12. Wu, Y.; He, Y.; Moya, I.A.; Qian, X.; Luo, Y. Crystal structure of archaeal recombinase RadA: A snapshot of its extended conformation. Mol. Cell 2004, 15, 423–435. [Google Scholar]
  13. Ariza, A.; Richard, D.J.; White, M.F.; Bond, C.S. Conformational flexibility revealed by the crystal structure of a crenarchaeal RadA. Nucleic Acids Res 2005, 33, 1465–1473. [Google Scholar]
  14. Yang, S.; Yu, X.; Seitz, E.M.; Kowalczykowski, S.C.; Egelman, E.H. Archaeal RadA protein binds DNA as both helical filaments and octameric rings. J. Mol. Biol 2001, 314, 1077–1085. [Google Scholar]
  15. McIlwraith, M.J.; Hall, D.R.; Stasiak, A.Z.; Stasiak, A.; Wigley, D.B.; West, S.C. RadA protein from Archaeoglobus fulgidus forms rings, nucleoprotein filaments and catalyses homologous recombination. Nucleic Acids Res 2001, 29, 4509–4517. [Google Scholar]
  16. Lee, M.-H.; Leng, C.-H.; Chang, Y.-C.; Chou, C.-C.; Chen, Y.-K.; Hsu, F.-F.; Chang, C.-S.; Wang, A.H.J.; Wang, T.-F. Self-polymerization of archaeal RadA protein into long and fine helical filaments. Biochem. Biophys. Res. Commun 2004, 323, 845–851. [Google Scholar]
  17. Chen, L.-T.; Ko, T.-P.; Chang, Y.-C.; Lin, K.-A.; Chang, C.-S.; Wang, A.H.J.; Wang, T.-F. Crystal structure of the left-handed archaeal RadA helical filament: Identification of a functional motif for controlling quaternary structures and enzymatic functions of RecA family proteins. Nucleic Acids Res 2007, 35, 1787–1801. [Google Scholar]
  18. Schlitter, J.; Engels, M.; Krüger, P.; Jacoby, E.; Wollmer, A. Targeted molecular dynamics simulation of conformational change: Application to the T↔R transition in insulin. Mol. Simul 1993, 10, 291–308. [Google Scholar]
  19. Ma, J.; Karplus, M. Molecular switch in signal transduction: Reaction paths of the conformational changes in ras p21. Proc. Natl. Acad. Sci. USA 1997, 94, 11905–11910. [Google Scholar]
  20. Ferrara, P.; Apostolakis, J.; Caflisch, A. Targeted molecular dynamics simulations of protein unfolding. J. Phys. Chem 2000, 104, 4511–4518. [Google Scholar]
  21. Wang, Y.-T.; Liao, J.-M.; Chen, C.-L.; Su, Z.-Y.; Chen, C.-H.; Hu, J.-J. Potential of mean force for human lysozyme-camelid vhh hl6 antibody interaction studies. Chem. Phys. Lett 2008, 455, 284–288. [Google Scholar]
  22. Wang, Y.-T.; Su, Z.-Y.; Chen, C.-L. Potential of mean force of the hepatitis C virus core protein-monoclonal 19D9D6 antibody interaction. Biophys. Chem 2009, 145, 86–90. [Google Scholar]
  23. Wang, Y.-T.; Su, Z.-Y.; Liao, J.-M.; Chen, C.-L. Potential of mean force for Syrian hamster prion epitope protein—Monoclonal fab 3f4 antibody interaction studies. Eur. J. Med. Chem 2009, 44, 3504–3508. [Google Scholar]
  24. Liu, J.; Bradley, R.; Eckmann, D.M.; Ayyaswamy, P.S.; Radhakrishnan, R. Multiscale modeling of functionalized nanocarriers in targeted drug delivery. Curr. Nanosci 2011, 7, 727–735. [Google Scholar]
  25. Liu, J.; Weller, G.E.R.; Zern, B.; Ayyaswamy, P.S.; Eckmann, D.M.; Muzykantov, V.R.; Radhakrishnan, R. Computational model for nanocarrier binding to endothelium validated using in vivo, in vitro, and atomic force microscopy experiments. Proc. Natl. Acad. Sci. USA 2010, 107, 16530–16535. [Google Scholar]
  26. Bahar, I.; Erman, B.; Monnerie, L. Kinematics of polymer chains with freely rotating bonds in a restrictive environment. 2. Conformational and orientational correlations. Macromolecules 1992, 25, 6315–6321. [Google Scholar]
  27. Smith, C.A.; Kortemme, T. backrub-like backbone simulation recapitulates natural protein conformational variability and improves mutant side-chain prediction. J. Mol. Biol 2008, 380, 742–756. [Google Scholar]
  28. Gunasekaran, K.; Gomathi, L.; Ramakrishnan, C.; Chandrasekhar, J.; Balaram, P. Conformational interconversions in peptide β-turns: Analysis of turns in proteins and computational estimates of barriers. J. Mol. Biol 1998, 284, 1505–1516. [Google Scholar]
  29. Ariza, A.; Richard, D.J.; White, M.F.; Bond, C.S. Conformational flexibility revealed by the crystal structure of a crenarchaeal RadA. Nucleic Acids Res 2005, 33, 1465–1473. [Google Scholar]
  30. Shin, D.S.; Pellegrini, L.; Daniels, D.S.; Yelent, B.; Craig, L.; Bates, D.; Yu, D.S.; Shivji, M.K.; Hitomi, C.; Arvai, A.S.; et al. Full-length archaeal Rad51 structure and mutants: Mechanisms for RAD51 assembly and control by BRCA2. EMBO J 2003, 22, 4566–4576. [Google Scholar]
  31. Thompson, J.D.; Gibson, T.J.; Plewniak, F.; Jeanmougin, F.; Higgins, D.G. The CLUSTAL_X windows interface: Flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res 1997, 25, 4876–4882. [Google Scholar]
  32. Fiser, A.; Do, R.K.G.; Šali, A. Modeling of loops in protein structures. Protein Sci 2000, 9, 1753–1773. [Google Scholar]
  33. Wu, Y.; Qian, X.; He, Y.; Moya, I.A.; Luo, Y. Crystal structure of an ATPase-active form of Rad51 homolog from methanococcus voltae. J. Biol. Chem 2005, 280, 722–728. [Google Scholar]
  34. Case, D.A.; Cheatham, T.E.; Darden, T.; Gohlke, H.; Luo, R.; Merz, K.M.; Onufriev, A.; Simmerling, C.; Wang, B.; Woods, R.J. The Amber biomolecular simulation programs. J. Comput. Chem 2005, 26, 1668–1688. [Google Scholar]
  35. Laskowski, R.A.; Macarthur, M.W.; Moss, D.S.; Thornton, J.M. PROCHECK: A program to check the stereochemical quality of protein structures. J. Appl. Cryst 1993, 26, 283–291. [Google Scholar]
  36. Wiederstein, M.; Sippl, M.J. ProSA-web: Interactive web service for the recognition of errors in three-dimensional structures of proteins. Nucleic Acids Res 2007, 35, W407–W410. [Google Scholar]
  37. Jiang, W.; Hardy, D.J.; Phillips, J.C.; MacKerell, A.D.; Schulten, K.; Roux, B.T. High-performance scalable molecular dynamics simulations of a polarizable force field based on classical drude oscillators in NAMD. J. Phys. Chem. Lett 1997, 2, 87–92. [Google Scholar]
Figure 1. The RadA rotary motor model is from left-hand RadA filament (green) to closed-ring RadA filament (red), then to MvRadA filament (pink), then to right-hand RadA filament (orange) and finally to left-hand RadA filament (green). The N-terminal domain (NTD) is indicated in gray.
Figure 1. The RadA rotary motor model is from left-hand RadA filament (green) to closed-ring RadA filament (red), then to MvRadA filament (pink), then to right-hand RadA filament (orange) and finally to left-hand RadA filament (green). The N-terminal domain (NTD) is indicated in gray.
Ijms 13 07138f1
Figure 2. The multiple sequence alignment for the four RadA protein filaments. 2DFL is the left-hand RadA filament, 2BKE is the right-hand RadA filament. 1PZN is the closed-ring RadA filament, and 1T4G is the MvRadA filament.
Figure 2. The multiple sequence alignment for the four RadA protein filaments. 2DFL is the left-hand RadA filament, 2BKE is the right-hand RadA filament. 1PZN is the closed-ring RadA filament, and 1T4G is the MvRadA filament.
Ijms 13 07138f2
Figure 3. The total MD simulation energy profiles of the four RadA filaments.
Figure 3. The total MD simulation energy profiles of the four RadA filaments.
Ijms 13 07138f3
Figure 4. The energy profiles of the RadA filament simulations process from TMD calculations. Ijms 13 07138f7, Left-hand RadA filaments; Ijms 13 07138f8, Close-ring RadA filaments; Ijms 13 07138f9, MvRadA filaments; Ijms 13 07138f10, Right-hand RadA filaments.
Figure 4. The energy profiles of the RadA filament simulations process from TMD calculations. Ijms 13 07138f7, Left-hand RadA filaments; Ijms 13 07138f8, Close-ring RadA filaments; Ijms 13 07138f9, MvRadA filaments; Ijms 13 07138f10, Right-hand RadA filaments.
Ijms 13 07138f4
Figure 5. The RMSD profiles of the 40 ns TMD simulations. (A) Left-hand RadA is indicated in brown; (B) Closed-ring RadA is indicated in orange; (C) MvRadA is indicated in light green; (D) Right-hand RadA is indicated in dark green.
Figure 5. The RMSD profiles of the 40 ns TMD simulations. (A) Left-hand RadA is indicated in brown; (B) Closed-ring RadA is indicated in orange; (C) MvRadA is indicated in light green; (D) Right-hand RadA is indicated in dark green.
Ijms 13 07138f5
Figure 6. The profile of the rotational residues’ (I71-K88) cumulate changed dihedral angles (phi and psi).
Figure 6. The profile of the rotational residues’ (I71-K88) cumulate changed dihedral angles (phi and psi).
Ijms 13 07138f6
Table 1. The cumulate changed dihedral angles (I71-K88) of the rotational residues.
Table 1. The cumulate changed dihedral angles (I71-K88) of the rotational residues.
Residue NumberCumulate Changed Dihedral Angles

Phi (Φ)Psi (Ψ)
I7100
R720360
F733600
K7400
T75360360
A760360
L770360
E780360
V790360
K803600
K81360360
E8200
R8300
M843600
N850360
V86360360
K870720
K8800

Share and Cite

MDPI and ACS Style

Su, Z.-Y.; Lee, W.-J.; Su, W.-S.; Wang, Y.-T. Target Molecular Simulations of RecA Family Protein Filaments. Int. J. Mol. Sci. 2012, 13, 7138-7148. https://doi.org/10.3390/ijms13067138

AMA Style

Su Z-Y, Lee W-J, Su W-S, Wang Y-T. Target Molecular Simulations of RecA Family Protein Filaments. International Journal of Molecular Sciences. 2012; 13(6):7138-7148. https://doi.org/10.3390/ijms13067138

Chicago/Turabian Style

Su, Zhi-Yuan, Wen-Jay Lee, Wan-Sheng Su, and Yeng-Tseng Wang. 2012. "Target Molecular Simulations of RecA Family Protein Filaments" International Journal of Molecular Sciences 13, no. 6: 7138-7148. https://doi.org/10.3390/ijms13067138

APA Style

Su, Z. -Y., Lee, W. -J., Su, W. -S., & Wang, Y. -T. (2012). Target Molecular Simulations of RecA Family Protein Filaments. International Journal of Molecular Sciences, 13(6), 7138-7148. https://doi.org/10.3390/ijms13067138

Article Metrics

Back to TopTop