Steffen Lindert

Accounting for receptor flexibility and enhanced sampling methods in computer-aided drug design.

Protein flexibility plays a major role in biomolecular recognition. In many cases, it is not obvious how molecular structure will change upon association with other molecules. In proteins, these changes can be major, with large deviations in overall backbone structure, or they can be more subtle as in a side‐chain rotation. Either way the algorithms that predict the favorability of biomolecular association require relatively accurate predictions of the bound structure to give an accurate assessment of the energy involved in association. Here, we review a number of techniques that have been proposed to accommodate receptor flexibility in the simulation of small molecules binding to protein receptors. We investigate modifications to standard rigid receptor docking algorithms and also explore enhanced sampling techniques, and the combination of free energy calculations and enhanced sampling techniques. The understanding and allowance for receptor flexibility are helping to make computer simulations of ligand protein binding more accurate. These developments may help improve the efficiency of drug discovery and development. Efficiency will be essential as we begin to see personalized medicine tailored to individual patients, which means specific drugs are needed for each patient’s genetic makeup.

Insight into the mechanisms of adenovirus capsid disassembly from studies of defensin neutralization.

Defensins are effectors of the innate immune response with potent antibacterial activity. Their role in antiviral immunity, particularly for non-enveloped viruses, is poorly understood. We recently found that human alpha-defensins inhibit human adenovirus (HAdV) by preventing virus uncoating and release of the endosomalytic protein VI during cell entry. Consequently, AdV remains trapped in the endosomal/lysosomal pathway rather than trafficking to the nucleus. To gain insight into the mechanism of defensin-mediated neutralization, we analyzed the specificity of the AdV-defensin interaction. Sensitivity to alpha-defensin neutralization is a common feature of HAdV species A, B1, B2, C, and E, whereas species D and F are resistant. Thousands of defensin molecules bind with low micromolar affinity to a sensitive serotype, but only a low level of binding is observed to resistant serotypes. Neutralization is dependent upon a correctly folded defensin molecule, suggesting that specific molecular interactions occur with the virion. CryoEM structural studies and protein sequence analysis led to a hypothesis that neutralization determinants are located in a region spanning the fiber and penton base proteins. This model was supported by infectivity studies using virus chimeras comprised of capsid proteins from sensitive and resistant serotypes. These findings suggest a mechanism in which defensin binding to critical sites on the AdV capsid prevents vertex removal and thereby blocks subsequent steps in uncoating that are required for release of protein VI and endosomalysis during infection. In addition to informing the mechanism of defensin-mediated neutralization of a non-enveloped virus, these studies provide insight into the mechanism of AdV uncoating and suggest new strategies to disrupt this process and inhibit infection.

Exploring the Role of Receptor Flexibility in Structure-Based Drug Discovery

The proper understanding of biomolecular recognition mechanisms that take place in a drug target is of paramount importance to improve the efficiency of drug discovery and development. The intrinsic dynamic character of proteins has a strong influence on biomolecular recognition mechanisms and models such as conformational selection have been widely used to account for this dynamic association process. However, conformational changes occurring in the receptor prior and upon association with other molecules are diverse and not obvious to predict when only a few structures of the receptor are available. In view of the prominent role of protein flexibility in ligand binding and its implications for drug discovery, it is of great interest to identify receptor conformations that play a major role in biomolecular recognition before starting rational drug design efforts. In this review, we discuss a number of recent advances in computer-aided drug discovery techniques that have been proposed to incorporate receptor flexibility into structure-based drug design. The allowance for receptor flexibility provided by computational techniques such as molecular dynamics simulations or enhanced sampling techniques helps to improve the accuracy of methods used to estimate binding affinities and, thus, such methods can contribute to the discovery of novel drug leads.

Antibacterial drug leads targeting isoprenoid biosynthesis.

With the rise in resistance to antibiotics such as methicillin, there is a need for new drugs. We report here the discovery and X-ray crystallographic structures of 10 chemically diverse compounds (benzoic, diketo, and phosphonic acids, as well as a bisamidine and a bisamine) that inhibit bacterial undecaprenyl diphosphate synthase, an essential enzyme involved in cell wall biosynthesis. The inhibitors bind to one or more of the four undecaprenyl diphosphate synthase inhibitor binding sites identified previously, with the most active leads binding to site 4, outside the catalytic center. The most potent leads are active against Staphylococcus aureus [minimal inhibitory concentration (MIC)90 ∼0.25 µg/mL], and one potently synergizes with methicillin (fractional inhibitory concentration index = 0.25) and is protective in a mouse infection model. These results provide numerous leads for antibacterial development and open up the possibility of restoring sensitivity to drugs such as methicillin, using combination therapies.

Cryo-Electron Microscopy Structure of Adenovirus Type 2 Temperature-Sensitive Mutant 1 Reveals Insight into the Cell Entry Defect

ABSTRACT The structure of the adenovirus type 2 temperature-sensitive mutant 1 (Ad2ts1) was determined to a resolution of 10 Å by cryo-electron microscopy single-particle reconstruction. Ad2ts1 was prepared at a nonpermissive temperature and contains the precursor forms of the capsid proteins IIIa, VI, and VIII; the core proteins VII, X (mu), and terminal protein (TP); and the L1-52K protein. Cell entry studies have shown that although Ad2ts1 can bind the coxsackievirus and Ad receptor and undergo internalization via αv integrins, this mutant does not escape from the early endosome and is targeted for degradation. Comparison of the Ad2ts1 structure to that of mature Ad indicates that Ad2ts1 has a different core architecture. The Ad2ts1 core is closely associated with the icosahedral capsid, a connection which may be mediated by preproteins IIIa and VI. Density within hexon cavities is assigned to preprotein VI, and membrane disruption assays show that hexon shields the lytic activity of both the mature and precursor forms of protein VI. The internal surface of the penton base in Ad2ts1 appears to be anchored to the core by interactions with preprotein IIIa. Our structural analyses suggest that these connections to the core inhibit the release of the vertex proteins and lead to the cell entry defect of Ad2ts1.

Cryo-Electron Microscopy Structure of an Adenovirus-Integrin Complex Indicates Conformational Changes in both Penton Base and Integrin

ABSTRACT A structure of adenovirus type 12 (HAdV12) complexed with a soluble form of integrin αvβ5 was determined by cryo-electron microscopy (cryoEM) image reconstruction. Subnanometer resolution (8 Å) was achieved for the icosahedral capsid with moderate resolution (27 Å) for integrin density above each penton base. Modeling with αvβ3 and αIIbβ3 crystal structures indicates that a maximum of four integrins fit over the pentameric penton base. The close spacing (∼60 Å) of the RGD protrusions on penton base precludes integrin binding in the same orientation to neighboring RGD sites. Flexible penton-base RGD loops and incoherent averaging of bound integrin molecules explain the moderate resolution observed for the integrin density. A model with four integrins bound to a penton base suggests that integrin might extend one RGD-loop in the direction that could induce a conformational change in the penton base involving clockwise untwisting of the pentamer. A global conformational change in penton base could be one step on the way to the release of Ad vertex proteins during cell entry. Comparison of the cryoEM structure with bent and extended models for the integrin ectodomain reveals that integrin adopts an extended conformation when bound to the Ad penton base, a multivalent viral ligand. These findings shed further light on the structural basis of integrin binding to biologically relevant ligands, as well as on the molecular events leading to HAdV cell entry.

Antiinfectives targeting enzymes and the proton motive force

Significance Uncoupling agents might be expected to be generally cytotoxic, but many US Food and Drug Administration (FDA)-approved drugs do have activity as uncouplers, in addition to targeting enzymes. There is therefore interest in the discovery of antibiotics that have such multitarget activity. Here, we show that some FDA-approved drugs, such as clofazimine, clomiphene, and bedaquiline, with antiinfective activity act as uncouplers. Using molecular dynamics-based in silico screening, we also discovered that the brain cancer drug lead vacquinol is an uncoupler that inhibits an enzyme involved in the formation of tuberculosis (TB) virulence factors, in addition to killing TB bacteria. Our results indicate strong drug–membrane interactions, and that screening for combined enzyme inhibition plus uncoupler activity will lead to new antibiotic leads. There is a growing need for new antibiotics. Compounds that target the proton motive force (PMF), uncouplers, represent one possible class of compounds that might be developed because they are already used to treat parasitic infections, and there is interest in their use for the treatment of other diseases, such as diabetes. Here, we tested a series of compounds, most with known antiinfective activity, for uncoupler activity. Many cationic amphiphiles tested positive, and some targeted isoprenoid biosynthesis or affected lipid bilayer structure. As an example, we found that clomiphene, a recently discovered undecaprenyl diphosphate synthase inhibitor active against Staphylococcus aureus, is an uncoupler. Using in silico screening, we then found that the anti-glioblastoma multiforme drug lead vacquinol is an inhibitor of Mycobacterium tuberculosis tuberculosinyl adenosine synthase, as well as being an uncoupler. Because vacquinol is also an inhibitor of M. tuberculosis cell growth, we used similarity searches based on the vacquinol structure, finding analogs with potent (∼0.5–2 μg/mL) activity against M. tuberculosis and S. aureus. Our results give a logical explanation of the observation that most new tuberculosis drug leads discovered by phenotypic screens and genome sequencing are highly lipophilic (logP ∼5.7) bases with membrane targets because such species are expected to partition into hydrophobic membranes, inhibiting membrane proteins, in addition to collapsing the PMF. This multiple targeting is expected to be of importance in overcoming the development of drug resistance because targeting membrane physical properties is expected to be less susceptible to the development of resistance.

Computational methods in drug discovery

The process for drug discovery and development is challenging, time consuming and expensive. Computer-aided drug discovery (CADD) tools can act as a virtual shortcut, assisting in the expedition of this long process and potentially reducing the cost of research and development. Today CADD has become an effective and indispensable tool in therapeutic development. The human genome project has made available a substantial amount of sequence data that can be used in various drug discovery projects. Additionally, increasing knowledge of biological structures, as well as increasing computer power have made it possible to use computational methods effectively in various phases of the drug discovery and development pipeline. The importance of in silico tools is greater than ever before and has advanced pharmaceutical research. Here we present an overview of computational methods used in different facets of drug discovery and highlight some of the recent successes. In this review, both structure-based and ligand-based drug discovery methods are discussed. Advances in virtual high-throughput screening, protein structure prediction methods, protein–ligand docking, pharmacophore modeling and QSAR techniques are reviewed.

Dynamics and calcium association to the N-terminal regulatory domain of human cardiac troponin C: a multiscale computational study.

Troponin C (TnC) is an important regulatory molecule in cardiomyocytes. Calcium binding to site II in TnC initiates a series of molecular events that result in muscle contraction. The most direct change upon Ca2+ binding is an opening motion of the molecule that exposes a hydrophobic patch on the surface allowing for Troponin I to bind. Molecular dynamics simulations were used to elucidate the dynamics of this crucial protein in three different states: apo, Ca2+-bound, and Ca2+-TnI-bound. Dynamics between the states are compared, and the Ca2+-bound system is investigated for opening motions. On the basis of the simulations, NMR chemical shifts and order parameters are calculated and compared with experimental observables. Agreement indicates that the simulations sample the relevant dynamics of the system. Brownian dynamics simulations are used to investigate the calcium association of TnC. We find that calcium binding gives rise to correlative motions involving the EF hand and collective motions conducive of formation of the TnI-binding interface. We furthermore indicate the essential role of electrostatic steering in facilitating diffusion-limited binding of Ca2+.

Hybrid Approaches: Applying Computational Methods in cryo-Electron Microscopy

Recent advances in cryo-electron microscopy have led to an increasing number of high (3-5A) to medium (5-10A) resolution cryoEM density maps. These density maps contain valuable information about the protein structure but frequently require computational algorithms to aid their structural interpretation. It is these hybrid approaches between cryoEM and computational protein structure prediction algorithms that will shape protein structure elucidation from density maps.