Human coronavirus OC43 3CL protease and the potential of ML188 as a broad-spectrum lead compound: Homology modelling and molecular dynamic studies
© Berry et al.; licensee BioMed Central. 2015
Received: 17 November 2014
Accepted: 2 April 2015
Published: 28 April 2015
The coronavirus 3 chymotrypsin-like protease (3CLpro) is a validated target in the design of potential anticoronavirus inhibitors. The high degree of homology within the protease’s active site and substrate conservation supports the identification of broad spectrum lead compounds. A previous study identified the compound ML188, also termed 16R, as an inhibitor of the Severe Acute Respiratory Syndrome coronavirus (SARS-CoV) 3CLpro. This study will detail the generation of a homology model of the 3CLpro of the human coronavirus OC43 and determine the potential of 16R to form a broad-spectrum lead compound. MODELLER was used to generate a suitable three-dimensional model of the OC43 3CLpro and the Prime module of Schrӧdinger predicted the binding conformation and free energy of binding of 16R within the 3CLpro active site. Molecular dynamics further confirmed ligand stability and hydrogen bonding networks.
A high quality homology model of the OC43 3CLpro was successfully generated in an active conformation. Further studies reproduced the binding pose of 16R within the active site of the generated model, where its free energy of binding was shown to equal that of the 3CLpro of SARS-CoV, a receptor it is experimentally proven to inhibit. The stability of the ligand was subsequently confirmed by molecular dynamics.
The lead compound 16R may represent a broad-spectrum inhibitor of the 3CLpro of OC43 and potentially other coronaviruses. This study provides an atomistic structure of the 3CLpro of OC43 and supports further experimental validation of the inhibitory effects of 16R. These findings further confirm that the 3CLpro of coronaviruses can be inhibited by broad spectrum lead compounds.
KeywordsHuman coronavirus OC43 3CLpro Homology modelling Molecular dynamics
Human coronaviruses have a worldwide distribution and are frequently associated with self-limiting upper respiratory tract disease or “the common cold”. They can, however, also present with high morbidity outcomes of the lower respiratory tract including bronchiolitis, pneumonia, [1-3], asthmatic exacerbations  and acute exacerbations of chronic obstructive pulmonary disease (COPD) , where human coronavirus OC43 has been shown to be the prevalent strain in several regions . Given the high burden of coronaviruses to human health, and their potential for genetic recombination to give rise to the emergence of completely novel strains, the presence of antiviral strategies is paramount. However, there is currently no commercially available molecular entity which is capable of inhibiting infection with human coronaviruses . The 3CLpro has been validated as an effective drug target in several studies and has even been termed “the Achilles’ heel of coronaviruses”  making it an ideal target for the identification of novel lead compounds.
During coronavirus replication the large open reading frame (ORF) 1a/1ab genes, located at the 5’ end of the genome, are responsible for expressing two large replicase polyproteins (pp). These are co- or post-translationally cleaved by the virally encoded 3CLpro and Papain-like protease to yield 16 non-structural proteins responsible for viral replication . The 3CLpro is also termed the main protease as it cleaves a total of 11 cleavage sites within pp1a and pp1ab , in comparison to only three cleavage sites predicted for the papain-like protease .
The 3CLpro has three distinct domains. Domains I and II are largely composed of several anti-parallel β barrels  and is connected to domain III by a long loop region. Domain III is composed of several globular α-helices and plays a role in protein dimerization, along with the N-terminal N-finger (residues 1–7). The active site is located in a chymotrypsin-like fold between domains I and II and contains a catalytic dyad of His41, which acts as a proton acceptor, and Cys144/5, which undergoes nucleophilic attack on the carbonyl carbon of the substrate [8,13]. The 3CLpro is so named for its close structural and sequence homology to the 3 chymotrypsin protease (3Cpro) of rhinoviruses, which contains a catalytic triad composed of His, Cys and Glu or Asp. Superimposition of the structures of 3CLpro and 3Cpro shows that the His and Cys of both proteases is almost perfectly aligned, which may explain the similar substrate specificity and catalytic mechanism . The position of the Glu/Asp in 3Cpro is, however, replaced by a water molecule in 3CLpro, which shares three hydrogen bonds with surrounding residues including the catalytic dyad member His41. It is assumed that this water molecule is responsible for stabilization of the histidine in the intermediate state during enzymatic cleavage .
All coronavirus 3CLpro share a high sequence homology, as well as main chain architecture and substrate conservation [15,16]. The substrate binding site of the 3CLpro has two deeply buried S1 and S2 subsites, with shallow S1’, S3 and S4 subsites with varying degrees of solvent exposure. Substrate specificity of coronavirus 3CLpro is mainly determined by the P1, P2 and P1’ positions . The P1 position has an absolute specificity for glutamine which stabilizes the S1 subsite via a hydrogen bond with the imidazole Nε2 of His162/3 and van der Waals interactions with surrounding residues of the S1 pocket. The P2 site has a preference for leucine or methionine to fill the hydrophobic S2 pocket. The sidechains of the S3 site are solvent exposed and therefore this site is expected to tolerate a wide range of functionality, but shows a preference for basic residues . Sidechains and backbones of residues surrounding the S4 site create a highly congested pocket which favors a small, hydrophobic residue in the P4 position, either Ser, Thr, Val or Pro [17-19]. The S1’ and S2’ subsites also accommodate small residues in the P1’ and P2’ positions, which may include Ser, Ala or Gly [18,20]. A typical cleavage recognition site is therefore (Ser, Ala)-(Val, Thr)-Leu-Glu ↓ (Ser, Ala, Gly), which is conserved among all coronavirus 3CLpro . Theoretically, since these viruses share a highly homologous binding site and substrate conservation, it is reasonable to anticipate that broad spectrum inhibitor strategies might be successful. .
A noncovalent, furyl amide ligand, 16R, which was identified by a high throughput screen, with subsequent lead optimization based on structure-activity relationships within the S1, S1’ and S2 binding pockets was shown to inhibit the 3CLpro of SARS-CoV at an IC50 of 1.5 ± 0.3 μM . This study will detail the generation of a theoretical model of the OC43 3CLpro based on a high sequence homology with the solved 3CLpro structure of HKU1 and assess the potential of the inhibitor, 16R, to inhibit this complex in a broad-spectrum manner. This study will therefore provide further evidence to support the identification of broad spectrum 3CLpro inhibitors.
Results and discussion
Homology modelling of the OC43 3CLpro
For the latter purposes of this model, in characterization of the inhibitory potential of 16R, it is essential to confirm that the receptor in its active conformation. To ensure the modelled 3CLpro of OC43 is in its active conformation, several features can be used to differentiate between the active and inactive states. These features include the maintenance of a loop conformation of the loop connecting domain II and III and the oxyanion loop. If these features assume an alpha helical conformation it directly or indirectly leads to the collapse of the oxyanion hole and thereby forms an inactive state . An alpha helix conformation of residues 186–190 in the loop region between domain II and III causes a further collapse of the S2 and S4 subsites . The orientation of the His163 sidechain in the S1 pocket is also vital for substrate binding . This residue is largely kept in place by a stacking interaction with the benzene ring of Phe140 and imidazole ring of His163 . The orientation of His172 is also stabilized by a hydrogen bond with the sidechain of Glu166. This prevents steric interactions between His172 and His163 allowing it to maintain its orientation in the S1 pocket. A hydrogen bond between His163 and Tyr161 has also been implicated in stabilizing the sidechain of His163 , however we could not observe the formation of this bond in the crystal structure of SARS-CoV 3CLpro (3V3M), raising the question of its importance to maintain the orientation of His162/3.
Assessing the binding conformation and free energy of binding of 16R
Analysis of the crystallographic structure of the SARS-CoV 3CLpro, with the bound inhibitor 16R, illustrates a pyridine occupying the S1 pocket and forming a hydrogen bond with the imidazole Nε2 of His163. A further two hydrogen bonds are formed between ligand carbonyls and Gly143 and Glu166 with three methyl groups inserting into the deep hydrophobic S2 pocket . The conformation of 16R is almost completely conserved when bound to the 3CLpro of OC43. Notable exceptions include the shifting of the hydrogen bond formed with Gly143 to a furan ring oxygen as opposed to the carbonyl oxygen seen in SARS-CoV 3CLpro. The hydrogen bond formed with Glu166 is also absent in OC43 3CLpro, however the distance between the ligand carbonyl and backbone hydrogen is 2.75 Å and therefore this may be capable of mediating hydrogen bond formation in a dynamic system (Figure 8). Free energy of binding between 16R and SARS-CoV and OC43 3CLpro, as assessed by MM-GBSA, is -85 kcal/mol for both receptors. These results suggest that the inhibitor 16R may be capable of inhibiting both complexes.
Molecular dynamic simulation of OC43 and 16R
To further confirm the findings described in the previous section, and to assess the potential of Glu166 to form a hydrogen bond in a dynamic system, molecular dynamic simulations were used. Molecular dynamics of the SARS-CoV 3CLpro-16R complex are in strong agreement with the crystal data of 3V3M . The simulation predicts an average of 3.17 hydrogen bonds formed over the 10 ns trajectory, with these bonds predominantly forming with Glu166, His163 and Gly143. The bond formed at Gly143 does however display the potential to form a bifurcated interaction with Asn142. As predicted by the previous results, detailed in the section describing the binding conformation, the 3CLpro-16R complex of OC43 displayed a high occupancy of 100% and 85% at both His163 and Gly143 positions respectively. These results however suggested that the bond formed with Gly143 was mediated via the furan ring oxygen, where molecular dynamic simulations indicated that this bond was formed between the carbonyl oxygen and Gly143, as seen with the crystal structure of SARS-CoV. Molecular dynamics also confirmed that the hydrogen bond at the Glu166 does not stably form in a dynamic environment and only displays the potential to form at 0-2 ns and 4-5 ns, with an overall occupancy of 21.62%.
The stability of the bound ligand within the active pocket was also assessed via RMSD, radius of gyration and changes in solvent accessible surface area. A RMS deviation below 2 Å is an indication of a ligand which is stably bound to its receptor , where an increase in radius of gyration from the starting reference would suggest ligand instability. An increase in solvent accessible surface area of the ligand may indicate the ligand is dissociating into surrounding solvent . Ligand RMSD within the active site was approximately 0.75 Å and radius of gyration and solvent accessible surface area remained within a stable range of 3.9 Å to 4.0 Å and 4.8 nm2 to 5.4 nm2 respectively (Figure 9). These parameters suggest that the ligand is highly stable within the active site.
The three-dimensional structure of the OC43 3CLpro was successfully solved by homology modelling with MODELLER. Analysis of various side chains and loop conformations within and surrounding the substrate-binding site is indicative of an active conformation of the enzyme. The solved structure has been made publicly available in the Protein Model Database (PMDB; https://bioinformatics.cineca.it/PMDB/) [PMDB ID: PM0079872]. Further analysis of a previously identified lead compound, 16R, suggests that it makes extensive and stable interactions with the 3CLpro binding site of OC43. Additional analysis of the free energy of binding by MM-GBSA suggests that the ligand binds with the same affinity to the 3CLpro of OC43 as it does to SARS-CoV, a receptor it is experimentally proven to inhibit. These results are strongly suggestive that the lead compound, 16R, not only displays the potential to inhibit the 3CLpro of SARS-CoV but also that of the highly homologous OC43. In additional studies (data not shown) 16R was shown to remain stably bound to the 3CLpro of 229E and NL63 in molecular dynamic simulations and may therefore even display the potential to inhibit the 3CLpro of additional coronaviruses. 16R is therefore an excellent lead compound in the pursuit of true broad spectrum inhibitors of all coronaviruses.
Homology modelling of the OC43 3CLpro
The amino acid sequence of the OC43 3CLpro was obtained from the NCBI database [GenBank: AEN19363]. Fold assignment with the profile.build module of MODELLER (9v10) identified the most suitable template structure from PDB with an 82.3% sequence identity and e-value of 0. Template-target alignment was conducted with the align.2d module, which also incorporates structural information from the template. Models of the OC43 3CLpro were built by satisfaction of homology derived spatial restraints on distances and dihedral angles , stereochemical restraints obtained from the CHARMM22 force field  and statistical preferences on dihedral angles and non-bonded distances obtained from known protein structures [22,31]. Models were built by the automodel class by minimizing the violations on all restraints and united-atom PDB models were generated.
DOPE scores are assigned by considering the positions of all non-hydrogen atoms, where the lowest DOPE score corresponds to the model that is predicted to be most accurate . The evaluate_model.py script was further utilised to calculate per residue DOPE scores, which were superimposed over the template structure DOPE scores. A general conservation in the pattern of respective DOPE scores was used as an indication of an accurate model. A low RMSD, as assessed by PyMOL , indicated structural homology between solved structure and template. The GA341 score was used to analyse fold-assignment. Z-scores, obtained from ProSA-web [33,34] were used to analyse overall model quality and to assure template and query models were in a range acceptable for structures of their respective size. A Z-score outside a range characteristic for native proteins will suggest the production of an erroneous model . Stereochemical analysis of dihedral angles was then further conducted using PROCHECK [35,36] on models corresponding to the lowest DOPE score with GA341 scores amounting to 1.0. Stereochemical analysis by PROCHECK assessed residue geometry, where more than 90% of residues should be located in the “most favoured region” to indicate a stereochemically favourable structure.
Once an appropriate model accuracy was achieved based on these previous observations, the active site of the enzyme was analysed to ensure it was in its active state. Several features are able to differentiate an active from inactive state of the 3CLpro. These include the conformation of the oxyanion loop and loop connecting domain II and III [19,26] and the sidechain orientation of residues, including His163, Phe140, His172, Tyr161 and Glu166, within the S1 pocket [26,27].
Assessing the binding conformation and free energy of binding of 16R
Protein structures of 3V3M (3CLpro of SARS-CoV) and the homology modelling derived structure of OC43 3CLpro were prepared with the Schrӧdinger Protein Preparation Wizard. Histidine protonation was assigned to the epsilon nitrogen in accordance with crystal data. A restrained minimization was then performed using the OPLS2005 force field [37,38], during which heavy atoms were restrained to remain within 0.3 Å of the original structure. To assess the conformation of the known inhibitor, 16R, when binding to the OC43 3CLpro, the experimental structure of 3V3M and OC43 3CLpro were aligned in Maestro, the interface for all Schrӧdinger’s software. The inhibitor, 16R, was subsequently superimposed over the OC43 3CLpro active site using coordinates present in the 3V3M crystal structure. Side chains of residues in the active pocket assuming orientations that resulted in a steric clash with the ligand were refined and the ligand was minimized with Prime [39,40] to better fit the pocket. Prime utilizes both empirically and theoretically derived constraints to predict and replicate the dynamic motion of protein sidechains, allowing for conformational changes that arise through the “induced-fit” phenomenon, which is often neglected in molecular docking screens.
The Prime/MM-GBSA method was used to calculate the binding free energy of 16R within the respective receptor. The free energy of binding is the calculated difference between the minimized receptor-inhibitor complex and the sum of the energies of the minimized unbound inhibitor and receptor. Ligand poses were minimized using the local optimization feature in Prime, where energies of the complex were calculated with the OPLS-2005 force field and GBSA continuum solvent model.
Molecular dynamic simulation of the 3CLpro-16R complex
The CHARMM27 all atom force field was assigned to receptor structures using the three point TIP3P water model. All input hydrogens in the coordinate file were ignored and were assigned according to the force field. Histidine protonation states were assigned to the epsilon nitrogen in the neutral form in accordance with crystallographic data. The ligand, in the previously predicted pose, was prepared by SwisParam . The protein-ligand complex was then placed in the centre of a solvated, cubic box, 1 nm from the box edge. The system was solvated with the spc216 water model and sodium counter ions were added to neutralize the system. The system was energy minimized by a brief dynamics simulation using Steepest Descent Minimization until steepest descent converged to Fmax < 1000. Long range electrostatic interactions were calculated with the Particle Mesh Ewald (PME) method  with a cutoff of 1 nm and periodic boundary conditions were applied. Both NVT and NPT equilibrations were run for 50 000 steps or 100 ps using a 2 fs time step and a leap-frog integrator. All bonds were constrained by the lincs algorithm. PME was again used for long range electrostatics with a 0.16 fourier spacing. Short range electrostatic and van der Waals cutoffs were 1.0 nm and 1.4 nm respectively. Coordinates, velocities and energies were saved every 0.2 ps. During NVT equilibration temperature coupling was achieved by V-rescale algorithm with a target temperature of 300 K. Pressure coupling during NPT equilibration was achieved via the Parrinello-Rahman algorithm and is incorporated in the NPT equilibration once temperature is stable to ensure the proper density is reached (approximately 1000 kg/m3). Following NVT and NPT system equilibration an unrestrained, 10 ns simulation (5 000 000 steps) was run with the leap frog integrator, saving coordinates, velocities and energies every 2 ps. All other parameters were maintained from the NVT and NPT equilibration including both temperature and pressure coupling.
This research was supported by grants provided by the Council for Scientific and Industrial Research (CSIR, South Africa), Poliomyelitis Research Foundation (PRF) and Ernst and Ethel Erikson Trust. The authors thank Novartis (Basel, Switzerland) and H3D (University of the Cape Town) for significant intellectual input and access to resources.
- Pene F, Merlat A, Vabret A, Rozenberg F, Buzyn A, Dreyfus F, et al. Coronavirus 229E-related pneumonia in immunocompromised patients. Clin Infect Dis. 2003;37(7):929–32.View ArticlePubMedGoogle Scholar
- Woo PC, Lau SK, Tsoi H-w, Huang Y, Poon RW, Chu C-m, et al. Clinical and molecular epidemiological features of coronavirus HKU1–associated community-acquired pneumonia. J Infect Dis. 2005;192(11):1898–907.View ArticlePubMedGoogle Scholar
- Vabret A, Mourez T, Gouarin S, Petitjean J, Freymuth F. An outbreak of coronavirus OC43 respiratory infection in Normandy, France. Clin Infect Dis. 2003;36(8):985–9.View ArticlePubMedGoogle Scholar
- Lau SKP, Woo PCY, Yip CCY, Tse H, Tsoi H, Cheng VCC, et al. Coronavirus HKU1 and Other Coronavirus Infections in Hong Kong. J Clin Microbiol. 2006;44(6):2063–71.View ArticlePubMed CentralPubMedGoogle Scholar
- Gorse GJ, O’Connor TZ, Hall SL, Vitale JN, Nichol KL. Human coronavirus and acute respiratory illness in older adults with chronic obstructive pulmonary disease. J Infect Dis. 2009;199(6):847–57.View ArticlePubMedGoogle Scholar
- Jean A, Quach C, Yung A, Semret M. Severity and outcome associated with human coronavirus OC43 infections among children. Pediatr Infect Dis J. 2013;32(4):325–9.View ArticlePubMedGoogle Scholar
- Ramajayam R, Tan KP, Liu HG, Liang PH. Synthesis, docking studies, and evaluation of pyrimidines as inhibitors of SARS-CoV 3CL protease. Bioorg Med Chem Lett. 2010;20(12):3569–72.View ArticlePubMedGoogle Scholar
- Yang H, Bartlam M, Rao Z. Drug Design targeting the main protease, the Achilles heel of Coronaviruses. Curr Pharm Des. 2006;12(35):4573–90.View ArticlePubMedGoogle Scholar
- Pyrc K, Berkhout B, van der Hoek L. The novel human coronaviruses NL63 and HKU1. J Virol. 2007;81(7):3051–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Ziebuhr J: The coronavirus replicase. Coronavirus replication and reverse genetics. Springer Berlin Heidelberg, 2005. 57-94.Google Scholar
- Gao F, Ou H-Y, Chen L-L, Zheng W-X, Zhang C-T. Prediction of proteinase cleavage sites in polyproteins of coronaviruses and its applications in analyzing SARS-CoV genomes. FEBS Lett. 2003;553(3):451–6.View ArticlePubMedGoogle Scholar
- Jacobs J, Grum-Tokars V, Zhou Y, Turlington M, Saldanha SA, Chase P, et al. Discovery, synthesis, and structure-based optimization of a series of N-(tert-butyl)-2-(N-arylamido)-2-(pyridin-3-yl) acetamides (ML188) as potent noncovalent small molecule inhibitors of the severe acute respiratory syndrome coronavirus (SARS-CoV) 3CL protease. J Med Chem. 2013;56(2):534–46.View ArticlePubMed CentralPubMedGoogle Scholar
- Yang H, Xie W, Xue X, Yang K, Ma J, Liang W, et al. Design of wide-spectrum inhibitors targeting coronavirus main proteases. PLoS Biol. 2005;3(10), e324.View ArticlePubMed CentralPubMedGoogle Scholar
- Lee CC, Kuo CJ, Ko TP, Hsu MF, Tsui YC, Chang SC, et al. Structural basis of inhibition specificities of 3C and 3C-like proteases by zinc-coordinating and peptidomimetic compounds. J Biol Chem. 2009;284(12):7646–55.View ArticlePubMed CentralPubMedGoogle Scholar
- Zhao Q, Li S, Xue F, Zou Y, Chen C, Bartlam M, et al. Structure of the main protease from a global infectious human coronavirus, HCoV-HKU1. J Virol. 2008;82(17):8647–55.View ArticlePubMed CentralPubMedGoogle Scholar
- Hegyi A, Ziebuhr J. Conservation of substrate specificities among coronavirus main proteases. J Gen Virol. 2002;83(3):595–9.PubMedGoogle Scholar
- Chuck C-P, Chow H-F, Wan DC-C, Wong K-B. Profiling of substrate specificities of 3C-like proteases from group 1, 2a, 2b, and 3 coronaviruses. PLoS One. 2011;6(11), e27228.View ArticlePubMed CentralPubMedGoogle Scholar
- Anand K, Ziebuhr J, Wadhwani P, Mesters JR, Hilgenfeld R. Coronavirus main proteinase (3CLpro) structure: basis for design of anti-SARS drugs. Science. 2003;300(5626):1763–7.View ArticlePubMedGoogle Scholar
- Xue X, Yu H, Yang H, Xue F, Wu Z, Shen W, et al. Structures of two coronavirus main proteases: implications for substrate binding and antiviral drug design. J Virol. 2008;82(5):2515–27.View ArticlePubMed CentralPubMedGoogle Scholar
- Hsu MF, Kuo CJ, Chang KT, Chang HC, Chou CC, Ko TP, et al. Mechanism of the maturation process of SARS-CoV 3CL protease. J Biol Chem. 2005;280(35):31257–66.View ArticlePubMedGoogle Scholar
- Ziebuhr J, Snijder EJ, Gorbalenya AE. Virus-encoded proteinases and proteolytic processing in the Nidovirales. J Gen Virol. 2000;81(4):853–79.PubMedGoogle Scholar
- Šali A, Overington JP. Derivation of rules for comparative protein modeling from a database of protein structure alignments. Protein Sci. 1994;3(9):1582–96.View ArticlePubMed CentralPubMedGoogle Scholar
- Šali A, Blundell TL. Comparative protein modelling by satisfaction of spatial restraints. J Mol Biol. 1993;234(3):779–815.View ArticlePubMedGoogle Scholar
- Arnold K, Bordoli L, Kopp J, Schwede T. The SWISS-MODEL workspace: a web-based environment for protein structure homology modelling. Bioinformatics. 2006;22(2):195–201.View ArticlePubMedGoogle Scholar
- Costanzi S. On the applicability of GPCR homology models to computer-aided drug discovery: a comparison between in silico and crystal structures of the β2-adrenergic receptor. J Med Chem. 2008;51(10):2907–14.View ArticlePubMed CentralPubMedGoogle Scholar
- Wu CG, Cheng SC, Chen SC, Li JY, Fang YH, Chen YH, et al. Mechanism for controlling the monomer-dimer conversion of SARS coronavirus main protease. Acta Crystallogr D Biol Crystallogr. 2013;69(Pt 5):747–55.View ArticlePubMedGoogle Scholar
- Zheng K, Ma G, Zhou J, Zen M, Zhao W, Jiang Y, et al. Insight into the activity of SARS main protease: Molecular dynamics study of dimeric and monomeric form of enzyme. Proteins. 2007;66(2):467–79.View ArticlePubMedGoogle Scholar
- Ni H, Sotriffer CA, McCammon JA. Ordered water and ligand mobility in the HIV-1 integrase-5CITEP complex: a molecular dynamics study. J Med Chem. 2001;44(19):3043–7.View ArticlePubMedGoogle Scholar
- Eisenhaber F, Lijnzaad P, Argos P, Sander C, Scharf M. The double cubic lattice method: efficient approaches to numerical integration of surface area and volume and to dot surface contouring of molecular assemblies. J Comput Chem. 1995;16(3):273–84.View ArticleGoogle Scholar
- MacKerell AD, Bashford D, Bellott M, Dunbrack R, Evanseck J, Field MJ, et al. All-atom empirical potential for molecular modeling and dynamics studies of proteins. J Phys Chem B. 1998;102(18):3586–616.View ArticlePubMedGoogle Scholar
- My S, Šali A. Statistical potential for assessment and prediction of protein structures. Protein Sci. 2006;15(11):2507–24.View ArticleGoogle Scholar
- DeLano WL. The PyMOL molecular graphics system. 2002.Google Scholar
- Wiederstein M, Sippl MJ. ProSA-web: interactive web service for the recognition of errors in three-dimensional structures of proteins. Nucleic Acids Res. 2007;35 suppl 2:W407–10.View ArticlePubMed CentralPubMedGoogle Scholar
- Sippl MJ. Recognition of errors in three‐dimensional structures of proteins. Proteins: Struct, Funct, Bioinf. 1993;17(4):355–62.View ArticleGoogle Scholar
- Laskowski RA, MacArthur MW, Moss DS, Thornton JM. PROCHECK: a program to check the stereochemical quality of protein structures. J Appl Crystallogr. 1993;26(2):283–91.View ArticleGoogle Scholar
- Laskowski R, MacArthur M, Moss D, Thornton J. PROCHECK: A program to produce both detailed and schematic plots of proteins. J Appl Crystallogr. 1993;24:946–56.Google Scholar
- Jorgensen WL, Maxwell DS, Tirado-Rives J. Development and testing of the OPLS all-atom force field on conformational energetics and properties of organic liquids. J Am Chem Soc. 1996;118(45):11225–36.View ArticleGoogle Scholar
- Shivakumar D, Williams J, Wu Y, Damm W, Shelley J, Sherman W. Prediction of absolute solvation free energies using molecular dynamics free energy perturbation and the OPLS force field. J Chem Theory Computation. 2010;6(5):1509–19.View ArticleGoogle Scholar
- Jacobson MP, Friesner RA, Xiang Z, Honig B. On the role of the crystal environment in determining protein side-chain conformations. J Mol Biol. 2002;320(3):597–608.View ArticlePubMedGoogle Scholar
- Jacobson MP, Pincus DL, Rapp CS, Day TJ, Honig B, Shaw DE, et al. A hierarchical approach to all‐atom protein loop prediction. Proteins: Struct, Funct, Bioinf. 2004;55(2):351–67.View ArticleGoogle Scholar
- Zoete V, Cuendet MA, Grosdidier A, Michielin O. SwissParam: a fast force field generation tool for small organic molecules. J Comput Chem. 2011;32(11):2359–68.View ArticlePubMedGoogle Scholar
- Essmann U, Perera L, Berkowitz ML, Darden T, Lee H, Pedersen LG. A smooth particle mesh Ewald method. J Chem Physics. 1995;103(19):8577–93.View ArticleGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.