Resolution of discordant HIV-1 protease resistance rankings using molecular dynamics simulations

The emergence of drug resistance is a major challenge for the effective treatment of HIV. In this article, we explore the application of atomistic molecular dynamics simulations to quantify the level of resistance of a patient-derived HIV-1 protease sequence to the inhibitor lopinavir. A comparative drug ranking methodology was developed to compare drug resistance rankings produced by the Stanford HIVdb, ANRS, and RegaDB clinical decision support systems. The methodology was used to identify a patient sequence for which the three rival online tools produced differing resistance rankings. Mutations at only three positions ( L10I , A71IV, and L90M ) influenced the resistance level assigned to the sequence. We use ensemble molecular dynamics simulations to elucidate the origin of these discrepancies and the mechanism of resistance. By simulating not only the full patient sequences but also systems containing the constituent mutations, we gain insight into why resistance estimates vary and the interactions between the various mutations. In the same way, we also gain valuable knowledge of the mechanistic causes of resistance. In particular, we identify changes in the relative conformation of the two beta sheets that form the protease dimer interface which suggest an explanation of the relative frequency of different amino acids observed in patients at residue 71.     

Analyzing the molecular basis of enzyme stability in ethanol/water mixtures using molecular dynamics simulations

One of the drawbacks of nonaqueous enzymology is the fact that enzymes tend to be less stable in organic solvents than in water. There are, however, some enzymes that display very high stabilities in nonaqueous media. In order to take full advantage of the use of nonaqueous solvents in enzyme catalysis, it is essential to elucidate the molecular basis of enzyme stability in these media. Toward this end, we performed μs-long molecular dynamics simulations using two homologous proteases, pseudolysin, and thermolysin, which are known to have considerably different stabilities in solutions containing ethanol. The analysis of the simulations indicates that pseudolysin is more stable than thermolysin in ethanol/water mixtures and that the disulfide bridge between C30 and C58 is important for the stability of the former enzyme, which is consistent with previous experimental observations. Our results indicate that thermolysin has a higher tendency to interact with ethanol molecules (especially through van der Waals contacts) than pseudolysin, which can lead to the disruption of intraprotein hydrophobic interactions and ultimately result in protein unfolding. In the absence of the C30-C58 disulfide bridge, pseudolysin undergoes larger conformational changes, becoming more open and more permeable to ethanol molecules which accumulate in its interior and form hydrophobic interactions with the enzyme, destroying its structure. Our observations are not only in good agreement with several previous experimental findings on the stability of the enzymes studied in ethanol/water mixtures but also give an insight on the molecular determinants of this stability. Our findings may, therefore, be useful in the rational development of enzymes with increased stability in these media.     

Characterizing the dynamics and ligand-specific interactions in the human leukocyte elastase through molecular dynamics simulations

The human leukocyte elastase (HLE), a neutrophil serine protease of the chymotrypsin superfamily, is a major therapeutic target for a number of inflammatory diseases, such as chronic obstructive pulmonary disease (COPD). In this work, we present a comparative explicit water molecular dynamics (MD) study on the free and inhibitor-bound HLE. Knowledge of the flexibility and conformational changes induced by this irreversible inhibitor, whether in a prebound stage or covalently bound at the enzyme binding site, encases fundamental biological interest and is particularly relevant to ongoing structure-based drug design studies. Our results suggest that HLE operates by an induced-fit mechanism with direct intervention of a surface loop which is open toward the solvent in the free enzyme and closed while in the presence of the ligand. MM-PBSA free energy calculations furthermore elucidate the energetic contributions to the distinct conformations adopted by this loop. Additionally, a survey of the major contributions to the inhibitor binding free energies was attained. Our findings enforce the need to account for HLE flexibility, whether through the use of MD-generated ensembles of HLE conformations as targets for molecular docking or via sophisticated flexible-docking algorithms. We anticipate that inclusion of the observed HLE dynamic behavior into future drug design methodologies will have a relevant impact in the development of novel, more efficient, inhibitors.     

Understanding the nonfouling mechanism of surfaces through molecular simulations of sugar-based self-assembled monolayers

This paper presents a molecular simulation study of the interactions of a protein (lysozyme) with self-assembled monolayers (SAMs) of mannitol and sorbitol terminated alkanethiols in the presence of explicit water molecules and ions. The all-atom simulations were performed to calculate the force generated on the protein as a function of its distance above the SAM surfaces. The structural and dynamic properties of water molecules both above the SAM surfaces and around the SAM head groups were analyzed to provide a better understanding of the nonfouling behavior of the sugar-based SAM surfaces. Results from this work suggest that both mannitol and sorbitol SAMs generate a tightly bound, structured water layer around the SAM chains. This hydration layer creates a repulsive force on the protein when it approaches the surface, resulting in a nonfouling surface despite the presence of hydrogen-bond donor groups. This work demonstrates the importance of strong surface-water interactions for surface resistance to nonspecific protein adsorption.     

Sommelet-hauser rearrangement of an ammonium ylide derived from the HIV-1 reverse transcriptase inhibitor nevirapine

Functionalization at the 3-position of the dipyridodiazepinone nevirapine ( 1 ) has been accomplished by Sommelet-Hauser rearrangement of an ylide derived from 1 . Treatment of N-cyanomethylpyrrolidinium salt 4 with potassium tert-butoxide in a mixture of dimethylsulfoxide and tetrahydrofuran at ?10°, followed by acid hydrolysis, afforded a mixture of compounds 5 and 6 in a ratio of 1:1.8. Upon treatment of 4 with sodium amide in liquid ammonia, 5 and 6 were obtained in a ratio of 1.5:1 and a combined yield of 83%. Compound 5 is the desired product resulting from Sommelet-Hauser rearrangement of 4 , whereas 6 derives from competing Stevens rearrangement and intramolecular cyclization of the aldehyde produced upon hydrolysis. Baeyer-Villiger oxidation of 5 afforded the 3-hydroxy derivative 2 , a recently identified metabolite of nevirapine.     

The effect of the counterion on water mobility in reverse micelles studied by molecular dynamics simulations

In this study, mobility and structure of water molecules in Aerosol OT (bis(2-ethylhexyl) sulfosuccinate, AOT) reverse micelles with water content w0 = 5 and Na+, K+, Cs+ counterions have been explored with molecular dynamics (MD) simulations. Using the Faeder/Ladanyi model (J. Phys. Chem. B, 2000, 104, 1033) of the reverse micelle interior, MD simulations were performed to calculate the self-intermediate scattering function, FS(Q,t), for water hydrogen atoms that could be measured in a quasielastic neutron scattering experiment. Separate intermediate scattering functions FRS(Q,t) and FCMS(Q,t) were determined for rotational and translational motion. We find that the decay of FCMS(Q,t) is nonexponential and our analysis of the MD data indicates that this behavior arises from decreased water mobility for molecules close to the interface and from confinement-induced restrictions on the range of translational displacements. Rotational relaxation also exhibits nonexponential decay, which is consistent with relatively rapid restricted rotation and slower rotational relaxation over the full angular range. Rotational relaxation is anisotropic, with the O-H bond short-time rotational mobility considerably higher than that of the molecular dipole. This behavior is related to the decreased density of water-water hydrogen bonds in the vicinity of the interface compared to core or bulk water. We find that the interfacial mobility of water molecules is quite different for the three counterion types, but that the core mobility exhibits weak counterion dependence. Differences in interfacial mobility are strongly correlated with structural features, especially ion-water coordination, and the extent of disruption by the counterions of the water hydrogen bond network.     

Structure, dynamics and solvation of HIV-1 protease/saquinavir complex in aqueous solution and their contributions to drug resistance: molecular dynamic simulations

As it is known that the understanding of the basic properties of the enzyme/inhibitor complex leads directly to enhancing the capability in drug designing and drug discovery. Molecular dynamics simulations have been performed to examine detailed information on the structure and dynamical properties of the HIV-1 PR complexed with saquinavir in the three protonated states, monoprotonates at Asp25 (Mono-25) and Asp25'(Mono-25') and diprotonate (Di-Pro) at both Asp25 and Asp25'. The obtained results support clinical data which reveal that Ile84 and Gly48 are two of the most frequent residues where mutation toward a protease inhibitor takes place. In contrast to the Ile84 mutation due to high displacement of Ile84 in the presence of saquinavir, source of the Gly48 mutation was observed to be due to the limited space in the HIV-1 PR pocket. The Gly48 was, on one side, found to form strong hydrogen bonds with saquinavir, while on the other side this residue was repelled by the hydrophobic Phe53 residue. In terms of inhibitor/enzyme binding, interactions between saquinavir and a catalytic triad of the HIV-1 PR were calculated using the ab initio method. The results show an order of the binding energy of Mono25相似文献