Cystic fibrosis (CF) is a recessive genetic disease that is caused by mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) protein. The recent development of a class of drugs ...called “correctors”, which repair the structure and function of mutant CFTR, has greatly enhanced the life expectancy of CF patients. These correctors target the most common disease causing CFTR mutant F508del and are exemplified by the FDA-approved VX-809. While one binding site of VX-809 to CFTR was recently elucidated by cryo-electron microscopy, four additional binding sites have been proposed in the literature and it has been theorized that VX-809 and structurally similar correctors may engage multiple CFTR binding sites. To explore these five binding sites, ensemble docking was performed on wild-type CFTR and the F508del mutant using a large library of structurally similar corrector drugs, including VX-809 (lumacaftor), VX-661 (tezacaftor), ABBV-2222 (galicaftor), and a host of other structurally related molecules. For wild-type CFTR, we find that only one site, located in membrane spanning domain 1 (MSD1), binds favorably to our ligand library. While this MSD1 site also binds our ligand library for F508del-CFTR, the F508del mutation also opens a binding site in nucleotide binding domain 1 (NBD1), which enables strong binding of our ligand library to this site. This NBD1 site in F508del-CFTR exhibits the strongest overall binding affinity for our library of corrector drugs. This data may serve to better understand the structural changes induced by mutation of CFTR and how correctors bind to the protein. Additionally, it may aid in the design of new, more effective CFTR corrector drugs.
Carbohydrate recognition by proteins, such as lectins and other (bio)molecules, can be essential for many biological functions. Recently, interest has arisen due to potential protein and drug design ...and future bioengineering applications. A quantitative measurement of carbohydrate-protein interaction is thus important for the full characterization of sugar recognition. We focus on the aspect of utilizing computer simulations and biophysical models to evaluate the strength and specificity of carbohydrate recognition in this review. With increasing computational resources, better algorithms and refined modeling parameters, using state-of-the-art supercomputers to calculate the strength of the interaction between molecules has become increasingly mainstream. We review the current state of this technique and its successful applications for studying protein-sugar interactions in recent years.
Understanding allosteric mechanisms is essential for the physical control of molecular switches and downstream cellular responses. However, it is difficult to decode essential allosteric motions in a ...high-throughput scheme. A general two-pronged approach to performing automatic data reduction of simulation trajectories is presented here. The first step involves coarse-graining and identifying the most dynamic residue–residue contacts. The second step is performing principal component analysis of these contacts and extracting the large-scale collective motions expressed via these residue–residue contacts. We demonstrated the method using a protein complex of nuclear receptors. Using atomistic modeling and simulation, we examined the protein complex and a set of 18 glycine point mutations of residues that constitute the binding pocket of the ligand effector. The important motions that are responsible for the allostery are reported. In contrast to conventional induced-fit and lock-and-key binding mechanisms, a novel “frustrated-fit” binding mechanism of RXR for allosteric control was revealed.
The promiscuous protein retinoid X receptor (RXR) displays essential allosteric regulation of several members in the nuclear hormone receptor superfamily via heterodimerization and (anti)cooperative ...binding of cognate ligands. Here, the structural basis of the positive allostery of RXR and constitutive androstane receptor (CAR) is revealed. In contrast, a similar computational approach had previously revealed the mechanism for negative allostery in the complex of RXR and thyroid receptor (TR). By comparing the positive and negative allostery of RXR complexed with CAR and TR respectively, we reported the promiscuous allosteric control involving RXR. We characterize the allosteric mechanism by expressing the correlated dynamics of selected residue–residue contacts which was extracted from atomistic molecular dynamics simulation and statistical analysis. While the same set of residues in the binding pocket of RXR may initiate the residue–residue interaction network, RXR uses largely different sets of contacts (only about one-third identical) and allosteric modes to regulate TR and CAR. The promiscuity of RXR control may originate from multiple factors, including (1) the frustrated fit of cognate ligand 9c to the RXR binding pocket and (2) the different ligand-binding features of TR (loose) versus CAR (tight) to their corresponding cognate ligands.
The amphiphilic peptide of the triacylglycerol lipase derived from Pseudomonas aeruginosa plays a critical role in guarding the gate for ligand access. Conformations of this peptide at several ...water–oil interfaces and in protein environments were compared using atomistic simulations with explicit solvents. In oil-containing solvents, this peptide is able to retain a folded structure. Interestingly, when the peptide is immersed in a low-polarity solvent environment, it exhibits a “coalesced” helix structure, which has both α- and 310-helix components. The observation that the 310-helical conformation is populated in a highly nonpolar environment is consistent with a previous report on polymethylalanine. Frequent interconversions of the secondary structure (between α-helix and 310-helix) of the peptide are also observed. We further studied how this solvent-induced structural transition may be connected to the trigger mechanism of lipase gating and how the lipase senses the hydrophobic–hydrophilic interface.