Molecular Mechanisms of Drug Action: X-ray Crystallography at the Basis of Structure-based and Ligand-based Drug Design

Biological systems are recognized for their complexity and diversity and yet we sometimes manage to cure disease via the administration of small chemical drug molecules. At first, active ingredients were found accidentally and at that time there did not seem a need to understand the molecular mechanism of drug functioning. However, the urge to develop new drugs, the discovery of multipurpose characteristics of some drugs, and the necessity to remove unwanted secondary drug effects, incited the pharmaceutical sector to rationalize drug design. This did not deliver success in the years directly following its conception, but it drove the evolution of biochemical and biophysical techniques to enable the characterization of molecular mechanisms of drug action. Functional and structural data generated by biochemists and structural biologists became a valuable input for computational biologists, chemists and bioinformaticians who could extrapolate in silico, based on variations in the structural aspects of the drug molecules and their target. This opened up new avenues with much improved predictive power because of a clearer perception of the role and impact of structural elements in the intrinsic affinity and specificity of the drug for its target. In this chapter, we review how crystal structures can initiate structure-based drug design in general.

Balance between Coiled-Coil Stability and Dynamics Regulates Activity of BvgS Sensor Kinase in Bordetella

The two-component system BvgAS controls the expression of the virulence regulon of Bordetella pertussis. BvgS is a prototype of bacterial sensor kinases with extracytoplasmic Venus flytrap perception domains. Following its transmembrane segment, BvgS harbors a cytoplasmic Per-Arnt-Sim (PAS) domain and then a predicted 2-helix coiled coil that precede the dimerization-histidine-phosphotransfer domain of the kinase. BvgS homologs have a similar domain organization, or they harbor only a predicted coiled coil between the transmembrane and the dimerization-histidine-phosphotransfer domains. Here, we show that the 2-helix coiled coil of BvgS regulates the enzymatic activity in a mechanical manner. Its marginally stable hydrophobic interface enables a switch between a state of great rotational dynamics in the kinase mode and a more rigid conformation in the phosphatase mode in response to signal perception by the periplasmic domains. We further show that the activity of BvgS is controlled in the same manner if its PAS domain is replaced with the natural α-helical sequences of PAS-less homologs. Clamshell motions of the Venus flytrap domains trigger the shift of the coiled coil's dynamics. Thus, we have uncovered a general mechanism of regulation for the BvgS family of Venus flytrap-containing two-component sensor kinases.

IMPORTANCE

The two-component system BvgAS of the whooping cough agent Bordetella pertussis regulates the virulence factors necessary for infection in a coordinated manner. BvgS is the prototype of a family of sensor kinase proteins found in major bacterial pathogens. When BvgS functions as a kinase, B. pertussis is virulent, and the bacterium shifts to an avirulent phase after BvgS senses chemicals that make it switch to phosphatase. Our goal is to decipher the signaling mechanisms of BvgS in order to understand virulence regulation in Bordetella, which may lead to new antimicrobial treatments targeting those two-component systems. We discovered that the activity of BvgS is regulated in a mechanical manner. A short region of the protein that precedes the enzymatic domain switches between two states in response to signal perception by other BvgS domains. This switch region is conserved among BvgS homologs, and thus, the regulation uncovered here will likely be relevant for the family.

Oligoarginine vectors for intracellular delivery: role of arginine side-chain orientation in chain length-dependent destabilization of lipid membranes

Arginine-rich peptides receive increased attention due to their capacity to cross different types of membranes and to transport cargo molecules inside cells. Even though peptide-induced destabilization has been investigated extensively, little is known about the peptide side-chain and backbone orientation with respect to the bilayer that may contribute to a molecular understanding of the peptide-induced membrane perturbations. The main objective of this work is to provide a detailed description of the orientation of arginine peptides in the lipid bilayer of PC and negatively charged PG liposomes using ATR-IR spectroscopy and molecular modeling, and to relate these orientational preferences to lipid bilayer destabilization. Molecular modeling showed that above the transition temperature arginine side-chains are preferentially solvent-directed at the PC/water interface whereas several arginine side-chains are pointing towards the PG hydrophobic core. IR dichroic spectra confirmed the orientation of the arginine side chains perpendicular to the lipid-water interface. IR spectra shows an randomly distributed backbone that seems essential to optimize interactions with the lipid membrane. The observed increase of permeation to a fluorescent dye is related to the peptide induced-formation of gauche bonds in the acyl chains. In the absence of hydrophobic residues, insertion of side-chains that favors phosphate/guanidium interaction is another mechanism of membrane permeabilization that has not been further analyzed so far.

Recognition-induced conformational changes in protein-protein docking

The ability to predict the three-dimensional structure of a protein complex starting from the isolated binding partners is becoming increasingly relevant. As our understanding of the molecular mechanisms behind protein-protein binding improves, so do the docking methods, however, it remains a challenge to adequately predict the unbound to bound transition. Side-chain flexibility is routinely handled and most docking methods allow for a certain degree of backbone flexibility, but systems undergoing moderate to large conformational changes can at present not correctly be modeled. The docking community is therefore putting an increased effort in the treatment of protein flexibility. Here we present a survey of the existing computational techniques to model protein flexibility in the context of protein-protein docking.

Penetratin-membrane association: W48/R52/W56 shield the peptide from the aqueous phase

Using molecular dynamics simulations, we studied the mode of association of the cell-penetrating peptide penetratin with both a neutral and a charged bilayer. The results show that the initial peptide-lipid association is a fast process driven by electrostatic interactions. The homogeneous distribution of positively charged residues along the axis of the helical peptide, and especially residues K46, R53, and K57, contribute to the association of the peptide with lipids. The bilayer enhances the stability of the penetratin helix. Oriented parallel to the lipid-water interface, the subsequent insertion of the peptide through the bilayer headgroups is significantly slower. The presence of negatively charged lipids considerably enhances peptide binding. Lateral side-chain motion creates an opening for the helix into the hydrophobic core of the membrane. The peptide aromatic residues form a pi-stacking cluster through W48/R52/W56 and F49/R53, protecting the peptide from the water phase. Interaction with the penetratin peptide has only limited effect on the overall membrane structure, as it affects mainly the conformation of the lipids which interact directly with the peptide. Charge matching locally increases the concentration of negatively charged lipids, lateral lipid diffusion locally decreases. Lipid disorder increases, through decreased order parameters of the lipids interacting with the penetratin side chains. Penetratin molecules at the membrane surface do not seem to aggregate.

A conserved arginine plays a role in the catalytic cycle of the protein disulphide isomerases

The pK(a) values of the CXXC active-site cysteine residues play a critical role in determining the physiological function of the thioredoxin superfamily. To act as an efficient thiol-disulphide oxidant the thiolate state of the N-terminal cysteine must be stabilised and the thiolate state of the C-terminal cysteine residue destabilised. While increasing the pK(a) value of the C-terminal cysteine residue promotes oxidation of substrates, it has an inhibitory effect on the reoxidation of the enzyme, which is promoted by the formation of a thiolate at this position. Since reoxidation is essential to complete the catalytic cycle, the differential requirement for a high and a low pK(a) value for the C-terminal cysteine residue for different steps in the reaction presents us with a paradox. Here, we report the identification of a conserved arginine residue, located in the loop between beta5 and alpha4 of the catalytic domains of the human protein disulphide isomerase (PDI) family, which is critical for the catalytic function of PDI, ERp57, ERp72 and P5, specifically for reoxidation. An examination of the published NMR structure for the a domain of PDI combined with molecular dynamic studies suggest that the side-chain of this arginine residue moves into and out of the active-site locale and that this has a very marked effect on the pK(a) value of the active-site cysteine residues. This intra-domain motion resolves the apparent dichotomy of the pK(a) requirements for the C-terminal active-site cysteine.

Response of SCP-2L domain of human MFE-2 to ligand removal: binding site closure and burial of peroxisomal targeting signal

In the study of the structure and function relationship of human MFE-2, we have investigated the dynamics of human MFE-2SCP-2L (hSCP-2L) and its response to ligand removal. A comparison was made with homologous rabbit SCP-2. Breathing and a closing motion are found, identifiable with an adjustment in size and a closing off of the binding pocket. Crucial residues for structural integrity have been identified. Particularly mobile areas of the protein are loop 1 that is connecting helices A and C in space, and helix D, next to the entrance of the pocket. In hSCP-2L, the binding pocket gets occupied by Phe93, which is making a tight hydrophobic contact with Trp36. In addition, it is found that the C-terminal peroxisomal targeting signal (PTS1) that is solvent exposed in the complexed structure becomes buried when no ligand is present. Moreover, an anti-correlation exists between burial of PTS1 and the size of the binding pocket. The results are in accordance with plant nsLTPs, where a similar accommodation of binding pocket size was found after ligand binding/removal. Furthermore, the calculations support the suggestion of a ligand-assisted targeting mechanism.