Adhesive water networks facilitate binding of protein interfaces

Ligand Gaussian Accelerated Molecular Dynamics 3 (LiGaMD3): Improved Calculations of Binding Thermodynamics and Kinetics of Both Small Molecules and Flexible Peptides

Binding thermodynamics and kinetics play critical roles in drug design. However, it has proven challenging to efficiently predict ligand binding thermodynamics and kinetics of small molecules and flexible peptides using conventional molecular dynamics (cMD), due to limited simulation time scales. Based on our previously developed ligand Gaussian accelerated molecular dynamics (LiGaMD) method, we present a new approach, termed "LiGaMD3″, in which we introduce triple boosts into three individual energy terms that play important roles in small-molecule/peptide dissociation, rebinding, and system conformational changes to improve the sampling efficiency of small-molecule/peptide interactions with target proteins. To validate the performance of LiGaMD3, MDM2 bound by a small molecule (Nutlin 3) and two highly flexible peptides (PMI and P53) were chosen as the model systems. LiGaMD3 could efficiently capture repetitive small-molecule/peptide dissociation and binding events within 2 μs simulations. The predicted binding kinetic constant rates and free energies from LiGaMD3 were in agreement with the available experimental values and previous simulation results. Therefore, LiGaMD3 provides a more general and efficient approach to capture dissociation and binding of both small-molecule ligands and flexible peptides, allowing for accurate prediction of their binding thermodynamics and kinetics.

Ligand Gaussian accelerated Molecular Dynamics 3 (LiGaMD3): Improved Calculations of Binding Thermodynamics and Kinetics of Both Small Molecules and Flexible Peptides

Binding thermodynamics and kinetics play critical roles in drug design. However, it has proven challenging to efficiently predict ligand binding thermodynamics and kinetics of small molecules and flexible peptides using conventional Molecular Dynamics (cMD), due to limited simulation timescales. Based on our previously developed Ligand Gaussian accelerated Molecular Dynamics (LiGaMD) method, we present a new approach, termed "LiGaMD3", in which we introduce triple boosts into three individual energy terms that play important roles in small-molecule/peptide dissociation, rebinding and system conformational changes to improve the sampling efficiency of small-molecule/peptide interactions with target proteins. To validate the performance of LiGaMD3, MDM2 bound by a small molecule (Nutlin 3) and two highly flexible peptides (PMI and P53) were chosen as model systems. LiGaMD3 could efficiently capture repetitive small-molecule/peptide dissociation and binding events within 2 microsecond simulations. The predicted binding kinetic constant rates and free energies from LiGaMD3 agreed with available experimental values and previous simulation results. Therefore, LiGaMD3 provides a more general and efficient approach to capture dissociation and binding of both small-molecule ligand and flexible peptides, allowing for accurate prediction of their binding thermodynamics and kinetics.

Protein-based bioactive coatings: from nanoarchitectonics to applications

Protein-based bioactive coatings have emerged as a versatile and promising strategy for enhancing the performance and biocompatibility of diverse biomedical materials and devices. Through surface modification, these coatings confer novel biofunctional attributes, rendering the material highly bioactive. Their widespread adoption across various domains in recent years underscores their importance. This review systematically elucidates the behavior of protein-based bioactive coatings in organisms and expounds on their underlying mechanisms. Furthermore, it highlights notable advancements in artificial synthesis methodologies and their functional applications in vitro. A focal point is the delineation of assembly strategies employed in crafting protein-based bioactive coatings, which provides a guide for their expansion and sustained implementation. Finally, the current trends, challenges, and future directions of protein-based bioactive coatings are discussed.

The Advances and Limitations of the Determination and Applications of Water Structure in Molecular Engineering

Water is a key actor of various processes of nature and, therefore, molecular engineering has to take the structural and energetic consequences of hydration into account. While the present review focuses on the target-ligand interactions in drug design, with a focus on biomolecules, these methods and applications can be easily adapted to other fields of the molecular engineering of molecular complexes, including solid hydrates. The review starts with the problems and solutions of the determination of water structures. The experimental approaches and theoretical calculations are summarized, including conceptual classifications. The implementations and applications of water models are featured for the calculation of the binding thermodynamics and computational ligand docking. It is concluded that theoretical approaches not only reproduce or complete experimental water structures, but also provide key information on the contribution of individual water molecules and are indispensable tools in molecular engineering.

Interfacial water molecules contribute to antibody binding to the receptor-binding domain of SARS-CoV-2 spike protein

Antibodies that recognize the spike protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), especially the neutralizing antibodies, carry great hope in the treatment and final elimination of COVID-19. Driven by a synchronized global effort, thousands of antibodies against the spike protein have been identified during the past two years, with the structural information available at atomistic detail for hundreds of these antibodies. We developed an improved molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) method including explicitly treated interfacial water to calculate the binding free energy between representative antibodies and the receptor binding domain (RBD) domain of SARS-COV-2 spike proteins. We discovered that explicit treatment of water molecules located at the interface between RBD and antibody effectively improves the results for the WT and variants of concern (VOC) systems. Interfacial water molecules, together with surface and internal water molecules, behave drastically from bulk water and exert peculiar impacts on protein dynamics and energy, and thus warrant explicit treatment to complement implicit solvent models. Our results illustrate the importance of including interfacial water molecules to approach efficient and reliable prediction of binding free energy.Communicated by Ramaswamy H. Sarma.

Total free energy analysis of fully hydrated proteins

The total free energy of a hydrated biomolecule and its corresponding decomposition of energy and entropy provides detailed information about regions of thermodynamic stability or instability. The free energies of four hydrated globular proteins with different net charges are calculated from a molecular dynamics simulation, with the energy coming from the system Hamiltonian and entropy using multiscale cell correlation. Water is found to be most stable around anionic residues, intermediate around cationic and polar residues, and least stable near hydrophobic residues, especially when more buried, with stability displaying moderate entropy-enthalpy compensation. Conversely, anionic residues in the proteins are energetically destabilized relative to singly solvated amino acids, while trends for other residues are less clear-cut. Almost all residues lose intraresidue entropy when in the protein, enthalpy changes are negative on average but may be positive or negative, and the resulting overall stability is moderate for some proteins and negligible for others. The free energy of water around single amino acids is found to closely match existing hydrophobicity scales. Regarding the effect of secondary structure, water is slightly more stable around loops, of intermediate stability around β strands and turns, and least stable around helices. An interesting asymmetry observed is that cationic residues stabilize a residue when bonded to its N-terminal side but destabilize it when on the C-terminal side, with a weaker reversed trend for anionic residues.

Antibody CDR amino acids underlying the functionality of antibody repertoires in recognizing diverse protein antigens

Antibodies recognize protein antigens with exquisite specificity in a complex aqueous environment, where interfacial waters are an integral part of the antibody–protein complex interfaces. In this work, we elucidate, with computational analyses, the principles governing the antibodies’ specificity and affinity towards their cognate protein antigens in the presence of explicit interfacial waters. Experimentally, in four model antibody–protein complexes, we compared the contributions of the interaction types in antibody–protein antigen complex interfaces with the antibody variants selected from phage-displayed synthetic antibody libraries. Evidently, the specific interactions involving a subset of aromatic CDR (complementarity determining region) residues largely form the predominant determinant underlying the specificity of the antibody–protein complexes in nature. The interfacial direct/water-mediated hydrogen bonds accompanying the CDR aromatic interactions are optimized locally but contribute little in determining the epitope location. The results provide insights into the phenomenon that natural antibodies with limited sequence and structural variations in an antibody repertoire can recognize seemingly unlimited protein antigens. Our work suggests guidelines in designing functional artificial antibody repertoires with practical applications in developing novel antibody-based therapeutics and diagnostics for treating and preventing human diseases.

A dynamical view of protein-protein complexes: Studies by molecular dynamics simulations

Protein-protein interactions are at the basis of many protein functions, and the knowledge of 3D structures of protein-protein complexes provides structural, mechanical and dynamical pieces of information essential to understand these functions. Protein-protein interfaces can be seen as stable, organized regions where residues from different partners form non-covalent interactions that are responsible for interaction specificity and strength. They are commonly described as a peripheral region, whose role is to protect the core region that concentrates the most contributing interactions, from the solvent. To get insights into the dynamics of protein-protein complexes, we carried out all-atom molecular dynamics simulations in explicit solvent on eight different protein-protein complexes of different functional class and interface size by taking into account the bound and unbound forms. On the one hand, we characterized structural changes upon binding of the proteins, and on the other hand we extensively analyzed the interfaces and the structural waters involved in the binding. Based on our analysis, in 6 cases out of 8, the interfaces rearranged during the simulation time, in stable and long-lived substates with alternative residue-residue contacts. These rearrangements are not restricted to side-chain fluctuations in the periphery but also affect the core interface. Finally, the analysis of the waters at the interface and involved in the binding pointed out the importance to take into account their role in the estimation of the interaction strength.

Protein-Protein Interaction-Gaussian Accelerated Molecular Dynamics (PPI-GaMD): Characterization of Protein Binding Thermodynamics and Kinetics

Protein-protein interactions (PPIs) play key roles in many fundamental biological processes such as cellular signaling and immune responses. However, it has proven challenging to simulate repetitive protein association and dissociation in order to calculate binding free energies and kinetics of PPIs due to long biological timescales and complex protein dynamics. To address this challenge, we have developed a new computational approach to all-atom simulations of PPIs based on a robust Gaussian accelerated molecular dynamics (GaMD) technique. The method, termed "PPI-GaMD", selectively boosts interaction potential energy between protein partners to facilitate their slow dissociation. Meanwhile, another boost potential is applied to the remaining potential energy of the entire system to effectively model the protein's flexibility and rebinding. PPI-GaMD has been demonstrated on a model system of the ribonuclease barnase interactions with its inhibitor barstar. Six independent 2 μs PPI-GaMD simulations have captured repetitive barstar dissociation and rebinding events, which enable calculations of the protein binding thermodynamics and kinetics simultaneously. The calculated binding free energies and kinetic rate constants agree well with the experimental data. Furthermore, PPI-GaMD simulations have provided mechanistic insights into barstar binding to barnase, which involves long-range electrostatic interactions and multiple binding pathways, being consistent with previous experimental and computational findings of this model system. In summary, PPI-GaMD provides a highly efficient and easy-to-use approach for binding free energy and kinetics calculations of PPIs.

Probing the Bovine Hemoglobin Adsorption Process and its Influence on Interfacial Water Structure at the Air-Water Interface

*These authors contributed equally to this work.The molecular-level insight of protein adsorption and its kinetics at interfaces is crucial because of its multifold role in diverse fundamental biological processes and applications. In the present study, the sum frequency generation (SFG) vibrational spectroscopy has been employed to demonstrate the adsorption process of bovine hemoglobin (BHb) protein molecules at the air-water interface at interfacial isoelectric point of the protein. It has been observed that surface coverage of BHb molecules significantly influences the arrangement of the protein molecules at the interface. The time-dependent SFG studies at two different frequencies in the fingerprint region elucidate the kinetics of protein denaturation process and its influence on the hydrogen-bonding network of interfacial water molecules at the air-water interface. The initial growth kinetics suggests the synchronized behavior of protein adsorption process with the structural changes in the interfacial water molecules. Interestingly, both the events carry similar characteristic time constants. However, the conformational changes in the protein structure due to the denaturation process stay for a long time, whereas the changes in water structure reconcile quickly. It is revealed that the protein denaturation process is followed by the advent of strongly hydrogen-bonded water molecules at the interface. In addition, we have also carried out the surface tension kinetics measurements to complement the findings of our SFG spectroscopic results.

Native or Non-Native Protein-Protein Docking Models? Molecular Dynamics to the Rescue

Molecular docking excels at creating a plethora of potential models of protein-protein complexes. To correctly distinguish the favorable, native-like models from the remaining ones remains, however, a challenge. We assessed here if a protocol based on molecular dynamics (MD) simulations would allow distinguishing native from non-native models to complement scoring functions used in docking. To this end, the first models for 25 protein-protein complexes were generated using HADDOCK. Next, MD simulations complemented with machine learning were used to discriminate between native and non-native complexes based on a combination of metrics reporting on the stability of the initial models. Native models showed higher stability in almost all measured properties, including the key ones used for scoring in the Critical Assessment of PRedicted Interaction (CAPRI) competition, namely the positional root mean square deviations and fraction of native contacts from the initial docked model. A random forest classifier was trained, reaching a 0.85 accuracy in correctly distinguishing native from non-native complexes. Reasonably modest simulation lengths of the order of 50-100 ns are sufficient to reach this accuracy, which makes this approach applicable in practice.

Structure of autoinhibited Akt1 reveals mechanism of PIP3-mediated activation

The protein kinase Akt is one of the primary effectors of growth factor signaling in the cell. Akt responds specifically to the lipid second messengers phosphatidylinositol-3,4,5-trisphosphate [PI(3,4,5)P3] and phosphatidylinositol-3,4-bisphosphate [PI(3,4)P2] via its PH domain, leading to phosphorylation of its activation loop and the hydrophobic motif of its kinase domain, which are critical for activity. We have now determined the crystal structure of Akt1, revealing an autoinhibitory interface between the PH and kinase domains that is often mutated in cancer and overgrowth disorders. This interface persists even after stoichiometric phosphorylation, thereby restricting maximum Akt activity to PI(3,4,5)P3- or PI(3,4)P2-containing membranes. Our work helps to resolve the roles of lipids and phosphorylation in the activation of Akt and has wide implications for the spatiotemporal control of Akt and potentially lipid-activated kinase signaling in general.

Residue interaction networks of K-Ras protein with water molecules identifies the potential role of switch II and P-loop

The mutant K-Ras with aberrant signaling is the primary cause of several cancers. The proposed study investigated the influence of water molecules in K-Ras crystal structure, where they have a significant function by understanding their residue interaction networks (RINs). We analyzed the RINs of K-Ras with and without water molecules and determined their interaction properties. RINs were developed with the help of StructureViz2 and RINspector; further, the changes in K-Ras backbone flexibility were predicted with the DynaMine. We found that the residues K42, I142, and L159 are the hotspots from water, including the K-Ras-GTP complex with the highest residue centrality analysis (RCA) Z-score. The DynaMine prediction calculated the NMR S² value for the frequently mutated positions G12, G13, and Q61 showing a minor shift in flexibility, which make up the P-Loop and switch II of the K-Ras protein. This flexibility shift can account for changes in conformational activity and the protein's GTPase activity, making it difficult to recognize by the effectors and exchange factors. Taken together, our study helps in understanding the functional importance of the water molecules in K-Ras protein and the impact of mutation that modulate the conformational state of the protein.

Nonspecific Binding-Fundamental Concepts and Consequences for Biosensing Applications

Nature achieves differentiation of specific and nonspecific binding in molecular interactions through precise control of biomolecules in space and time. Artificial systems such as biosensors that rely on distinguishing specific molecular binding events in a sea of nonspecific interactions have struggled to overcome this issue. Despite the numerous technological advancements in biosensor technologies, nonspecific binding has remained a critical bottleneck due to the lack of a fundamental understanding of the phenomenon. To date, the identity, cause, and influence of nonspecific binding remain topics of debate within the scientific community. In this review, we discuss the evolution of the concept of nonspecific binding over the past five decades based upon the thermodynamic, intermolecular, and structural perspectives to provide classification frameworks for biomolecular interactions. Further, we introduce various theoretical models that predict the expected behavior of biosensors in physiologically relevant environments to calculate the theoretical detection limit and to optimize sensor performance. We conclude by discussing existing practical approaches to tackle the nonspecific binding challenge in vitro for biosensing platforms and how we can both address and harness nonspecific interactions for in vivo systems.

Coulomb interactions between dipolar quantum fluctuations in van der Waals bound molecules and materials

Mutual Coulomb interactions between electrons lead to a plethora of interesting physical and chemical effects, especially if those interactions involve many fluctuating electrons over large spatial scales. Here, we identify and study in detail the Coulomb interaction between dipolar quantum fluctuations in the context of van der Waals complexes and materials. Up to now, the interaction arising from the modification of the electron density due to quantum van der Waals interactions was considered to be vanishingly small. We demonstrate that in supramolecular systems and for molecules embedded in nanostructures, such contributions can amount to up to 6 kJ/mol and can even lead to qualitative changes in the long-range van der Waals interaction. Taking into account these broad implications, we advocate for the systematic assessment of so-called Dipole-Correlated Coulomb Singles in large molecular systems and discuss their relevance for explaining several recent puzzling experimental observations of collective behavior in nanostructured materials.

Confinement-Induced Enhancement of Parallel Dielectric Permittivity: Super Permittivity Under Extreme Confinement

Enhancement of parallel (x-y plane) dielectric permittivity of confined fluids has been shown previously. However, a theoretical model that explains this enhancement is lacking thus far. In this study, using statistical-mechanical theories and molecular dynamics simulations, we show an explicit relation between the parallel dielectric permittivity, density variations, and dipolar correlations for protic and aprotic fluids confined in slit-like channels. We analyze the importance of dipolar correlations on enhancement of parallel dielectric permittivity inside large channels and extreme confinements. In large channels, beyond the interfacial region, dipolar correlations exhibit bulk-like behavior. Under extreme confinement, the correlations become stronger to the extent that they give rise to a giant increase in the parallel dielectric permittivity. This sudden increase in dielectric permittivity can be a signature of a liquid transition into higher-ordered structures and has important consequences for understanding ion transport, molecular dissociation, and chemical reactions inside nanoconfined environments.

Computational design of small molecular modulators of protein-protein interactions with a novel thermodynamic cycle: Allosteric inhibitors of HIV-1 integrase

Targeting protein-protein interactions for therapeutic development involves designing small molecules to either disrupt or enhance a known PPI. For this purpose, it is necessary to compute reliably the effect of chemical modifications of small molecules on the protein-protein association free energy. Here we present results obtained using a novel thermodynamic free energy cycle, for the rational design of allosteric inhibitors of HIV-1 integrase (ALLINI) that act specifically in the early stage of the infection cycle. The new compounds can serve as molecular probes to dissect the multifunctional mechanisms of ALLINIs to inform the discovery of new allosteric inhibitors. The free energy protocol developed here can be more broadly applied to study quantitatively the effects of small molecules on modulating the strengths of protein-protein interactions.

Exploring the Free-Energy Landscape and Thermodynamics of Protein-Protein Association

We report the free-energy landscape and thermodynamics of the protein-protein association responsible for the drug-induced multimerization of HIV-1 integrase (IN). Allosteric HIV-1 integrase inhibitors promote aberrant IN multimerization by bridging IN-IN intermolecular interactions. However, the thermodynamic driving forces and kinetics of the multimerization remain largely unknown. Here, we explore the early steps in the IN multimerization by using umbrella sampling and unbiased molecular dynamics simulations in explicit solvent. In direct simulations, the two initially separated dimers spontaneously associate to form near-native complexes that resemble the crystal structure of the aberrant tetramer. Most strikingly, the effective interaction of the protein-protein association is very short-ranged: the two dimers associate rapidly within tens of nanoseconds when their binding surfaces are separated by d ≤ 4.3 Å (less than two water diameters). Beyond this distance, the oligomerization kinetics appears to be diffusion controlled with a much longer association time. The free-energy profile also captured the crucial role of allosteric IN inhibitors in promoting multimerization and explained why several C-terminal domain mutations are remarkably resistant to the drug-induced multimerization. The results also show that at small separation, the protein-protein binding process contains two consecutive phases with distinct thermodynamic signatures. First, interprotein water molecules are expelled to the bulk, resulting in a small increase in entropy, as the solvent entropy gain from the water release is nearly cancelled by the loss of side-chain entropies as the two proteins approach each other. At shorter distances, the two dry binding surfaces adapt to each other to optimize their interaction energy at the expense of further protein configurational entropy loss. Although the binding interfaces feature clusters of hydrophobic residues, overall, the protein-protein association in this system is driven by enthalpy and opposed by entropy.

Regulation of Thrombin Activity with Ultrasmall Nanoparticles: Effects of Surface Chemistry

Nanomaterials displaying well-tailored sizes and surface chemistries can provide novel ways with which to modulate the structure and function of enzymes. Recently, we showed that gold nanoparticles (AuNPs) in the ultrasmall size regime could perform as allosteric effectors inducing partial inhibition of thrombin activity. We now find that the nature of the AuNP surface chemistry controls the interactions to the anion-binding exosites 1 and 2 on the surface of thrombin, the allosterically induced changes to the active-site conformation, and, by extension, the enzymatic activity. Ultrasmall AuNPs passivated with p-mercaptobenzoic acid ligands (AuMBA) and a peptide-based (Ac-ECYN) biomimetic coat (AuECYN) were utilized in our investigations. Remarkably, we found that while AuMBA binds to exosites 1 and 2, AuECYN interacts primarily with exosite 2. It was further established that AuMBA behaves as a "mild denaturant" of thrombin leading to catalytic dysfunction over time. Conversely, AuECYN resembles a proper allosteric effector leading to partial and reversible inhibition of the activity. Collectively, our findings reveal how the distinct binding modes of different AuNP types may uniquely influence thrombin structure and catalysis. The present study further contributes to our understanding of how synthetic nanomaterials could be exploited in the allosteric regulation of enzymes.

Stretching Single Collagen Fibrils Reveals Nonlinear Mechanical Behavior

The mechanical properties of collagen fibrils play an important role in cell-matrix interactions and are a manifestation of their molecular structure. Using a, to our knowledge, novel combination of uniaxial, longitudinal straining and radial nanoindentation, we found that type I collagen fibrils show a pronounced nonlinear behavior in the form of strain stiffening at strains from 0 to 15%, followed by strain softening at strains from 15 to 25%. At the molecular scale, this surprising phenomenon can be explained by the combination of unfolding of disordered domains and breaking of native cross-links at different stages of strain. Fibrils cross-linked artificially by glutaraldehyde do not show such a behavior, and nanoindentation allowed us to measure the mechanics of the overlap and gap regions in the D-banding individually. The results could have consequences for our understanding of matrix mechanics and the influence of excessive glycation, which has been linked with age-related diseases such as diabetes. Furthermore, the simplicity of the straining method could be attractive in other areas of biophysics at the nanometer scale because it does not require any bespoke instrumentation and is easy to use.

Explicit Characterization of the Free Energy Landscape of pKID-KIX Coupled Folding and Binding

The most fundamental aspect of the free energy landscape of proteins is that it is globally funneled such that protein folding is energetically biased. Then, what are the distinctive characteristics of the landscape of intrinsically disordered proteins, apparently lacking such energetic bias, that nevertheless fold upon binding? Here, we address this fundamental issue through the explicit characterization of the free energy landscape of the paradigmatic pKID-KIX system (pKID, phosphorylated kinase-inducible domain; KIX, kinase interacting domain). This is done based on unguided, fully atomistic, explicit-water molecular dynamics simulations with an aggregated simulation time of >30 μs and on the computation of the free energy that defines the landscape. We find that, while the landscape of pKID before binding is considerably shallower than the one for a protein that autonomously folds, it becomes progressively more funneled as the binding of pKID with KIX proceeds. This explains why pKID is disordered in a free state, and the binding of pKID with KIX is a prerequisite for pKID's folding. In addition, we observe that the key event in completing the pKID-KIX coupled folding and binding is the directed self-assembly where pKID is docked upon the KIX surface to maximize the surface electrostatic complementarity, which, in turn, require pKID to adopt the correct folded structure. This key process shows up as the free energy barrier in the pKID landscape separating the intermediate nonspecific complex state and the specific complex state. The present work not only provides a detailed molecular picture of the coupled folding and binding of pKID but also expands the funneled landscape perspective to intrinsically disordered proteins.

Nanosecond pulsed electric signals can affect electrostatic environment of proteins below the threshold of conformational effects: The case study of SOD1 with a molecular simulation study

Electric fields can be a powerful tool to interact with enzymes or proteins, with an intriguing perspective to allow protein manipulation. Recently, researchers have focused the interest on intracellular enzyme modifications triggered by the application of nanosecond pulsed electric fields. These findings were also supported by theoretical predictions from molecular dynamics simulations focussing on significant variations in protein secondary structures. In this work, a theoretical study utilizing molecular dynamics simulations is proposed to explore effects of electric fields of high intensity and very short nanosecond duration applied to the superoxide dismutase (Cu/Zn-SOD or SOD-1), an important enzyme involved in the cellular antioxidant defence mechanism. The effects of 100-nanosecond pulsed electric fields, with intensities ranging from 108 to 7x108 V/m, on a single SOD1 enzyme are presented. We demonstrated that the lowest intensity of 108 V/m, although not inducing structural changes, can produce electrostatic modifications on the reaction centre of the enzyme, as apparent from the dipolar response and the electric field distribution of the protein active site. Electric pulses above 5x108 V/m produced a fast transition between the folded and a partially denatured state, as inferred by the secondary structures analysis. Finally, for the highest field intensity used (7x108 V/m), a not reversible transition toward an unfolded state was observed.

Intramolecular Domain Movements of Free and Bound pMHC and TCR Proteins: A Molecular Dynamics Simulation Study

The interaction of antigenic peptides (p) and major histocompatibility complexes (pMHC) with T-cell receptors (TCR) is one of the most important steps during the immune response. Here we present a molecular dynamics simulation study of bound and unbound TCR and pMHC proteins of the LC13-HLA-B*44:05-pEEYLQAFTY complex to monitor differences in relative orientations and movements of domains between bound and unbound states of TCR-pMHC. We generated local coordinate systems for MHC α1- and MHC α2-helices and the variable T-cell receptor regions TCR V_α and TCR V_β and monitored changes in the distances and mutual orientations of these domains. In comparison to unbound states, we found decreased inter-domain movements in the simulations of bound states. Moreover, increased conformational flexibility was observed for the MHC α2-helix, the peptide, and for the complementary determining regions of the TCR in TCR-unbound states as compared to TCR-bound states.

Applications of in Silico Methods for Design and Development of Drugs Targeting Protein-Protein Interactions

BACKGROUND

Identification of Protein-Protein Interactions (PPIs) is a major challenge in modern molecular biology and biochemistry research, due to the unquestionable role of proteins in cells, biological process and pathological states. Over the past decade, the PPIs have evolved from being considered a highly challenging field of research to being investigated and examined as targets for pharmacological intervention.

OBJECTIVE

Comprehension of protein interactions is crucial to known how proteins come together to build signalling pathways, to carry out their functions, or to cause diseases, when deregulated. Multiplicity and great amount of PPIs structures offer a huge number of new and potential targets for the treatment of different diseases.

METHODS

Computational techniques are becoming predominant in PPIs studies for their effectiveness, flexibility, accuracy and cost. As a matter of fact, there are effective in silico approaches which are able to identify PPIs and PPI site. Such methods for computational target prediction have been developed through molecular descriptors and data-mining procedures.

RESULTS

In this review, we present different types of interactions between protein-protein and the application of in silico methods for design and development of drugs targeting PPIs. We described computational approaches for the identification of possible targets on protein surface and to detect of stimulator/ inhibitor molecules.

CONCLUSION

A deeper study of the most recent bioinformatics methodologies for PPIs studies is vital for a better understanding of protein complexes and for discover new potential PPI modulators in therapeutic intervention.

Using Graphs of Dynamic Hydrogen-Bond Networks To Dissect Conformational Coupling in a Protein Motor

DExD/H-box proteins are soluble enzymes that couple binding and hydrolysis of adenosine triphosphate (ATP) with reactions involving RNA metabolism or bind and push newly synthesized proteins across bacterial cell membranes. Knowledge of the reaction mechanism of these enzymes could help the development of new therapeutics. In order to explore the mechanism of long-distance conformational coupling in SecA, the DEAD-box motor of the Sec protein secretion in bacteria, we implemented algorithms that provide simplified graph representations of the protein's dynamic hydrogen-bond networks. We find that mutations near the nucleotide-binding site or changes of the nucleotide-binding state of SecA associate with altered dynamics at the preprotein binding domain and identify extended networks of hydrogen bonds that connect the active site of SecA to the region where SecA binds newly synthesized secretory proteins. Water molecules participate in hydrogen-bonded water chains that bridge functional domains of SecA and could contribute to long-distance conformational coupling.

Assembly of peptides in mica-graphene nanocapillaries controlled by confined water

Water in nanoscale-confined geometries has unique physicochemical properties in contrast to bulk water, and is believed to play important roles in biological processes although there is less direct information available in the literature. Here, we report the self-assembly behaviors of a neurodegenerative disease related peptide termed GAV-9 encapsulated in mica-graphene nanocapillaries interacting with water nanofilms condensed under ambient conditions, based on atomic force microscopy (AFM) imaging and molecular dynamics (MD) simulations. The results revealed that, upon increase in the humidity, the GAV-9 peptide monomers adsorbed the confined water molecules and transitioned to unexpected hydrogel-like structures. Our MD simulations also suggested that in the confined mica-graphene nanocapillaries, the GAV-9 peptide monomers would indeed form water-rich hydrogel structures instead of highly ordered nanofilaments. The interfacial water confined in the mica-graphene nanocapillary is found to be crucial for such a transition. Moreover, the distribution of confined water layers largely depended on the locations of the preformed peptide nanofilaments, and the peptide nanofilaments further assembled into nanosheets with the water layer increasing, but depolymerized to amorphous peptide assemblies with the water layer decreasing. The polymerization and depolymerization of the peptide nanofilaments could be controlled in a reversible manner. Our results have supplied a simplified model system to uncover the effects of the confined interfacial water on the biological process at the molecular level.

Structural mechanism for Bruton's tyrosine kinase activation at the cell membrane

Bruton’s tyrosine kinase (Btk) activation on the cell membrane is critical for B cell proliferation and development, and Btk inhibition is a promising treatment for several hematologic cancers and autoimmune diseases. Here, we examine Btk activation using the results of long-timescale molecular dynamics simulations. In our simulations, Btk lipid-binding modules dimerized on the membrane in a single predominant conformation. We observed that the phospholipid PIP₃—in addition to its expected role of recruiting Btk to the membrane—allosterically mediated dimer formation and stability by binding at two novel sites. Our results provide strong evidence that PIP₃-mediated dimerization of Btk at the cell membrane is a critical step in Btk activation and suggest a potential approach to allosteric Btk inhibitor development.

Bruton’s tyrosine kinase (Btk) is critical for B cell proliferation and activation, and the development of Btk inhibitors is a vigorously pursued strategy for the treatment of various B cell malignancies. A detailed mechanistic understanding of Btk activation has, however, been lacking. Here, inspired by a previous suggestion that Btk activation might depend on dimerization of its lipid-binding PH–TH module on the cell membrane, we performed long-timescale molecular dynamics simulations of membrane-bound PH–TH modules and observed that they dimerized into a single predominant conformation. We found that the phospholipid PIP₃ stabilized the dimer allosterically by binding at multiple sites, and that the effects of PH–TH mutations on dimer stability were consistent with their known effects on Btk activity. Taken together, our simulation results strongly suggest that PIP₃-mediated dimerization of Btk at the cell membrane is a critical step in Btk activation.

Thick Two-Dimensional Water Film Confined between the Atomically Thin Mica Nanosheet and Hydrophilic Substrate

The interesting properties of water molecules confined in a two-dimensional (2D) environment have aroused great attention. However, the study of 2D-confined water at the hydrophilic-hydrophilic interface is largely unexplored due to the lack of appropriate system. In this work, the behavior of water molecules confined between an atomically thin mica nanosheet and a hydrophilic SiO₂/Si substrate was investigated using an atomic force microscope in detail at ambient conditions. The confined water molecules aggregated as droplets when the relative humidity (RH) of the environment was 11%. A large-area 2D water film with a uniform thickness of ∼2 nm was observed when the mica flake was incubated at 33% RH for 1 h before being mechanically exfoliated on a SiO₂/Si substrate. Interestingly, the water film showed ordered edges with a predominant angle of 120°, which was the same with the lattice orientation of the mica nanosheet on top of it. The water film showed a fluidic behavior at the early stage and reached a stable state after 48 h under ambient conditions. The surface properties of the upper mica nanosheet and the underlying substrate played a crucial role in manipulating the behavior of confined water molecules. When the surface of the upper mica nanosheet was modified by Na⁺, Ni²⁺_, and aminopropyltriethoxysilane (APS), only some small water droplets were observed instead of a water film. The surface of the underlying SiO₂/Si substrate was functionalized by hydrophilic APS and hydrophobic octadecyltrimethoxysiliane (OTS). The small water droplets were imaged on a hydrophobic OTS-SiO₂/Si substrate, while the water film with regular edges was maintained on a hydrophilic APS-SiO₂/Si substrate. Our results might provide an alternative molecular view for investigating structures and properties of confined water molecules in 2D environments.

Collective Transformation of Water between Hyperactive Antifreeze Proteins: RiAFPs

We demonstrate, by molecular dynamics simulations, that water confined between a pair of insect hyperactive antifreeze proteins from the longhorn beetle Rhagium inquisitor is discontinuously expelled as the two proteins approach each other at a certain distance. The extensive striped hydrophobic–hydrophilic pattern on the surface, comprising arrays of threonine residues, enables water to form three independent ice channels through the assistance of hydroxyl groups, even at 300 K. The transformation is reminiscent of a freezing–melting transition rather than a drying transition and governs the stable protein–protein separation in the evaluation of the potential of mean force. The collectivity of water penetration or expulsion and the hysteresis in the time scale of ten nanoseconds predict a potential first-order phase transition at the limit of infinite size and provide a new framework for the water-mediated interaction between solutes.

Unbiased Atomistic Insight into the Mechanisms and Solvent Role for Globular Protein Dimer Dissociation

Association and dissociation of proteins are fundamental processes in nature. Although simple to understand conceptually, the details of the underlying mechanisms and role of the solvent are poorly understood. Here, we investigate the dissociation of the hydrophilic β-lactoglobulin dimer by employing transition path sampling. Analysis of the sampled path ensembles reveals a variety of mechanisms: (1) a direct aligned dissociation (2) a hopping and rebinding transition followed by unbinding, and (3) a sliding transition before unbinding. Reaction coordinate and transition-state analysis predicts that, besides native contact and neighboring salt-bridge interactions, solvent degrees of freedom play an important role in the dissociation process. Bridging waters, hydrogen-bonded to both proteins, support contacts in the native state and nearby lying transition-state regions, whereas they exhibit faster dynamics in further lying transition-state regions, rendering the proteins more mobile and assisting in rebinding. Analysis of the structure and dynamics of the solvent molecules reveals that the dry native interface induces enhanced populations of both disordered hydration water near hydrophilic residues and tetrahedrally ordered hydration water nearby hydrophobic residues. Although not exhaustive, our sampling of rare unbiased reactive molecular dynamics trajectories enhances the understanding of protein dissociation via complex pathways including (multiple) rebinding events.

Atomic-level characterization of protein-protein association

Most proteins associate with other proteins to function, forming complexes that are central to almost all physiological processes. Determining the structures of these complexes and understanding how they associate are problems of fundamental importance. Using long-timescale molecular dynamics simulations, some performed using a new enhanced sampling method, we observed spontaneous association and dissociation of five protein–protein systems to and from their experimentally determined native complexes. By analyzing the simulations of these five systems, which include members of diverse structural and functional classes, we are able to draw general mechanistic conclusions about protein association.

Despite the biological importance of protein–protein complexes, determining their structures and association mechanisms remains an outstanding challenge. Here, we report the results of atomic-level simulations in which we observed five protein–protein pairs repeatedly associate to, and dissociate from, their experimentally determined native complexes using a molecular dynamics (MD)–based sampling approach that does not make use of any prior structural information about the complexes. To study association mechanisms, we performed additional, conventional MD simulations, in which we observed numerous spontaneous association events. A shared feature of native association for these five structurally and functionally diverse protein systems was that if the proteins made contact far from the native interface, the native state was reached by dissociation and eventual reassociation near the native interface, rather than by extensive interfacial exploration while the proteins remained in contact. At the transition state (the conformational ensemble from which association to the native complex and dissociation are equally likely), the protein–protein interfaces were still highly hydrated, and no more than 20% of native contacts had formed.

Hydrophobic Surface Enhances Electrostatic Interaction in Water

A high dielectric constant is one of the peculiar properties of liquid water, indicating that the electrostatic interaction between charged substances is largely reduced in water. We show by molecular dynamics simulation that the dielectric constant of water is decreased near the hydrophobic surface. We further show that the decrease in the dielectric constant is due to both the decreased water density and the reduced water dipole correlation in the direction perpendicular to the surface. We finally demonstrate that electrostatic interaction in water is actually strengthened near the hydrophobic surface.

Molecular dynamics study on the role of solvation water in the adsorption of hyperactive AFP to the ice surface

Using computer simulations, the early stages of the adsorption of the CfAFP molecule to the ice surface were analyzed. We found that the ice and the protein interact at least as early as when the protein is about 1 nm away from the ice surface. These interactions are mediated by interfacial solvation water and are possible thanks to the structural ordering of the solvent. This ordering leads to positional preference of the protein relative to the ice crystal before the final attachment to the ice surface takes place, accompanied by the solidification of the interfacial water. It is possible because the solvation water of the ice-binding plane of CfAFP is susceptible to the overlapping with the solvation water of ice and is mostly changeable into ice itself. These remote interactions significantly increase efficacy of the adsorption process by facilitating the geometric adjustment of the active region of the CfAFP molecule to the ice surface. Because of the ordered nature of the water molecules at the ice-binding plane, the energy of their interactions with the ice-binding surface of the protein does not change upon the ongoing solidification of solvation water. However, the structure of the solvation water is not strictly ice-like and the growth of ice in the interfacial water is not initiated at the side of the protein. On the contrary, we find that solvation water of CfAFP solidifies slower than solvation water of ice - the solidification of interfacial water starts at the surface of ice. Moreover, in the presence of the CfAFP molecule, also solvation water of ice solidifies slower compared to the situation when the protein is not present next to the ice surface. Additionally, the presence of the protein molecule shifts the ratio of cubic to hexagonal ice that spontaneously forms at the ice surface, by introducing another layer of ordered water molecules - opposite to the ice lattice, at the other side of the crystallizing layer of water.

The Inclusion of Water Molecules in Residue Interaction Networks Identifies Additional Central Residues

The relevance of water molecules for the recognition and the interaction of biomolecules is widely appreciated. In this paper we address the role that water molecules associated to protein complexes play for the functional relevance of residues by considering their residue interaction networks (RINs). These are commonly defined on the basis of the amino acid composition of the proteins themselves, disregarding the solvation state of the protein. We determine properties of the RINs of two protein complexes, colicin E2/Im2 and barnase/barstar, with and without associated water molecules, using a previously developed methodology and its associated application RINspector. We find that the inclusion of water molecules in RINs leads to an increase in the number of central residues which adds a novel mechanism to the relevance of water molecules for protein function.

Protein-Assisted Formation and Stabilization of Catalytically Active Polyoxometalate Species

The effect of the protein environment on the formation and stabilization of an elusive catalytically active polyoxometalate (POM) species, K₆ [Hf(α₂ -P₂ W₁₇ O₆₁ )] (1), is reported. In the co-crystal of hen egg-white lysozyme (HEWL) with 1, the catalytically active monomeric species is observed, originating from the dimeric 1:2 POM form, while it is intrinsically unstable under physiological pH conditions. The protein-assisted dissociation of the dimeric POM was rationalized by means of DFT calculations. The dissociation process is unfavorable in bulk water, but becomes favorable in the protein-POM complex due to the low dielectric response at the protein surface. The crystal structure shows that the monomeric form is stabilized by electrostatic and water-mediated hydrogen bonding interactions with the protein. It interacts at three distinct sites, close to the aspartate-containing hydrolysis sites, demonstrating high selectivity towards peptide bonds containing this residue.

An Efficient Implementation of the Nwat-MMGBSA Method to Rescore Docking Results in Medium-Throughput Virtual Screenings

Nwat-MMGBSA is a variant of MM-PB/GBSA based on the inclusion of a number of explicit water molecules that are the closest to the ligand in each frame of a molecular dynamics trajectory. This method demonstrated improved correlations between calculated and experimental binding energies in both protein-protein interactions and ligand-receptor complexes, in comparison to the standard MM-GBSA. A protocol optimization, aimed to maximize efficacy and efficiency, is discussed here considering penicillopepsin, HIV1-protease, and BCL-XL as test cases. Calculations were performed in triplicates on both classic HPC environments and on standard workstations equipped by a GPU card, evidencing no statistical differences in the results. No relevant differences in correlation to experiments were also observed when performing Nwat-MMGBSA calculations on 4 or 1 ns long trajectories. A fully automatic workflow for structure-based virtual screening, performing from library set-up to docking and Nwat-MMGBSA rescoring, has then been developed. The protocol has been tested against no rescoring or standard MM-GBSA rescoring within a retrospective virtual screening of inhibitors of AmpC β-lactamase and of the Rac1-Tiam1 protein-protein interaction. In both cases, Nwat-MMGBSA rescoring provided a statistically significant increase in the ROC AUCs of between 20 and 30%, compared to docking scoring or to standard MM-GBSA rescoring.

Triggering molecular assembly at the mesoscale for advanced Raman detection of proteins in liquid

An advanced optofluidic system for protein detection based on Raman signal amplification via dewetting and molecular gathering within temporary mesoscale assemblies is presented. The evaporation of a microliter volume of protein solution deposited in a circular microwell precisely follows an outward-receding geometry. Herein the combination of liquid withdrawal with intermolecular interactions induces the formation of self-assembled molecular domains at the solid-liquid interface. Through proper control of the evaporation rate, amplitude of the assemblies and time for spectral collection at the liquid edge are extensively raised, resulting in a local enhancement and refinement of the Raman response, respectively. Further signal amplification is obtained by taking advantage of the intense local electromagnetic fields generated upon adding a plasmonic coating to the microwell. Major advantages of this optofluidic method lie in the obtainment of high-quality, high-sensitivity Raman spectra with detection limit down to sub-micromolar values. Peculiarly, the assembled proteins in the liquid edge region maintain their native-like state without displaying spectral changes usually occurring when dried drop deposits are considered.

Enhancing Structure Prediction and Design of Soluble and Membrane Proteins with Explicit Solvent-Protein Interactions

Solvent molecules interact intimately with proteins and can profoundly regulate their structure and function. However, accurately and efficiently modeling protein solvation effects at the molecular level has been challenging. Here, we present a method that improves the atomic-level modeling of soluble and membrane protein structures and binding by efficiently predicting de novo protein-solvent molecule interactions. The method predicted with unprecedented accuracy buried water molecule positions, solvated protein conformations, and challenging mutational effects on protein binding. When applied to homology modeling, solvent-bound membrane protein structures, pockets, and cavities were recapitulated with near-atomic precision even from distant homologs. Blindly refined atomic-level structures of evolutionary distant G protein-coupled receptors imply strikingly different functional roles of buried solvent between receptor classes. The method should prove useful for refining low-resolution protein structures, accurately modeling drug-binding sites in structurally uncharacterized receptors, and designing solvent-mediated protein catalysis, recognition, ligand binding, and membrane protein signaling.

The Role of Water Distribution Controlled by Transmembrane Potentials in the Cytochrome c-Cardiolipin Interaction: Revealing from Surface-Enhanced Infrared Absorption Spectroscopy

The interaction of cytochrome c (cyt c) with cardiolipin (CL) plays a crucial role in apoptotic functions, however, the changes of the transmembrane potential in governing the protein behavior at the membrane-water interface have not been studied due to the difficulties in simultaneously monitoring the interaction and regulating the electric field. Herein, surface-enhanced infrared absorption (SEIRA) spectroelectrochemistry is employed to study the mechanism of how the transmembrane potentials control the interaction of cyt c with CL membranes by regulating the electrode potentials of an Au film. When the transmembrane potential decreases, the water content at the interface of the membranes can be increased to slow down protein adsorption through decreasing the hydrogen-bond and hydrophobic interactions, but regulates the redox behavior of CL-bound cyt c through a possible water-facilitated proton-coupled electron transfer process. Our results suggest that the potential drop-induced restructure of the CL conformation and the hydration state could modify the structure and function of CL-bound cyt c on the lipid membrane.

Signatures of Solvation Thermodynamics in Spectra of Intermolecular Vibrations

This study explores the thermodynamic and vibrational properties of water in the three-dimensional environment of solvated ions and small molecules using molecular simulations. The spectrum of intermolecular vibrations in liquid solvents provides detailed information on the shape of the local potential energy surface, which in turn determines local thermodynamic properties such as the entropy. Here, we extract this information using a spatially resolved extension of the two-phase thermodynamics method to estimate hydration water entropies based on the local vibrational density of states (3D-2PT). Combined with an analysis of solute–water and water–water interaction energies, this allows us to resolve local contributions to the solvation enthalpy, entropy, and free energy. We use this approach to study effects of ions on their surrounding water hydrogen bond network, its spectrum of intermolecular vibrations, and resulting thermodynamic properties. In the three-dimensional environment of polar and nonpolar functional groups of molecular solutes, we identify distinct hydration water species and classify them by their characteristic vibrational density of states and molecular entropies. In each case, we are able to assign variations in local hydration water entropies to specific changes in the spectrum of intermolecular vibrations. This provides an important link for the thermodynamic interpretation of vibrational spectra that are accessible to far-infrared absorption and Raman spectroscopy experiments. Our analysis provides unique microscopic details regarding the hydration of hydrophobic and hydrophilic functional groups, which enable us to identify interactions and molecular degrees of freedom that determine relevant contributions to the solvation entropy and consequently the free energy.

Dispersion Interactions and the Stability of Amine Dimers

Weak, intermolecular interactions in amine dimers were studied by using the combination of a dispersionless density functional and a function that describes the dispersion contribution to the interaction energy. The validity of this method was shown by comparison of structural and energetic properties with data obtained with a conventional density functional and the coupled cluster method. The stability of amine dimers was shown to depend on the size, the shape, and the relative orientation of the alkyl substituents, and it was shown that the stabilization energy for large substituents is dominated by dispersion interactions. In contrast to traditional chemical explanations that attribute stability and condensed matter properties solely to hydrogen bonding and, thus, to the properties of the atoms forming the hydrogen bridge, we show that without dispersion interactions not even the stability and structure of the ammonia dimer can be correctly described. The stability of amine dimers depends crucially on the interaction between the non-polar alkyl groups, which is dominated by dispersion interactions. This interaction is also responsible for the energetic part of the free energy interaction used to describe hydrophobic interactions in liquid alkanes. The entropic part has its origin in the high degeneracy of the interaction energy for complexes of alkane molecules, which exist in a great variety of conformers, having their origin in internal rotations of the alkane chains.

In Silico Studies of Small Molecule Interactions with Enzymes Reveal Aspects of Catalytic Function

Small molecules, such as solvent, substrate, and cofactor molecules, are key players in enzyme catalysis. Computational methods are powerful tools for exploring the dynamics and thermodynamics of these small molecules as they participate in or contribute to enzymatic processes. In-depth knowledge of how small molecule interactions and dynamics influence protein conformational dynamics and function is critical for progress in the field of enzyme catalysis. Although numerous computational studies have focused on enzyme-substrate complexes to gain insight into catalytic mechanisms, transition states and reaction rates, the dynamics of solvents, substrates, and cofactors are generally less well studied. Also, solvent dynamics within the biomolecular solvation layer play an important part in enzyme catalysis, but a full understanding of its role is hampered by its complexity. Moreover, passive substrate transport has been identified in certain enzymes, and the underlying principles of molecular recognition are an area of active investigation. Enzymes are highly dynamic entities that undergo different conformational changes, which range from side chain rearrangement of a residue to larger-scale conformational dynamics involving domains. These events may happen nearby or far away from the catalytic site, and may occur on different time scales, yet many are related to biological and catalytic function. Computational studies, primarily molecular dynamics (MD) simulations, provide atomistic-level insight and site-specific information on small molecule interactions, and their role in conformational pre-reorganization and dynamics in enzyme catalysis. The review is focused on MD simulation studies of small molecule interactions and dynamics to characterize and comprehend protein dynamics and function in catalyzed reactions. Experimental and theoretical methods available to complement and expand insight from MD simulations are discussed briefly.

Complete protein-protein association kinetics in atomic detail revealed by molecular dynamics simulations and Markov modelling

Protein-protein association is fundamental to many life processes. However, a microscopic model describing the structures and kinetics during association and dissociation is lacking on account of the long lifetimes of associated states, which have prevented efficient sampling by direct molecular dynamics (MD) simulations. Here we demonstrate protein-protein association and dissociation in atomistic resolution for the ribonuclease barnase and its inhibitor barstar by combining adaptive high-throughput MD simulations and hidden Markov modelling. The model reveals experimentally consistent intermediate structures, energetics and kinetics on timescales from microseconds to hours. A variety of flexibly attached intermediates and misbound states funnel down to a transition state and a native basin consisting of the loosely bound near-native state and the tightly bound crystallographic state. These results offer a deeper level of insight into macromolecular recognition and our approach opens the door for understanding and manipulating a wide range of macromolecular association processes.

Perturbation of hydration layer in solvated proteins by external electric and electromagnetic fields: Insights from non-equilibrium molecular dynamics

Given the fundamental role of water in governing the biochemistry of enzymes, and in regulating their wider biological activity (e.g., by local water concentration surrounding biomolecules), the influence of extraneous electric and electromagnetic (e/m) fields thereon is of central relevance to biophysics and, more widely, biology. With the increase in levels of local and atmospheric microwave-frequency radiation present in modern life, as well as other electric-field exposure, the impact upon hydration-water layers surrounding proteins, and biomolecules generally, becomes a particularly pertinent issue. Here, we present a (non-equilibrium) molecular-dynamics-simulation study on a model protein (hen egg-white lysozyme) hydrated in water, in which we determine, inter alia, translational self-diffusivities for both hen egg-white lysozyme and its hydration layer together with relaxation dynamics of the hydrogen-bond network between the protein and its hydration-layer water molecules on a residue-per-residue basis. Crucially, we perform this analysis both above and below the dynamical-transition temperature (at ∼220 K), at 300 and 200 K, respectively, and we compare the effects of external static-electric and e/m fields with linear-response-régime (r.m.s.) intensities of 0.02 V/Å. It was found that the translational self-diffusivity of hen egg-white lysozyme and its hydration-water layer are increased substantially in static fields, primarily due to the induced electrophoretic motion, whilst the water-protein hydrogen-bond-network-rearrangement kinetics can also undergo rather striking accelerations, primarily due to the enhancement of a larger-amplitude local translational and rotational motion by charged and dipolar residues, which serves to promote hydrogen-bond breakage and re-formation kinetics. These external-field effects are particularly evident at 200 K, where they serve to induce the protein- and solvation-layer-response effects redolent of dynamical transition at a lower temperature (∼200 K) vis-à-vis the zero-field case (∼220 K).

Mechanism of the Association Pathways for a Pair of Fast and Slow Binding Ligands of HIV-1 Protease

Equilibrium constants, together with kinetic rate constants of binding, are key factors in the efficacy and safety of drug compounds, informing drug design. However, the association pathways of protein-ligand binding, which contribute to their kinetic behaviors, are little understood. In this work, we used unbiased all-atom molecular dynamics (MD) simulations with an explicit solvent model to study the association processes of protein-ligand binding. Using the HIV protease (HIVp)-xk263 and HIVp-ritonavir protein-ligand systems as cases, we observed that ligand association is a multistep process involving diffusion, localization, and conformational rearrangements of the protein, ligand, and water molecules. Moreover, these two ligands preferred different routes of binding, which reflect two well-known binding mechanisms: induced-fit and conformation selection models. Our study shows that xk263 has a stronger capacity for desolvating surrounding water molecules, thereby inducing a semiopen conformation of the HIVp flaps (induced-fit model). In contrast, the slow dehydration characteristic of ritonavir allows for gradual association with the binding pocket of HIVp when the protein's flap conformation is fully open (conformation selection model). By studying the mechanism of ligand association and understanding the role of solvent molecules during the binding event, we can obtain a different perspective on the mechanism of macromolecule recognition, providing insights into drug discovery.