The TRAnsient Pockets in Proteins (TRAPP) webserver provides an automated workflow that allows users to explore the dynamics of a protein binding site and to detect pockets or sub-pockets that may transiently open due to protein internal motion. These transient or cryptic sub-pockets may be of interest in the design and optimization of small molecular inhibitors for a protein target of interest. The TRAPP workflow consists of the following three modules: (i) TRAPP structure— generation of an ensemble of structures using one or more of four possible molecular simulation methods; (ii) TRAPP analysis—superposition and clustering of the binding site conformations either in an ensemble of structures generated in step (i) or in PDB structures or trajectories uploaded by the user; and (iii) TRAPP pocket—detection, analysis, and visualization of the binding pocket dynamics and characteristics, such as volume, solvent-exposed area or properties of surrounding residues. A standard sequence conservation score per residue or a differential score per residue, for comparing on- and off-targets, can be calculated and displayed on the binding pocket for an uploaded multiple sequence alignment file, and known protein sequence annotations can be displayed simultaneously. The TRAPP webserver is freely available at http://trapp.h-its.org.

The optimization of compounds with multiple targets is a difficult multidimensional problem in the drug discovery cycle. Here, we present a systematic, multidisciplinary approach to the development of selective antiparasitic compounds. Computational fragment-based design of novel pteridine derivatives along with iterations of crystallographic structure determination allowed for the derivation of a structure-activity relationship for multitarget inhibition. The approach yielded compounds showing apparent picomolar inhibition of T. brucei pteridine reductase 1 (PTR1), nanomolar inhibition of L. major PTR1, and selective submicromolar inhibition of parasite dihydrofolate reductase (DHFR) versus human DHFR. Moreover, by combining design for polypharmacology with a property-based on-parasite optimization, we found three compounds that exhibited micromolar EC₅₀ values against T. brucei brucei while retaining their target inhibition. Our results provide a basis for the further development of pteridine-based compounds, and we expect our multitarget approach to be generally applicable to the design and optimization of anti-infective agents.

Pteridine reductase 1 (PTR1) is a folate and pterin pathway enzyme unique for pathogenic trypanosomatids. As a validated drug target, PTR1 has been the focus of recent research efforts aimed at finding more effective treatments against human parasitic diseases such as leishmaniasis or sleeping sickness. Previous PTR1-centered structural studies highlighted the enzyme characteristics, such as flexible regions around the active site, highly conserved structural waters, and species-specific differences in pocket properties and dynamics, which likely impacts the binding of natural substrates and inhibitors. Furthermore, several aspects of the PTR1 function, such as the substrate inhibition phenomenon and the level of ligand binding cooperativity in the enzyme homotetramer, likely related to the global enzyme dynamics, are poorly known at the molecular level. We postulate that future drug design efforts could greatly benefit from a better understanding of these phenomena through studying both the local and global PTR1 dynamics. This review highlights the key aspects of the PTR1 structure and dynamics relevant to structure-based drug design that could be effectively investigated by modeling approaches. Particular emphasis is given to the perspective of molecular dynamics, what has been accomplished in this area to date, and how modeling could impact the PTR1-targeted drug design in the future.

The conformational properties of the aminoacyl-tRNA binding site (A-site), and its surroundings in the Escherichia coli 30S ribosomal subunit, are of great relevance in designing antibacterial agents. The 30S subunit A-site is near ribosomal protein S12, which neighbors helices h27 and H69; this latter helix, of the 50S subunit, is a functionally important component of an intersubunit bridge. Experimental work has shown that specific point mutations in S12 (K42A, R53A) yield hyper-accurate ribosomes, which in turn confers resistance to the antibiotic ‘paromomycin’ (even when this aminoglycoside is bound to the A-site). Suspecting that these effects can be elucidated in terms of the local atomic interactions and detailed dynamics of this region of the bacterial ribosome, we have used molecular dynamics simulations to explore the motion of a fragment of the E. coli ribosome, including the A-site. We found that the ribosomal regions surrounding the A-site modify the conformational space of the flexible A-site adenines 1492/93. Specifically, we found that A-site mobility is affected by stacking interactions between adenines A1493 and A1913, and by contacts between A1492 and a flexible side-chain (K43) from the S12 protein. In addition, our simulations reveal possible indirect pathways by which the R53A and K42A mutations in S12 are coupled to the dynamical properties of the A-site. Our work extends what is known about the atomistic dynamics of the A-site, and suggests possible links between the biological effects of hyper-accurate mutations in the S12 protein and conformational properties of the ribosome; the implications for S12 dynamics help elucidate how the miscoding effects of paromomycin may be evaded in antibiotic-resistant mutants of the bacterial ribosome.