Martiniano Bello

Automated docking for novel drug discovery

Introduction: The volume of three-dimensional structural information of macromolecules and the number of computational tools to predict binding modes and affinities of molecular complexes are increasing daily. Molecular docking is a rational structural approach employed to predict thermodynamic parameters based on molecular recognition between two or more molecules. In addition, docking studies have become very important for therapeutic applications in modern structure-based drug design because this computational tool uses few economic resources. However, they omit many biological conditions that critically influence small and macromolecular structural motions. To mimic physiological conditions, it is necessary to consider other environmental factors, such as the presence of water molecules and the flexibility of ligands and side chain residues of proteins. Furthermore, molecular dynamics simulations have been coupled with docking procedures to expand the boundaries and obtain more reliable information. Areas covered: In this article, we review current advances in protein-small molecule docking and possible future directions. Expert opinion: Docking studies include many conformations to predict binding free energies (scoring functions) and to search (scoring sampling) for the most representative binding conformations. Therefore, several biological properties, from side chain residues to complete protein motions, have been included in docking studies to improve theoretical predictions.

Energetics of protein homodimerization: Effects of water sequestering on the formation of β-lactoglobulin dimer

Transient protein–protein interactions are functionally relevant as a control mechanism in a variety of biological processes. Analysis of the 3D structure of protein–protein complexes indicates that water molecules trapped at the interface are very common; however, their role in the stability and specificity of protein homodimer interactions has been not addressed yet. To provide new insights into the energetic bases that govern the formation of highly hydrated interfaces, the dissociation process of bovine βlg variant A at a neutral pH was characterized here thermodynamically by conducting dilution experiments with an isothermal titration calorimeter. Association was enthalpically driven throughout the temperature range spanned. ΔH and ΔCp were significantly more negative than estimates based on surface area changes, suggesting the occurrence of effects additional to the dehydration of the contact surfaces between subunits. Near‐UV CD spectra proved to be independent of protein concentration, indicating a rigid body‐like association. Furthermore, the process proved not to be coupled to significant changes in the protonation state of ionizable groups or counterion exchange. In contrast, both osmotic stress experiments and a computational analysis of the dimers 3D structure indicated that a large number of water molecules are incorporated into the interface upon association. Numerical estimates considering the contributions of interface area desolvation and water immobilization accounted satisfactorily for the experimental ΔCp. Thus, our study highlights the importance of explicitly considering the effects of water sequestering to perform a proper quantitative analysis of the formation of homodimers with highly hydrated interfaces. Proteins 2008.

Ligand binding and self-association cooperativity of β-lactoglobulin.

Unlike most small globular proteins, lipocalins lack a compact hydrophobic core. Instead, they present a large central cavity that functions as the primary binding site for hydrophobic molecules. Not surprisingly, these proteins typically exhibit complex structural dynamics in solution, which is intricately modified by intermolecular recognition events. Although many lipocalins are monomeric, an increasing number of them have been proven to form oligomers. The coupling effects between self‐association and ligand binding in these proteins are largely unknown. To address this issue, we have calorimetrically characterized the recognition of dodecyl sulfate by bovine β‐lactoglobulin, which forms weak homodimers at neutral pH. A thermodynamic analysis based on coupled‐equilibria revealed that dimerization exerts disparate effects on the ligand‐binding capacity of β‐lactoglobulin. Protein dimerization decreases ligand affinity (or, reciprocally, ligand binding promotes dimer dissociation). The two subunits in the dimer exhibit a positive, entropically driven cooperativity. To investigate the structural determinants of the interaction, the crystal structure of β‐lactoglobulin bound to dodecyl sulfate was solved at 1.64 Å resolution. Copyright

Structure and dynamics of β-lactoglobulin in complex with dodecyl sulfate and laurate: A molecular dynamics study

Bovine β-lactoglobulin (βlg) is able to recognize a wide variety of hydrophobic ligands. Although binding promiscuity is characteristic of highly hydrophobic interactions, the structural plasticity of the βlg binding cavity entrance seems to be crucial for the interaction with polar moieties of different ligands. On the other hand, thermodynamic studies have shown that βlg can associate to cognate ligands with distinctly different binding energetics, as in the case of the closely related molecules lauric acid (LA) and dodecyl sulfate (DS). In the recognition of LA, βlg shows a classical hydrophobic signature (entropically driven), whereas the interaction of βlg with DS exhibits a nonclassical hydrophobic signature (enthalpically driven). To gain insights into these opposed binding behaviors, MD simulations were carried out on βlg in apo-form and bound to DS or LA. Overall, the results suggested that the distinct energetic signatures of these ligands come from distinct optimizations of both hydrophilic and hydrophobic contacts with the protein.

Molecular Dynamics of a Thermostable Multicopper Oxidase from Thermus thermophilus HB27: Structural Differences between the Apo and Holo Forms

Molecular dynamic (MD) simulations have been performed on Tth-MCO, a hyperthermophilic multicopper oxidase from thermus thermophilus HB27, in the apo as well as the holo form, with the aim of exploring the structural dynamic properties common to the two conformational states. According to structural comparison between this enzyme and other MCOs, the substrate in process to electron transfer in an outer-sphere event seems to transiently occupy a shallow and overall hydrophobic cavity near the Cu type 1 (T1Cu). The linker connecting the β-strands 21 and 24 of the second domain (loop (β21–β24)D2) has the same conformation in both states, forming a flexible lid at the entrance of the electron-transfer cavity. Loop (β21–β24)D2 has been tentatively assigned a role occluding the access to the electron-transfer site. The dynamic of the loop (β21–β24)D2 has been investigated by MD simulation, and results show that the structures of both species have the same secondary and tertiary structure during almost all the MD simulations. In the simulation, loop (β21–β24)D2 of the holo form undergoes a higher mobility than in the apo form. In fact, loop (β21–β24)D2 of the holo form experiences a conformational change which enables exposure to the electron-transfer site (open conformation), while in the apo form the opposite effect takes place (closed conformation). To confirm the hypothesis that the open conformation might facilitate the transient electron-donor molecule occupation of the site, the simulation was extended another 40 ns with the electron-donor molecule docked into the protein cavity. Upon electron-donor molecule stabilization, loops near the cavity reduce their mobility. These findings show that coordination between the copper and the protein might play an important role in the general mobility of the enzyme, and that the open conformation seems to be required for the electron transfer process to T1Cu.

Binding free energy calculations between bovine β-lactoglobulin and four fatty acids using the MMGBSA method.

The bovine dairy protein β‐lactoglobulin (βlg) is a promiscuous protein that has the ability to bind several hydrophobic ligands. In this study, based on known experimental data, the dynamic interaction mechanism between bovine βlg and four fatty acids was investigated by a protocol combining molecular dynamics (MD) simulations and molecular mechanics generalized Born surface area (MMGBSA) binding free energy calculations. Energetic analyses revealed binding free energy trends that corroborated known experimental findings; larger ligand size corresponded to greater binding affinity. Finally, binding free energy decomposition provided detailed information about the key residues stabilizing the complex.

Molecular Dynamics Simulations to Provide Insights into Epitopes Coupled to the Soluble and Membrane-Bound MHC-II Complexes

Epitope recognition by major histocompatibility complex II (MHC-II) is essential for the activation of immunological responses to infectious diseases. Several studies have demonstrated that this molecular event takes place in the MHC-II peptide-binding groove constituted by the α and β light chains of the heterodimer. This MHC-II peptide-binding groove has several pockets (P1-P11) involved in peptide recognition and complex stabilization that have been probed through crystallographic experiments and in silico calculations. However, most of these theoretical calculations have been performed without taking into consideration the heavy chains, which could generate misleading information about conformational mobility both in water and in the membrane environment. Therefore, in absence of structural information about the difference in the conformational changes between the peptide-free and peptide-bound states (pMHC-II) when the system is soluble in an aqueous environment or non-covalently bound to a cell membrane, as the physiological environment for MHC-II is. In this study, we explored the mechanistic basis of these MHC-II components using molecular dynamics (MD) simulations in which MHC-II was previously co-crystallized with a small epitope (P7) or coupled by docking procedures to a large (P22) epitope. These MD simulations were performed at 310 K over 100 ns for the water-soluble (MHC-IIw, MHC-II-P7w, and MHC-II-P22w) and 150 ns for the membrane-bound species (MHC-IIm, MHC-II-P7m, and MHC-II-P22m). Our results reveal that despite the different epitope sizes and MD simulation environments, both peptides are stabilized primarily by residues lining P1, P4, and P6-7, and similar noncovalent intermolecular energies were observed for the soluble and membrane-bound complexes. However, there were remarkably differences in the conformational mobility and intramolecular energies upon complex formation, causing some differences with respect to how the two peptides are stabilized in the peptide-binding groove.

Ligand entry into the calyx of β‐lactoglobulin

Although the thermodynamic principles that control the binding of drug molecules to their protein targets are well understood, the detailed process of how a ligand reaches a protein binding site has been an intriguing question over decades. The short time interval between the encounter between a ligand and its receptor to the formation of the stable complex has prevented experimental observations. Bovine β‐lactoglobulin (βlg) is a lipocalin member that carries fatty acids (FAs) and other lipids in the cellular environment. Βlg accommodates a FA molecule in its highly hydrophobic cavity and exhibits the capability of recognizing a wide variety of hydrophobic ligands. To elucidate the ligand entry process on βlg, we report molecular dynamics simulations of the encounter between palmitate (PA) or laurate (LA) and βlg. Our results show that residues localized in loops at the cavity entrance play an important role in the ligand penetration process. Analysis of the short‐term interaction energies show that the forces operating on the systems lead to average conformations very close to the crystallographic holo‐forms. Whereas the binding free energy analysis using the molecular mechanics Generalized Born surface area method shows that these conformations were thermodynamically favorable.

Predicting peptide vaccine candidates against H1N1 influenza virus through theoretical approaches

Identification of potential epitopes that might activate the immune system has been facilitated by the employment of algorithms that use experimental data as templates. However, in order to prove the affinity and the map of interactions between the receptor (major histocompatibility complex, MHC, or T-cell receptor) and the potential epitope, further computational studies are required. Docking and molecular dynamics (MDs) simulations have been an effective source of generating structural information at molecular level in immunology. Herein, in order to provide a detailed understanding of the origins of epitope recognition and to select the best peptide candidate to develop an epitope-based vaccine, docking and MDs simulations in combination with MMGBSA free energy calculations and per-residue free energy decomposition were performed, taking as starting complexes those formed between four designed epitopes (P1–P4) from hemagglutinin (HA) of the H1N1 influenza virus and MHC-II anchored in POPC membrane. Our results revealed that the energetic contributions of individual amino acids within the pMHC-II complexes are mainly dictated by van der Waals interactions and the nonpolar part of solvation energy, whereas the electrostatic interactions corresponding to hydrogen bonds and salt bridges determine the binding specificity, being the most favorable interactions formed between p4 and MHC-II. Then, P1–P4 epitopes were synthesized and tested experimentally to compare theoretical and experimental results. Experimental results show that P4 elicited the highest strong humoral immune response to HA of the H1N1 and may induce antibodies that are cross-reactive to other influenza subtypes, suggesting that it could be a good candidate for the development of a peptide-based vaccine.

Theoretical analysis of the neuraminidase epitope of the Mexican A H1N1 influenza strain, and experimental studies on its interaction with rabbit and human hosts

The neuraminidase (NA) epitope from the Mexican AH1N1 influenza virus was identified by using sequences registered at the GenBank during the peak of a pandemic (from April 2009 to October 2010). First, NA protein sequences were submitted for multiple alignment analysis, and their three-dimensional models (3-D) were then built by using homology modeling. The most common sequence (denominated wild-type) and its mutants were submitted to linear and nonlinear epitope predictors, which included the major histocompatibility complex type II (MHC II) and B-cell peptides. The epitope prediction was in accordance with evolutionary behavior and some protein structural properties. The latter included a low NA mutation rate, NA 3-D surface exposure, and the presence of high hindrance side chain residues. After selecting the epitope, docking studies and molecular dynamics (MD) simulations were used to explore interactions between the epitope and MHC II. Afterward, several experimental assays were performed to validate the theoretical study by using antibodies from humans (infected by pandemic H1N1) and rabbits (epitope vaccination). The results show 119 complete sequences that were grouped into 28 protein sequences according to their identity (one wild-type and 27 representative mutants (1–5 mutations)). The predictors yielded several epitopes, with the best fit being the one located in the C-terminal region. Theoretical methods demonstrated that the selected epitope reached the P4, P6, P7, and P9 pockets of MHC II, whereas the experimental evidence indicates that the epitope is recognized by human antibodies and also by rabbit antibodies immunized with the peptide.