Electronic Structure, Binding Energy, and Solvation Structure of the Streptavidin−Biotin Supramolecular Complex: ONIOM and 3D-RISM Study

Binding Site Programmable Self-Assembly of 3D Hierarchical DNA Origami Nanostructures

DNA nanotechnology has broad applications in biomedical drug delivery and programmable materials. Characterization of the self-assembly of DNA origami and quantum dots (QDs) is necessary for the development of new DNA-based nanostructures. We use computation and experiment to show that the self-assembly of 3D hierarchical nanostructures can be controlled by programming the binding site number and their positions on DNA origami. Using biotinylated pentagonal pyramid wireframe DNA origamis and streptavidin capped QDs, we demonstrate that DNA origami with 1 binding site at the outer vertex can assemble multimeric origamis with up to 6 DNA origamis on 1 QD, and DNA origami with 1 binding site at the inner center can only assemble monomeric and dimeric origamis. Meanwhile, the yield percentages of different multimeric origamis are controlled by the QD:DNA-origami stoichiometric mixing ratio. DNA origamis with 2 binding sites at the αγ positions (of the pentagon) make larger nanostructures than those with binding sites at the αβ positions. In general, increasing the number of binding sites leads to increases in the nanostructure size. At high DNA origami concentration, the QD number in each cluster becomes the limiting factor for the growth of nanostructures. We find that reducing the QD size can also affect the self-assembly because of the reduced access to the binding sites from more densely packed origamis.

Reactivity and Selectivity of the Diels-Alder Reaction of Anthracene in [Pd6L4]12+ Supramolecular Cages: A Computational Study

The field of supramolecular metal-organic cage catalysis has grown rapidly in recent years. However, theoretical studies regarding the reaction mechanism and reactivity and selectivity controlling factors for supramolecular catalysis are still underdeveloped. Herein, we demonstrate a detailed density functional theory study on the mechanism, catalytic efficiency, and regioselectivity of the Diels-Alder reaction in bulk solution and within two [Pd₆L₄]¹²⁺ supramolecular cages. Our calculations are consistent with experiments. The origins of the catalytic efficiency of the bowl-shaped cage 1 have been elucidated to be the host-guest stabilization of the transition states and the favorable entropy effect. The reasons for the switch of the regioselectivity from 9,10-addition to 1,4-addition within the octahedral cage 2 were attributed to the confinement effect and the noncovalent interactions. This work would shed light on the understanding of [Pd₆L₄]¹²⁺ metallocage-catalyzed reactions and provide a detailed mechanistic profile otherwise difficult to obtain from experiments. The findings of this study could also aid to the improvement and development of more efficient and selective supramolecular catalysis.

Computational methods for calculation of protein-ligand binding affinities in structure-based drug design

Drug design is an expensive and time-consuming process. Any method that allows reducing the time the costs of the drug development project can have great practical value for the pharmaceutical industry. In structure-based drug design, affinity prediction methods are of great importance. The majority of methods used to predict binding free energy in protein-ligand complexes use molecular mechanics methods. However, many limitations of these methods in describing interactions exist. An attempt to go beyond these limits is the application of quantum-mechanical description for all or only part of the analyzed system. However, the extensive use of quantum mechanical (QM) approaches in drug discovery is still a demanding challenge. This chapter briefly reviews selected methods used to calculate protein-ligand binding affinity applied in virtual screening (VS), rescoring of docked poses, and lead optimization stage, including QM methods based on molecular simulations.

Integration of In Silico Strategies for Drug Repositioning towards P38α Mitogen-Activated Protein Kinase (MAPK) at the Allosteric Site

P38α mitogen-activated protein kinase (p38α MAPK), one of the p38 MAPK isoforms participating in a signaling cascade, has been identified for its pivotal role in the regulation of physiological processes such as cell proliferation, differentiation, survival, and death. Herein, by shedding light on docking- and 100-ns dynamic-based screening from 3210 FDA-approved drugs, we found that lomitapide (a lipid-lowering agent) and nilotinib (a Bcr-Abl fusion protein inhibitor) could alternatively inhibit phosphorylation of p38α MAPK at the allosteric site. All-atom molecular dynamics simulations and free energy calculations including end-point and QM-based ONIOM methods revealed that the binding affinity of the two screened drugs exhibited a comparable level as the known p38α MAPK inhibitor (BIRB796), suggesting the high potential of being a novel p38α MAPK inhibitor. In addition, noncovalent contacts and the number of hydrogen bonds were found to be corresponding with the great binding recognition. Key influential amino acids were mostly hydrophobic residues, while the two charged residues including E71 and D168 were considered crucial ones due to their ability to form very strong H-bonds with the focused drugs. Altogether, our contributions obtained here could be theoretical guidance for further conducting experimental-based preclinical studies necessary for developing therapeutic agents targeting p38α MAPK.

Combined Machine Learning and Molecular Modelling Workflow for the Recognition of Potentially Novel Fungicides

Novel machine learning and molecular modelling filtering procedures for drug repurposing have been carried out for the recognition of the novel fungicide targets of Cyp51 and Erg2. Classification and regression approaches on molecular descriptors have been performed using stepwise multilinear regression (FS-MLR), uninformative-variable elimination partial-least square regression, and a non-linear method called Forward Stepwise Limited Correlation Random Forest (FS-LM-RF). Altogether, 112 prediction models from two different approaches have been built for the descriptor recognition of fungicide hit compounds. Aiming at the fungal targets of sterol biosynthesis in membranes, antifungal hit compounds have been selected for docking experiments from the Drugbank database using the Autodock4 molecular docking program. The results were verified by Gold Protein-Ligand Docking Software. The best-docked conformation, for each high-scored ligand considered, was submitted to quantum mechanics/molecular mechanics (QM/MM) gradient optimization with final single point calculations taking into account both the basis set superposition error and thermal corrections (with frequency calculations). Finally, seven Drugbank lead compounds were selected based on their high QM/MM scores for the Cyp51 target, and three were selected for the Erg2 target. These lead compounds could be recommended for further in vitro studies.

Water in protein hydration and ligand recognition

This review describes selected basics of water in biomolecular recognition. We focus on a qualitative understanding of the most important physical aspects, how these change in magnitude between bulk water and protein environment, and how the roles that water plays for proteins arise from them. These roles include mechanical support, thermal coupling, dielectric screening, mass and charge transport, and the competition with a ligand for the occupation of a binding site. The presence or absence of water has ramifications that range from the thermodynamic binding signature of a single ligand up to cellular survival. The large inhomogeneity in water density, polarity and mobility around a solute is hard to assess in experiment. This is a source of many difficulties in the solvation of protein models and computational studies that attempt to elucidate or predict ligand recognition. The influence of water in a protein binding site on the experimental enthalpic and entropic signature of ligand binding is still a point of much debate. The strong water‐water interaction in enthalpic terms is counteracted by a water molecule's high mobility in entropic terms. The complete arrest of a water molecule's mobility sets a limit on the entropic contribution of a water displacement process, while the solvent environment sets limits on ligand reactivity.

Molecular dynamics simulation and free energy analysis of the interaction of platinum-based anti-cancer drugs with DNA

Cisplatin and oxaliplatin are two widely-used anti-cancer drugs which covalently bind to a same location in DNA strands. Platinum agents make intrastrand and interstrand cross-links with the N7 atoms of guanine nucleotides which prevent DNA from polymerization by causing a distortion in the double helix. Molecular dynamics simulations and free energy calculations were carried out to investigate the binding of two platinum-based anti-cancer drugs with DNA. We compared the binding of these drugs which differ in their carrier ligands, and hence their potential interactions with DNA. When a platinum agent binds to nucleotides, it causes a high amount of deformation in DNA structure. To find the extent of deformation, torsion angles and base pair and groove parameters of DNA were considered. These parameters were compared with normal B-DNA which was considered as the undamaged DNA. The formation of hydrogen bonds between drugs and DNA nucleotides was examined in solution. It was shown that oxaliplatin forms more hydrogen bonds than cisplatin. Our results confirm that the structure of the platinated DNA rearranges significantly and cisplatin tries to deform DNA more than oxaliplatin. The binding free energies were also investigated to understand the affinities, types and the contributions of interactions between drugs and DNA. It was concluded that oxaliplatin tendency for binding to DNA is more than cisplatin in solvent environment. The binding free energy was calculated based on the MM/PBSA and MM/GBSA methods and the results of QM/MM calculations verified them.

MTS-MD of Biomolecules Steered with 3D-RISM-KH Mean Solvation Forces Accelerated with Generalized Solvation Force Extrapolation

We developed a generalized solvation force extrapolation (GSFE) approach to speed up multiple time step molecular dynamics (MTS-MD) of biomolecules steered with mean solvation forces obtained from the 3D-RISM-KH molecular theory of solvation (three-dimensional reference interaction site model with the Kovalenko-Hirata closure). GSFE is based on a set of techniques including the non-Eckart-like transformation of coordinate space separately for each solute atom, extension of the force-coordinate pair basis set followed by selection of the best subset, balancing the normal equations by modified least-squares minimization of deviations, and incremental increase of outer time step in motion integration. Mean solvation forces acting on the biomolecule atoms in conformations at successive inner time steps are extrapolated using a relatively small number of best (closest) solute atomic coordinates and corresponding mean solvation forces obtained at previous outer time steps by converging the 3D-RISM-KH integral equations. The MTS-MD evolution steered with GSFE of 3D-RISM-KH mean solvation forces is efficiently stabilized with our optimized isokinetic Nosé-Hoover chain (OIN) thermostat. We validated the hybrid MTS-MD/OIN/GSFE/3D-RISM-KH integrator on solvated organic and biomolecules of different stiffness and complexity: asphaltene dimer in toluene solvent, hydrated alanine dipeptide, miniprotein 1L2Y, and protein G. The GSFE accuracy and the OIN efficiency allowed us to enlarge outer time steps up to huge values of 1-4 ps while accurately reproducing conformational properties. Quasidynamics steered with 3D-RISM-KH mean solvation forces achieves time scale compression of conformational changes coupled with solvent exchange, resulting in further significant acceleration of protein conformational sampling with respect to real time dynamics. Overall, this provided a 50- to 1000-fold effective speedup of conformational sampling for these systems, compared to conventional MD with explicit solvent. We have been able to fold the miniprotein from a fully denatured, extended state in about 60 ns of quasidynamics steered with 3D-RISM-KH mean solvation forces, compared to the average physical folding time of 4-9 μs observed in experiment.

A Review of Membrane-Based Biosensors for Pathogen Detection

Biosensors are of increasing interest for the detection of bacterial pathogens in many applications such as human, animal and plant health, as well as food and water safety. Membranes and membrane-like structures have been integral part of several pathogen detection platforms. Such structures may serve as simple mechanical support, function as a part of the transduction mechanism, may be used to filter out or concentrate pathogens, and may be engineered to specifically house active proteins. This review focuses on membrane materials, their associated biosensing applications, chemical linking procedures, and transduction mechanisms. The sensitivity of membrane biosensors is discussed, and the state of the field is evaluated and summarized.

Characterization of selective binding of biologically relevant inorganic ions with the proline zwitterion by 3D-RISM theory

The features of selective binding of several biologically relevant mono- and divalent inorganic ions with the proline zwitterion were studied over a wide range of electrolyte concentrations.

Efficient asymmetric synthesis of 1-cyano-tetrahydroisoquinolines from lipase dual activity and opposite enantioselectivities in α-Aminonitrile resolution

Dual promiscuous racemization/amidation activities of lipases leading to efficient dynamic kinetic resolution protocols of racemic α-aminonitrile compounds are described. α-Amidonitrile products of high enantiomeric purity could be formed in high yields. Several lipases from different sources were shown to exhibit the dual catalytic activities, where opposite enantioselectivities could be recorded for certain substrates.

Weak sulfur-sulfur interactions between chemically-identical atoms

Abstract Experimentally-known sulfur-sulfur distances shorter than the sum of van der Waals radii and involving two chemically-identical sulfur atoms are examined at several levels of theory (BP86/6-31G** to CCSD(T)/6-311+G**). None of the theoretical methods predict an attractive interaction from an energetic point of view, even though molecular orbitals stretching between the two sulfur atoms have been identified. Most likely, if there is indeed an attractive interaction force between chemically identical sulfur atoms, its value is comparable to the accuracy of the methods employed here — implying an attractive interaction below 1 kcal/mol. The investigation includes some simple models of 1,6,12,17-tetrathiacyclodocosa-2,4,13,15-tetrayne which was previously shown to have an S—S interaction involving two chemically-identical atoms. Attractive interactions calculated for these latter models are shown to arise from S—HC weak bonding, with the S—S interaction being again repulsive. Graphical abstract

Dissociation kinetics of the streptavidin-biotin interaction measured using direct electrospray ionization mass spectrometry analysis

Dissociation rate constants (k (off)) for the model high affinity interaction between biotin (B) and the homotetramer of natural core streptavidin (S(4)) were measured at pH 7 and temperatures ranging from 15 to 45 °C using electrospray ionization mass spectrometry (ESI-MS). Two different approaches to data analysis were employed, one based on the initial rate of dissociation of the (S(4) + 4B) complex, the other involving nonlinear fitting of the time-dependent relative abundances of the (S(4) + iB) species. The two methods were found to yield k (off) values that are in good agreement, within a factor of two. The Arrhenius parameters for the dissociation of the biotin-streptavidin interaction in solution were established from the k (off) values determined by ESI-MS and compared with values measured using a radiolabeled biotin assay. Importantly, the dissociation activation energies determined by ESI-MS agree, within 1 kcal mol(-1), with the reported value. In addition to providing a quantitative measure of k (off), the results of the ESI-MS measurements revealed that the apparent cooperative distribution of (S(4) + iB) species observed at short reaction times is of kinetic origin and that sequential binding of B to S(4) occurs in a noncooperative fashion with the four ligand binding sites being kinetically and thermodynamically equivalent and independent.

Streptavidin and its biotin complex at atomic resolution

Atomic resolution crystallographic studies of streptavidin and its biotin complex have been carried out at 1.03 and 0.95 Å, respectively. The wild-type protein crystallized with a tetramer in the asymmetric unit, while the crystals of the biotin complex contained two subunits in the asymmetric unit. Comparison of the six subunits shows the various ways in which the protein accommodates ligand binding and different crystal-packing environments. Conformational variation is found in each of the polypeptide loops connecting the eight strands in the β-sandwich subunit, but the largest differences are found in the flexible binding loop (residues 45-52). In three of the unliganded subunits the loop is in an `open' conformation, while in the two subunits binding biotin, as well as in one of the unliganded subunits, this loop `closes' over the biotin-binding site. The `closed' loop contributes to the protein's high affinity for biotin. Analysis of the anisotropic displacement parameters included in the crystallographic models is consistent with the variation found in the loop structures and the view that the dynamic nature of the protein structure contributes to the ability of the protein to bind biotin so tightly.

How the biotin-streptavidin interaction was made even stronger: investigation via crystallography and a chimaeric tetramer

The interaction between SA (streptavidin) and biotin is one of the strongest non-covalent interactions in Nature. SA is a widely used tool and a paradigm for protein–ligand interactions. We previously developed a SA mutant, termed Tr (traptavidin), possessing a 10-fold lower off-rate for biotin, with increased mechanical and thermal stability. In the present study, we determined the crystal structures of apo-Tr and biotin–Tr at 1.5 Å resolution. In apo-SA the loop (L3/4), near biotin's valeryl tail, is typically disordered and open, but closes upon biotin binding. In contrast, L3/4 was shut in both apo-Tr and biotin–Tr. The reduced flexibility of L3/4 and decreased conformational change on biotin binding provide an explanation for Tr's reduced biotin off- and on-rates. L3/4 includes Ser⁴⁵, which forms a hydrogen bond to biotin consistently in Tr, but erratically in SA. Reduced breakage of the biotin–Ser⁴⁵ hydrogen bond in Tr is likely to inhibit the initiating event in biotin's dissociation pathway. We generated a Tr with a single biotin-binding site rather than four, which showed a simi-larly low off-rate, demonstrating that Tr's low off-rate was governed by intrasubunit effects. Understanding the structural features of this tenacious interaction may assist the design of even stronger affinity tags and inhibitors.

Computational and experimental study of the structure, binding preferences, and spectroscopy of nickel(II) and vanadyl porphyrins in petroleum

We present a computational exploration of five- and six-coordinate Ni(II) and vanadyl porphyrins, including prediction of UV-vis spectroscopic behavior and metalloporphyrin structure as well as determination of a binding energy threshold between strongly bound complexes that have been isolated as single crystals and weakly bound ones that we detect by visible absorption spectroscopy. The excited states are calculated using the tandem of the time-dependent density functional theory (TD-DFT) and the conductor-like polarizable continuum model (CPCM). The excited-state energies in chloroform solvent obtained by using two density functionals are found to correlate linearly with the experimental Soret and alpha-band energies for a known series of five-coordinate vanadyl porphyrins. The established linear correction allows simulation of the excited states for labile octahedral vanadyl porphyrins that have not been isolated and yields Soret and alpha-band bathochromic shifts that are in agreement with our UV-vis spectroscopic results. The PBE0 and PW91 functionals in combination with DNP basis set perform best for both structure and binding energy prediction. The reactivity preferences of Ni(II) and vanadyl porphyrins toward aromatic fragments of large petroleum molecules are explored by using the density functional theory (DFT). Analysis of electrostatic potentials and Fukui functions matching shows that axial coordination and hydrogen bonding are the preferred aggregation modes between vanadyl porphyrins and nitrogen-containing heterocycle fragments. This investigation improves our understanding on the cause for broadening of the Ni and V porphyrin Soret band in heavy oils. Our findings can be useful for the development of metals removal methods for heavy oil upgrading.