Evaluating the accuracy of the AMBER protein force fields in modeling dihydrofolate reductase structures: misbalance in the conformational arrangements of the flexible loop domains.

Protein flexible loop regions were once thought to be simple linkers between other more functional secondary structural elements. However, as it becomes clearer that these loop domains are critical players in a plethora of biological processes, accurate conformational sampling of 3D loop structures is vital to the advancement of drug design techniques and the overall growth of knowledge surrounding molecular systems. While experimental techniques provide a wealth of structural information, the resolution of flexible loop domains is sometimes low or entirely absent due to their complex and dynamic nature. This highlights an opportunity for de novo structure prediction using in silico methods with molecular dynamics (MDs). This study evaluates some of the AMBER protein force field's (ffs) ability to accurately model dihydrofolate reductase (DHFR) conformations, a protein complex characterized by specific arrangements and interactions of multiple flexible loops whose conformations are determined by the presence or absence of bound ligands and cofactors. Although the AMBER ffs, including ff19SB, studied well model most protein structures with rich secondary structure, results obtained here suggest the inability to significantly sample the expected DHFR loop-loop conformations - of the six distinct protein-ligand systems simulated, a majority lacked consistent stabilization of experimentally derived metrics definitive the three enzyme conformations. Although under-sampling and the chosen ff parameter combinations could be the cause, given past successes with these MD approaches for many protein systems, this suggests a potential misbalance in available ff parameters required to accurately predict the structure of multiple flexible loop regions present in proteins.Communicated by Ramaswamy H. Sarma.

Computational Modeling of Stapled Peptides toward a Treatment Strategy for CML and Broader Implications in the Design of Lengthy Peptide Therapeutics.

The oncogenic gene product Bcr-Abl is the principal cause of chronic myeloid leukemia, and although several therapies exist to curb the aberrant kinase activity of Bcr-Abl through targeting of the Abl kinase domain, these therapies are rendered ineffective by frequent mutations in the corresponding gene. It has been demonstrated that a designed protein, known as CCmut3, is able to produce a dominant negative inactivating effect on Bcr-Abl kinase by preferentially oligomerizing with the N-terminal coiled-coil oligomerization domain of Bcr-Abl (Bcr-CC) to effectively reduce the oncogenic potential of Bcr-Abl. However, the sheer length of the CCmut3 peptide introduces a high degree of conformational variability and opportunity for targeting by intracellular proteolytic mechanisms. Here, we have examined the effects of introducing one or two molecular staples, or cross-links, spanning i, i + 7 backbone residues of the CCmut3 construct, which have been suggested to reinforce α-helical conformation, enhance cellular internalization, and increase resistance to proteolytic degradation, leading to enhanced pharmacokinetic properties. The importance of optimizing staple location along a highly tuned biological construct such as CCmut3 has been widely emphasized and, as such, we have employed in silico techniques to swiftly build, relax, and characterize a large number of candidates. This approach effectively allowed exploring each and every possible staple location along the peptide backbone so that every possible candidate is considered. Although many of the stapled candidate peptides displayed enhanced binding characteristics for Bcr-CC and improved conformational stability in the (Bcr-CC) bound form, simulations of the stapled peptides in the unbound form revealed widespread conformational variability among stapled candidates dependent on staple type and location, implicating the molecular replacement of helix-stabilizing residues with staple-containing residues in disrupting the native α-helical conformation of CCmut3, further highlighting a need for careful optimization of the CCmut3 construct. A candidate set has been assembled, which retains the native backbone α-helical integrity in both the bound and unbound forms while providing enhanced binding affinity for the Bcr-CC target, as research disseminated in this manuscript is intended to guide the development of a next-generation CCmut3 inhibitor peptide in an experimental setting.

Application of Thiol-yne/Thiol-ene Reactions for Peptide and Protein Macrocyclizations.

The application of thiol-yne/thiol-ene reactions to synthesize mono- and bicyclic-stapled peptides and proteins is reported. First, a thiol-ene-based peptide-stapling method in aqueous conditions was developed. This method enabled the efficient stapling of recombinantly expressed coil-coiled proteins. The resulting stapled protein demonstrated higher stability in its secondary structure than the unstapled version. Furthermore, a thiol-yne coupling was performed by using an α,ω-diyne to react with two cysteine residues to synthesize a stapled peptide with two vinyl sulfide groups. The stapled peptide could further react with another biscysteine peptide to yield a bicyclic stapled peptide with enhanced properties. For example, the cell permeability of a stapled peptide was further increased by appending an oligoarginine cell-penetrating peptide. The robustness and versatility of thiol-yne/thiol-ene reactions that can be applied to both synthetic and expressed peptides and proteins were demonstrated.

Re-engineered p53 chimera with enhanced homo-oligomerization that maintains tumor suppressor activity.

The use of the tumor suppressor p53 for gene therapy of cancer is limited by the dominant negative inactivating effect of mutant endogenous p53 in cancer cells. We have shown previously that swapping the tetramerization domain (TD) of p53 with the coiled-coil (CC) from Bcr allows for our chimeric p53 (p53-CC) to evade hetero-oligomerization with endogenous mutant p53. This enhances the utility of this construct, p53-CC, for cancer gene therapy. Because domain swapping to create p53-CC could result in p53-CC interacting with endogenous Bcr, which is ubiquitous in cells, modifications on the CC domain are necessary to minimize potential interactions with Bcr. Hence, we investigated the possible design of mutations that will improve homodimerization of CC mutants and disfavor hetero-oligomerization with wild-type CC (CCwt), with the goal of minimizing potential interactions with endogenous Bcr in cells. This involved integrated computational and experimental approaches to rationally design an enhanced version of our chimeric p53-CC tumor suppressor. Indeed, the resulting lead candidate p53-CCmutE34K-R55E avoids binding to endogenous Bcr and retains p53 tumor suppressor activity. Specifically, p53-CCmutE34K-R55E exhibits potent apoptotic activity in a variety of cancer cell lines, regardless of p53 status (in cells with mutant p53, wild-type p53, or p53-null cells). This construct overcomes the dominant negative effect limitation of wt p53 and has high significance for future gene therapy for treatment of cancers characterized by p53 dysfunction, which represent over half of all human cancers.