Decoding drug resistant mechanism of V32I, I50V and I84V mutations of HIV-1 protease on amprenavir binding by using molecular dynamics simulations and MM-GBSA calculations

Mutations V32I, I50V and I84V in the HIV-1 protease (PR) induce drug resistance towards drug amprenavir (APV). Multiple short molecular dynamics (MSMD) simulations and molecular mechanics generalized Born surface area (MM-GBSA) method were utilized to investigate drug-resistant mechanism of V32I, I50V and I84V towards APV. Dynamic information arising from MSMD simulations suggest that V32I, I50V and I84V highly affect structural flexibility, motion modes and conformational behaviours of two flaps in the PR. Binding free energies calculated by MM-GBSA method suggest that the decrease in binding enthalpy and the increase in binding entropy induced by mutations V32I, I50V and I84V are responsible for drug resistance of the mutated PRs on APV. The energetic contributions of separate residues on binding of APV to the PR show that V32I, I50V and I84V highly disturb the interactions of two flaps with APV and mostly drive the decrease in binding ability of APV to the PR. Thus, the conformational changes of two flaps in the PR caused by V32I, I50V and I84V play key roles in drug resistance of three mutated PR towards APV. This study can provide useful dynamics information for the design of potent inhibitors relieving drug resistance.

Can the solvent enhance the rate of chemical reactions through C-H/π interactions? insights from theory

The current computational study with density functional theory (DFT) shows that the rate of chemical reactions can be influenced through non-covalent C-H/π interactions between substrates and the solvent. It is shown that intramolecular carbon-carbon interaction and CO2 activation by a low valent silicon complex are both favourably affected by the explicit presence of the solvent toluene, due to C-H/π interactions between toluene and the silicon complex. Furthermore, ab initio molecular dynamics (AIMD) simulations demonstrate that even if the C-H/π interacting solvent molecule is displaced from the complex, another would quickly take its place, thus maintaining the interaction. Hence, the current work shows how non-covalent interactions between solvent and substrate can enhance the rate of the reaction and expands our understanding of the role and influence of the solvent in effecting important chemical transformations.

An insight into Csp–H⋯π hydrogen bonds and stability of complexes formed by acetylene and its substituted derivatives with benzene and borazine

Theoretical calculations at the MP2/aug-cc-pVDZ level are used to investigate the C_sp–H⋯π interactions of C₂HX (X = H, F, Cl, Br, CH₃, NH₂) with C₆H₆ and B₃N₃H₆ molecules.

A unified set of experimental organometallic data used to evaluate modern theoretical methods

We applied a test set of ligand dissociation enthalpies derived entirely from a unified experimental approach to evaluate the efficacy of various methods for modeling organometallic chemistry.

CH-π hydrogen bonds in biological macromolecules

This is a sequel to the previous Perspective "The CH-π hydrogen bond in chemistry. Conformation, supramolecules, optical resolution and interactions involving carbohydrates", which featured in a PCCP themed issue on "Weak Hydrogen Bonds - Strong Effects?": Phys. Chem. Chem. Phys., 2011, 13, 13873-13900. Evidence that weak hydrogen bonds play an enormously important role in chemistry and biochemistry has now accumulated to an extent that the rigid classical concept of hydrogen bonds formulated by Pauling needs to be seriously revised and extended. The concept of a more generalized hydrogen bond definition is indispensable for understanding the folding mechanisms of proteins. The CH-π hydrogen bond, a weak molecular force occurring between a soft acid CH and a soft base π-electron system, among all is one of the most important and plays a functional role in defining the conformation and stability of 3D structures as well as in many molecular recognition events. This concept is also valuable in structure-based drug design efforts. Despite their frequent occurrence in organic molecules and bio-molecules, the importance of CH-π hydrogen bonds is still largely unknown to many chemists and biochemists. Here we present a review that deals with the evidence, nature, characteristics and consequences of the CH-π hydrogen bond in biological macromolecules (proteins, nucleic acids, lipids and polysaccharides). It is hoped that the present Perspective will show the importance of CH-π hydrogen bonds and stimulate interest in the interactions of biological macromolecules, one of the most fascinating fields in bioorganic chemistry. Implication of this concept is enormous and valuable in the scientific community.

Enhanced semiempirical QM methods for biomolecular interactions

Recent successes and failures of the application of 'enhanced' semiempirical QM (SQM) methods are reviewed in the light of the benefits and backdraws of adding dispersion (D) and hydrogen-bond (H) correction terms. We find that the accuracy of SQM-DH methods for non-covalent interactions is very often reported to be comparable to dispersion-corrected density functional theory (DFT-D), while computation times are about three orders of magnitude lower. SQM-DH methods thus open up a possibility to simulate realistically large model systems for problems both in life and materials science with comparably high accuracy.

C-H···π interactions and the nature of the donor carbon atom

The influence of multiple substituents (F, CH3, NO2, CN, Cl, OH and NH2) on the C-H···π interaction in benzene-ethylene complex was investigated using the estimated CCSD(T) method and complete basis set limit. The results were compared with our earlier reported complexes of benzene-acetylene and benzene-methane, thus completing the sp, sp(2) and sp(3) series of C-H donors. The stabilization energy values for multiple fluoro-substituted benzene-ethylene complexes are found to be very close to those of the multiple fluoro-substituted benzene-methane complexes. Expectedly, the stabilization energies for the multiple methyl-substituted benzene-ethylene complexes lie between those of the multiple methyl-substituted benzene-methane and benzene-acetylene complexes. Energy decomposition analysis using the DFT-SAPT method predicts the dispersion energy to be dominant, similar to the benzene-methane complexes. For the symmetrically disubstituted complexes (-OH, -Cl, -NH2, -CN and -NO2), additional C-H···X interaction was observed, possibly due to the angular orientation of the ethylene molecule. Multidimensional correlation analysis between the electrostatic, dispersion and exchange-repulsion with the C-H···π interaction distance (r), Hammett constant (σ) and the molar refractivity (MR) revealed strong correlation between dispersion energy and the C-H···π interaction distance (r) as well as molar refractivity (MR).

Toward a more complete understanding of noncovalent interactions involving aromatic rings

Noncovalent interactions involving aromatic rings, which include π-stacking interactions, anion-π interactions, and XH-π interactions, among others, are ubiquitous in chemical and biochemical systems. Despite dramatic advances in our understanding of these interactions over the past decade, many aspects of these noncovalent interactions have only recently been uncovered, with many questions remaining. We summarize our computational studies aimed at understanding the impact of substituents and heteroatoms on these noncovalent interactions. In particular, we discuss our local, direct interaction model of substituent effects in π-stacking interactions. In this model, substituent effects are dominated by electrostatic interactions of the local dipoles associated with the substituents and the electric field of the other ring. The implications of the local nature of substituent effects on π-stacking interactions in larger systems are discussed, with examples given for complexes with carbon nanotubes and a small graphene model, as well as model stacked discotic systems. We also discuss related issues involving the interpretation of electrostatic potential (ESP) maps. Although ESP maps are widely used in discussions of noncovalent interactions, they are often misinterpreted. Next, we provide an alternative explanation for the origin of anion-π interactions involving substituted benzenes and N-heterocycles, and show that these interactions are well-described by simple models based solely on charge-dipole interactions. Finally, we summarize our recent work on the physical nature of substituent effects in XH-π interactions. Together, these results paint a more complete picture of noncovalent interactions involving aromatic rings and provide a firm conceptual foundation for the rational exploitation of these interactions in a myriad of chemical contexts.

Localized frustration and binding-induced conformational change in recognition of 5S RNA by TFIIIA zinc finger

Protein TFIIIA is composed of nine tandemly arranged Cys2His2 zinc fingers. It can bind either to the 5S RNA gene as a transcription factor or to the 5S RNA transcript as a chaperone. Although structural and biochemical data provided valuable information on the recognition between the TFIIIIA and the 5S DNA/RNA, the involved conformational motions and energetic factors contributing to the binding affinity and specificity remain unclear. In this work, we conducted MD simulations and MM/GBSA calculations to investigate the binding-induced conformational changes in the recognition of the 5S RNA by the central three zinc fingers of TFIIIA and the energetic factors that influence the binding affinity and specificity at an atomistic level. Our results revealed drastic interdomain conformational changes between these three zinc fingers, involving the exposure/burial of several crucial DNA/RNA binding residues, which can be related to the competition between DNA and RNA for the binding of TFIIIA. We also showed that the specific recognition between finger 4/finger 6 and the 5S RNA introduces frustrations to the nonspecific interactions between finger 5 and the 5S RNA, which may be important to achieve optimal binding affinity and specificity.

Understanding the molecular mechanism of the broad and potent neutralization of HIV-1 by antibody VRC01 from the perspective of molecular dynamics simulation and binding free energy calculations

VRC01 is one of the most broadly and potently neutralizing HIV-1 antibodies known-it has been shown to neutralize 91 % of the tested primary isolate Env pseudoviruses by recognizing the viral envelope glycoprotein gp120. To explore the mechanism of HIV-1 neutralization by VRC01 and thus obtain valuable information for vaccine design, we performed molecular dynamics simulations and binding free energy calculations for apo-VRC01, apo-gp120, and the gp120-VRC01 complex. For gp120, residue energy decomposition analysis showed that the hotspot residues Asn280, Lys282, Asp368, Ile371, and Asp457 are located in three primary loops, including the CD4-binding loop, loop D, and loop V5. For VRC01, the hotspot residues Trp47, Trp50, Asn58, Arg61, Gln64, Trp100, and Tyr91 mainly come from CDR2 of the heavy chain. By decomposing the binding free energy into different components, intermolecular van der Waals interactions and nonpolar solvation were found to dominate the binding process. Principal component analysis of loops D and V5, which are related to neutralization resistance, indicated that these two areas have a larger conformational space in apo-gp120 compared to bound gp120. A comparison of three representative structures from the cluster analysis of loops D and V5 indicated that changes primarily occur at the tip of loop V5, and are caused by fluctuations in the terminal Glu1 residue of the antibody. This information can be used to guide the design of vaccines and small molecule inhibitors.

Substituent effects on non-covalent interactions with aromatic rings: insights from computational chemistry

Non-covalent interactions with aromatic rings pervade modern chemical research. The strength and orientation of these interactions can be tuned and controlled through substituent effects. Computational studies of model complexes have provided a detailed understanding of the origin and nature of these substituent effects, and pinpointed flaws in entrenched models of these interactions in the literature. Here, we provide a brief review of efforts over the last decade to unravel the origin of substituent effects in π-stacking, XH/π, and ion/π interactions through detailed computational studies. We highlight recent progress that has been made, while also uncovering areas where future studies are warranted.