Coarse-grained (CG) modelling with the Martini force field has come of age. By combining a variety of bead types and sizes with a new mapping approach, the newest version of the model is able to accurately simulate large biomolecular complexes at millisecond timescales. In this perspective, we discuss possible applications of the Martini 3 model in drug discovery and development pipelines and highlight areas for future development. Owing to its high simulation efficiency and extended chemical space, Martini 3 has great potential in the area of drug design and delivery. However, several aspects of the model should be improved before Martini 3 CG simulations can be routinely employed in academic and industrial settings. These include the development of automatic parameterisation protocols for a variety of molecule types, the improvement of backmapping procedures, the description of protein flexibility and the development of methodologies enabling efficient sampling. We illustrate our view with examples on key areas where Martini could give important contributions such as drugs targeting membrane proteins, cryptic pockets and protein–protein interactions and the development of soft drug delivery systems.

The aim of rational drug design is to develop small molecules using a quantitative approach to optimize affinity. This should enhance the development of chemical compounds that would specifically, selectively, reversibly, and with high affinity interact with a target protein. It is not yet possible to develop such compounds using computational (i.e., in silico) approach and instead the lead molecules are discovered in high-throughput screening searches of large compound libraries. The main reason why in silico methods are not capable to deliver is our poor understanding of the compound structure–thermodynamics and structure–kinetics correlations. There is a need for databases of intrinsic binding parameters (e.g., the change upon binding in standard Gibbs energy (ΔGint), enthalpy (ΔHint), entropy (ΔSint), volume (ΔVintr), heat capacity (ΔCp,int), association rate (ka,int), and dissociation rate (kd,int)) between a series of closely related proteins and a chemically diverse, but pharmacophoric group-guided library of compounds together with the co-crystal structures that could help explain the structure–energetics correlations and rationally design novel compounds. Assembly of these data will facilitate attempts to provide correlations and train data for modeling of compound binding. Here, we report large datasets of the intrinsic thermodynamic and kinetic data including over 400 primary sulfonamide compound binding to a family of 12 catalytically active human carbonic anhydrases (CA). Thermodynamic parameters have been determined by the fluorescent thermal shift assay, isothermal titration calorimetry, and by the stopped-flow assay of the inhibition of enzymatic activity. Kinetic measurements were performed using surface plasmon resonance. Intrinsic thermodynamic and kinetic parameters of binding were determined by dissecting the binding-linked protonation reactions of the protein and sulfonamide. The compound structure–thermodynamics and kinetics correlations reported here helped to discover compounds that exhibited picomolar affinities, hour-long residence times, and million-fold selectivities over non-target CA isoforms. Drug-lead compounds are suggested for anticancer target CA IX and CA XII, antiglaucoma CA IV, antiobesity CA VA and CA VB, and other isoforms. Together with 85 X-ray crystallographic structures of 60 compounds bound to six CA isoforms, the database should be of help to continue developing the principles of rational target-based drug design.

The apicomplexan protozoan parasites include the causative agents of animal and human diseases ranging from malaria (Plasmodium spp.) to toxoplasmosis (Toxoplasma gondii). The complex life cycle of T. gondii is regulated by a unique family of calcium-dependent protein kinases (CDPKs) that have become the target of intensive efforts to develop new therapeutics. In this review, we will summarize structure-based strategies, recent successes and future directions in the pursuit of specific and selective inhibitors of T. gondii CDPK1.

In recent decades, in silico absorption, distribution, metabolism, excretion (ADME), and toxicity (T) modelling as a tool for rational drug design has received considerable attention from pharmaceutical scientists, and various ADME/T-related prediction models have been reported. The high-throughput and low-cost nature of these models permits a more streamlined drug development process in which the identification of hits or their structural optimization can be guided based on a parallel investigation of bioavailability and safety, along with activity. However, the effectiveness of these tools is highly dependent on their capacity to cope with needs at different stages, e.g. their use in candidate selection has been limited due to their lack of the required predictability. For some events or endpoints involving more complex mechanisms, the current in silico approaches still need further improvement. In this review, we will briefly introduce the development of in silico models for some physicochemical parameters, ADME properties and toxicity evaluation, with an emphasis on the modelling approaches thereof, their application in drug discovery, and the potential merits or deficiencies of these models. Finally, the outlook for future ADME/T modelling based on big data analysis and systems sciences will be discussed.

The need for new drugs to treat microbial infections is pressing. The great progress made in the middle part of the twentieth Century was followed by a period of relative inactivity as the medical needs relating to infectious disease in the wealthier nations receded. Growing realisation that anti-infectives are needed in many parts of the world, to treat neglected diseases as well as to combat the burgeoning risk of resistance to existing drugs, has galvanised a new wave of research into anti-microbial drugs. The transfer of knowledge from the Pharmaceutical industry relating to the importance of understanding how to target drugs successfully within the body, and improved understanding of how pathogens interact with their hosts, is driving a series of new paradigms in anti-infective drug design. Here we provide an overview of those processes as an introduction to a series of articles from experts in this area that emerged from a meeting entitled “Emerging Paradigms in Anti-Infective Drug Design” held in London on the 17th and 18th September 2012. The symposium was organised jointly by British Society for Parasitology (BSP) and the Biological & Medicinal Chemistry sector of the Royal Society of Chemistry (RSC) and held at the London School of Hygiene & Tropical Medicine. The symposium set out to cover all aspects of the identification of new therapeutic modalities for the treatment of neglected and tropical diseases. We aimed to bring together leading scientists from all the disciplines working in this field and cover the pharmacology, medicinal chemistry and drug delivery of potential new medicines. Sessions were held on: “Target diseases and targets for drugs”, “Target based medicinal chemistry”, “Bioavailability and chemistry”, “Targeting intracellular microbes”, “Alternative approaches and models”, and “New anti-infectives – how do we get there?”

This symposium was organised by Simon Croft (LSHTM) and Mike Barrett (University of Glasgow) for the BSP, and David Alker (David Alker Associates) and Andrew Stachulski (University of Liverpool) for the Biological & Medicinal Chemistry sector of the RSC.

There is an alarming rate of human African trypanosomiasis recrudescence in many parts of sub-Saharan Africa. Yet, the disease has no successful chemotherapy. Trypanosoma lacks the enzymatic machinery for the de novo synthesis of purine nucleotides, and is critically dependent on salvage mechanisms. Inosine 5′-monophosphate dehydrogenase (IMPDH) is responsible for the rate-limiting step in guanine nucleotide metabolism. Here, we characterize recombinant Trypanosoma brucei IMPDH (TbIMPDH) to investigate the enzymatic differences between TbIMPDH and host IMPDH. Size-exclusion chromatography and analytical ultracentrifugation sedimentation velocity experiments reveal that TbIMPDH forms a heptamer, different from type 1 and 2 mammalian tetrameric IMPDHs. Kinetic analysis reveals calculated Km values of 30 and 1300 μm for IMP and NAD, respectively. The obtained Km value of TbIMPDH for NAD is approximately 20–200-fold higher than that of mammalian enzymes and indicative of a different NAD binding mode between trypanosomal and mammalian IMPDHs. Inhibition studies show Ki values of 3·2 μm, 21 nM and 3·3 nM for ribavirin 5′-monophosphate, mycophenolic acid and mizoribine 5′-monophosphate, respectively. Our results show that TbIMPDH is different from its mammalian counterpart and thus may be a good target for further studies on anti-trypanosomal drugs.

Infection with parasitic helminths takes a heavy toll on the health and well-being of humans and their domestic livestock, concomitantly resulting in major economic losses. Analyses have consistently revealed bioactive molecules in extracts of helminths or in their excretory/secretory products that modulate the immune response of the host. It is our view that parasitic helminths are an untapped source of immunomodulatory substances that, in pure form, could become new drugs (or models for drug design) to treat disease. Here, we illustrate the range of immunomodulatory molecules in selected parasitic trematodes, cestodes and nematodes, their impact on the immune cells in the host and how the host may recognize these molecules. There are many examples of the partial characterization of helminth-derived immunomodulatory molecules, but these have not yet translated into new drugs, reflecting the difficulty of isolating and fully characterizing proteins, glycoproteins and lipid-based molecules from small amounts of parasite material. However, this should not deter the investigator, since analytical techniques are now being used to accrue considerable structural information on parasite-derived molecules, even when only minute quantities of tissue are available. With the introduction of methodologies to purify and structurally-characterize molecules from small amounts of tissue and the application of high throughput immunological assays, one would predict that an assessment of parasitic helminths will yield a variety of novel drug candidates in the coming years.

Dihydrofolate reductase-thymidylate synthase (DHFR-TS) from Plasmodium falciparum, a validated target for antifolate antimalarials, is a dimeric enzyme with interdomain interactions significantly mediated by the junction region as well as the Plasmodium-specific additional sequences (inserts) in the DHFR domain. The X-ray structures of both the wild-type and mutant enzymes associated with drug resistance, in complex with either a drug which lost, or which still retains, effectiveness for the mutants, reveal features which explain the basis of drug resistance resulting from mutations around the active site. Binding of rigid inhibitors like pyrimethamine and cycloguanil to the enzyme active site is affected by steric conflict with the side-chains of mutated residues 108 and 16, as well as by changes in the main chain configuration. The role of important residues on binding of inhibitors and substrates was further elucidated by site-directed and random mutagenesis studies. Guided by the active site structure and modes of inhibitor binding, new inhibitors with high affinity against both wild-type and mutant enzymes have been designed and synthesized, some of which have very potent antimalarial activities against drug-resistant P. falciparum bearing the mutant enzymes.

The aim of this paper is to provide an outline of the current developments in drug discovery. However, in doing so, it is important to put these developments into an historical perspective and to acknowledge the impact that the previous research effort has had on research in progress.

Subdomain IIId from the hepatitis C virus (HCV) internal ribosome entry site (IRES) has been shown to be essential for cap-independent translation. We have conducted a structural study of a 27-nt fragment, identical in sequence to IIId, to explore the structural features of this subdomain. The proposed secondary structure of IIId is comprised of two 3 bp helical regions separated by an internal loop and closed at one end by a 6-nt terminal loop. NMR and molecular modeling were used interactively to formulate a validated model of the three-dimensional structure of IIId. We found that this fragment contains several noncanonical structural motifs and non-Watson–Crick base pairs, some of which are common to other RNAs. In particular, a motif characteristic of the rRNA α-sarcin/ricin loop was located in the internal loop. The terminal loop, 5′-UUGGGU, was found to fold to form a trinucleotide loop closed by a trans-wobble U˚G base pair. The sixth nucleotide was bulged out to allow stacking of this U˚G pair on the adjacent helical region. In vivo mutational analysis in the context of the full IRES confirmed the importance of each structural motif within IIId for IRES function. These findings may provide clues as to host cellular proteins that play a role in IRES-directed translation and, in particular, the mechanism through which host ribosomes are sequestered for viral function.

The monoclonal antibody 1696, directed against the HIV-1 protease, displays strong inhibitory effects toward the catalytic activity of the enzyme of both the HIV-1 and HIV-2 isolates. This antibody cross-reacts with peptides that include the N-terminus of the enzyme, a region that is well conserved in sequence among different viral strains and which, furthermore, is crucial for homodimerization to the active enzymatic form. This observation, as well as antigen-binding studies in the presence of an active site inhibitor, suggest that 1696 inhibits the HIV protease by destabilizing its active homodimeric form. To characterize further how the antibody 1696 inhibits the HIV-1 and HIV-2 proteases, we have solved the crystal structure of its Fab fragment by molecular replacement and refined it at 3.0 Å resolution. The antigen binding site has a deep cavity at its center, which is lined mainly by acidic and hydrophobic residues, and is large enough to accommodate several antigen residues. The structure of the Fab 1696 could form a starting basis for the design of alternative HIV protease-inhibiting molecules of broad specificity.

The oligopeptide-binding protein OppA provides a useful model system for studying the physical chemistry underlying noncovalent interactions since it binds a variety of readily synthesized ligands. We have studied the binding of eight closely related tripeptides of the type Lysine-X-Lysine, where X is an abnormal amino acid, by isothermal titration calorimetry (ITC) and X-ray crystallography. The tripeptides fall into three series of ligands, which have been designed to examine the effects of small changes to the central side chain. Three ligands have a primary amine as the second side chain, two have a straight alkane chain, and three have ring systems. The results have revealed a definite preference for the binding of hydrophobic residues over the positively charged side chains, the latter binding only weakly due to unfavorable enthalpic effects. Within the series of positively charged groups, a point of lowest affinity has been identified and this is proposed to arise from unfavorable electrostatic interactions in the pocket, including the disruption of a key salt bridge. Marked entropy-enthalpy compensation is found across the series, and some of the difficulties in designing tightly binding ligands have been highlighted.

Effective inhibitors of matrix metalloproteinases (MMPs), a family of connective tissue-degrading enzymes, could be useful for the treatment of diseases such as cancer, multiple sclerosis, and arthritis. Many of the known MMP inhibitors are derived from peptide substrates, with high potency in vitro but little selectivity among MMPs and poor bioavailability. We have discovered nonpeptidic MMP inhibitors with improved properties, and report here the crystal structures of human stromelysin-1 catalytic domain (SCD) complexed with four of these inhibitors. The structures were determined and refined at resolutions ranging from 1.64 to 2.0 Å. Each inhibitor binds in the active site of SCD such that a bulky diphenyl piperidine moiety penetrates a deep, predominantly hydrophobic S′1 pocket. The active site structure of the SCD is similar in all four inhibitor complexes, but differs substantially from the peptide hydroxamate complex, which has a smaller side chain bound in the S′1 pocket. The largest differences occur in the loop forming the “top” of this pocket. The occupation of these nonpeptidic inhibitors in the S′1 pocket provides a structural basis to explain their selectivity among MMPs. An analysis of the unique binding mode predicts structural modifications to design improved MMP inhibitors.