Gated Binding of Ligands to HIV-1 Protease: Brownian Dynamics Simulations in a Coarse-Grained Model

ABSTRACT

The internal motions of proteins may subserve as a "gate" in a certain quantity of systems, which controls ligand-protein association. This research applies Brownian dynamics simulations in a coarse-grained archetype to study the gated association rate constants of HIV-1 proteases and medicines The computed gated association rate constants of three protease mutants, G48V/V82A/I84V/L90M, G48V and L90M with three physics amprenavir, indinavir, and saquinavir, yield serviceable agreements with experiments. The work exhibits that the flap dynamics leads to "slow gating". The simulations allude to that the flap flexibility and the opening frequent occurrence of the wild-type, the G48V and L90M mutants are similar, on the contrary the flaps of the variant G48V/V82A/I84V/L90M render free of access less frequently, resulting in a lower gated rate constant. The make knowned methodology is fast and provides an efficient way to predict the gated association rate constants for various protease mutants and ligands.

INTRODUCTION

The important first pace in many biological processes is the clash of protein-protein or protein-ligand indivisible particles For example, in the combination of parts to form a whole of human immunodeficiency virus emblem 1 (HIV-1), the viral polyprotein has to bind to the active site of the HIV-1 protease to bring forward active structural and replicative proteins (1) HIV-1 protease remedys also have to associate with the protease and typically cope with the polyprotein binding. The binding begins with the rencounter of the two molecules through random walks, and the association presumably appears at a rate that approaches the diffusioncontrolled limit (2-4) However, in the HIV-1 protease combination of parts to form a whole mentioned above, the experimentally measured association rate constant is a not many orders of magnitude smaller than that wait fored if the molecules were uniformly reactive spheres (56) Deviation of this limiting rate may spring from geometric constraints of the binding sites, interaction potentials, hydrodynamic interactions, and the binding site accessibility owed to protein internal motions (7-11) brawny electrostatic steering may enhance the rate of association, and the requirement of high steric specificity of the two molecules may reduce the association rate (312) In a certain quantity of cases, the association rate constants may be further lowered by means of the dynamic nature of a protein, which modulates the binding site accessibility and can be viewed as a "gated" binding site (1314)

A number of experimental works exhibit that some enzymes have bights over their active sites whose conformational changes can "gate" ligand binding (15-17) For example, the active site bight of triosephosphate isomerase (TIM) exhibits a hinged-lid motion, which alternates between the sum of two units well-defined open and closed conformations (18) The fluctuations of the protein chain of myoglobin be subservient to as a gate that render free of accesss and closes the reaction binding site, and permit small ligands (O2 and CO) to enterthe heme "pocket" (1920) The general [i]or[/i] abstract notion of gating has been applied to research association rate constants coupled to conformational changes in the active sites of enzyme one as well as the other in equilibrium and out of equilibrium (21-26) Earlier theoretical works use the kinetics schemes of gate opening and closing, together with the characteristic diffusional relaxation time of the ligand-protein a whole to estimate gated association rate constants (14) Computer simulations have also been applied directly to estimate the gating events on protein-ligand association, i.e., with explicit modeling of the gate motion (2728) notwithstanding that such studies are challenging to be paid to the limited contigurational sampling that is possible with typical simulation techniques.

This close attention focuses on the fluctuations of HIV-1 protease flaps, which modulate the binding site accessibility. We studied wild-type protease and three mutants, G48V/V82A/ I84V/L90M, L90M and G48V with three clinically used remedys amprenavir, indinavir, and saquinavir, where experimental data are available (29) The work uses a coarsegrained original to simulate substantial protein conformational changes (30) and Brownian dynamics to perform microsecond (/as) timescale simulations (3132) The use of a flexible force field allows inquiry of the opening/closing of the flaps of the active site. Applying of that kind a flexible force field is necessary, as simpler Go-like prototypes cannot describe the complicated flap dynamics (30) The comput association rate constants of HIV-I variants yield profitable agreements with experimental values.

Theoretical background

However, if the ligand-protein interactions include realistic descriptions of the irregular surface topography and the molecular interactions, the above equation is not accurate, and individual must rely on computer simulations, as described in previous publications (734-36)

Simulation manners

Coarse-grained model of HIV-1 protease

This is essentially the Stokes-Einstein equation, and the hydrodynamic radius is approximated by the agency of the effective radius of each bead, r^sub i^ plus the water radius, 14 [Angstrom], used in this application of mind (31,32). The viscosity of water ?· is 1 cp (T = 293 K)