Our interpretation of ligand-protein interactions is often informed by high-resolution structures, which represent the cornerstone of structure-based drug design. However, visual inspection and molecular mechanics approaches cannot explain the full complexity of molecular interactions. Quantum Mechanics approaches are often too computationally expensive, but one method, Fragment Molecular Orbital (FMO), offers an excellent compromise and has the potential to reveal key interactions that would otherwise be hard to detect. To illustrate this, we have applied the FMO method to 18 Class A GPCR-ligand crystal structures, representing different branches of the GPCR genome. Our work reveals key interactions that are often omitted from structure-based descriptions, including hydrophobic interactions, nonclassical hydrogen bonds, and the involvement of backbone atoms. This approach provides a more comprehensive picture of receptor-ligand interactions than is currently used and should prove useful for evaluation of the chemical nature of ligand binding and to support structure-based drug design.

Carbapenems are the most potent β-lactam antibiotics and key drugs for treating infections by Gram-negative bacteria. In such organisms, β-lactam resistance arises principally from β-lactamase production. Although carbapenems escape the activity of most β-lactamases, due in the class A enzymes to slow deacylation of the covalent acylenzyme intermediate, carbapenem-hydrolyzing class A β-lactamases are now disseminating in clinically relevant bacteria. The reasons why carbapenems are substrates for these enzymes, but inhibit other class A β-lactamases, remain to be fully established. Here, we present crystal structures of the class A carbapenemase SFC-1 from Serratia fonticola and of complexes of its Ser70 Ala (Michaelis) and Glu166 Ala (acylenzyme) mutants with the carbapenem meropenem. These are the first crystal structures of carbapenem complexes of a class A carbapenemase. Our data reveal that, in the SFC-1 acylenzyme complex, the meropenem 6α-1R-hydroxyethyl group interacts with Asn132, but not with the deacylating water molecule. Molecular dynamics simulations indicate that this mode of binding occurs in both the Michaelis and acylenzyme complexes of wild-type SFC-1. In carbapenem-inhibited class A β-lactamases, it is proposed that the deacylating water molecule is deactivated by interaction with the carbapenem 6α-1R-hydroxyethyl substituent. Structural comparisons with such enzymes suggest that in SFC-1 subtle repositioning of key residues (Ser70, Ser130, Asn132 and Asn170) enlarges the active site, permitting rotation of the carbapenem 6α-1R-hydroxyethyl group and abolishing this contact. Our data show that SFC-1, and by implication other such carbapenem-hydrolyzing enzymes, uses Asn132 to orient bound carbapenems for efficient deacylation and prevent their interaction with the deacylating water molecule.

Carbapenems, 'last resort' antibiotics for many bacterial infections, can now be broken down by several class A β-lactamases (i.e. carbapenemases). Here, carbapenemase activity is predicted through QM/MM dynamics simulations of acyl-enzyme deacylation, requiring only the 3D structure of the apo-enzyme. This may assist in anticipating resistance and future antibiotic design.

Inhibition of inducible T-cell kinase (ITK), a nonreceptor tyrosine kinase, may represent a novel treatment for allergic asthma. In our previous reports, we described the discovery of sulfonylpyridine (SAP), benzothiazole (BZT), indazole (IND), and tetrahydroindazole (THI) series as novel ITK inhibitors and how computational tools such as dihedral scans and docking were used to support this process. X-ray crystallography and modeling were applied to provide essential insight into ITK-ligand interactions. However, "visual inspection" traditionally used for the rationalization of protein-ligand affinity cannot always explain the full complexity of the molecular interactions. The fragment molecular orbital (FMO) quantum-mechanical (QM) method provides a complete list of the interactions formed between the ligand and protein that are often omitted from traditional structure-based descriptions. FMO methodology was successfully used as part of a rational structure-based drug design effort to improve the ITK potency of high-throughput screening hits, ultimately delivering ligands with potency in the subnanomolar range.

Simulation details including system setup, all steps of the efficient computational assay, WHAM analysis and energies for the alternative proton transfer reaction.

The understanding of binding interactions between any protein and a small molecule plays a key role in the rationalization of affinity and selectivity and is essential for an efficient structure-based drug discovery (SBDD) process. Clearly, to begin SBDD, a structure is needed, and although there has been fantastic progress in solving G-protein-coupled receptor (GPCR) crystal structures, the process remains quite slow and is not currently feasible for every GPCR or GPCR–ligand complex. This situation significantly limits the ability of X-ray crystallography to impact the drug discovery process for GPCR targets in ‘real-time’ and hence there is still a need for other practical and cost-efficient alternatives. We present here an approach that integrates our previously described hierarchical GPCR modelling protocol (HGMP) and the fragment molecular orbital (FMO) quantum mechanics (QM) method to explore the interactions and selectivity of the human orexin-2 receptor (OX2R) and its recently discovered nonpeptidic agonists. HGMP generates a 3D model of GPCR structures and its complexes with small molecules by applying a set of computational methods. FMO allows ab initio approaches to be applied to systems that conventional QM methods would find challenging. The key advantage of FMO is that it can reveal information on the individual contribution and chemical nature of each residue and water molecule to the ligand binding that normally would be difficult to detect without QM. We illustrate how the combination of both techniques provides a practical and efficient approach that can be used to analyse the existing structure–function relationships (SAR) and to drive forward SBDD in a real-world example for which there is no crystal structure of the complex available.

QM/MM modelling of FAAH inactivation by O-biphenyl-3-yl carbamates identifies the deprotonation of Ser241 as the key reaction step, explaining why FAAH is insensitive to the electrondonor effect of conjugated substituents; this may aid design of new inhibitors with improved selectivity and in vivo potency.Carbamate-based compounds are widely used as covalent inhibitors of serine hydrolases of therapeutic interest, including fatty acid amide hydrolase (FAAH). 1 This enzyme is characterized by an uncommon Ser-Ser-Lys catalytic triad. It is the main enzyme responsible for the hydrolysis of the endocannabinoid anandamide. Despite its unusual catalytic mechanism, 2-4 FAAH is inhibited by classical serine hydrolase inhibitors, 1 and by O-biphenyl-3-yl carbamates, which are promising clinical candidates for the treatment of central nervous system and peripheral disorders. 5 These carbamate inhibitors, e.g. URB524 (Fig. 1), can be docked in two possible orientations (called orientations I and II) within the FAAH catalytic site. 6,7 Recently, hybrid quantum mechanical/molecular mechanics (QM/MM) modelling, 8 using the B3LYP/6-31G(d)//PM3-CHARMM22 potential, showed that the inhibitory process is energetically preferred in orientation II. 9 This orientation allows the catalytic nucleophile, Ser241, to efficiently attack the carbonyl group of URB524 (Fig. 2), yielding a carbamoylated enzyme. 10

Self-consistent charge density functional tight binding (SCC-DFTB) is a promising method for hybrid quantum mechanics/molecular mechanics (QM/MM) simulations of enzyme-catalyzed reactions. The acylation reaction of fatty acid amide hydrolase (FAAH), a promising drug target, was investigated by applying a SCC-DFTB/CHARMM27 scheme. Calculated potential energy barriers resulted in reasonable agreement with experiments for oleamide (OA) and oleoylmethyl ester (OME) substrates, outperforming previous calculations performed at the PM3/CHARMM22 level. Furthermore, the experimental preference of FAAH in hydrolyzing OA faster than OME was adequately reproduced by calculations. All these findings indicate that the SCC-DFTB/CHARMM27 approach can be successfully applied to mechanistic investigations of FAAH-catalyzed reactions.